Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	71,644	T→C	L98L (CTT→CTC)	yabI →	DedA family protein YabI

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	71,644	0	T	C	97.8%	151.4 / NA	46	L98L (CTT→CTC)	yabI	DedA family protein YabI
Reads supporting (aligned to +/- strand):  ref base T (0/0);  major base C (22/23);  minor base A (1/0);  total (23/23)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.54e-01

ATTGATGGCGGTGGTGTTGGTGGCCTTTCTGGAGTCGCTGGCGCTGGTCGGTTTGATTCTACCCGGTACGGTGCTGATGGCGGGGCTGGGAGCGCTGATTGGCAGCGGCGAGTTAAGTTTCTGGCACGCCTGGCTGGCAGGGATTATTGGCTGCTTGATGGGCG‑ACTGGATTTCTTTCTGGCTGGGTTGGCGTTTTAAAAAGCCGTTGCATCGCTGGTCATTTCTG‑AAGAAAAACAAAGCACTACTTGATAAAACTGAACATGCGTTGCATCAACACAGCATGTTCACCATTCTGGTCGGTCGTTTTGTTGGCCCGACGCGTCCGCTGGTGCCAATGGTGGCGGGAATGCTGGATCT‑GCCGGTGGCTAAATTTATTACGCCGAATATTATCGGCTGCCTGCTGTGGCCGCCGTTTTACTTCCTGCCAGGGATTCTGGCGGGCGCGGCGATCGATATTCCTGCCGGAATGCAGAGCGGT  >  NC_000913/71398‑71875                                                                                                                                                                                                                                                         |                                                                                                                                                                                                                                         aTTGATGGCGGTGGTGTTGGTGGCCTTTCTGGAGTCGCTGGCTCTGGTCGGTTTGATTCTACCCGGTACGGTGCTGATGGCGGGGCTGGGAGCGCTGATTGGCAGCGGCGAGTTAAGTTTCTGGCACGCCTGGCTGGCAGGGATTGTCGGTTGTCTGCTGGGCG‑ACTGGATCTCCTTCTGGCTGGGCTGGCGTTTTAAAAAGCCGTTGCATCGCTGGTCATTTCTG‑AAGAAAAACAAAGCATTACTCGATa                                                                                                                                                                                                                                      <  1:375760/251‑1 (MQ=255)   ggATGGCGGTGGTGTTGGTGGCCTTTCTGGAGTCGCTGGCGCTGGTCGGTTTGATTCTACCCGGTACGGTGCTGATGGCGGGGCTGGGAGCGCTGATTGGCAGCGGCGAGTTAAGTTTCTGGCACGCCTGGCTGGCAGGGATTGTCGGTTGTCTGCTGGGCG‑ACTGGATCTCCTTCTTGCTGGGCTGGCGTTTTAAAAAGCCGTTGCATCGCTGGTCATTTCTG‑AATAAAAACAAAGCATTACTCGATaaa                                                                                                                                                                                                                                    <  1:465011/250‑1 (MQ=255)         cGGTGGTGTTGGTGGCCTTTCTGGAGGCGCTGGCGCTGGTCGGTTTGATTCTACCCGGTACGGTGCTGATGGCGGGGCTGGGAGCGCTGATTGGCAGCGGCGAGTTAAGTTTCTGGCACGCCTGGCTGGCAGGGATTGTCGGTTGTCTGCTGGGCG‑ACTGGATCTCCTTCTGGCTGGGCTGGCGTTTTAAAAAGCCGTTGCATCGCTGGTCATTTCTG‑AAGAAAAACAAAGCATTACTCGATAAAACTGaa                                                                                                                                                                                                                              <  1:692147/251‑1 (MQ=255)                    gtggCCTTTCTGGAGTCGCTGGCGCTGGTCGGTTTGATTCTACCCGGTACGGTGCTGATGGCGGGGCTGGGAGCGCTGATTGGCAGCGGCGAGTTAAGTTTCTGGCACGCCTGGCTGGCAGGGATTGTCGGTTGTCTGCTGGGCG‑ACTGGATCTCCTTCTGGCTGGGCTGGCGTTTTAAAAAGCCGTTGCATCGCTGGTCATTTCTG‑AAGAAAAACAAAGCATTACTCGATAAAACTGAACATGCGTTGCa                                                                                                                                                                                                                   <  1:132109/251‑1 (MQ=255)                              tGGAGTCGCTGGCGCTGGTCGGTTTGATTCTACCCGGTACGGTGCTGATGGCGGGGCTGGGAGCGCTGATTGGCAGCGGCGAGTTAAGTTTCTGGCACGCCTGGCTGGCAGGGGTTGTCGGTTGTCTGCTGGGCG‑GCTGGGTCTCCTTCTGGCTGGGCTGGCGTTTTTAAAAGCCGGTGCATCGCTGGTCATTTCTG‑AAGAAAAAAAAAGCATTACTCGATAAAACTGAACATGCCGTGCATCAACACAGc                                                                                                                                                                                                         >  1:427604/1‑251 (MQ=255)                              tGGAGTCGCTGGCGCTGGCCGGTTTGATTCTCCCCGGTACGGTGCTGATGGCGGGGCTGGGAGCGCTGTTTGGCAGCGGCGGGTTAAGTTTCTGGCACGCCTGGCTGGCGGGGATTGTCGGTTGTCTGCTGGGGG‑ACTGGATTTCCTTCTGGCGGGGCTGGCGTTTTAAAAAGCCGTTGCATCGCTGGTCATTTCTG‑AAGAAAAACAAAGCATTACTCGATAAAACTGAACATGCGTTGCTTCAaccacgc                                                                                                                                                                                                         >  1:686860/1‑246 (MQ=255)                                gAGTCGCTGGCGCTGGTCGGTTTGATTCTACACGGTACGGTGCCGATGGCGGGGATGGGGACGCTTATTGGCAGCGGCGAGGTAAGGTTTTGGCACGCCTGGCTGGCAGGGATTGTCGGGTGTCTGGTGGGCGGAC‑GGATCTACTTCTGGCTGGGCTGGAGTTTTAAAAAACCGTTGGATCGCTTGGCATTTCTGCAGGAAAAACAA‑GAATTTCTCGATAAAACTGCACAAGCGATGCATCaaaccaacat                                                                                                                                                                                                       >  1:629131/1‑243 (MQ=255)                                 ggTCGCTGGCGCTGGTCGGTTTGATGCTACCCGGTACGGTGCTGATGGCGGGGCTGGGGGCGCTGATGGGCAGCGGCGAGTTAAGTTTCTGGCACGCCTGGCTGGCAGGGATTGTCGGTTGTCTGCTGGGCG‑ACTGGATCTCCTTCTGGCTGGGCTGGCGTTTTAAAAAGCCGTTGCATCGCTGGTCATTTCTG‑AAGAAAAACAAAGCATTACTCGATAAAACTGAACATGCGTTGCATCAACACAGCATg                                                                                                                                                                                                      <  1:434101/250‑1 (MQ=255)                                    cgctggcgctggTCGGTTTGATTCTGCCCGGGACGGTGCTGATGGCGGGGCTGGGAGCGCTGATTGGCAGCGGCGAGTTAAGTTTCTGGCACGCCTGGCTGGCAGGGATTGTCGGTTGTCTGCTGGGCG‑ACTGGATCTCCTTCTGGCTGGGCTGGCGTTTTAAAAAGCCGTTGCATCGCTGGTCATTTCTG‑AAGAAAAACAAAGCATTACTCGATAAATCTGAACATGCGTTGCATCAACACAGCATGtt                                                                                                                                                                                                    <  1:530345/250‑1 (MQ=255)                                      ctggcgctggTCGGTTTGATTCTACCCGGTACGGTGCTGATGGCGGGGCTGGGAGCGCTGATTGGCAGCGGCGAGTTAAGTTTCTGGCACGCCTGGCTGGGAGGGATTGTCGGGTGGCTGCTGGGCG‑ACTGGATCTCCTTCTGGCTGGGGTGGCGTTTTAAAAAGCCGTTGCATCGCTGGTCATTTCTG‑AAGAAAAAAAAAGCATTAATCGATAAAACTGAACATGCGTTGCATCAACACAGCATGTTCAc                                                                                                                                                                                                 >  1:126273/1‑251 (MQ=255)                                               gTCGGTTTGATTCTACCCGGGACGGTGCTGATGGCGGGGCTGGGAGCGCTGATTGGCAGCGGCGAGGTAAGTTTCTGGCACGCCTGGCTGGCAGGGATTGTCGGTTGTCTGCTGGGCG‑ACTGGATCTCCTTCTGGCTGGGCTGGCGTTTTAAAAAGCCGTTGCATCGCTGGTCATTTCTG‑AAGAAAAACAAAGCATTACGCGATAAAACGGAACATGCGTTGCATCAACACAGCATGTTCACCAGTCTggg                                                                                                                                                                                        <  2:261483/251‑2 (MQ=255)                                                    tttGATTCTACCCGGTACGGTGCTGATGGCGGGGCTGGGAGCGCTGATTGGCAGCGGCGAGTTAAGGTTCTGGCACGCCTGGCTGGCAGGGGTTGTCGGTTGTCTGCTGGGCG‑ACTGGATCTCCTTTTGGCTGGGCTGGCGTTTTTAAAAGCCGTTGCCTCGCTGGTCATTTCTG‑GAGAAAAACAAAGCATTACTCGGTAAAAATGAAAATACGTTGCATCAACACAGCAAGTTCCCCCTTCTggtcggtc                                                                                                                                                                                   >  1:104046/1‑251 (MQ=255)                                                      tGATTCTTCCCGGTACGGTGCTGATGGGGGGGCTGGGGGCGCTGATTGGCAGCGGCGAGTTAAGTTTCGGGCACGCCTGGCTGGCAGGGATTGTCGGTTGTCTGCTGGGCG‑ACTGGATCTCCTTCTGGCTGGGCTGGCGTTTTAAAAAGCCGTTGCATCGCTGGTCATTTCTG‑AAGAAAAACAAAGCATTACTCGATAAAACTGAACATGCGTTGCATCAACACAGCATGTGCACCATTCTGgtcggtcg                                                                                                                                                                                  <  2:384902/250‑1 (MQ=255)                                                                                 gggggCTGGGAGCGCTTGTTGGCTGCGGCGAGTTAAGTTTCTGGGACGCCTGGCTGGCAGGGATTGTCGGTTGTCTGCTGGGCG‑ACTGGATCTCCTGCTGGCTGGGCTGGCGTTTTAAAAAGCCGTTGCATCGCTGGTCATTTCTG‑AAGAAAAACAAAGCATTACTCGATAAAACTGAACATGCGTTGCATCAACACAGCATGTTCACCATTCTGGTCGGTCGCTTTGTTGGCCCGACGCGTCCGCTGGTg                                                                                                                                                      <  2:289168/250‑1 (MQ=255)                                                                                 cGGGGCTGGGAGCGCTGATTGGCAGCGGCGAGTTAAGTTTCTGGCACGCCTGGCTGGCAGGGATTGTCGGTTGTCTGCTGGGCG‑ACTGGATCTCCTTCTGGCTGGGCTGGCGTTTTAAAAAGCCGTTGCATCGCTGGTCATTTCTG‑AAGAAAAACAAAGCATTACTCGATAAAACTGAACATGCGTTGCATCAACACAGCATGTTCACCATTCTGGTCGGTCGCTTTGTTGGCCCGACGCGTCCGCTGGTg                                                                                                                                                      <  1:392368/251‑1 (MQ=255)                                                                                                    tGGCAGCGGCGAGATAAGTTTCTGGCACGCATGGCTGGCAGGGATTGTCGGTTGTCTGCTGGGGG‑ACTGGATCTCCTTCTGGCTGGGCTGGCGTTTTAAAAAGCCGTTGTATCGCTGGTCATTTCTG‑AAGAAAAACAAAGCATTACTCGATAAAACTGAACATGCGTTGCATCAACACAGCATGTTCACCATACTTGGCGGGCGCTTTGGTGGCCCGACGCGACCGCTGGTGCCAATGGGGGCGGGGATgg                                                                                                                                   >  1:483648/1‑250 (MQ=255)                                                                                                        aGCGGCGAGTTAAGTTTCTGGCACGCCTGGCTGGCAGGGATTGTCGGTTGTCTGCTGGGCG‑ACTGGATCTCCTTCTGGCTGGGCTGGCGTTTTAAAAAGCCGTTGCATCGCTGGTCATTTCTG‑AAGAAAAACAAAGCATTACTCGATAAAACTGAACATGCGTTGCATCAACACAGCATGTTCACCATTc                                                                                                                                                                                            >  1:403397/1‑190 (MQ=255)                                                                                                        aGCGGCGAGTTAAGTTTCTGGCACGCCTGGCTGGCAGGGATTGTCGGTTGTCTGCTGGGCG‑ACTGGATCTCCTTCTGGCTGGGCTGGCGTTTTAAAAAGCCGTTGCATCGCTGGTCATTTCTG‑AAGAAAAACAAAGCATTACTCGATAAAACTGAACATGCGTTGCATCAACACAGCATGTTCACCATTc                                                                                                                                                                                            <  2:403397/190‑1 (MQ=255)                                                                                                           ggcgAGTTAAGTTGCGGGCACGCCTGGCTGGCAGGGATTGTCGGTTGTCTGGTGGGCG‑ACTGGATCTCCTTCTGGCTGGGCTGGCGTTTTAAAAAGCCGTTGCATCGCTGGTCATTTCTG‑AAGAAAAACAAAGCATTACTCGATAAAACTGAACATGCGTTGCATCAACACAGCATGTTCACCATTCTGGTCGGTCGCTTTGTTGGCCCGACGCGTCCGCTGGTGCCAATGGTGGCGGGAATGCTGGATCt                                                                                                                            <  2:340538/251‑1 (MQ=255)                                                                                                                           ggCACGCCTGGCTGGCAGGGATTGTCGGTTGTCTGCTGGGCG‑ACTGGATCTCCTTCTGGCTGGGCTGGCGTTTTAAAAAGCCGTTGCATCGCTGGTCATTTCTG‑AAGAAAAACAAAGCATTACTCGATAAAACTGAACATGCGTTGCATCAACACAGCATGTTCACCATTCTGGTCGGTCGCTTTGTTGGCCCGACGCGTCCGCTGGTGCCAATGGTGGCGGGAATGCTGGATCT‑GCCGGTGGCAAAAttc                                                                                                           <  2:126273/251‑2 (MQ=255)                                                                                                                            gCACGCCTGGCTGGCAGGGATTGTCGGTTGTCTGCTAGGCG‑ACTGGATCTCCTTCTGGCTGGGCTGGCGTTTTAAAAAGCCGTTGCATCGCTGGTCATTTCTG‑AAGAAAAACAAAGCATTACTCGATAAAACTGAACATGCGTTGCATCAACACAGCATGTTCACCATTCTGGTCGGTCGCTTTGTTGGCCCGACGCGTCCGCTGGTGCCAATGGTGGCGGGAATGCTGGATCT‑GCCGGTGGCAAAAttc                                                                                                           >  1:267865/1‑249 (MQ=255)                                                                                                                                       tggcagggattgtcggttgtctgctGGGCG‑ACTGGATCTCCTTCTGGCTGGGCTGGCGTTTTAAAAAGCCGTTGCATCGCTGGTCATTTCTG‑AAGAAAAACAAAGCATTACTCGATAAAACTGAACATGCGTTGCATCAACACAGCATGTTCACCATTCTGGTCGGTCGCTTTGTTGGCCCGACGCGTCCGCTGGTGCCAATGGTGGCGGGAATGCTGGGTCT‑TCCGGTGGCAAAATTCATTACGCCGAat                                                                                               >  1:177059/25‑251 (MQ=255)                                                                                                                                                 ggtcggttgtctgctGGGCG‑ACTGGATCTCCTTCTGGCTGGGCTGGCGTTTTAAAAAGCCGTTGCATCGCTGGTCATTTCTG‑AAGAAAAACAAAGCATTACGCGATAAAACTGAACATGCGTTGCATCAACACAGCATGTTCACCATTCGGGTCGGTCGCTTTGTTGGCCCGACGCGTCCGCTGGTGCCAATGGTGGCGGGAATGCTGGATCT‑GCCGGTGGCAAAATTCATTACGCCGAATATTATCGGc                                                                                      <  1:319090/236‑1 (MQ=255)                                                                                                                                                                gggCG‑ACTGGATCTCCTTCTGGCTGGGCTGGCGTTTTAAAAAGCCGTTGCATCGCTGGTCATTTCTG‑AAGAAAAACAAAGCATTACTCGATAAAACTGAACATGCGTTGCATCAACACAGCATGTTCACCATTCTGGTCGGTCGCTTTGTTGGCCCGACGCGTCCGCTGGTGCCAATGGTGGCGGGAATGCTGGATCT‑GCCGGTGGCAAAATTCATTACGCCGAATATTATCGGCTGCCTGCTGTGgccac                                                                      >  1:248788/1‑249 (MQ=255)                                                                                                                                                                 ggCG‑ACTGGATCTCCTTCTGGCTGGGCTGGCGTTTTAAAAAGCCGTTGCATCGCTGGTCATTTCTG‑AAGAAAAACAAAGCATTACTCGATAAAACTGAACATGCGTTGCATCAACACAGCATGTTCACCATTCTGGTCGGTCGCTTTGTTGGCCCGACGCGTCCGCTGGTGCCAATGGTGGCGGGAATGCTGGATCT‑GCCGGTGGCAAAATTCATTACGCCGAATATTATCGGCTGCCTGCTGTGgccacc                                                                     <  1:432338/251‑4 (MQ=255)                                                                                                                                                                  ggcgactggttctccgtatggctggGCTGGCGTTTTAAAAAGCCGTTGCATCGGTGGTCATTGCGG‑AAGAAAAACAAAGCATTACTCGATAAATCTGAACATGCGTTGCATCAACACAGCATGGTCACCATTCTGGTCGGTCGCTTGGTTGGCCCGACGCGTCCGCTGGTGCCAATGGTGGCGGGAATGCTGGATCT‑GCCGGTGGCAAAATTCATTACGCCGAATATTATCGGCTGCCTGCTGTGgcaacc                                                                     <  2:174077/233‑5 (MQ=255)                                                                                                                                                                     gactgaaactccttctGGATGGGCTTGCGTTTTAAAAAGCCGCTGCATCGCAGGTCATTTCTG‑AAGAAAATCAAAGCATTACTAGATAAAACTGAACATGCGTTGCATCAACACAGCATGTTCACTATTCTGGTCGGTCTCTTTGTTGGCACGACGCGTACGCTGGTGCCAATGGTGGCGGGGATGCTGGATCTTGCCGG‑GGCGAAATTCATTACCCAGAATATTTTCGGCGGCctgctgctgccaccctt                                                                  >  2:307781/13‑240 (MQ=255)                                                                                                                                                                         ggatctccttctGGCTGGTGTGGCGTTTTAAAAAGCCGTTGCATCGCTGGTCATTCCTG‑AAGAAAAACAAAGCATTACTCGATAAAACTGAACATGCGTTGAATCAACACAGCATGTTCACCATCCTGGTCGGTCGCTTTGTTGGCCCGATGCGGCCGCTGGTGCCAATGGTGGTGGGAATGCTGGATCT‑TCCGGTTGCAAAATTCATTACGCCGAATATTATCGGCTGCCTGCTGTGGCCACCGTTTTa                                                               <  1:408613/242‑1 (MQ=255)                                                                                                                                                                         ggatctccttctGGCTGGGCTGGCGTTTTAAAAAGCCGTTGCATCGCTGGTCATTTCTG‑AAGAAAAACAAAGCATTACTCGATAAAACTGAACATGCGTTGCATCAACACAGCATGTTCACCATTCTGGTCGGTCGCTTTGTTGGCCCGACGCGTCCGCTGGTGCCAATGGTGGCGGGAATGCTGGATCT‑GCCGGTGGCAAAATTCATTACGCCGAATATTATCGGCTGCCTGCTGTGGCCACCGtttt                                                                <  1:384339/241‑1 (MQ=255)                                                                                                                                                                                           gctGGCGTTTTAAAAAGCCGTTGCATCGCTGGTCATTTCTG‑AAGAAAAACAAAGCATTGCTCGATAAAACTGAACATGCGTTGCATCAACACAGCATGTTCACCATTCTGGTCGGTCGCTTTGTTGGCCCGACGCGTCCGCTGGTGCCAATGGTGGCGGGAATGCTGGATCT‑GCCGGTGGCAAAATTCATTACGCCGAATATTATCg                                                                                        <  1:31735/205‑1 (MQ=255)                                                                                                                                                                                             tGGCGTTTTAAAAAGCCGATGCATCGCTGGTCATTTCTG‑AAGAAAAACAAAGCATTACTCGATAAAACTGAACATGCGTTGCATCAACACAGCATGTTCACCATTCTGGTCGGTCGCTTTGTTGGCCCGACGAGTCCGCTGGTGCCAATGGTGGCGGGGATGGTGGATCT‑GCCGGTGGCAAAATTCAATACGCCGTCTGTTTTCGGCTGTCTGCTGTGGGCACTGTTTTACTTGCTgactggggattgggg                                          >  2:306171/1‑237 (MQ=255)                                                                                                                                                                                                 gTTTTAAAAAGCCGTTGCATCGCTGGTCATTTCTG‑AAGAAAAACAAAGCATTACTCGATAAAACTGAACATGCGTTGCATCAACACAGCATGTTCACCATTCTGGTCGGTCGCTTTGTTGGCCCGACGCGTCCGCTGGTGCCAATGGTGGCGGGAATGCTGGATCT‑GCCGGTGGCAAAATTCATTACGCCGAATATTATCGGCTGCCTGCTGTGGCCACCGTTTTACTTCCTGCCTGGGATTCTggcgggg                                      >  1:543892/1‑250 (MQ=255)                                                                                                                                                                                                   tttAAAAAGCCGTTGCATCGATGGTCATTTCTG‑AAGAAAAACAAAACATTACTCGATAAAACTGAACATGCGTTGCATCAACACAGCATGTTCACCATTCTGGTCGGTCGCTTTGTTGGCCCGACGCGTCCGCTGGTGCCAATGGTGGCGGGAATGCTGGATCT‑TCCGGTGGCAAAATTCATTACGCCGAATATTAACGGCTGCCTGCTGTGGACACCGTTTTTCTTCCTGCCTGGGCATCTggcgggggg                                    >  2:396293/1‑248 (MQ=255)                                                                                                                                                                                                    ttAAAAAGCCGTTGCATCGCTGGTCATTTCTG‑AAGAAAAACAAAGCATTACTCGATAAAACTGAACATGCGTTGCATCAACACAGCATGTTCACCATTCTGGTCGGTCGCTTTGTTGGCCCGACGCGTCCGCTGGTGCCAATGGTGGCGGGAATGCTGGATCT‑GCCGGTGGCAAAATTCATTACGCCGAATATTATCGGCTGCCTGCTGTGGCCACCGTTTTACTTCCTGCCTGGGGTTCTGGCGGGcgcg                                   >  1:326382/1‑251 (MQ=255)                                                                                                                                                                                                     tAAAAAGCCGTTGCATCGCTGGTCATTTCTG‑AAGAAAAACAAAGCATTACTCGATAAAACTGAACATGCGTTGCATCAACACAGCATGTTCACCATTCTGGTCGGTCGCTTTGTTGGCCCGACGCGTCCGCTGGTGCCAATGGTGGCGGGAATGCTGGATCT‑GCCGGTGGCAAAATTCATTACGCCGAATATTATCGGCTGCCTGCTGTGGCCACCGTTTTACTTCCTGCCTGGGATTCTGGCGGGcgcg                                   <  1:396293/250‑1 (MQ=255)                                                                                                                                                                                                     tAAAAAGCCGTTGCATCGCTGGTCATTTCTG‑AAGAAAAACAAAGCATTACTCGATAAAACTGAACATGCGTTGCATCAACACAGCATGTTCACCATTCTGGTCGGTCGCTTTGTTGGCCCGACGCGTCCGCTGGTGCCAATGGTGGCGGGAATGCTGGATCT‑GACGGGGGCAAAATTCATTACGCCGAATATTATCGGCTGCCTGCTGTGGCCACCGTTTTACTTCCTGACTGGGATTCTGGCGGGCgcgg                                  >  1:396221/1‑251 (MQ=255)                                                                                                                                                                                                                   aTCGCTGGTCATTTCTG‑AAGAAAAACAAAGCATTACTCGATAAAACTGAACATGCGTTGCATCAACACAGCATGTTCACCATTCTGGTCGGTCGCTTTGTTGGCCCGACGCGTCCGCTGGTGCCAATGGTGGCGGGAATGCTGGATCT‑GCCGGTGGCAAAATTCATTACGCCGAATATTATCGGCTGCCTGCTGTGGCCACCGTTTTACTTCCTGCCTGGGATTCTGGCGGGCGCGGCGGTCGATATTcc                     >  1:321187/1‑250 (MQ=255)                                                                                                                                                                                                                       cTGGTCATTTCTG‑AAGAAAAACAAAGCATTACTCGATAAAACTGAACATGCGTTGCATCAACACAGCATGTTCACCATTCTGGTCGGTCGCTTTGTTGGCCCGACGCGTCCGCTGGTGCCAATGGTGGCGGGAATGCTGGATCT‑GCCGGTGGCAAAATTCATTACGCCGAATATTATCGGctgcct                                                                                 <  1:52392/186‑1 (MQ=255)                                                                                                                                                                                                                       cTGGTCATTTCTG‑AAGAAAAACAAAGCATTACTCGATAAAACTGAACATGCGTTGCATCAACACAGCATGTTCACCATTCTGGTCGGTCGCTTTGTTGGCACGACGCGTCCGCTGGTGCCAATGGTGGCGGGAATTCTGGATCT‑GCCGGTGGCAAAATTCATTACACCGAATATTATCGGctgcct                                                                                 >  2:52392/1‑186 (MQ=255)                                                                                                                                                                                                                            aagttCTG‑AAGAAAACCAAAGCATTACTCAGTAAATCGGAACAGCCGTGTCATCAACATAGCATGTCCACCAGTCTGTGCGGTCGTTTGGTTGTGCCGAGGCGTCCGCTGGTGCCAATGGTGGCGGGAATGCTGGATCG‑GCCGGTGGCAAAATTCCTGACGCCGACTATTATCGCCTGCCTGCTGTGGCCACCGGTTTACTTCCTGCCTTGGATCCGGGCGGGCGCGTCGATCGATATTCCAGCCGGaag            <  2:396719/247‑2 (MQ=255)                                                                                                                                                                                                                             aTTTCTG‑AAGAAAAACAAAGCATTACTCGATAAAACTGAACATGCGTTGCATCAACACAGCATGTTCACCATTCTGGTCGGTCGCTTTGTTGGCCCGACGCGTCCGCTGGTGCCAATGGTGGCGGGAATGCTGGATCT‑GCCGGTGGCAAAATTCATTACGCCGAATATTATCGGCTGCCTGCTGTGGCCACCGTTTTACTTCCTGCCTGGGATTCTGGCGGGCGCGGCGGTCGATATTCCAGCCGGAATGc          >  1:684186/1‑251 (MQ=255)                                                                                                                                                                                                                                 cTG‑AAGAAAAACAAAGCATTACTCGATAAACCTGAACATGCGTTGCATCAACACAGCATGTTCACCATTCTGGTCGGTCGCTTTGTTGGCCCGACGCGTCCGCTGGTGCCAATGGTGGCGGGACTGCTGGATCT‑GCCGGTGGCAAAATTCATTACGCCGAATCTTATCGGCTGCCTGCTGTGGCAACCGTTTTAATTCCTGCCTGGGATTTGGGCGGGGGCGGCGAGCCATATTACCGCCGGGATGCag        >  2:153498/1‑249 (MQ=255)                                                                                                                                                                                                                                  tG‑AAGAAAAACAAAGCATTACTCGATAAAACTGAACATGCGTTGCATCAACACAGCATGTTCACCATTCTGGTCGGTCGCTTTGTTGGCCCGACGCGTCCGCTGGTGCCAATGGTGGCGGGAATGCTGGATCT‑GCCGGTGGCAAAATTCATTACGCCGAATATTATCGGCTGCCTGCTGTGGCCACCGTTTTACTTCCTGCCTGGGATTCTGGCGGGCGCGGCGATCGATATTCCCGCCGGAATGCAGAGc     >  1:13277/1‑251 (MQ=255)                                                                                                                                                                                                                                   g‑aagaaAAACAAAGCATTACTCGATAACACTGAACATGCGTTGCATCAACACAGCATGTTCACCATTCTGGTCGGTCGCTTTGTTGGCCCGACGCGTCCGCTGGTGCCAATGGTGGCGGGAATGCTGGATCT‑GCCGGGGGCAAAACTCATTACGCCGAATATTATCGGCTGCCTGCTGTGGCCACCGATTTACTTCCTGCCTGGGATTCTGGCGGGGGCGGCGAGCGGTATTCCAGCCGGAATGCAGAGCg    >  1:276772/1‑251 (MQ=255)                                                                                                                                                                                                                                    aaagaaAAACAAAGGATTACGCTATAAAACTGAACATGCGTTGCATCAACACAGCATGTTCCCCATTCGGGTCGGTCGCTTTGTTGGCCCGAGGCGTCCGCTGGTGCCAATGGTGGCGGGAATGCTGGATCT‑GCCGGTGGCAAAATTCATTACGCCGAATATTATCGGCTGCCTGCTGTGGCCACCGTTTTACTTCCTGCCTGGGATTCTGGCGGGCGCGGCGATCGATAGTCCAGCCGGAATGCAGAGCg    <  2:686909/250‑1 (MQ=255)                                                                                                                                                                                                                                      agcaaaaCAAAGCATTACTCGATAAAACTGAACATGCGTTGCATCAACACAGCATGTTCACCATTCTGGTCGGTCGCTTTGTTGGCCCGACGCGTCCGCTGGTGCCAATGGTGGCGGGGATGCTGGATCT‑GCCGGTGGCAAAATTCATTAAGCCGAATATTATCGGCTGCCTGCTGTGGCCACCGTTTTACTTCCTGCCTGGGATTCTGGCGGGGGCGGCGAGCGATATTcaatccgggatgaagagcggg  >  2:120214/4‑231 (MQ=255)                                                                                                                                                                                                                                                         |                                                                                                                                                                                                                                         ATTGATGGCGGTGGTGTTGGTGGCCTTTCTGGAGTCGCTGGCGCTGGTCGGTTTGATTCTACCCGGTACGGTGCTGATGGCGGGGCTGGGAGCGCTGATTGGCAGCGGCGAGTTAAGTTTCTGGCACGCCTGGCTGGCAGGGATTATTGGCTGCTTGATGGGCG‑ACTGGATTTCTTTCTGGCTGGGTTGGCGTTTTAAAAAGCCGTTGCATCGCTGGTCATTTCTG‑AAGAAAAACAAAGCACTACTTGATAAAACTGAACATGCGTTGCATCAACACAGCATGTTCACCATTCTGGTCGGTCGTTTTGTTGGCCCGACGCGTCCGCTGGTGCCAATGGTGGCGGGAATGCTGGATCT‑GCCGGTGGCTAAATTTATTACGCCGAATATTATCGGCTGCCTGCTGTGGCCGCCGTTTTACTTCCTGCCAGGGATTCTGGCGGGCGCGGCGATCGATATTCCTGCCGGAATGCAGAGCGGT  >  NC_000913/71398‑71875

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 26 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑