Edit main.py

main.py

@@ -30,9 +30,11 @@ if __name__ == "__main__":
     qc_object = Concat(output)
     qc_sequence = Sequence(qc_object.get_sequence())
 
-    with open(args.out, "w") as output_file:
-        output_file.write("Read_ID\tScore\tQuery_Align\tGenome_Align\n")  # En-tête du fichier TSV
+    resultat = []
+    count = 0
 
+    with open(args.out, "w") as output_file:
+        output_file.write("ALIGNMENT : \n\n")
         kmers = qc_sequence.kmers(args.taille_kmers)
         for n in range(genome.__len__()):
             read = genome.get_read(n)
@@ -40,8 +42,16 @@ if __name__ == "__main__":
             result = alignement.align_sequences()
 
             if result is not None:
+                count += 1
                 score, aligned_query, aligned_genome = result
-                output_file.write(f"{score}\t{aligned_query}\t{aligned_genome}\n")
+
+                resultat.append(f"score = {score}")
+                resultat.append(f"genome number = {n}\n")
+                resultat.append(f"Query      {aligned_query}")
+                resultat.append(f"Genome     {aligned_genome}\n\n\n" )
+
+            output_file.write(f"total alignment number = {count}\n")
+            output_file.write(f"{resultat}")
 
         print(f"Alignment done ! Results saved at {args.out}")