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_0009133,173,4310TA97.3% 257.6 / NA 73I97N (ATC→AAC) ygiNputative quinol monooxygenase YgiN
Reads supporting (aligned to +/- strand):  ref base T (0/0);  major base A (36/35);  minor base C (2/0);  total (38/35)
Fisher's exact test for biased strand distribution p-value = 4.94e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.97e-01

CATCACCGTCAGGCGGTATTGGATCAGTTTGCTAAAATCGTTCCAACCGTACTGAAAGAAGAAGGTTGCCACGGCTATGCGCCAATGGTGGATTGCGCAGCTGGCGTGAGTTTCCAGTCTATGGCACCGGATTCTATCGTGATGATTGAGCAGTGGGAAAGCATCGCGCATCTTGAAGCGCATCTGCAAACCCCGCACATGAAGGCGTATAGCGAAGCCGTAAAAGGTGACGTGCTGGAGATGAATATCCGTATTCTGCAGCCAGGGATTTAATCCTGCCTTGTTTGCCCGGCCATCCTGACCGGGCAATGTTCTTTCCTTTAAACCTCAATCTCCGCCATGTCGCCTTTCTCTTGCAACCAGTTGCGGCGATCTTCCGAGCGTTTCTTCGCCAGCAGCATATCCATCATCGCGTCAGTACGCTGATCGTCTTCATCATCGATAGTCAACTGCACCAGACGGCGAGTGTTCGGATCAAGCGTGGTTTCGCG  >  NC_000913/3173184‑3173674
                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                   
aataaacggcaacgcacggtccatgatcacgtacatggagtagtttaagtaggcgttttccgtaaattcatgtagcgcaaggcgctctgccatatcgctcattaattctgattcctcaacttattcgccagccttcaaaggggcaatattgccggagatactacattatctgccgcgtcgagtcacaggagatcaacgtgaatagcgccgccggacgtgctccgacggcgataaaggcatcaGAATAACCg                                                                                                                                                                                                                                                  >  1:270198‑M2/243‑251 (MQ=255)
   aaacggcaacgcacggtccatgatcacgtacatggagtagtttaagtaggcgttttccgtaaattcatgtagcgcaaggcgctctgccatatcgctcattaattctgattcctcaacttattcgccagccttcaaaggggcaatattgccggagatactacattatctgccgcgtcgagtcacaggagatcaacgtgaatagcgccgccggacgtgctccgacggcgataaaggcatcaGAATAACCGTAt                                                                                                                                                                                                                                               >  1:306609‑M2/240‑251 (MQ=255)
         caacgcacggtccatgatcacgtacatggagtagtttaagtaggcgttttccgtaaattcatgtagcgcaaggcgctctgccatatcgctcattaattctgattcctcaacttattcgccagccttcaaaggggcaatattgccggagatactacattatctgccgcgtcgcgtcacaggagatcaacgtgaatagagccgccggacgtgctccgacggcgataaaggcatcaGAATAACCGTATTCTGc                                                                                                                                                                                                                                          >  2:273557‑M2/234‑250 (MQ=255)
              cacggtccatgatcacgtacatggagtagtttaagtaggcgttttccgtaaattcatgtagcgcaaggcgctctgccatatcgctcattaattctgattcctcaacttattcgccagccttcaaaggggcaatattgccggagatactacattatctgccgcgtcgagtcacaggagatcaacgtgaatagcgccgccggacgtgctccgacggcgataaaggcatcaGAATAACCGTATTCTGCagcccg                                                                                                                                                                                                                                    >  2:449670‑M2/229‑249 (MQ=255)
               acggtccatgatcacgtacatggagtagtttaagtaggcgttttccgtaaattcatgtagcgcaaggcgctctgccatatcgctcattaattctgattcctcaacttattcgccagccttcaaaggggcaatattgccggagatactacattatctgccgcgtcgagtcacaggagatcaacgtgaatagcgccgccggacgtgctccgacggcgataaaggcatcaGAATAACCGTATTCTGCAGCCAgg                                                                                                                                                                                                                                   >  1:280308‑M2/228‑251 (MQ=255)
                 ggtccatgatcacgtacatggagtagtttaagtaggcgttttccgtaaattcatgtagcgcaaggcgctctgccatatcgctcattaattctgattcctcaacttattcgccagccttcaaaggggcaatattgccggagatactacattatctgccgcgtcgagtcacaggagatcaacgtgaatagcgccgccggacgtgctccgacagcgataaagncatcaGAATAACCGTATTCTGCAGCCAggg                                                                                                                                                                                                                                  >  2:232401‑M2/226‑250 (MQ=255)
                               tacatggagtagtttaagtaggcgttttccgtaaattcatgtagcgcaaggcgctctgccatatcgctcattaattctgattcctcaacttattcgccagccttcaaaggggcaatattgccggagatactacattatctgccgcgtcgagtcacaggagatcaacgtgaatagcgccgccggacgtgctccgacggcgataaaggcatcaGAATAACCGTATTCTGCAGCCAGGGATTTAATCCTGCCtt                                                                                                                                                                                                                   >  2:87058‑M2/212‑251 (MQ=255)
                                acatggagtagtttaagtaggcgttttccgtaaattcatgtagcgcaaggcgctctgccatatcgctcattaattctgattcctcaacttattcgccagccttcaaaggggcaatattgccggagatactacattatctgccgcgtcgagtcacaggagatcaacgtgaatagcgccgccggacgtgctccgacggcgataaaggcatcaGAATAACCGTATTCTGCAGCCAGGGATTTAATCCTGCCtt                                                                                                                                                                                                                   <  2:270198‑M2/40‑1 (MQ=255)
                                acatggagtagcttatgtaggcgttttccgtaaatgcaggtagcgcaaggcgctctgccatatcgctcattaattcggatgcctcaacttattcgccagccttcaaaggcgcaagattgccggagatactacattatctgctgcgtcgagtcacaggagatcaactggaatagcgccgccggacttgctccgacggcgataaaggcatcaTAATAACCGTATTCTGCAGCCAGGTATTTAATCCTGCCtt                                                                                                                                                                                                                   <  1:486186‑M2/40‑1 (MQ=255)
                                   gggagtagtttaagtaggcgttttccgtaaattcatgtagcgcaaggcgctctgccatatcgctcattaattctgattcctcaacttattcgccagccttcaaaggggcaatattgccggagatactacattatctgccgcgtcgagtcacaggagatcaacgtgaatagcgccgccggacgtgctccgacggcgataaaggcatcaGAATAACCGTATTCTGCAGCCAGGGATTTAATCCTGCCTTGttt                                                                                                                                                                                                               <  2:473199‑M2/44‑1 (MQ=255)
                                          gtttaagtaggcgttttccgtaaattcatgtagcgcaaggcgctctgccatatcgctcattaattctgattcctcaacttattcgccagccttcaaaggggcaatattgccggagatactacattatctgccgcgtcgagtcacaggagatcaacgtgaatagcgccgccggacgtgctccgacggcgataaaggcatcaGAATAACCGTATTCTGCAGCCAGGGATTTAATCCTGCCTTGTTTGCCCGGc                                                                                                                                                                                                        >  1:369477‑M2/201‑251 (MQ=255)
                                                   ggcttttgccgtaaagtcgtgtagcgcaaggcgctctgccatatcgctcattaattctgattcctcaacttattcgccagccttcaaaggggcaatattgccggggatactacattatcttccgtggcgagtcacagtagatcaacgtgaatagcgccgccggacgtgctccgacggcgataaaggcatcaGAATAACCGTATTCTGCAGCCAGGGATTTAATCCTGCCTTGTTTGCCCGGCCATCCTGAc                                                                                                                                                                                               <  2:377946‑M2/60‑1 (MQ=255)
                                                             gtaaattcatgtagcgcaaggcgctctgccatatcgctcattaattctgattcctcaacttattcgccagccttcaaaggggcaatattgccggagatactacattatctgccgcgtcgagtcacaggagatcaacgtgaatagcgccgccggacgtgctccgacggcgataaaggcatcaGAATAACCGTATTCTGCAGCCAGGGATTTAATCCTGCCTTGTTTGCCCGGCCATCCTGACCGGGCAATGt                                                                                                                                                                                     >  2:28850‑M2/182‑251 (MQ=255)
                                                               aaattcatgtagcgcaaggcgctctgccatatcgctcattaattctgattcctcaacttattcgccagccttcaaaggggcaatattgccggagatactacattatctgccgcgtcgagtcacaggagatcaacgtgaatagcgccgccggacgtgctccgacggcgataaaggcatcaGAATAACCGTATTCTGCAGCCAGGGATTTAATCCTGCCTTGTTTGCCCGGCCATCCTGACCGGGCAATGTTc                                                                                                                                                                                   <  2:23960‑M2/72‑1 (MQ=255)
                                                               aaatgcaggaagcccaaggcgctctgccatatcgctcattaattctgtttcctcaagttattcgccagccttcaaaggggcgatattgccggagattctacattatctgccgcgtcgagtcacaggagatcaacgtgaacagtgccgccggacgtgctccgacggcgataaaggcttcaGAATAACCGTATTCTGCAGCCAGGGATTTAATCCTGCCTTGTTTGCCCTGCCATCCTGACCGGGCATTGTTc                                                                                                                                                                                   <  2:377015‑M2/72‑1 (MQ=255)
                                                                    catgtagcgcaagtcgctctgccatatcgctcattaattctgattcctcaacttattcgccagccttcaaaggggcaatattgccggagatactacattatctgccgcgtcgagtcacaggagatcaacgtgaatagcgccgccggacgtgctccgacggcgataaaggcatcaGAATAACCGTATTCTGCAGCCAGGGATTTAATCCTGCCTTGTTTGCCCGGCCATCCTGACCGGGCAATGTTCTTTc                                                                                                                                                                               <  1:411650‑M2/76‑1 (MQ=255)
                                                                     atgtagcgcaaggcgctctgccatatcgctcattaattctgattcctcaacttattcgccagccttccaaggggcaatattgccggagatactacattatctgccgcgtcgagtcacaggagatcaacgtgaatagcgccgccggacgtgctccgacggcgctaaaggcatcaGAATACCCGTATTCTGCATCCAGGGATTTAATCCTGCCTTGTTTGCCCGGCCATCCTGCCCGGGCACTGTTCTTTCCt                                                                                                                                                                             >  2:209943‑M2/174‑251 (MQ=255)
                                                                       gtagcgcaaggcgctctgccatatcgctcattaattctgattcctcaacttattcgccagccttcaaaggggcaatattgccggagatactacattatctgccgcgtcgagtcacaggagatcaacgtgaatagcgccgccggacgtgctccgacggcgataaaggcatcaGAATAACCGTATTc                                                                                                                                                                                                                                             <  1:226168‑M2/14‑1 (MQ=255)
                                                                       gtagcgcaaggcgctctgccatatcgctcattaattctgattcctcaacttattcgccagccttcaaaggggcaatattgccggagatactacattatctgccgcgtcgagtcacaggagatcaacgtgaatagcgccgccggacgtgctccgacggcgataaaggcatcaGAATAACCGTATTc                                                                                                                                                                                                                                             >  2:226168‑M2/172‑185 (MQ=255)
                                                                       gtagcgcaaggcgctctgccatatcgctcattaattctgattcctcaacttattcgccagccttcaaaggggcaatattgccggagatactacattatcggccgcgtcgagtcacaggagatcaacgtgaatagcgccgccggacgtgctccgacggcgataaaggcatcaGAATAACCGTATTCTGCAGCCAGGGATGTAATCCTGCCTTGTTTGCCCGGCCATCCTGACCGGGCAATGTTCTTTCCttt                                                                                                                                                                           <  1:34391‑M2/80‑1 (MQ=255)
                                                                          gcgcaaggcgctctgccatatcgctcattaattctgattcctcaacttattcgccagccttcaaaggggcaatattgccggagatactacattatctgccgcgtcgagtcacaggagatcaacgtgaatagcgccgccggacgtgctccgacggcgataaaggcatcaGAATAACCGTATTCTGCAGCCAGGGATTTAATCCTGCCTTGTTTGCCCGGACATCCTGACCGGGCAATGGTCTTTCCTTTaa                                                                                                                                                                         >  1:197304‑M2/169‑250 (MQ=255)
                                                                           cgcaaggcgctctgccatatcgctcattaattctgattcctcaacttattcgccagccttcaaaggggcaatattgccggagatactacattatctgccgcgtcgagtcacaggagatcaacgtgaatagcgccgccggacgtgctccgacggcgataaaggcatcaGAATAACCGTATTCTGCAGCCAGGGATTTAATCCTGCCTTGTTTGCCCGGCCATCCTGACCGGGCAATGTTCTTTCCTTTAAAc                                                                                                                                                                       >  1:418752‑M2/168‑251 (MQ=255)
                                                                              aaggcgctctgccatatcgctcattaattctgattcctcaacttattcgccagccttcaaaggggcaatattgccggagatactacattatctgccgcgtcgagtcacaggagatcaacgtgaatagcgccgccggacgtgctccgacggcgataaaggcatcaGAATACCCGTATTCTGCAGCCAGGGATTTAATCCTGCCTTGTTTGCCCGGCCATCCTGACCGGGCAATGTTCTTTCCTTTAATCCTc                                                                                                                                                                    >  2:211486‑M2/165‑251 (MQ=255)
                                                                                      ctgccttatcgctcagtaattctgattcctcaacttattcgccagccttcaaaggggcaatattgccggagatactacattatctgccgcgtcgagtcacaggagatcaacgtgaatagcgccgccggacgtgctccgacggcgataaaggcatcaGAATAACCGTATTCTGCAGCCAGGGATTTAATCCTGCCTTGTTTGCCCGGCCATCCTGACCGGGCAATGTTCTTTCCTTTAAACCTCAAtctc                                                                                                                                                              <  2:480728‑M2/93‑1 (MQ=255)
                                                                                          catatcgctcattaattctgattcctcaacttattcgccagccttcaaaggggcaatattgccggagatactacattatctgccgcgtcgagtcacaggcgatcaacgtgaatagcgccgccggacgtgctccgacggcgataaaggcatcaGAATAACCGTATTCTGCAGCCAGGGATTTAATCCTGCCTTGTTTGCCCGGCCATCCTGACCGGGCAATGTTCTTTCCTTTAAACCTCAATCTCCGCCAt                                                                                                                                                        >  1:498012‑M2/153‑251 (MQ=255)
                                                                                             atcgctcattaattctgattcctcaacttatgcgccagccttcaaaggggcaatattgccggagatactacattatctgccgcgtcgagtcacaggagatcaacgtgaatagcgccgccggacgtgctccgacggcgataaaggcatcaGAATAACCGTATTCTGCAGCCAGGGATTTAATCCTGCCTTGTTTGCCCGGCCATCCTGACCGGGCAATGTTCTTTCCTTTAAACCTCAATCTCCGCCATGt                                                                                                                                                      <  1:433013‑M2/101‑1 (MQ=255)
                                                                                                gctcattaattctgattcctcaacttattcgccagccttcaaaggggcaatattgccggagatactacattatctgccgcgtcgagtcacaggagatcaacgtgaatagcgccgccggacgtgctccgacggcgataaaggcatcaGAATAACCGTATTCTGCAGCCAGGGATTTAATCCTGCCTTGTTTGCCCGGCCATCCTGACCGGGCAATGTTCTTTCCTTTAAACCTCAATCTCCGCCATGTCGcc                                                                                                                                                  >  1:217277‑M2/147‑251 (MQ=255)
                                                                                                     ttaattctgattcctcaacttattcgccagccttcaaaggggcaatattgccggagatactacattatctgccgcgtcgagtcacaggagatcaacgtgaatagcgccgccggacgtgctccgacggcgataaaggcatcaGAATAACCGTATTCTGCAGCCAGGGATTTAATCCTGCCTTGTTTGCCCGGCCATCCTGACCGGGCAATGTTCTTTCCTTTAAACCTCAATCTCCGCCATGTCGCCTTtct                                                                                                                                             >  1:333483‑M2/142‑251 (MQ=255)
                                                                                                        attctgatcccgcaacttattcgccagccttnaaaggggcaagatggccggagatactacattatctgccgcgtcgagtcacaggagatcaacgtgaatagcgacgccggacgtgctccgacggcgataaaggcatcaGAATAACCGTATTCTGCAGCCAGGGATTTAATCCTGCCTTGTTTGCCCGGCCACCCTGACCGGGCAATGTTCTTTCCTTTAAACCTCAATCTCCGCCATGTCGCCTTTCTCtt                                                                                                                                          <  2:256435‑M2/113‑1 (MQ=255)
                                                                                                             gattcctccacttattcgccagccttcaaaggggcaatattgccggagatactacattatctgccgcgtcgagtcacaggagatcaacgtgaatagcgccgccggacgtgctccgacggcgataaaggcatcaGAATAACCGTATTCTGCAGCCAGGGATTTAATCCTGCCTTGTTTGCCCGGCCATCCTGACCGGGCAATGTTCTTTCCTTTAAACCTCAATATCCGCCATGTCGCCTTTCTCTTGCAAc                                                                                                                                     >  1:88267‑M2/134‑251 (MQ=255)
                                                                                                             gattcctcaacttattcgccagccttcaaaggggcaatattgccggagatactacattatctgccgcgtcgagtcacaggagatcaacgtgaatagcgccgccggacgtgctccgacggcgataaaggcatcaGAATAACCGTATTCTGCAGCCAgg                                                                                                                                                                                                                                   <  2:388903‑M2/24‑1 (MQ=255)
                                                                                                             gattcctcaacttattcgccagccttcaaaggggcaatattgccggagatactacattatctgccgcgtcgagtcacaggagatcaacgtgaatagcgccgccggacgtgctccgacggcgataaaggcatcaGAATAACCGTATTCTGCAGCCAgg                                                                                                                                                                                                                                   >  1:388903‑M2/134‑157 (MQ=255)
                                                                                                                         tattcgccagccttcaaaggggcaatattgccggagatactacattatctgccgcgtcgagtcacaggagatcaacgtgaatagcgccgccggacgtgctccgacggcgataaaggcatcaGAATAACCGTATTCTGCAGCCAGGGATTTAATCCTGCCTTGTTTGCCCGGCCATCCTGACCGGGCAATGTTCTTTCCTTTAAACCTCAATCTCCGCCANGTCGCCTTTCTCTTGCAACCAGTTgcggcg                                                                                                                          >  2:225811‑M2/122‑250 (MQ=255)
                                                                                                                          attcgccagccttcaaaggggcaatattgccggagatactacattatctgccgcgtcgagtcacaggagatcaacgtgaatagcgccgccggacgtgctccgacggcgataaaggcatcaGAATAACCGTATTCTGCAGCCAGGGATTTAATCCTGCCTTGTTTGCCCGGCCATCCTGACCGGGCAATGTTCTTTCCTTTAAACCTCAATCTCCGCCATGTCGCCTTTCTCTTGCAACCAGTTGCGGCGa                                                                                                                         <  1:28850‑M2/130‑1 (MQ=255)
                                                                                                                                cagcattaaaaggggcaatattgccggagatactacattatctgccgcgtcgaggcacaggagatcaacgtgaatagcgccgccggacgtgctccgacggcgataaaggcatcaGAATAACCGTATTCTGCAGCCAGGGATTTAATCCTGCCTTGTTTGCCCGGCCATCCTGACCGGGCAATGTTCTTTCCTTTAAACCTCAATCTCCGCCATGTCGCCTTTCTCTTGCAACCAGTTGCGGCGATCTTCCg                                                                                                                  <  2:6961‑M2/137‑1 (MQ=255)
                                                                                                                                           ggggcaatattgccggagatactacattatctgccgcgtcgagtcacaggagatcaacgtgaatagcgccgccggacgttctccgacggcgataaaggcatcaGAATAACCGTATTCTGCAGCCAGGGATTTAATCCTGCCTTGTTTGCCCGGCCATCCTGACCGGGCAATGTTCTTTCCTTTAAACCTCAATCTCCGCCATGTCGCCTTTCTCTTGCAACCAGTTGCGGCGATCTTCCGAGCGTttcttc                                                                                                       <  2:306609‑M2/148‑1 (MQ=255)
                                                                                                                                                   attgccggagatactacattatctgccgcgtcgagtcacaggagatcaacgtgaatagcgccgccggacgtgctccgacggcgataaaggcatcaGAATAACCGTATTCTGCAGCCAGGGATTTAATCCTGCCTTGTTTGCCCGGCCATCCTGACCGGGCAATGTTCTTTCCTTTAAACCTCAATCTCCGCCATGTCGCCTTTCTCTTGCAACCAGTTGCGGCGATCTTCCGAGCGTTTCTTCGCcagcag                                                                                               <  1:87058‑M2/156‑1 (MQ=255)
                                                                                                                                                   attgccggagatactacattatctgccgcgtcgagtcacaggagatcaacgtgaatagcgccgccggacgtgctccgacggcgataaaggcatcaGAATAACCGTATTCTGCAGCCAGGGATTTAATCCTGCCTTGTTTGCCCGGCCATCCTGACCGGGCAATGTTCTTTCCTTTAAACCTCAATCTCCGCCATGTCGCCTTTCTCTTGCAACCAGTTGCGGCGATCTTCCGAGCGTTTCTTCGCcagca                                                                                                <  1:146353‑M2/155‑1 (MQ=255)
                                                                                                                                                    ttgccggagatactacattatctgccgcgtcgagtcacaggagatcaacgtgaatagcgccgccggacgtgctccgacggcgataaaggcatcaGAATAACCGTATTCTGCAGCCAGGGATTTAATCCTGCCTTGTTTGCCCGGCCATCCTGACCGGGCAATGTTCTTTCCTTTAAACCTCAATCTCCGCCATGTCGCCTTTCTCTTGCAACCAGTTGCGGCGATCTTCCGCGCGTTTCTTCGCCatcagc                                                                                              >  2:140845‑M2/95‑251 (MQ=255)
                                                                                                                                                             atactacattatctgccgcgtcgagtcacaggagatcaacgtgaatagcgccgccggacgtgctccgacggcgataaaggcatcaGAATAACCGTATTCTGCAGCCAGGGATTTAATCCTGCCTTGTTTGCCCGGCCATCCTGACCGGGCAATGTTCTTTCCTTTAAACCTCAATCTCCGCCATGTCGCCTTTCTCTTGCAACCAGTTGCGGCGATCTTCCGAGCGTTTCTTCGCCAGCAGCATATCcatc                                                                                     >  2:489058‑M2/86‑251 (MQ=255)
                                                                                                                                                             atactacattatctgccgcgtcgagtcacaggagatcaacgtgaatagcgccgccggacgtgctccgacagcgataaaggcatcaGAATAACCGTATTCTGCAGCCAGGGATTTAATCCTGCCTTGTTTGCCCGGCCATCCTGACCGGGCAATGTTCTTTCCTTTAAACCTCAATCTCCGCCATGTCGCCTTTCTCTTGCAACCAGTTGCGGCGATCTTCCGAGCGTTTCTTCGCCAGCAGCATATCcatc                                                                                     <  1:232401‑M2/166‑1 (MQ=255)
                                                                                                                                                               actacattatctgccgcgtcgagtcacaggagatcaacgtgaatagcgccgccggacgtgctccgacggcgataaaggcatcaGAATAACCGTATTCTGCAGCCAGGGATTTAATCCTGCCTTGTTTGCCCGGCCATCCTGACCGGGCAATGTTCTTTcc                                                                                                                                                                              <  1:387202‑M2/77‑1 (MQ=255)
                                                                                                                                                               actacattatctgccgcgtcgagtcacaggagatcaacgtgaatagcgccgccggacgtgctccgacggcgataaaggcatcaGAATAACCGTATTCTGCAGCCAGGGATTTAATCCTGCCTTGTTTGCCCGGCCATCCTGACCGGGCAATGTTCTTTcc                                                                                                                                                                              >  2:387202‑M2/84‑160 (MQ=255)
                                                                                                                                                                  acattatctgccgcgtcgcgtcacaggagatcaacgtgaatagcgccgcccgacgtgctccgacggcgataaaggcatcaGAATAACCGTATTCTGCAGCCAGGGATTTAATCCTGCCTTGTTTGCCCGGCAATCCTGACCGGGCAATGTTCTTTCCTTTAACCCTCAATCTCCGCCATGTCGCCTTTCTCTTGCAACCAGTTGCGGCGATCGTCCGAGCGTTTCTTCGCCAGCAGCATATCCATCCTcgc                                                                                >  2:119625‑M2/81‑251 (MQ=255)
                                                                                                                                                                  acattatctgccgcgtcgagtcacaggagatcaacgtgaatagcgccgccggacgtgctccgacggcgataaaggcatcaGAATAACCGTATTCTGCAGCCAGGGATTTAATCCTGCCTTGTTTGCCCGGCCATCCTGACCGGGCAATGTTCTTTCCTTTAAACCTCAATCTCCGCCATGTCGCCTTTCTCTTGCAACCAGTTGCGGCGATCTTCCGAGCGTTTCTTCGCCCGCAGCATATCCATCATcgc                                                                                >  2:429358‑M2/81‑251 (MQ=255)
                                                                                                                                                                  acattatctgccgcgtcgagtcacaggagatcaacgtgaatagcgccgccggacgtgctccgacggcgataaaggcatcaGAATAACCGTATTCTGCAGCCAGGGATTTAATCCTGCCTTGTTTGCCCGGCCATCCTGACCGGGCAATGTTCTTTCCTTTAAACCTCAATCTCCGCCATGTCGCCTTTCTCTTGCAACCAGTTGCGGCGATCTTCCGAGCGTTTCTTCGCCAGCAGCATATCCATCATcg                                                                                 >  2:470614‑M2/81‑250 (MQ=255)
                                                                                                                                                                      tatctgccgcgtcgagtcacaggagatcaacgtgaatagcgccgccggacgtgctccgacggcgataaaggcatcaGAATAACCGTATTCTGCAGCCAGGGATTTAATCCTGCCTTGTTTGCCCGGCCATCCTGACCGGGCAATGTTCTTTCCTTTAAACCTCAATCTCCGCCATGTCGCCTTTCTCTTGCAACCAGTTGCGGCGATCTTCCGAGCGTTTCTTCGCCAGCAGCATATCCATCATCGCGTc                                                                             >  1:119192‑M2/77‑250 (MQ=255)
                                                                                                                                                                       atctgccgcgtcgagtcacaggagatcaacgtgaatagcgccgccggacgtgctccgacggcgataaaggcatcaGAATAACCGTATTCTGCAGCCAGGGATTTAATCCTGCCTTGTTTGCCCGGCCATCCTGACCGGGCAATGTTCTTTCCTTTAAACCTCAATCTCCGCCATGTCGCCTTTCTCTTGCAACCAGTTGCGGCGATCTTCCGAGCGTTTCTTCGCCAGCAGCATATCCATCATCGCGTc                                                                             >  1:325098‑M2/76‑249 (MQ=255)
                                                                                                                                                                         ctgacgcgtcgagtcacaggagatcaacgtgaatagcgacgccgtaagtgctccgacggcgataaaggcatcaGAATAACCGTATTCTGCAGCCAGGGATTTAATCCTGCCTTGTTTGCCCGGCCATCCTGACCGGTCAATGTTCTTTCCTTTCAACCTCAATCTCCGCCATGTCGCCTTTCTCTTGCAACCAGTTGCTGCGCTCTTCCGAGCGCTTCTTCGCTAGCAGCGAACCCATCATCGCGTCAGt                                                                          >  2:75967‑M2/74‑250 (MQ=255)
                                                                                                                                                                          tgccgcgtcgagtcacaggagatcaacgtgaatagcgccgccggacgtgctccgacggcgataaaggcatcaGAATAACCGTATTCTGCAGCCAGGGATTTAATCCTGCCTTGTTTGCCCGGCCATCCTGACCGGGCAATGTTCTTTCCTTTAAACCTCAATCTCCGCCATGTCGCCTTTCTCTTGCAACCAGTTGCGGCGATCTTCCGAGCGTTTCTTCGCCAGCAGCATATCCATCATCGCGTCAGTAc                                                                        >  2:497737‑M2/73‑251 (MQ=255)
                                                                                                                                                                           gccgcgtcgagtcacaggagatcaacgtgaatagcgccgccggacgtgctccgacggcgataaaggcatcaGAATAACCGTATTCTGCAGCCAGGGATTTAATCCTGCCTTGTTTGCCCGGCCATCCTGACCGGGCAATGTTCTTTCCTTTAAACCTCAATCTCCGCCATGTCGCCTTTCTCTTGCAACCAGTTGCGGCGATCTTCCGAGCGTTTCTTCGCCAGCAGCATATCCATCATCGCGTCAGt                                                                          <  2:369477‑M2/177‑1 (MQ=255)
                                                                                                                                                                           gccgcggcgagtcacaggagatcaacgtgaatagcgccgccggacgtgctccgacggcgataaaggcatcaGAATAACCGTATTCTGCAGCCAGGGATTTAATCCTGCCTTGTTTGCCCGGCCATCCTGACCGGGCAATGTTCTTTCCTTTAAACCTCAATCGCCGCCATGTCGCCTTTCTCTTGCAACCAGTTGCGGCGATCTTCCGAGCGTTTCTTCGCCAGCAGCATATCCATCATCGCGTCAGt                                                                          <  1:281142‑M2/177‑1 (MQ=255)
                                                                                                                                                                             cgcgtcgagtcacaggagatcaacgtgaatagcgccgccggacgtgctccgacggcgataaaggcatcaGAATAACCGTATTCTGCAGCCAGGGATTTAATCCTGCCTTGTTTGCCCGGCCATCCTGACCGGGCAATGTTCTTTCCTTTAAACCTCAATCTCCGCCATGTCGCCTTTCTCTTGCAACCAGTTGCGGCGATCTTCCGAGCGTTTCTTCGCCAGCAGCATATCCATCATCGCGTCAGTACGCt                                                                     <  1:200953‑M2/182‑1 (MQ=255)
                                                                                                                                                                                gtcgagtcacaggagatcaacgtgaatagcgccgccggacgtgctccgacggcgataaaggcatcaGAATAACCGTATTCTGCAGCCAGGGATTTAATCCTGCCTTGTTTGCCCGGCCATCCTGACCGGGCAATGTTCTTTCCtt                                                                                                                                                                            <  2:175483‑M2/79‑1 (MQ=255)
                                                                                                                                                                                gtcgagtcacaggagatcaacgtgaatagcgccgccggacgtgctccgacggcgataaaggcataaGAATAACCGTATTCTGCAGCCAGGGATTTAATCCTGCCTTGTTTGCCCGGCCATCCTGACCGGGCAATGTTCTTTCCtt                                                                                                                                                                            >  1:175483‑M2/67‑145 (MQ=255)
                                                                                                                                                                                       cacaggagatcaacgtgaatagcgccgccggacgtgctccgacggcgataaaggcatcaGAATAACCGTATTCTGCAGCCAGGGATTTAATCCTGCCTTGTTTGCCCGGCCATCCTGACCGGGCAATGTTCTTTCCTTTAAACCTCAATCTCCGCCATGTCGCCTTTCTCTTGCAACCAGTTGCGGCGATCTTCCGAGCGTTTCTTCGCCAGCAGCATATCCATCATCGCGTCAGTACGCTGATCGTCtt                                                            <  1:17904‑M2/191‑1 (MQ=255)
                                                                                                                                                                                        acaggagatcaacgtgaatagcgccgccggacgtgctccgacggcgataaaggcatcaGAATAACCGTATTCTGCAGCCAGGGATTTAATCCTGCCTTGTTTGCCCGGCCATCCTGACCGGGCAATGTTCTTTCCTTTAAACCTCAATCTCCGCCATGTCGCCTTTCTCTTGCAACCAGTTGCGGCGATCTTCCGAGCGTTTCTTCGCCAGCAGCATATCCATCATCGCGTCAGTACGCTGATCGTCTtca                                                          <  1:388441‑M2/193‑1 (MQ=255)
                                                                                                                                                                                             agatcaacgtgaatagcgccgcctgacgtgccccgacggcgataaaggcatcaGAATAACCGTATTCTGCAGCCAGGGATTTAATCCTGCCTTGTTTGCCCGGCCATCCTGACCGGGCAATGTTCTTTCCTTTAAACCTCAATCTCCGCCATGTCGCCTTTCTCTTGCAACCAGTTGCGGCGATCTTCCGAGCGTTTCTTCGCCAGCAGCATATCCATCATCGCGTCAGTACGCTGATCGTCTTCagcatc                                                     <  2:127231‑M2/198‑1 (MQ=255)
                                                                                                                                                                                                 caacgtgaatagcgccgccggacgtgctccgacggcgataaaggcatcaGAATAACCGTATTCTGCAGCCAGGGATTTAATCCTGCCTTGTTTGCCCGGCCATCCTGACCGGGCAATGTTCTTTCCTTTAAACCTCAATCTCCGCCATGTCGCCTTTCTCTTGCAACCAGTTGCGGCGATCTTCCGAGCGTTTCTTCGCCAGCAGCATATCCATCATCGCGTCAGTACGCTGTTCGTCTTCATCATCGATa                                                 <  2:418752‑M2/202‑1 (MQ=255)
                                                                                                                                                                                                  aacgtgaatagcgccgccggacgtgctccgacggcgataaaggcatcaGAATAACCGTATTCTGCAGCCAGGGATTTAATCCTGCCTTGTTTGCCCGGCCATCCTGACCGGGCAATGTTCTTTCCTTTAAACCTCAATCTCCGCCATGTCGCCTTTCTCTTGCAACCAGTTGCGGCGATCTTCCGAGCGTTTCTTCGCCAGCAGCATATCCATCATCGCGTCAGTACGCTGATCGTCTTCATCATCGAt                                                  >  1:494309‑M2/49‑249 (MQ=255)
                                                                                                                                                                                                        aatagctccgccggacgtgctccgacggcgataaaggcatcaGAATAACCGTATTCTGCAGCCAGGGATTTAATCCTGCCTTGTTTGCCCGGCCATCCTGACCGGGCAATGTTCTTTCCTTTAAACCTCAATCTCCGCCATGTCGCCTTTCTCTTGCAACCAGTTGCGGCGATCTTCCGAGCGTTTCTTCGCCAGCAGCATATCCATCATCGCGTCAGTACGCTGATCGTCTTCATCATCGATAGTCAACt                                          <  1:397036‑M2/209‑1 (MQ=255)
                                                                                                                                                                                                               ccgccggacgtgctccgacggcgataaaggcatcaGAATAACCGTATTCTGCAGCCAGGGATTTAATCCTGCCTTGTTTGCCCGGCCATCCTGACCGGGCAATGTTCTTTCCTTTAAACCTCAATCTCCGCCATGTCGCCTTTCTCTTGCAACCAGTTGCGGCGATCTTCCGAGCGTTTCTTCGCCAGCAGCATATCCATCATCGCGTCAGTACGCTGATCGTCTTCATCATCGATAGTCAACTGCACCAg                                   <  2:454052‑M2/216‑1 (MQ=255)
                                                                                                                                                                                                                    ggacgtgctccgacggcgataaaggcatcaGAATAACCGTATTCTGCAGCCAGGGATTTAATCCTGCCTTGTTTGCCCGGCCATCCTGACCGGGCAATGTTCTTTCCTTTAAACCTCAATCTCCGCCATGTCGCCTTTCTCTTGCAACCAGTTGCGGCGATCTTCCGAGCGTTTCTTCGCCAGCAGCATATCCATCATCGCGTCAGTACGCTGATCGTCTTCATCATCGATAGTCAACTGCACCAGACGGc                              >  1:281453‑M2/31‑251 (MQ=255)
                                                                                                                                                                                                                            cccgacggcgataaaggcatcaGAATAACCGTATTCTGCAGCCAGGGATTTAATCCTGCCTTGTTTGCCCGGCCATCCTGACCGGGCAATGTTCTTTCCTTTAAACCTCAATCTCCGCCATGTCGCCTTTCTCTTGCAACCAGTTGCGGCGATCTTCCGAGCGTTTCTTCGCCAGCAGCATATCCATCATCGCGTCAGTACGCTGATCGTCTTCATCATCGATAGTCAACTGCACCAGACGGCGAGTGTTc                      <  2:374292‑M2/229‑1 (MQ=255)
                                                                                                                                                                                                                            cccgacggcgataaaggcatcaGAATAACCGTATTCTGCAGCCAGGGATTTAATCCTGCCTTGTTTGCCCGGCCATCCTGACCGGGCAATGTTCTTTCCTTTAAACCTCAATCTCCGCCATGTCGCCTTTCTCTTGCAACCAGTTGCGGCGATCTTCCGAGCGTTTCTTCGCCAGCAGCATATCCATCATCGCGTCAGTACGCTGATCGTCTTCATCATCGATAGTCAACTGCACCAGACGGCGAGTGTTc                      <  1:225811‑M2/229‑1 (MQ=255)
                                                                                                                                                                                                                             ccgacggcgataaaggcatcaGAATAACCGTATTCTGCAGCCAGGGATTTAATCCTGCCTTGTTTGCCCGGCCATCCTGACCGGGCAATGTTCTTTCCTTTAAACCTCAATCTCCGCCATGTCGCCTTTCTCTTGCAACCAGTTGCGGCGATCTTCCGAGCGTTTCTTCGCCAGCAGCATATCCATCATCGCGTCAGTACGCTGATCGTCTTCATCATCNATAGTCAACTGCACCAGACGGCGAGTGTTCg                     >  2:249320‑M2/22‑251 (MQ=255)
                                                                                                                                                                                                                             ccgacggcgataaaggcatcaGAATAACCGTATTCTGCAGCCAGGGATTTAATCCTGCCTTGTTTGCCCGGCCATCCTGACCGGGCAATGTTCTTTCCTTTAAACCTCAATCTCCGCCATGTCGCCTTTCTCTTGCAACCAGTTGCGGCGATCTTCCGAGCGTTTCTTCGCCAGCAGCATATCCATCATCGCGTCAGTACGCTGATCGTCTTCATCATCGATAGTCAACTGCACCAGACGGCGAGTGTTCg                     <  1:68420‑M2/230‑1 (MQ=255)
                                                                                                                                                                                                                                  ggcgataaaggcatcaGAATAACCGTATTCTGCAGCCAGGGATTTAATCCTGCCTTGTTTGCCCGGCCATACTGACCGGGCAATGTTCTTTCCTTTAAACCTCAATCTCCGCCATGTCGCCTTTCTCTTGCAACCCGTTGCGGCGATCTTCCGAGCGTTTCTTCGCCAGCAGCATATCCATCATCGCGTCAGTACGCTGATCGTCTTCATCATCGATAGTCAACTGCACCAGCCGGCGCGTGTTCGGATc                 >  2:479672‑M2/17‑250 (MQ=255)
                                                                                                                                                                                                                                   gcgataaaggcatcaGAATAACCGTATTCTGCAGCCAGGGATTTAATCCTGCCTTGTTTGCCCGGCCATCCTGACCGGGCAATGTTCTTTCCTTTAAACCTCAATCTCCGCCATGTCGCCTTTCTCTTGCAACCAGTTGCGGCGATCTTCCGAGCGTTTCTTCGCCAGCAGCATATCCATCATCGCGTCAGTACGCTGATCGTCTTCATCATCGATAGTCAACTGCACCAGACGGCGAGTGTTCGGATCa                >  1:460387‑M2/16‑250 (MQ=255)
                                                                                                                                                                                                                                   gcgataaaggcatcaGAATAACCGTATTCTGCAGCCAGGGATTTAATCCTGCCTTGTTTGCCCGGCCATCCTGACCGGGCAATGTTCTTTCCTTTAAACCTCAATCTCCGCCATGTCGCCTTTCTCTTGCAACCAGTTGCGGCGATCTTCCGAGCGTTTCTTCGCCAGCAGCATATCCATCATCGCGTCAGTACGCTGATCGTCTTCATCATCGATAGTCAACTGCACCAGACGGCGAGTGTTCGGATCa                >  1:132341‑M2/16‑250 (MQ=255)
                                                                                                                                                                                                                                        aaaggcatcaGAATAACCGTATTCTGCAGCCAGAGATTTAATCCTGCCTTGTTTGCCCGGCCATCCTGACCGGGCAATGTTCTTTCCTTTAAACCTCAATCTCCGCCATGTCGCCTTTCTCTTGCAACCAGTTGCGGCGATCTTCCGAGCGTTTCTTCGCCAGCAGCATATCCATCATCGCGTCAGTACGCTGATCGTCTTCATCATCGATAGTCAACTGCCCCAGACGGCGTGTGTTCGGATCAAGAGTg          >  2:11547‑M2/11‑251 (MQ=255)
                                                                                                                                                                                                                                             catcaGAATAACCGTATTCCGCAGCCAGGGATTTGATCCTGCTTTGTTTGCCCGGCCATCCTTCCCGGGCAATGTTCTTTCCTTTAAACCTCAATCTCCGCCATGTCGCCTTTCTCTTGCAACCAGTTGCGGCGGTCTTCCGAGCGTTTCTTCGCCAGCAGCATATCCATCATCGCGTCAGTACGCTGATCGTCTTCATCATCGATAGTCAACTGCACCAGACGTCGAGTGTTCGGATCAAGCGTGGTTTc     <  2:435313‑M2/246‑1 (MQ=255)
                                                                                                                                                                                                                                                caGAATAACCGTATTCTGCAGCCAGGGATTTAATCCTGCCTTGTTTGCCCGGCCATCCTGACCGGGCAATGTTCTTTCCTTTAAACCTCAATCTCCGCCATGTCGCCTTTCTCTTGCAACCAGTTGCGGCGATCTTCCGAGCGTTTCTTCGCCAGCAGCATATCCATCATCGCGTCAGTACGCTGATCGTCTTCATCATCGATAGTCAACTGCACCAGACGGCGAGTGTTCGGATCAAGCGTGGTTTcgcg  <  1:130259‑M2/249‑1 (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 < 16 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: