Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,173,431	T→A	I97N (ATC→AAC)	ygiN →	putative quinol monooxygenase YgiN

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,173,431	0	T	A	100.0%	246.1 / NA	68	I97N (ATC→AAC)	ygiN	putative quinol monooxygenase YgiN
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base A (42/26);  total (42/26)

CATCACCGTCAGGCGGTATTGGATCAGTTTGCTAAAATCGTTCCAACCGTACTGAAAGAAGAAGGTTGCCACGGCTATGCGCCAATGGTGGATTGCGCAGCTGGCGTGAGTTTCCAGTCTATGGCACCGGATTCTATCGTGATGATTGAGCAGTGGGAAAGCATCGCGCATCTTGAAGCGCATCTGCAAACCCCGCACATGAAGGCGTATAGCGAAGCCGTAAAAGGTGACGTGCTGGAGATGAATATCCGTATTCTGCAGCCAGGGATTTAATCCTGCCTTGTTTGCCCGGCCATCCTGACCGGGCAATGTTCTTTCCTTTAAACCTCAATCTCCGCCATGTCGCCTTTCTCTTGCAACCAGTTGCGGCGATCTTCCGAGCGTTTCTTCGCCAGCAGCATATCCATCATCGCGTCAGTACGCTGATCGTCTTCATCATCGATAGTCAACTGCACCAGACGGCGAGTGTTCGGATCAAGCGTGGTTTCGCG  >  NC_000913/3173184‑3173674                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                    aataaacggcaacgcacggtccatgatcacgtacatggagtagtttaagtaggcgttttccgtaaattcatgtagcgcaaggcgctctgccatatcgctcattaattctgattcctcaacttattcgccagccttcaaaggggcaatattgccggagatactacattatctgccgcgtcgagtcacaggagatcaacgtgaatagcgccgccggacgtgctccgacggcgataaaggcatcaGAATAACCg                                                                                                                                                                                                                                                  <  1:457676‑M2/9‑1 (MQ=255)  ataaacggcaacgcacggtccatgatcacgtacatggagtagtttaagtaggcgttttccgtaaattcatgtagcgcaaggcgctctgccatatcgctcattaattctgattcctcaacttattcgccagccttcaaaggggcaatattgccggagatactacattatctgccgcgtcgagtcacaggagatcaacgtgaatagcgccgccggacgtgctccgacggcgataaaggcatcaGAATAACCGt                                                                                                                                                                                                                                                 >  1:240963‑M2/242‑251 (MQ=255)   taaacggcaacgcacggtccatgatcacgtacatggagtagtttaagtaggcgttttccgtaaattcatgtagcgcaaggcgctctgccatatcgctcattaattctgattcctcaacttattcgccagccttcaaaggggcaatattgccggagatactacattatctgccgcgtcgagtcacaggagatcaacgtgaatagcgccgccggacgtgctccgacggcgataaaggcatcaGAATAACCGt                                                                                                                                                                                                                                                 >  1:357498‑M2/241‑250 (MQ=255)      actgactcgcacggcccatgatctcgtacgtggagtagtttaagtaggcgttttccgtaaattcatgtatcgcaaggcgctctgccatatcgctcattaattctgagtcctcaccttattcgccagcattcaaagggtcaataatgccggagatactacattatctgccgtgtcgagtcacaggagatcaacgtgaatagcgccgccggacgtgctccgacggcgataaaggcatcaGAATAACCGTATTc                                                                                                                                                                                                                                             <  2:437406‑M2/14‑1 (MQ=255)          caacgcacggtccatgatcacgtatatggagaagtttaagtaggcgttttccgaaaattcatgtagcgcaaggcgctctgccatatcgctcattaattctgattcctcaacttattcgccagccttcaaaggggcaatattgccggagatactacattatctgccgcgtcgagtcacaggagatcaacgtgaatagcgccgccggacgtgctccgacggcgataaaggcatcaGAAAAACCGGATTCTGc                                                                                                                                                                                                                                          >  2:352946‑M2/234‑250 (MQ=255)                   gtccatgatcacgtacatggagtagtttaagtaggcgttttccgtaaattcatgtagcgcaaggcgctctgccatatcgctcattaattctgattcctcaacttattcgccagccttcaaaggggcaatattgccggagatactacattatctgccgcgtcgagtcacaggagatcaacgtgaatagcgccgccggacgtgctccgacggcgataaaggcatcaGAATAACCGTATTCTGCAGCCAGGGAt                                                                                                                                                                                                                                >  2:112227‑M2/225‑251 (MQ=255)                     ccatgatcacgtacatggagtagtttaagtaggcgttttccgtaaattcatgtagcgcaaggcgctctgccatatcgctcattaattctgattcctcaacttattcgccagccttcaaaggggcaatattgccggagatactacattatctgccgcgtcgagtcacaggagatcaacgtgaatagcgccgccggacgtgctccgccggcgataaaggcatcaGAATAACCGTATTCTGCAGCCAGGGCttt                                                                                                                                                                                                                              >  2:550373‑M2/223‑251 (MQ=255)                     ccatgatcacgtacatggagtagtttaagtaggcgttttccgtaaattcatgtagcgcaaggcgctctgccatatcgctcattaattctgattcctcaacttattcgccagccttcaaaggggcaatattgccggagatactacattatctgccgcgtcgagtcacaggagatcaacgtgaatagcgccgccggacgtgctccgacggcgataaaggcatcaGAATAACCGTATTCTGCAGCCAGGGAttt                                                                                                                                                                                                                              >  1:2485‑M2/223‑251 (MQ=255)                         gatcacgtacatgtagtagtttaagtaggcgttttccgtaaattcatgtagcgcaaggcgctctgccatatcgctcattaattctgattcctcaacttattcgccagccttcaaaggggcaatattgccggagatactacattatctgccgcgtcgagtcacaggagatcaacgtgaatagcgccgccggacgtgctccgacggcgataaaggcatcaGAATAACCGTATTCTGCAGCCAGGGATTTAATc                                                                                                                                                                                                                          <  1:298465‑M2/33‑1 (MQ=255)                                   atggagtagtttaagtaggcgttttccgtaaattcatgtagcgcaaggcgctctgccatatcgctcattaattctgattcctcaacttattcgccagccttcaaaggggcaatattgccggagatactacattatctgccgcgtcgagtcacaggagatcaacgtgaatagcgccgccggacgtgctccgacggcgataaaggcatcaGAATAACCGTATTCTGCAGCCAGGGATTTAATCCTGCCTTGtt                                                                                                                                                                                                                >  1:105818‑M2/209‑251 (MQ=255)                                     ggagtagtttaagtaggcgttttccgtaaattcatgtagcgcaaggcgctctgccatatcgctcattaattctgattcctcaacttattcgccagccttcaaaggggcaatattgccggagatactacattatctgccgcgtcgagtcacaggagatcaacgtgaatagcgccgccggacgtgctccgacggcgataaaggcatcaGAATAACCGTATTCTGCAGCCAGGGATTTAATCCTGCCTTGttt                                                                                                                                                                                                               <  1:289635‑M2/44‑1 (MQ=255)                                      gagtagtttaagtaggcgttttccgtaaattcatgtagcgcaaggcgctctgccatatcgctcattaattctgattcctcaacttattcgccagccttcaaaggggcaatattgccggagatactacattatctgccgcgtcgagtcacaggagatcaacgtgaatagcgccgccggacgtgctccgacggcgataaaggcatcaGAATAACCGTATTCTGCAGCCAGGGATTTAATCCTGCCTTGTTTGc                                                                                                                                                                                                             >  1:50125‑M2/206‑251 (MQ=255)                                      gagtagtttaagtaggcgttttccgtaaattcatgtagcgcaaggcgctctgccatatcgctcattaattctgattcctcaacttattcgccagccttcaaaggggcaatattgccgcagatactacattatctgccgcgtctagtcacaggagatcaacatgaatcgcgccgccggacgtcctccgacggcgataaaggcatcaGAATAACCGTCTTCTTCATCCAGGGATTTACTCCTGCCTTGTTTGc                                                                                                                                                                                                             >  2:306453‑M2/206‑251 (MQ=255)                                          agtttaagtaggcgttttccgtaaattcatgtagcgcaaggcgctctgccatatcgctcattaattctgattcctcaacttattcgccagccttcaaaggggcaatattgccggagatactacattatctgccgcgtcgagtcacaggagatcaacgtgaatagcgccgccggacgtgctccgacggcgataaaggcatcaGAATAACCGTATTCTGCAGCCAGGGATTTAATCCTGCCTTGTTTGCCCg                                                                                                                                                                                                          <  1:154016‑M2/49‑1 (MQ=255)                                                      cgttttccgtaaattcatgtagcgcaaggcgctctgccatatcgctcattaattctgattcctcaacttattcgccagccttcaaaggggcaatattgccggagatactacattatctgccgcgtcgagtcacaggagatcaacgtgaatagcgccgccggacgtgctccgacggcgataaaggcatcaGAATAACCGTATTCTGCAGCCAGGGATTTAATCCTGCCTTGTTTGCCCGGCCATCCTGACCg                                                                                                                                                                                             <  1:119529‑M2/62‑1 (MQ=255)                                                       gttttccgtaaattcatgtagcgcaatgcgatctgccatatcgcttattaattctgattcctcaacttattcgccagccttcaaaggggcaatattgccggagatactacattatctgccgcgtcgagtcacaggagatccacgtgaatagcgccgccggacgtgctcagacgacgataaaggcatcaGACTAACGGTATTCTGCAGCCAGGGCGTTAATCCTACCTTGTTTGCACGGCCCTCCTGACCgg                                                                                                                                                                                            >  2:579161‑M2/189‑251 (MQ=255)                                                       gttttccgtaaaatcatgtagcgcaaggcgctctgccatatcgctcattaattctgattcctcaacttatccgccagccttcaaaggggcaatattgccggagatactacattatctgtcgcgtcgagtcacaggagatcaacgtgaatagcgccgccggacgttctccgacggcgataaaggcatcaGAATAACCGTATTCTGCAGCCAGGGATTTAATCCTGCCTTGTTTGCCCGGCCATCCTGACCgg                                                                                                                                                                                            <  1:436836‑M2/63‑1 (MQ=255)                                                          ttccgtaaattcatgtagcgcaaggcgctctgccatatcgatcattcattctgattcctcaccttattcgccagccttcaaaggggcaatattgccggagatactacattatctgccgcgtcgagtcacaggagatcaacgtgaatagcgccgccggacgtgctccgacggcgataaaggcatcaGAATAACCGTATTCTGCAGCCAGGGATTTAATCCTGCCTTGTTTGCCCGGCCATCCTGACCGGGc                                                                                                                                                                                          >  1:483832‑M2/186‑250 (MQ=255)                                                            ccgaaagttcagttcgcgcaaggcgctctgccatatcgcgcattaattctgattccacaacttattcgccagccttaaaaggggcaatattgccggagatatgtcattatctgcctcgtcgagtcacaggagatcaacgagaatagcgccgccggatgtgctccgacggcgataaaggcatcaGTATAACCGTATGCTTCAGCTAGGGATTTAATCCTGCCTTGTTTGCCCGGCCATCCTGACCTGGCAAt                                                                                                                                                                                       <  2:448251‑M2/68‑1 (MQ=255)                                                             cgtaaattcatgtagcgcaaggcgctctgccatatcgctcattaattctgattcctcaacttatgcgccagcctgcaaaggggcaatattgccggagatactacattatctgccgcgtcgagtcacaggagatcaacgtggagagcgccgccggacgtgctccgacggcgataaaggcatcaGAATAACCGTATGCTGCAGCCAGGGATTTAATCCTGCCTTGTTTGCCCGGCCATCCTGACCGGGCAATg                                                                                                                                                                                      <  1:273578‑M2/69‑1 (MQ=255)                                                                 aattcatgtagcgcaaggcgctctgccatatcgctcattaattctgattcctcaacttattcgccagccttcaaaggggcaatattgccggagatactacattatctgccgcgtcgagtcacaggagatcaacgtgaatagcgccgccggacgtgctccgacggcgataaaggcatcaGAATAACCGTATTCTGCAGCCAGGGATTTAATCCTGCCTTGTTTGCCCGGCCATCCTGACCGGGCAATGTTCt                                                                                                                                                                                  >  1:32499‑M2/179‑251 (MQ=255)                                                                 aattcatgtagcgcaaggcgctctgccatatcgctcattaattctgattcctcaacttattcgccagccttcaaaggggcaatattgccggagatactacattatctgccgcgtcgagtcacaggagatcaacgtgaatagcgccgccggacgtgctccgacggcgataaaggcatcaGAATAACCGTATTCTGCAGCCAGGGATTTAATCCTGCCTTGTTTGCCCGGCCATCCTGACCGGGCAATGTTCt                                                                                                                                                                                  <  1:48269‑M2/73‑1 (MQ=255)                                                                      atgtagcgcaaggcgctctgccatatcgctcattaattctgattcctcaacttattcgccagccttcaaaggggcaatattgccggagatactacattatctgccgcgtcgagtcacaggagatcaacgtgaatagcgccgccggacgtgctccgacggcgataaaggcatcaGAATAACCGTATTCTGCAGc                                                                                                                                                                                                                                       >  2:566466‑M2/174‑193 (MQ=255)                                                                      atgtagcgcaaggcgctctgccatatcgctcattaattctgattcctcaacttattcgccagccttcaaaggggcaatattgccggagatactacattatctgccgcgtcgagtcacaggagatcaacgtgaatagcgccgccggacgtgctccgacggcgataaaggcatcaGAATAACCGTATTCTGCAGc                                                                                                                                                                                                                                       <  1:566466‑M2/20‑1 (MQ=255)                                                                      atgtagcgcaaggcgctctgccatatcgctcattaattctgattcctcaacttattcgccagccttcaaaggggcaatattgccggagatactacattatctgccgcgtcgagtcacaggagatcaacgtgaatagcgccgccggacgtgctccgacggcgataaaggcatcaGAATAACCGTATTCTGCAGCCAGGGATTTAATCCTGCCTTGTTTGCCCGGCCATCCTGACCGGGCAATGTTCTTTCCt                                                                                                                                                                             >  1:420084‑M2/174‑251 (MQ=255)                                                                           gcgcaaggcgctctgccatatcgctcattaattctgattcctcaacttattcgccagccttcaaaggggcaatattgccggagatactacattatctgccgcgtcgagtcacaggagatcaacgtgaatagcgccgccggacgtgctccgacggcgataaaggcatcaGAATAACCGTATTCTGCAGCCAGGGATTTAATCCTGCCTTGTTTGCCCGGCCATCCTGACCGGGCAATGTTCTTTCCTTTaaa                                                                                                                                                                        <  2:106766‑M2/83‑1 (MQ=255)                                                                           gcgcaaggcgctctgccatatcgctcattaattctgattcctcaacttattcgccagccttcaaaggggcaatattgccggagatactacattatctgccgcgtcgagtcacaggagatcaacgtgaatagcgccgccggacgtgctccgacggcgataaaggcatcaGAATAACCGTATTCTGCAGCCAGGGATTTAATCCTGCCTTGTTTGCCCGGCCATCCTGACCGGGCAATGTTCTTTCCTTTaaa                                                                                                                                                                        <  1:73159‑M2/83‑1 (MQ=255)                                                                           gcgcaaggcgctctgccatatcgctcattaattctgattcctcaacttattcgccagccttcaaaggggcaatattgccggagatactacattatctgccgcgtcgagtcacaggagatcaacgtgaatagcgccgccggacgtgctccgacggcgataaaggcatcaGAATAACCGTATTCTGCAGCCAGGGATTTAATCCTGCCTTGTTTGCCCGGCCATCCTGACCGGGCAATGTTCTTTCCTTTaaa                                                                                                                                                                        <  1:475120‑M2/83‑1 (MQ=255)                                                                               aaggcgctctgccatatcgctcattaattctgattcctcaacttattcgccagccttcaaaggggcaatattgccggagatactacattatctgccgcgtcgagtcacaggagatcaacgtgaatagcgccgccggacgtgctccgacggcgataaaggcatcaGAATAACCGTATTCTGCAGCCAGGGATTTAATCCTGCCTTGTTTGCCCGGCCCTCCTGACCGGGCAATGTTCTTTCCTTTAAACCTc                                                                                                                                                                    >  2:574863‑M2/165‑251 (MQ=255)                                                                               aaggccctctgccatatcgctcattaattctntttcctcaacttattcgccagccttcaaaggggcaatattgccggagatgctacattatctgctgcgtcgagtcacaggagatcaacgtgaatagcgccgccggacgtgctccgacggcgataaaggcatcaGAATAACCGTATTCTGCAGCCAGGGATTTAATCCTGCCTTGTTTGCCCGGCCATCCTGACCGGGCAATGTTCTTTCCTTTAaaccac                                                                                                                                                                    <  2:254778‑M2/87‑3 (MQ=255)                                                                                 ggcgctctgccatatcgctcattaattctgattcctcaacttattcgccagccttcaaaggggcaatattgccggagatactacattatctgccgcgtcgagtcacaggagatcaacgtgaatagcgccgccggacgtgctccgacggcgataaaggcatcaGAATAACCGTATTCTGCAGCCAGGGATTTAATCCTGCCTTGTTTGCCCGGCCATCCTGACCGGGCAATGTTCTTTCCTTTAAACCTCa                                                                                                                                                                   >  1:486818‑M2/163‑250 (MQ=255)                                                                                       ctgccatatcgctcattaattctgattcctcaacttattcgccagccttcaaaggggcaatattgccggagatactacattatctgccgcgtcgagtcacaggagatcaacgtgaatagcgccgccggacgtgctccgacggcgataaaggcatcaGAATAACCGTATTCTGCAGCCAGGGATTTAATCCTGCCTTGTTTGCCCGGCCATCCTGACCGGGCAATGTTCTTTCCTTTAAACCTCAAtctc                                                                                                                                                              <  1:422757‑M2/93‑1 (MQ=255)                                                                                            atatcgctcattaattctgattcctcaacttattcgccagccttcaaaggggcaatattgccggagatactacattatctgccgcgtcgagtcacaggagatcaacgtgaatagcgccgccggacgtgctccgacggcgataaaggcatcaGAATAACCGTATTCTGCAGCCAGGGATTTAATCCTGCCTTGTTTGCCCGGCCATCCTGACCGGGCAATGTTCTTTCCTTTAAACCTCAATCTCCTCCATg                                                                                                                                                       <  2:160344‑M2/100‑1 (MQ=255)                                                                                                   tcattaattctgattcctcaacttattcgccagccttcaaaggggcaatattgccggagatactacattatctgccgcgtcgagtcacaggagatcaacgtgaatagcgccgccggacgtgctccgacggcgataaaggcatcaGAATAACCGTATTCTGCAGCCAGGGATTTAATCCTGCCTTGTTTGCCCGGCCATCCTGACCGGGCAATGTTCTTTCCTTTAAACCTCAATCTCCGCCATGTCGCCtt                                                                                                                                                >  1:363789‑M2/145‑251 (MQ=255)                                                                                                     attaattctgattcctcaacttattcgccagccttcaaaggggcaatattgccggagatactacattatctgccgcgtcgagtcacaggagatcaacgtgaatagcgccgccggacgtgctccgacggcgataaaggcatcaGAATAACCGTATTCTGCAGCCAGGGATTTAATCCTGCCTTGTTTGCCCGGCCATCCTGACCGGGCAATGTTCTTTCCTTTAAACCTCAATCTCCGCCATGTCGCCTTtc                                                                                                                                              >  1:479893‑M2/143‑251 (MQ=255)                                                                                                       taattctgattcctcaacttattcgccagccttcaaaggggcaatattaccggagatactacattatctgccgcgtcgagtcacaggagatcaacgtgaatagcgccgccggacgtgctccgacggcgataaaggcatcaGAATAACCGTATTCTGCAGCCAGGGATTTAATCCTGCCTTGTTTGCCCGGCCATCCTGACCGGGCAATGTTCTTTCCTTTAAACCTCAATCTCCGCCTCTTCGCCTTtctc                                                                                                                                            >  2:404103‑M2/141‑251 (MQ=255)                                                                                                              gattcctcaacttattcgccagccttcaaaggggcaatattgccggagatactacattatctgccgcgtcgagtcacaggagatcaacgtgaatagcgccgccggacgtgctccgacggcgataaaggcatcaGAATAACCGTATTCTGCAGCCAGGGATTTAATCCTGCCTTGTTTGCCCGGCCATCCTGACCGGGCAATGTTCTTTCCTTTAAACCTCAATCTCCGCCATGTCGCCTTTCTCTTGCAAc                                                                                                                                     <  1:416571‑M2/118‑1 (MQ=255)                                                                                                                 tcctcaacttattcgccaaccttcaaaggggcaatattgccggagatactacattatctgccgcgtcgagtcacaggagatcaacgtgaatagcgccgccggacgtgctccgacggcgataaaggcatcaGAATAACCGTATTCTGCAGCCAGGGATTTAATCCTGCCTTGTTTGCCCGGCCATCCTGACCGGGCAATGTTCTTTCCTTTAAACCTCAATCTCCGCCATGTCGCCTTTCTCTTGCAAcc                                                                                                                                    >  2:485199‑M2/131‑249 (MQ=255)                                                                                                                         ttattcgccagccttcaaaggggcaatattgccggagatactacattatctgccgcgtcgagtcacaggagatcaacgtgaatagcgccgccggacgtgctccgacggcgataaaggcatcaGAATAACCGTATTCTGCAGCCAGGGATTTAATCCTGCCTTGTTTGCCAGGCCATCCTGACCGGGCAATGTTCTTTCCTTTAAACCTCAATCTCCGCCATGTCGCCTTTCTCTTGCAACCAGTTgcggcg                                                                                                                          >  1:72165‑M2/123‑251 (MQ=255)                                                                                                                          tattcgccagccttcaaaggggcaatattgccggagatactacattatctgccgcgtcgagtcacaggagatcaacgtgaatagcgccgccggacgtgctccgacggcgataaaggcatcaGAATAACCGTATTCTGCAGCCAGGGATTTAATCCTGCCTTGTTTGCCCGGCCATCCTGACCGGGCAATGTTCTTTCCTTTAAACCTCAATCTCCGCCATGTCGCCTTTCTCTTGCAACCAGTTgcggc                                                                                                                           <  1:167694‑M2/128‑1 (MQ=255)                                                                                                                                ccagccttcaaaggggcaatattgccggagatactacattatctgccgcgtcgagtcacaggagatcaacgtgaatagcgccgccggacgtgctccgacggcgataaaggcatcaGAATAACCGTATTCTGCAGCCAGGGATTTAATCCTGCCTTGTTTGCCCGGCCATCCTGACCGGGCAATGTTCTTTCCTTTAAACCTCAATCTCCGCCATGTCGCCTTTCTCTTGCAACCAGTTGCGGCGATCTTcc                                                                                                                   <  1:384393‑M2/136‑1 (MQ=255)                                                                                                                                         aaaggggcaatattcccggagctactacattatctgccgcgtcgagtcacaggagctcaacgtgaatagcgccgccggacgtgctccgacggcgataaaggcatcaGAATAACCGTATTCTACAGCCAGGGATTTAATCCTGCCTTGTTTGACCCGCCATCCTGACCGGGCAATGTTCTTTCCTTTAAACCTCAATCTCCGCCATGTCGCCTTTATCTTGCAACCAGTTGCGGCGATCTTCCGAGCGTttc                                                                                                          >  1:84056‑M2/107‑251 (MQ=255)                                                                                                                                           aggggcaatattgccggagatactacattatctgccgcgtcgagtcacaggagctccacgtgaatcgcgccgccggacgtgctccgacggagataaaggcatcaGAATAACCGAATTCTACAGCCAGGGATTTAATCCTGCCTTGTTTGCCCGGCCATCCTGACCGGGCAATGTTCTTTCCTTTAAACCTCAATCTACGCCATGTCGCCTTTCTCTTGCAACCAGTTGCGGCGATCTTCCGAGCGTttctt                                                                                                        >  1:197666‑M2/105‑251 (MQ=255)                                                                                                                                                 aatattgccggagatactacattatctgccgcgtcgagtcacaggagatcaacgtgaatagcgccgccggacgtgctccgacggcgataaaggcatcaGAATAACCGTATTCTGCAGCCAGGGATTTAATCCTGCCTTGTTTGCCCGGCCATCCTGACCGGGCAATGTTCTTTCCTTTAAACCTCAATCTCCGCCATGTCGCCTTTCTCTTGCAACCAGTTGCGGCGATCTTCCGAGCGTTTCTTCGcc                                                                                                    >  1:73342‑M2/99‑249 (MQ=255)                                                                                                                                                   tattgccggagatactacattatctgccgcgtcgagtcacaggagatcaacgtgaatagcgccgccggacgtgctccgacggcgataaaggcatcaGAATAACCGTATTCTGCAGCCAGGGATTTAATCCTGCCTTGTTTGCCCGGCCATCCTGACCGGGCAATGTTCTTTCCTTTAAACCTCAATCTCCGCCATGTCGCCTTTCTCTTGCAACCAGTTGCGGCGATCTTCCGAGCGTTTCTTCGCcagc                                                                                                 >  1:67679‑M2/97‑250 (MQ=255)                                                                                                                                                    attgccggagatactacattatctgccgcgtcgagtcacaggagatcaacgtgaatagcgccgccggacgtgctccgacggcgataaaggcatcaGAATAACCGTATTCTGCAGCCAGGGATTTAATCCTGCCTTGTTTGCCCGGCCATCCTGACCGGGCAATGTTCTTTCCTTTAAACCTCAATCTCCGCCATGTCGCCTTTCTCTTGCAACCAGTTGCGGCGATCTTCCGAGCGTTTCTTCGCcagc                                                                                                 >  1:51507‑M2/96‑249 (MQ=255)                                                                                                                                                                actacattatctgccgcgtcgagtcacaggagatcaacgtgaatagcgccgccggacgtgctccgacggcgataaaggcatcaGAATAACCGTATTCTGCAGCCAGGGATTTAATCCTGCCTTGTTTGCCCGGCCATCCTGACCGGGCAATGTTCTTTCCTTTAAACCTCAATCTCCGCCATGTCGCCTTTCTCTTGCAACCAGTTGCGGCGATCTTCCGAGCGTTTCTTCGCCAGCAGCATATCcatcat                                                                                   <  1:574863‑M2/168‑1 (MQ=255)                                                                                                                                                                 ctacattatctgccgcgtcgaatcacaggagatcaacgtgaatagcgccgccggacgtgctccgacggcgataaaggcatcaGAATAACCGTATTCTGCAGCCAGGGATTTAATCCTGCCTTGTTTGCCCGGCCATCCTGACCGGGCAATGTTCTTTCATTTAAACCTCAATCTCCGCCATGTCGCCTTTCTCTTGCAACCAGTTGCGGCGATCTTCCGAGCGTTTATTCGCCAGCAGCATATCCatcatc                                                                                  >  1:465500‑M2/83‑251 (MQ=255)                                                                                                                                                                   ccattacctgccgcgtcgagtcacaggagatcaacgtgaatagcgccgccggacgggctccgacggcgataaaggcatcaGAATAACCGTATTCTGCAGCCAGGGATTTAATCCTGCCTTGTTTGCCCGGCCATCCTGACCGGGCAATGTTCTTTCCTTTAAACCTCAATCGCCGCCATGTCGCCTTTCTCTTGCACCCAGTTGCGGCGATCTTCCGAGCGTTTCTTCGCCAGCAGCATATCCATCATcg                                                                                 <  2:567320‑M2/170‑1 (MQ=255)                                                                                                                                                                   acattacctgccgcgccgagtcacaggagatcaacgtgaatagcgccgccggacgtgctccgacggcgataaaggcatcaGAATAACCGTATTCTGCAGCCAGGGATTTAATCCTGCCTTGTTTGCCCGGCCATCCTGCCCGGGCAATGTTCTTTCCTTTAAACCTCAATCTCCGCCATGTCGCCTTTCTCTTGCAACCAGTTGCGGCGATCTTCCGAGCGTTTCTTCGCCAGCAGCATATCCatcatc                                                                                  <  2:366836‑M2/169‑1 (MQ=255)                                                                                                                                                                     attatctgccgcgtcgagtcacaggagatcaacgtgaatagcgccgccggacgtgctccgacggcgataaaggcatcaGAATAACCGTATTCTGCAGCCAGGGATTTAATCCTGCCTTGTTTGCCCGGCCATCCTGACCGGGCAATGTTCTTTCCTTTAAACCTCAATCTCCGCCATGTCGCCTTTCTCTTGCAACCAGTTGCGGCGATCTTCCGAGCGTTTCTTCGCCAGCAGCATATCCATCATCGCGt                                                                              <  1:189850‑M2/173‑1 (MQ=255)                                                                                                                                                                        atctgccgcgtcgagtcacaggagatcaacgtgaatagcgccgccggacgtgctccgacggcgataaaggcatcaGAATAACCGTATTCTGCAGCCAGGGATTTAATCCTGCCTTGTTTGCCCGGCCATCCTGACCGGGCAATGTTCTTTCCTTTAAACCTCAATCTCCGCCATGTCGCCTTTCTCTTGCAACCAGTTGCGGCGATCTTCCGAGCGTTTCTTCGCCAGCAGCATATCCATCATCGCGTc                                                                             >  1:75437‑M2/76‑249 (MQ=255)                                                                                                                                                                           tgccgcgtcgcctcccaggagatcaacctgactagcgccgacggacgtgctccgtcgggggtaaaggcatcaGAATAACCGTATTCTGCAGCCAGGGATTTAATCCTGCCTTGTTTGCCCGGCCATCCTGACCGGTCAATGTTCTTTCCTTTAAACCTCAATCTCCGCCATGTCGCCTTTTTCTTGCAACCAGTTGCGGCGCTCTTCCGAGCGTTTTTTCGCCAGcagcatatccatcatcctgtacgcac                                                                        >  2:45068‑M2/73‑241 (MQ=255)                                                                                                                                                                                   cgagtcacaggagatcaacgtgaatagcgccgccggacgtgctccgacggcgataaaggcatcaGAATAACCGTATTCTGCAGCCAGGGATTTAATCCTGCCTTGTTTGCCCGGCCATCCTGACCGGGCAATGTTCTTTCCTTTAAACCTCAATCTCCGCCATGTCGCCTTTCTCTTGCAACCAGGTGCGGCGATCTTCCGAGCGTTTCTTCGCCAGCAGCATATCCATCATCGCGTCAGTACGCTGATCg                                                                >  2:92498‑M2/65‑251 (MQ=255)                                                                                                                                                                                     agtcacaggagatcaacgtgaatagcgccgccggacgtgctccgacggcgataaaggcatcaGAATAACCGTATTCTGCAGCCAGGGATTTAATCCTGCCTTGTTTGCCCGGCCATCCTGACCGGGCAATGTTCTTTCCTTTAAACCTCAATCTCCGCCATGTCGCCTTTCTCTTGCAACCAGTTGCGGCGATCTTCCGAGCGTTTCTTCGCCAGCAGCATATCCATCATCGCGTCAGTACGCTGATCGTc                                                              >  1:455956‑M2/63‑251 (MQ=255)                                                                                                                                                                                       tcacaggagatcaacgtgaatagcgccgccggacgtgctccgacggcgataaaggcatcaGAATAACCGTATTCTGCAGCCAGGGATTTAATCCTGCCTTGTTTGCCCGGCCATCCTGACCGGGCAATGTTCTTTCCTTTAAACCTCAATCTCCGCCATGTCGCCTTTCTCTTGCAACCAGTTGCGGCGATCTTCCGAGCGTTTCTTCGCCAGCAGCATATCCATCATCGCGTCAGTACGCTGCTCGTCtt                                                            >  1:185475‑M2/61‑251 (MQ=255)                                                                                                                                                                                          caggagatcaacgtgaatagcgccgccggacgtgctccgacggcgataaaggcatcaGAATAACCGAATTCTGCAGCCAGGGATTTAATCCTGCCTTGTTTGCCCGGCCATCCTGACCGGGCAATGTTCTTTCCTTTAAACCTCAATCTCCGCCATGTCGCCTTTCTCTTGCAACCAGTTGCGGCGATCTTCCGAGCGTTTCTTCGCCAGCAGCATATCCATCATCGCGTCAGTACGCTGATCGTCTtcat                                                         >  1:235102‑M2/58‑251 (MQ=255)                                                                                                                                                                                               gatcaacttgaatagcgccgccggacgtgctccgacggcgataaaggcatcaGAATAACCGTATTCTGCAGCCAGGGATTTAATCCTGCCTTGTTTGCCCGGCCATCCTGACCGGGCAATGTTCTTTCCTTTAAACCTCAAGCTCCGCCATGTCGCCTTTCTCTTGCAACCAGTTGCGGCGATCTTCCGAGCGTTTCTTCGCCAGCAGCATATCCATCATCGCGTCAGTACGCTGATCGTCTTCATCATCg                                                    <  1:138716‑M2/199‑1 (MQ=255)                                                                                                                                                                                                   aacgtgaatagcgccgccggacgtgctccgacggcgataaaggcatcaGAATAACCGTATTCTGCAGCCAGGGATTTAATCCTGCCTTGTTTGCCCGGCCATCCTGACCGGGCAATGTTCTTTCCTTTAAACCTCAATCTCCGCCATGTCGCCTTTCTCTTGCAACCAGTTGCGGCGATCTTCCGAGCGTTTCTTCGCCAGCAGCATATCCATCATCGCGTCAGTACGCTGATCGTCTTCATCATCGAt                                                  >  1:153337‑M2/49‑249 (MQ=255)                                                                                                                                                                                                              cgccgccggacgtgctccgacggcgataaaggcatcaGAATAACCGTATTCTGCAGCCAGGGATTTAATCCTGCCTTGTTTGCCCGGCCATCCTGACCGGGCAATGTTCTTTCCTTTAAACCTCAATCTCCGCCATGTCGCCTTTCTCTTGCAACCAGTTGCGGCGATCTTCCGAGCGTTTCTTCGCCAGCAGCATATCCATCATCGCGTCAGTACGCTGATCGTCTTCATCATCGATAGTCAACTGAAcc                                     >  1:298473‑M2/38‑251 (MQ=255)                                                                                                                                                                                                                  gccggacgtgctccgcccgcgataacggcatcaGAATAACCGTCTTCTGCAGCCAGGGATTTAATCCTGCCTTGTTTGCCCGGCCATCCTGACCGGGCCATGTTCTTTCCTTTAAACCTCAATCTCCGCCATGTCGCCATTCTCTTGCAACCAGTTGCGGCGATCTTCCGAGCGTTTCTTCGCCAGCAGCATATCCATCATAGCGTCAGTACGCTGATCGTCTTCATCATCGATAGTCAACTGCACCAGAc                                 >  2:246379‑M2/34‑251 (MQ=255)                                                                                                                                                                                                                     ggacgtgctccgtaggcgataaaggcatcaGAATAACCGTATTCTGCAGCCAGGGATTTACTCCTGCCTTGTTTGCCCGGCCATCCTGACCGGGCAATGTTCTTTCCTTTAACCCTCAATCTCCGCCATGTCGCCTTTCTCTTGCCACCAGTTGCGGCGATCTTCCGAGCGTTTCTTCGCCAGCCGCATATCCATCACCGCGCCAGTACGCTGATATTCTTCATCATCGATAGTCAACGCcaccagccgcc                              >  2:385278‑M2/31‑246 (MQ=255)                                                                                                                                                                                                                                 acggcgataatgtcatcagAATAACCGCATTCTGCAGCCAGCGATTTAATAATGCCTTGTTTGCACGGCCATCTTGACCGGGCAATGTTCTTTCCTTTAAACCTCAATATCCGCCATGTCGCCTTTCTCTTGCAACCAGTTGCGGCGATCTTCCCAGCGTTTCTTCGCCAGCATCATATCAATAATCGCGTCAGTCCGCTGATCGTCTACATCATCGACAGTCAACTGCACCAGACGGCGCGTGCTCGGAt                  >  2:494201/19‑251 (MQ=255)                                                                                                                                                                                                                                 acggcgataaaggcatcaGAATAACCGTATTCTGCAGCCAGGGATTTAATCCTGCCTTGTTTGCCCGGCCATCCTGACCGGGCAATGTTCTTTCCTTTAAACCTCAATCTCCGCCATGTCGCCTTTCTCTTGCAACCAGTTGCGGCGATCTTCCGAGCGTTTCTTCGCCAGCAGCATATCCATCATCGCGTCAGTACGCTGATCGTCTTCATCATCGATAGTCAACTGCACCAGACGGCGAGTGTTCGGAt                  <  1:56845‑M2/233‑1 (MQ=255)                                                                                                                                                                                                                                      gataaaggcatcaGAATAACCGTATTCTGCAGCCAGGGCTTTAATCCTGCCTTGTTTGCCCGGCCATCCTGACCGGGCAATGTTCTTTCCTTTAAACCTCAATCTCCGCCATGTCGCCTTTCTCTTGCAACCAGTTGCGGAGATCTTCCGAGCGTTTCTTCGCCAGCAGCATATCCATCATCGCGTCAGTACGCTGATCGGCTTCATCATCGATAGTCAACTGCACCAGACGGCAAGTGCTCGGACCAAGc             >  2:98806‑M2/14‑251 (MQ=255)                                                                                                                                                                                                                                       ataaaggcatcaGAATAACCGTATTCTGCAGCCAGGGATTTAATCCTGCCTTGTTTGCCCGGCCATCCTGACCGGGCAATGTTCTTTCCTTTAAACCTCAATCTCCGCCCTGTCGCCTTTCTCTTGCAACCAGTTGCGGCGATCTTCCGAGCGTTTCTTCGCCAGCAGCATATCCATCATCGCGTCAGTACGCTGATCGTCTTCATCATCGATAGTCAACTGCACCAGCCGGCGAGTGttcggataacgcg            >  1:581047‑M2/13‑245 (MQ=255)                                                                                                                                                                                                                                           aggcatcaGAATAACCGTATTCTGCAGCCAGGGATTTAATCCTGCCTTGTTTGCCCGGCCATCCTGACCGGGCAATGTTCTTTCCTTTAAACCTCAATCTCCGCCATGTCGCCTTTCTCTTGCAACCAGTTGCGGCGATCTTCCGAGCGTTTCTTCGCCAGCAGCATATCCATCATCGCGTCAGTACGCTGATCGTCTTCATCATCGATAGTCAACTGCACCAGACGGCGAGTGTTCGGATCAAGCGTGGt        >  1:169434‑M2/9‑251 (MQ=255)                                                                                                                                                                                                                                                 caGAATAACCGTATTCTGCAGCCAGGGATTTAATCCTGCCTTGTTTGCCCGGCCATCCTGACCGGGCAATGTTCTTTCCTTTAAACCTCAATCTCCGCCATGTCGCCTTTCTCTTGCAACCAGTTGCGGCGATCTTCCGAGCGTTTCTTCGCCAGCAGCATATCCATCATCGCGTCAGTACGCTGATCGTCTTCATCATCGATAGTCAACTGCACCAGACGGCGAGTGTTCGGATCAAGCGTGGTTTcgcg  >  1:165632‑M2/3‑251 (MQ=255)                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                    CATCACCGTCAGGCGGTATTGGATCAGTTTGCTAAAATCGTTCCAACCGTACTGAAAGAAGAAGGTTGCCACGGCTATGCGCCAATGGTGGATTGCGCAGCTGGCGTGAGTTTCCAGTCTATGGCACCGGATTCTATCGTGATGATTGAGCAGTGGGAAAGCATCGCGCATCTTGAAGCGCATCTGCAAACCCCGCACATGAAGGCGTATAGCGAAGCCGTAAAAGGTGACGTGCTGGAGATGAATATCCGTATTCTGCAGCCAGGGATTTAATCCTGCCTTGTTTGCCCGGCCATCCTGACCGGGCAATGTTCTTTCCTTTAAACCTCAATCTCCGCCATGTCGCCTTTCTCTTGCAACCAGTTGCGGCGATCTTCCGAGCGTTTCTTCGCCAGCAGCATATCCATCATCGCGTCAGTACGCTGATCGTCTTCATCATCGATAGTCAACTGCACCAGACGGCGAGTGTTCGGATCAAGCGTGGTTTCGCG  >  NC_000913/3173184‑3173674

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