Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,069,024	C→T	G272G (GGG→GGA)	rutG ←	pyrimidine:H(+) symporter

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,069,024	0	C	T	98.3%	174.7 / NA	58	G272G (GGG→GGA)	rutG	pyrimidine:H(+) symporter
Reads supporting (aligned to +/- strand):  ref base C (0/1);  new base T (35/22);  total (35/23)
Fisher's exact test for biased strand distribution p-value = 3.97e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ACCAATAACCGCGGCCGGAATGGTATGAATCAGCGCACCAAATTTCGGTGAAAAGCCCAACAGCATGGCGATGACGGCAGCAG‑CAACAAACACCAGCGTCGAGTAGACTTTGGTCACGGCCATCACACCGATATTTTCAGCATAGGTGG‑TCACGCCGCTACCGCCGACAGAGCCGGAAAGCATCGTTGCCAGACCATCGCCTACGAATGCCCGCCCCATATACGGG‑TCCATATTGCGTCCGGTCATCCCGGCGACTGCCTTGAGATGACCTAAGTTTTCCGCCACCAGAATCACCGCCACGGGCGCAATCAGCATCATTGCCTGACCATTAAAAGCAGGAGTGGAAAAATGTGGCAGACCGAACCAGGCAGCATGGCTGACGAGAGTAAAATCGACGGCTTTTCCCAGCCCTAAAACGTTGGTCATCACGCCATACAGCAGACAGGCGACAATTAATCCTACGAGAATCAATAACCGCTGGATCATGCCAC  >  NC_000913/1068778‑1069268                                                                                                                                                                                                                                                          |                                                                                                                                                                                                                                                     ccaaaTAACCGCGGCCGGAATGGTAGGAATCAGCGCCCCAAATTTCGGTGAAAAGCCAAACAGCAGGGCGATGACGGCAGCAG‑CAACAAACCCCAGCGTCGAGTAGACTTTGGTCACGGCCATCACACCGATATTTGCAGCATAGGTGG‑TCACGCCGCTACCGCCGACAGAGCCGGAAAGCATCGTTGCCAGACCATCGCCTACGAATGCCCGCCCCATATACGGG‑TCCATATTGCGTCCGGTCATTCCgg                                                                                                                                                                                                                                                  <  1:391319/248‑1 (MQ=255)   aaaTAACCGCGGCGGAAATGTAATTAACCAGCGCACAAAATTCCGGTGAAAAGCCCAACAGCATGGCGATGAGGGCAACAA‑CAACAAACCCCAGCGTCGAGTAGACTTTGGTCACGGCCATCACACCGATATTTTCAGCATAGGTGT‑TCACGCCGCTACCGCCGACAGAGCCGGAAAGCATCGTTGCCAGACCATCGCCTACGAATGCCCGCCCCATATACGGG‑TCCATATTGCGTCCGGTCATTCCGGCg                                                                                                                                                                                                                                                <  2:354608/250‑1 (MQ=255)     aTAACCGCGGCCGGAATGGAATGAATCAGCGCACCAAATTTCGATGAAAAGCCCAACAGCATGGCGATGACGGCAGCAG‑CAACAAACACCAGCGTCGAGTAGACTTTGGTCACGGCCATCACACCGATATTTTCAGCATAGGTGG‑TCACGCCGCTACCGCCGACAGAGCCGGAAAGCATCGTTGCCAGACCATCGCCTACGACTGCCCGCCCCATATACGGG‑TCCATATTGCGTCCGGTCATTCCGGCGAc                                                                                                                                                                                                                                              >  2:388736/1‑251 (MQ=255)                cGGAATGGTATAAATCAGCGCACAAAATTCCGGTAAAAAGCCCAACAGCATGGCAAGGACGGCAGCAG‑CAACAAACACCGGCTTCGAGTAGACTTTGGTCACGGCCATCACACCGATATTTTCAGCATGGGTTG‑TCACGCCGCTACCGCCGACAGAGCCGGAAAGCATCGTTGCCAGACCATCGCCTACGAATGCCCGCCCCATATACGGG‑TCCATATTGCGTCCGGTCATTCCGGCGACTGCCTTGAGAt                                                                                                                                                                                                                                   <  2:620955/251‑1 (MQ=255)                    aTGGTATGAATCAGCGCACCAAATTTCGGTGAAAAGCCCAACAGCATGGCGATGACGGCAGCAG‑CAACAAACACCAGCGTCGAGTAGACTTTGGTCACGGCCATCACACCGATATTTTCAGCATAGGTGG‑TCACGCCGCTACCGCCGACAGAGCCGGAAAGCATCGTTGCCAGACCATCGCCTACGAATGCCCGCCCCATATACGGG‑TCCATATTGCGTCCGGTCATTCCGGCGACTGCCTTGAGATGAcc                                                                                                                                                                                                                               >  1:87553/1‑251 (MQ=255)                           gAATCAGCGCACCAAATTTCGGTGAAAAGCCCAACAGCATGGCGATGACGGCAGCAG‑CAACAAACACCAGCGTCGAGTAGACTTTGGTCACGGCCATCACACCGATATTTTCAGCATAGGTGG‑TCACGCCGCTACCGCCGACAGAGCCGGAAAGCATCGTTGCCAGACCATCGCCTACGAATGCCCGCCCCATATACGGG‑TCCATATTGCGTCCGGTCATTCCGGCGACTGCCTTGAGATGACCTAAGttt                                                                                                                                                                                                                        >  1:650677/1‑251 (MQ=255)                            aaacAGCGCACCAAATTTCGGTGAAAAGCCCAACAGCAGGGCGATGACGGCAGCAG‑CAACAAACCCCAGCGTCGAGTAGACTTTGGTCACGGCCATCACACCGATATTTTCAGCATAGGTGG‑TCACGCCGCTACCGCCGACAGAGCCGGAAAGCATCGTTGCCAGACCATCGCCTACGAATGCCCGCCCCATATACGGG‑TCCATATTGCGTCCGGTCATTCCGGc                                                                                                                                                                                                                                                 <  2:180162/222‑1 (MQ=255)                            aaTCAGCGCACCAAATTTCGGTGAAAAGCCCAACAGCATGGCGATGACGGCAGCAG‑CAACAAACACCAGCGTCGAGTAGACTTTGGTCACGGCCATCACACCGATATTTTCAGCATAGGTGG‑TCACGCCGCTACCGCCGACAGAGCCGGAAAGCATCGTTGCCAGACCATCGCCTACGAATGCCCGCCCCATATACGGG‑TCCATATTGCGTCCGGTCATTCCGGc                                                                                                                                                                                                                                                 >  1:180162/1‑225 (MQ=255)                                         aaTTTCGGTGAAAAGCCCAACAGCATGGCGATGACGGCAGCAG‑CAACAAACACCAGCGTCGAGTAGACTTTGGTCACGGCCATCACACCGATATTTTCAGCATAGGTGG‑TCACGCCGCTACCGCCGACAGAGCCGGAAAGCATCGTTGCCAGACCATCGCCTACGAATGCCCGCCCCATATACGGG‑TCCATATTGCGTCCGGTCATTCCGGCGACTGCCTTGAGATGACCTAAGTTTTCCGCCACCAGAAt                                                                                                                                                                                                          >  2:414251/1‑251 (MQ=255)                                           tttCGGTGAAAAGCCCAACAGCATGGCGATGACGGCAGCAG‑CAACAAACACCAGCGTCGAGTAGACTTTGGTCACGGCCATCACACCGATATTTTCAGCATAGGTGG‑TCACGCCGCTACCGCCGACAGAGCCGGAAAGCATCGTTGCCAGACCATCGCCTACGAATGCCCGCCCCATATACGGG‑TCCATATTGCGTCCGGTCATTCCGGCGACTGCCTTTAGATGACCTAAGGTTTCCGCCACCAGAATc                                                                                                                                                                                                         >  1:636477/1‑250 (MQ=255)                                                 gaaaaaGCCCAACAGCGGGGCGATGCCGGCAGCAAACAA‑AAACCCCAGCGTGAAGTAGACTTGGGCCCCGCCCACCCCCCCGATATTTCCAGCAAGGGTGG‑TCACGCCGCTCCCGCCGACAGAGCCGGAAAGCATCGTTGCCAGACCATCGCCTACGAATGCCCGCCCCATATACGGG‑TCCATATTGCGTCCGGTCATTCCGGCGACTGCCTTGAGATGACCTAAGTTTTCCGCCACCAGAATCACCGCCa                                                                                                                                                                                                  <  2:446419/249‑1 (MQ=255)                                                   aaaaGCCCAACAGCATGGCGATGACGGCAGCAG‑CAACAAACACCAGCGTCGAGTAGCCTTTGGTCACGGCCATCACACCGATATTTTCAGCATAGGTGG‑TCACGCCGCTACCGCCGACAGAGCCGGAAAGCATCGTTGCCAGACCATCGCCTACGAATGCCCGCCCCATATACGGG‑TCCATATTGCGTCCGGTCATTCCGGCGACTGCCTTGAGATGACCTAACTTTTCCGCCACCAGAATCACCGCCACg                                                                                                                                                                                                >  2:362074/1‑251 (MQ=255)                                                        cccAACAGCATGGCGATGAAGGCAGCAG‑CAACAAACACCAGCGTCGAGTAGACTTTGGTCACGGCCATCACACCGATATTTTCAGCATAGGTGG‑TCACGCCGCTACCGCCGACAGAGCCGGAAAGCATCGTTGCCAGACCATCGCCTACGAATGCCCGCCCCATATACGGG‑TCCATATTGCGTCCGGTCATTCCGGCGACTGCCTTGAGATGACCTAAGTTTTCCGCCACCAGAATCACCGCCACGGgcgc                                                                                                                                                                                           >  1:577223/1‑251 (MQ=255)                                                              aGCATGGCGATGACGGCAGCAG‑CAACAAACACCAGCGTCGAGTAGACTTTGGTCACGGCCATCACACCGATATTTTCAGCATAGGTGG‑TCACGCCGCTACCGCCGACAGAGCCGGAAAGCATCGTTGCCAGACCATCGCCTACGAATGCCCGCCCCATATACGGG‑TCCATATTGCGTCCGGTCATTCCGGCGACTGCCTTGAGATGACCTAAGTTTTCCGCCACCAGAATCACCGCCACGGGCGCAATc                                                                                                                                                                                       >  2:126239/1‑249 (MQ=255)                                                                  tGGCGATGACGGCAGCAG‑CAACAAACACCAGCGTCGAGTAGACTTTGGTCACGGCCATCACACCGATATTTTCAGCATAGGTGG‑TCACGCCGCTACCGCCGACAGAGCCGGAAAGCATCGTTGCCAGACCATCGCCTACGAATGCCCGCCCCATATACGGG‑TCCATATTGCGTCCGGTCATTCCGGCGACTGCCTTGAGATGACCTAAGTTTTCCGCCACCAGAATCACCGCCACGGGCGCAATCAGcatc                                                                                                                                                                                 >  1:281527/1‑251 (MQ=255)                                                                   ggCGATGAGGGCAGCAA‑CAACAACCCCAAGCGTCGAGTAGACTTTGGTCACGGCCATACCACCGATATTTTCAGAATGGGTGG‑TCACGCCGCTCCCGCGGACAGAGCCGGAAAGCATCGTTCCCAGACCATCGCCTACGAATGCCCGCCCCATATCCGGG‑TCCATATTGCGTCCGGTCATTCCGGCGACTGCCTTGAGATGACCTAAGTTTTCCGCCCCCAGAATCACCGCCACGGGCGCAATCAGcatca                                                                                                                                                                                <  1:165177/251‑1 (MQ=255)                                                                       aTGAAGGCAGCAG‑CAACAAACCCCAGCGTCGAGTAGACTTTGGTCACGGACATCACACCGATATTTTCAGCATAGGCGG‑TCACGCCGCTACCGCCGACAGAGACGGAAAGCATCGTTGCCAGACCATCGCCTACGAATGCCCGCCCCATATACGGG‑TCCATATTGCGTCAGGTCATTCCGGCGCATGCCTTGAGATGACCTAAGTTTTCCGCCACCAGACTCAACGCCACGGGCGCAATCAGCATCATTGc                                                                                                                                                                            >  2:364384/1‑251 (MQ=255)                                                                            ggCAGCAG‑CAACAAACACCAGCGTCGAGTAGACTTTGGTCACGGCCATCACACCGATATTTTCAGCATAGGTGG‑TCACGCCGCTACCGCCGACAGAGCCGGAAAGCATCGTTGCCAGACCATCGCCTACGCATGCCCGCCCCATATACGGG‑TCCATATTGCGTCCGGTCATTCCGGCGACTGCCTTGAGATGACCTAAGTTTTCCGCCACCAGAATCACCGCCACGGGCGCAATCAGCATCATTGCCTGAc                                                                                                                                                                       >  1:91403/1‑251 (MQ=255)                                                                               agcag‑caACAAACACCAGCGTCGAGTAGACTTTGGTCACGGCCATCACACCGATATTTTCAGCATAGGTGG‑TCACGCCGCTACCGCCGACAGAGCCGGAAAGCATCGTTGCCAGACCATCGCCTACGAATGCCCGCCCCATATACGGG‑TCCATATTGCGTCCGGTCATTCCGGCGACTGCCTTGAGATGACCTAAGTTTTCCGCCACCAGAATCACCGCCACGGGCGCAATCAGCATCATTGCCTGACCAt                                                                                                                                                                    >  1:239135/1‑251 (MQ=255)                                                                               agcag‑caACAAACACCAGCGTCGAGTAGACTTTGGTCACGGCCATCACACCGATATTTTCAGCATAGGTGG‑TCACGCCGCTACCGCCGACAGAGCCGGAAAGCATCGTTGCCAGACCATCGCCTACGAATGCCCGCCCCATATACGGG‑TCCATATTGCGTCCGGTCATTCCGGCGACTGCCTTGAGATGACCTAAGTTTTCCGCCACCAGAATCACCGCCACGGGCGCAATCAGCATCATTGCCTGACCAt                                                                                                                                                                    >  2:442978/1‑251 (MQ=255)                                                                                         aaaccccggCGTCGAGTAGCTTTTGGCCCCGCCCATCACCCCGATATTTTCAGCAT‑GGTGGGTCACGCCGCTACCGCCGACAGAGCCGGAAAGCATCGTTGCCAGCCCATCGCCTACGAATGCCCGCCCCATATACGGG‑TCCATATTGCGTCCGGTCATTCCGGCGACTGCCTTGAGATGACCTAAGTTTTCCGCCACCAGAATCACCGCCACGGGCGCAATCAGCATCATTGCCTGACCATTAAAAGCAg                                                                                                                                                           <  1:426284/243‑1 (MQ=255)                                                                                         aaacaCCAGCGTCGAGTAGACTTTGGTCACGGCCATCACACCGATATTTTCAGCATAGGTGG‑TCACGCCGCTACCGCCGACAGAGCCGGAAAGCATCGTTGCCAGACCATCGCCTACGAATGCCCGCCCCATATACGGG‑TCCATATTGCGTCCGGTCATTCCGGCGACTGCCTTGAGATGACCTAAGTTTTCCGCCACCAGAATCACCGCCACGGGCGCAATCAGCATCATTGCCTGACCATTAAAAGc                                                                                                                                                             >  1:360462/1‑249 (MQ=255)                                                                                           acacCAGCGTCGAGTAGACTTTGGTCACGGCCATCACACCGATATTTTCAGCATAGGTGG‑TCACGCCGCTACCGCCGACAGAGCCGGAAAGCATCGTTGCCAGACCATCGCCTACGAATGCCCGCCCCATATACGGG‑TCCATATTGCGTCCGGTCATTCCGGCGACTGCCTTGAGATGACCTAAGTTTTCCGCCACCAGAATCACCGCCACGGGCGCAATCAGCATCATTGCCTGACCATTAAAAGCAggg                                                                                                                                                         >  1:535360/1‑250 (MQ=255)                                                                                                           gACTTTGGTCACGGCCATCACACCGATATTTTCAGCATAGGTGG‑TCACGCCGCTACCGCCGACAGAGCCGGAAAGCATCGTTGCCAGACCATCGCCTACGAATGCACGCCCCATATACGGG‑TCCATATTGCGTCCGGTCATTCCGGCGACTGCCTTGAGATGACCTAAGTTTTCCGCCACCAGAATCACCGCCACGGGCGCAATCAGCATCATTGCCTGACCATTAAAAGCAGGAGTGGAAAAATGTGGc                                                                                                                                          >  1:422605/1‑250 (MQ=255)                                                                                                                                 ccGATATTTTCAGCATAGGTGG‑TCACGCCGCTACCGCCGACAGAGCCGGAAAGCATCGTTGCCAGACCATCGCCTACGAATGCCCGCCCCATATACGGG‑TCCATATTGCGTCCGGTCATTCCGGCGACTGCCTTGAGATGACCTAAGTTTTCCGCCACCAGAATCACCGCCACGGGCGCAATCAGCATCATTGCCTGACCCTTAAAAGCAGGAGTGGAAAAATGTGGCAGACCGAACCAGGCAGCATGGCt                                                                                                                   >  1:479930/1‑251 (MQ=255)                                                                                                                                              cATAGGTGG‑TCACGCCGCTACCGCCGACAGAGCCGGAAAGCATCGTTGCCAGACCATCGCCTACGAATGCCCGCCCCATATACGGG‑TCCATATTGCGTCCGGTCATTCCGGCGACTGCCTTGAGATGACCTAAGTTTTCCGCCACCAGAATCACCGCCACGGGCGCAATCAGCATCATTGCCTGACCATTAAAAGCAGGAGTGGAAAAATGTGGCAGACCGAACCAGGCAGCATGGCTGACGAGAGTaaa                                                                                                       >  2:288356/1‑250 (MQ=255)                                                                                                                                               aatggggggccccccccctcccgccgACAGAGCCGGAAAGCATCGTTGCCAGCCCACCGCCTACGAATGCCCGCCCCATATCCGGG‑TCCATATTGCGCCCGGTCATTCCGGCGACTGCCTTGAGATGACCTAAGTTTTCCGCCACCAGAATCACCGCCACGGGCGCAATCAGCATCATTGCCTGACCATTAAAAGCAGGAGTGGAAAAATGTGGCAGACCGAACCAGGCAGCATGGCTGACGAGAGTaaaa                                                                                                      <  1:288356/231‑1 (MQ=255)                                                                                                                                                        ccccgccgcTCCCGCCGACAGAGCCGGAAAGCATCGTTGCCAGACCATCGCCTACGAATGCCCGCCCCATATACGGG‑TCCATATTGCGTCCGGTCATTCCGGCGACTGCCTGGAGATGACCTAAGTTTTCCGCCACCAGAACCACCGCCACGGGCGCAATCAGCATCATTGCCTGACCATTAAAAGCAGGAGTGGAAAAATGTGGCAGACCGAACCAGGCAGCATGGCTGACGAGAGTAAAATCGACGGCt                                                                                             <  1:37480/248‑1 (MQ=255)                                                                                                                                                           cgccgcTCCCGCCGACAGAGCCGAAAAGCATCGTGGCAAGACCATCGCCTACGAATGCCCGCCCAATATACGGG‑TCCATATTGCGTCCGGTCATTCCGGCGACTGCCTGGAGATGACCTAAGTTTTCCGCCACCAGAATCACCGCCACGGGCGCAATCAGCATCATTGCCTGACCATTAAAAGCAGGAGTGGAAAAATGTGGCAGACCGAACCAGGCAGCATGGCTGACGAGAGTAAAATCGACGGCtttt                                                                                          <  2:87553/251‑1 (MQ=255)                                                                                                                                                            gccgcTACCGCCGACAGAGCCGGAAAGCATCGTTGCCAGACCATCGCCTACGAATGCCCGCCCCATATACGGG‑TCCATATTGCGTCCGGTCATTCCGGCGACTGCCTTGAGATGACCTAAGTTTTCCGCCACCAGAATCACCGCCACGGGCGCAATCAGCATCATTGCCTGACCATTAAAAGCAGGAGTGGAAAAATGTGGCAGACCGAACCAGGCAGCATGGCTGACGAGAGTa                                                                                                         >  1:394561/1‑235 (MQ=255)                                                                                                                                                            gccgcTACCGCCGACAGAGCCGGAAAGCATCGTTGCCAGACCATCGCCTACGAATGCCCGCCCCATATACGGG‑TCCATATTGCGTCCGGTCATTCCGGCGACTGCCTTGAGATGACCTAAGTTTTCCGCCACCAGAATCACCGCCACGGGCGCAATCAGCATCATTGCCTGACCATTAAAAGCAGGAGTGGAAAAATGTGGCAGACCGAACCAGGCAGCATGGCTGACGAGAGTAAAATCGACGGCTNTTc                                                                                         >  2:586782/1‑251 (MQ=255)                                                                                                                                                                cTACCGCCGACAGAGCCGGAAAGCATCGATGCCAGACCATCGCCTACGAATGCCCGCCCCATATACGGG‑TCCATATTGCGTCCGGTCATTCCGGCGACTGCCTTGAGATGACCTAAGTTTTCCGCCACCAGAATCACCGCCACGGGCGCAATCAGCATCATTGCCTGACCATTAAAAGCAGGAGAGGACAAATGTGGCAGCCCGAACCAGGCAGCATGGCTGACGAGAGTAAAATCGACGGCTTTTCCCa                                                                                      >  2:36419/1‑250 (MQ=255)                                                                                                                                                                        gacggagCTGGAAAGCAGCGTGGCCAGACCAGCGCCTACAAATGCCCGCCCCATATCCGGG‑TCCATATTGCGCCCGGTCATTCCGGCGACTGCCTGGAGAGGACCTAAGTTTTCCGCCACCAGAAGCACCGCCACGGGCGCAATCAGCATCATTGCCTGACCATTAAAAGCAGGAGTGGAAAAATGTGGCAGACCGAACCAGGCAGCATGGCTGACGAGAGTAAAATCGACGGCTTTTCCCAGCCCTaaaa                                                                             <  2:257960/247‑1 (MQ=255)                                                                                                                                                                         aaanagCCGGAACGCATGGTTGCAGGCCCATGCCCTACGAATGCCAGCCCAAT‑TACGGGGTCCAGATTGCGTCCGGTCATTCCGGCGACTGCCTTGAGATGACCTAAGTTTTCCGCCCCAAGAAGCACCGCCACGGGCGCAATCAGCATCATTGCCTGACCATTAAAAGCAGGAGTGGAAAAATGTGGCAGACCGAACCAGGCAGCATGGCTGACGAGAGTAAAATCGTCGGCTTTTCCCAGCCCTAAAAc                                                                            <  2:606687/249‑1 (MQ=255)                                                                                                                                                                              gCCGGAAAGCATCGTTGCCAGACCATCGCCTACGAATGCCCCCCCCATATACGGG‑TCCATATTGCGTCCGGTCATTCCGGCGACTGCCTTGAGATGACCTAAGTTTTCCGCCCCCAGAATCACCGCCACGGGCGCAATCAGCATCATTGCCTGACCATTAAAAGCAGGAGTGGAAAAATGTGGCAGACCGAACCAGGCAGCATGGCTGACGAGAGTAAAATCGACGGCTTTTCCCAGCCCTAAAACGTTgg                                                                       <  1:414251/251‑1 (MQ=255)                                                                                                                                                                                cGGAAAGCATCGTTGCCACCCCAACGCCGACGAATCCCCGCCCCATATACGGG‑TCCATATTGCGTCCGGTCATTCCGGCGACTGCCTTGAGAGGACCTAAGTTTTCCGCCACCAGAATCACCGCCACGGGCGCAATCAGCATCATTGCCTGACCATTAAAAGCAGGAGTGGAAAAATGTGGCAGACCGAACCAGGCAGCATGGCTGACGAGAGTAAAATCGACGGCTTTTCCCAGCCCTAAAACGTTGGtc                                                                     <  2:535360/251‑1 (MQ=255)                                                                                                                                                                                    aaGCATCGTTGCCAGACCACCGCCTACGAACGCCCGCCCCATATACGGG‑TCCATATTGCGTCCGGTCATTCCGGCGACTGCCTTGAGATGACCTAAGTTTTCCGCCACCAGAATCACCGCCACGGGCGCAATCAGCATCATTGCCTGACCATTAAAAGCAGGAGTGGAAAAATGTGGCAGACCGAACCAGGCAGCATGGCTGACGAGAGTAAAATCGACGGCTTTTCCCAGCCCTAAAACGTTGGtcatca                                                                 <  1:523967/251‑1 (MQ=255)                                                                                                                                                                                              gCCAGACCATCGCCTACGAATGCCCGCCCCATATACGGG‑TCCATATTGCGTCCGGTCATTCCGGCGACTGCCTTGAGATGACCTAAGTTTTCCGCCACCAGAATCACCGCCACGGGCGCAATCAGCATCATTGCCTGACCATTAAAAGCAGGAGTGGAAAAATGTGGCAGACCGAACCAGGCAGCATGGCTGACGAGAGTAAAATCGACGGCTTTTCCCAGCCCTAAAACGTTGGTCATCACGCCATAc                                                         >  1:615891/1‑249 (MQ=255)                                                                                                                                                                                               ccAGACCATCGCCTACGAATGCCCGCCCCATATACGGG‑TCCATATTGCGTCCGGTCATTCCGGCGACTGCCTTGAGATGACCTAAGTTTTCCGCCACCAGAATCACCGCCACGGGCGCAATCAGCATCATTGCCTGACCATTAAAAGCAGGAGTGGAAAAATGTGGCAGACCGAACCAGGCAGCATGGCTGACGAGAGTAAAATCGACGGCTTTTCCCAGCCCTAAAACGTTGGTCATCACGCCATAcagc                                                      >  1:318926/1‑251 (MQ=255)                                                                                                                                                                                                 aGACCATCGCCTACGAATGCCCGCCCCATATACGGG‑TCCATATTGCGTCCGGTCATTCCGGCGACTGCCTTGAGATGACCTAAGTTTTCCGCCACCAGAATCACCGCCACGGGCGCAATCAGCATCATTGCCTGACCATTAAAAGCAGGAGTGGAAAAATGTGGCAGACCGAACCAGGCAGCATGGCTGACGAGAGTAAAATCGACGGCTTTTCCCAGCCCTAAAACGTTGGTCATCACGCCATAcagc                                                      >  2:337549/1‑249 (MQ=255)                                                                                                                                                                                                        gcactccaaatccccgcccCAAAGCCGGT‑TCCAGATTGCATCCGGTCATTCCGGCGACGCCCTTGAGATGACCTAGGTTTTCCGCACCCAGAATCGCCGCCACGGGCGCAATCAGCATCATTGCCTGACCATTAAAAGCAGGAGTGGAAAAATGTGGCAGACCGAACCAGGCAGCATGGCTGACGAGAGTAAAATCGACGGCTTTTCCCAGCCCTAAAACGTTGGTCAGCACGCCATACAGCAGACAGGCg                                             <  2:157980/239‑1 (MQ=255)                                                                                                                                                                                                         gCCTACGAATGCCCGCCCCATATACGGG‑TCCATATTGCGTCCGGTCATTCCGGCGACTGCCTTGAGATGACCTAAGTTTTCCGCCACCAGAATCACCGCCACGGGCGCAATCAGCATCATTGCCTGACCATTAAAAGCAGGAGTGGAAAAATGTGGCAGACCGAACCAGGCAGCATGGCTGACGAGAGTAAAATCGACGGCTTTTCCCAGCCCTAAAACGTTGGTCATCACGCCATACAGCAGACAGGCg                                             >  2:41879/1‑250 (MQ=255)                                                                                                                                                                                                         gCCTACGAATGCCCGCACCATATACGGG‑TCCATATTGCGTCCGGTCATTCCGGCGACTGCCTTGAGATGACCTAAGTTTTCCGCCACCAGAATCACCGCCACGGGCGCAATCAGCATCAtt                                                                                                                                                                              >  2:356650/1‑121 (MQ=255)                                                                                                                                                                                                                aagcccccccccATATAGGGG‑TCCTTATTGGGTCCGGTCATTCCGGATATTGCCTTGAGATGACCTTAGTTTTCCGCCCCCAAAATCACCGCCACGGGCGCAATCAGCATCATTGCCTGACCATTAAAAGCAGGAGTGGAAAAATGTGGCAGACCGAACCAGGCAGCATGGCTGACGAGAGTAAAATCGACGGCTTTTCCCAGCCCTAAAACGTTGGTCATCACGCCATACAGCAGACAGGCGACAATTaa                                     <  2:360462/243‑1 (MQ=255)                                                                                                                                                                                                                   gcccgcccCATATACGGG‑TCCATATTGCGTCCGGTCATTCCGGCGACTGCCTTGAGATGACCTAAGTTTTCCGCCACCAGAATCACCGCCACGGGCGCAATCAGCATCATTGCCTGACCATTAAAAGCAGGAGTGGAAAAATGAGGCAGACCGAACCAGGCAGCATGGCTGACGAGAGTAAAATCGACGGCTTTTCCCAGCCCTAAAACGTTGGTCATCACGCCATACAGCAGACAGGCGACAATTAATcc                                  >  2:507157/1‑251 (MQ=255)                                                                                                                                                                                                                   gcccgcccAATATACGGG‑CCCATATTGCGCCCGGTCATTCCGGGGACTGCCTTGAGATGCCCAAATTTTCCCGCCACCAGAATCACCGCCACGGGCGCAATCAGCATCATTGCCTGACCATTAAAAGCAGGAGTGGAAAAATGTGGCAGACCGAACCAGGCAGCATGGCTGACGAGAGTAAAATCGACGGCTTTTCCCAGCCCTAAAACGTTGGTCATCACGCCATACAGCAGACAGGCGACAATTAATcc                                  <  2:173005/251‑1 (MQ=255)                                                                                                                                                                                                                        ccccATATACGGG‑TCCATATTGCGTCCGGTCATTCCGGCGACTGCCTTGAGATGACCTAAGTTTTCCGCCACCAGAATCACCGCCACGGGCGCAATCAGCATCATTGCCTGACCATTAAAAGCAGGAGTGGAAAAATGTGGCAGACCGAACCAGGCAGCATGGCTGACGAGAGTAAAATCGACGGCTTTTCCCAGCCCTAAAACGTTGGTCATAACGCCATACAGCAGACAGGCGACAATTAATCCTACg                              >  2:630400/1‑250 (MQ=255)                                                                                                                                                                                                                         ccgatatACGGG‑TCCATATTGCGTCCGGTCATTCCGGCGACTGCCTTGAGATGACCTAAGTTTTCCGCCACCAGAATCACCGCCACGGTCGCAATCAGCATCATTGCCTGACCATTAAAAGCAGGAGTGGAAAAATGAGGCAGACCGAACCAGGCAGCACGGCTGACGAGCGTAAAATCGACTGCTTTTCCCAGCCCTAAAACGTTGGTCATCACGCCATCCAGCAGACAGGCGACAATTAATCCTACg                              >  2:526822/4‑249 (MQ=255)                                                                                                                                                                                                                            aaataCGGG‑TCATTTTTGCGTCCGGTCATTCCGGCGGCGGCCTGGAGATGACCTAAGTTTGCCGCCACCAGAATCACCGCCACGGGCGCAATCAGCATCATTGCCTGACCAGTAAAAGCAGGAGTGGAAAAATGTGGCAGACCGAACCAGGCAGCATGGCTGACGAGAGTAAAATCGACGGCTTTTCCCAGCCCTAAAACGTTGGTCATCACGCCATACAGCAGACAGGCGACAATTAATCCTACGAGAAt                         <  2:440163/249‑1 (MQ=255)                                                                                                                                                                                                                             tataCGGG‑TCCATATTGCGTCCGGTCATTCCGGCGACTGCCTTGAGATGACCTAAGTTTTCCGCCACCAGAATCACCGCCACGGGCGCAATCAGCATCATTGCCTGACCATTAAAAGCAGGAGTGGAAAAATGTGGCAGACCGAACCAGGCAGCATGGCTGACGAGAGTAAAATCGACGGCTTTTCCCAGCCCTAAAACGTTGGTCATCACGCCATACAGCAGACAGGCGACAATTAATCCTACGAGAATc                        >  2:497989/1‑251 (MQ=255)                                                                                                                                                                                                                             aataCGGG‑TTCATATTGCGTCCGGTCATTCCGGCGACTGCCTTAAGAAGACCTAATTTTACCGCCACCAAAACCACCGCCACGGGCGCAATCAGCATCATTGCCTGACCATTAAAAGCAGGAGTGGAAAAATGTGGCAGACCGAACCAGGCAGCATGGCTGACGAGAGTAAAATCGACGGCTTTTCCCAGCCCTAAAACGTTGGTCATCACGCCATACAGCAGACAGGCGACAATTAATCCTACGAGAATc                        <  2:654584/250‑1 (MQ=255)                                                                                                                                                                                                                                        cATATTCCGTCCGGCCATTCCGGCGACAGCCTTGAGATGACCTAAGTTTTCCGCCACCAGAATCACCGCCACGGGCGCAATCAGCATCATTGCCTGACCATTAAAAGCAGGAGTGGAAAAATGTGGCAGACCGAACCAGGCAGCATGGCTGACGAGAGTAAAATCGACGGCTTTTCCCAGCCCTAAAACGTTGGTCATCACGCCATACAGCAGACAGGCGACAATTAATCCTACGAGAATCAATACCCGCt              >  1:513703/1‑251 (MQ=255)                                                                                                                                                                                                                                           atTGCGTCCGGTCATTCCGGCGACTGCCTTGAGATGACCTAAGTTTTCCGCCACCAGAATCACCGCCACGGGCGCAATCAGCATCATTGCCTGACCATTAAAAGCAGGAGTGGAAAAATGTGGCAGACCGAACCAGGCAGCATGGCTGACGAGAGTAAAATCGACGGCTTTTCCCAGCCCTAAAACGTTGGTCATCACGCCATACAGCAGACAGGCGACAATTAATCCTACGAGAATCAATAACCGCTgg            >  2:515324/1‑250 (MQ=255)                                                                                                                                                                                                                                             tGCGTCCGGTCATCCCGGCGCCTGCCTTGAGAGGCCCTAAGTTTCCCCCCCCCAAAATCACCGCCACGGGCGCAATCAGCATCATTGCCTGACCATTAAAAGCAGGAGTGGAAAAATGTGGCAGACCGAACCAGGCAGCATGGCTGACGAGAGTAAAATCGACGGCTTTTCCCAGCCCTAAAACGTTGGTCATCACGCCATACAGCAGACAGGCGACAATTAATCCTACGAGAATCAATAACCGCTGGAt          <  2:318926/250‑1 (MQ=255)                                                                                                                                                                                                                                                gcccGGTCATTCCGGCGACGGCCTTGAAATGACTTATGTTTTCCGCCACCAGATCCACCGCCACGGGCGCAATCAGCATCATTGCCTGACCATTAAAAGCAGGAGTGGAAAAATGTGGCAGACCGAACCAGGCAGCATGGCTGACGAGAGTAAAATCGACGGCTTTTCCCAGCCCTAAAACGTTGGTCATCACGCCATACAGCAGACAGGCGACAATTAATCCTACGAGAATCAATAACCGCTGGATCATg      <  2:577223/249‑1 (MQ=255)                                                                                                                                                                                                                                                gTCCGGTCATTCCGGCGACTGCCTTGAGATGACCTAAGTTTTCCGCCACCAGAATCACCGCCACGGGCGCAATCAGCATCATTGCCTGACCATTAAAAGCAGGAGTGGAAAAATGTGGCAGACCGAACCAGGCAGCATGGCTGACGAGAGTAAAATCGACGGCTTTTCCCAGCCCTAAAACGTTGGTCATCACGCCATACAGCAGACAGGCGACAATTAATCCTACGAGAATCAATAACCGCTGGATCATg      >  1:167892/1‑251 (MQ=255)                                                                                                                                                                                                                                                    ggTCATTCCGGCTACTGCCTTGAGATGCCCTAAGTTTCCCGCCCCCAAAATCACCGCCACGGGCGCAATCAGCATCATTGCCTGACCATTAAAAGCAGGAGTGGAAAAATGTGGCAGACCGAACCAGGCAGCATGGCTGACGAGAGTAAAATCGACGGCTTTTCCCAGCCCTAAAACGTTGGTCATCACGCCATACAGCAGACAGGCGACAATTAATCCTACGAGAATCAATAACCGCTGGATCATGCCAc  <  1:362074/251‑1 (MQ=255)                                                                                                                                                                                                                                                    ggTCATTCCGGCGACTGCCTTGAGATGACCTAAGTTTTCCGCCACCAGAATCACCGCCACGGGCGCAATCAGCATCATTGCCTGACCATTAAAAGCAGGAGTGGAAAAATGTGGCAGACCGAACCAGGCAGCATGGCTGACGAGAGTAAAATCGACGGCTTTTCCCAGCCCTAAAACGTTGGTCATCACGCCATACAGCAGACAGGCGACAATTAATCCTACGAGAATCAATAACCGCTGGATCATGCCAc  >  1:132520/1‑251 (MQ=255)                                                                                                                                                                                                                                                          |                                                                                                                                                                                                                                                     ACCAATAACCGCGGCCGGAATGGTATGAATCAGCGCACCAAATTTCGGTGAAAAGCCCAACAGCATGGCGATGACGGCAGCAG‑CAACAAACACCAGCGTCGAGTAGACTTTGGTCACGGCCATCACACCGATATTTTCAGCATAGGTGG‑TCACGCCGCTACCGCCGACAGAGCCGGAAAGCATCGTTGCCAGACCATCGCCTACGAATGCCCGCCCCATATACGGG‑TCCATATTGCGTCCGGTCATCCCGGCGACTGCCTTGAGATGACCTAAGTTTTCCGCCACCAGAATCACCGCCACGGGCGCAATCAGCATCATTGCCTGACCATTAAAAGCAGGAGTGGAAAAATGTGGCAGACCGAACCAGGCAGCATGGCTGACGAGAGTAAAATCGACGGCTTTTCCCAGCCCTAAAACGTTGGTCATCACGCCATACAGCAGACAGGCGACAATTAATCCTACGAGAATCAATAACCGCTGGATCATGCCAC  >  NC_000913/1068778‑1069268

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