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,4310TA98.7% 286.4 / NA 77I97N (ATC→AAC) ygiNputative quinol monooxygenase YgiN
Reads supporting (aligned to +/- strand):  ref base T (0/1);  new base A (37/39);  total (37/40)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.97e-01

CACCGTCAGGCGGTATTGGATCAGTTTGCTAAAATCGTTCCAACCGTACTGAAAGAAGAAGGTTGCCACGGCTATGCGCCAATGGTGGATTGCGCAGCTGGCGTGAGTTTCCAGTCTATGGCACCGGATTCTATCGTGATGATTGAGCAGTGGGAAAGCATCGCGCATCTTGAAGCGCATCTGCAAACCCCGCACATGAAGGCGTATAGCGAAGCCGTAAAAGGTGACGTGCTGGAGATGAATATCCGTATTCTGCAGCCAGGGATTTAATCCTGCCTTGTTTGCCCGGCCATCCTGACCGGGCAATGTTCTTTCCTTTAAACCTCAATCTCCGCCATGTCGCCTTTCTCTTGCAACCAGTTGCGGCGATCTTCCGAGCGTTTCTTCGCCAGCAGCATATCCATCATCGCGTCAGTACGCTGATCGTCTTCATCATCGATAGTCAACTGCACCAGACGGCGAGTGTTCGGATCAAGCGTGGTTT  >  NC_000913/3173187‑3173670
                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                               
aaacggcaacgcacggtccatgatcacgtacatggcgtagtttaagtaggcgttttccgtaaattcatgtagcgcaaggcgctctgccatatcgctcattaattctgattcctcaacttattcgccagccttcaaaggggcaatattgccggagatactacattatctgccgcgtcgagtcacaggagatcaacgtgaatagcgccgccggacgtgatccgacggcgataaaggcatcaGAATAACcgtct                                                                                                                                                                                                                                           >  2:530139‑M2/240‑249 (MQ=255)
       aacgcacggtccatgatcacgtacatggagtagtttaagtaggcgttttccgtaaattcatgtagcgcaaggcgctctgccatatcgctcattaattctgattcctcaacttattcgccagccttcaaaggggcaatattgccggagatactacattatctgccgcgtcgagtcacaggagatcaacgtgaatagcgccgccggacgtgctccgacggcgataaaggcatcaGAATAACCGTATTCTGCAg                                                                                                                                                                                                                                    <  1:403541‑M2/19‑1 (MQ=255)
        acgcacggtccatgatcacgtacatggagtagtttaagtaggcgttttccgtaaattcatgtagcgcaaggcgctctgccatatcgctcattaattctgattcctcaacttattcgccagccttcaaaggggcaatattgccggagatactacattatctgccgcgtcgagtcacaggagatcaacgtgaatagcgccgccggacgtgctccgacggcgataaaggcatcaGAATAACCGTATTCTGCAGc                                                                                                                                                                                                                                   <  1:504580‑M2/20‑1 (MQ=255)
         cgcacggtccatgatcacgtacatggagtagtttaagtaggcgttttccgtcaattcatgtagcgcaaggcgctctgccatatcgctcattaattctgattcctcaacttattcgccagccttcaaaggggcaatattgccggcgatactacattatctgccgcgtcgagtcacaggagatcaacgtgaatagcgccgccggacgtgctccgacggcgataaaggcatcaGAATAACCGTATTCTGCAGcc                                                                                                                                                                                                                                  >  2:519041‑M2/231‑251 (MQ=255)
             cggtccatgatcacgtacatggagtagtttaagtaggcgttttccgtaaattcatgtagcgcaaggcgctctgccatatcgctcattaattctgattcctcaacttattcgccagccttcaaaggggcaatattgccggagatactacattatctgccgcgtcgagtcacaggagatcaacgtgaatagcgccgccggacgtgctccgacggcgataaaggcatcaGAATAACCGTATTCTGCAGCCAggg                                                                                                                                                                                                                              <  1:303118‑M2/25‑1 (MQ=255)
              ggtccatgatcacgtacatggagtagtttaagtaggcgttttccgtaaattcatgtagcgcaaggcgctctgccatatcgctcattaattctgattcctcaacttattcgccagccttcaaaggggcaatattgccggagatactacattatctgccgcgtcgagtcacaggagatcaacgtgaatagcgccgccggacgtgctccgacggcgataaaggcatcaGAATAACCGCATTCTGCAGCCAggg                                                                                                                                                                                                                              >  2:287795‑M2/226‑250 (MQ=255)
                  catgatcacgtacatggagtagtttaagtaggcgttttccgtaaattcatgtagcgcaaggcgctctgccatatcgctcattaattctgattcctcaacttattcgccagccttcaaaggggcaatattgccggagatactacattatctgccgcgtcgagtcacaggagatcaacgtgaatagcgccgccggacgtgctccgacggcgataaaggcatcaGAATAACCGTATTCTGCAGCCAGGGAttt                                                                                                                                                                                                                          >  2:26012‑M2/222‑250 (MQ=255)
                         acgtacatggagtagtttaagaaggcgttttccggaaatgcaggtagcgcaaggcgctctgccatatcgctcattaattctgattcctcaacttattcgccagccttcaaaggggcaatattgccggagatactacattatctgccgcgtcgagtcacaggagaacaacgtgaatagcgccgccggacgtgctccgacggcgataaaggcatcaGAATAACCGTATTCTGCAGCCAGTGATTTAATccttc                                                                                                                                                                                                                  <  2:370078‑M2/37‑3 (MQ=255)
                             acatggagtagtttaagtaggcgttttccgtaaattcatgtagcgcaaggcgctctgccatatcgctcattaattctgattcctcaacttattcgccagccttcaaaggggcaatattgccggagatactacattatctgccgcgtcgagtcacaggagatcaacgtgaatagcgccgccggacgtgctccgacggcgataaaggcatcaGAATAACCGTATTCTGCAGCCAGGGATTTAATCCTGCCTTg                                                                                                                                                                                                              <  1:320697‑M2/41‑1 (MQ=255)
                             acatggagtagtttaagtaggcgttttccgtaaattcatgtagcgcaaggcgctctgccatatcgctcattaattctgattcctcaacttattcgccagccttcaaaggggcaatattgccggagatactacattatctgccgcgtcgagtcacaggagatcaacgtgaatagcgccgccggacgtgctccgacggcgataaaggcatcaGAATAACCGTATTCTGCAGCCAGGGATTTAATCCTGCCTTg                                                                                                                                                                                                              >  1:257015‑M2/211‑251 (MQ=255)
                                 ggatgagtttaagtagtcgtgttccgtaattccatgtagggcaaggcgctctgccatatggctcattaagtctgatgccgctacttattcgccagccttcaaaggggcaatattgccggagatacaacattatcggccgcgtcgagtcacaggagatcaacgtgaatagctccgccgggtgtgctccgacggcgataaaggcatcaTAATATCCGTATTCTGCATCCAGGGATTTAATCCTGCCTTGttt                                                                                                                                                                                                           <  2:44302‑M2/44‑1 (MQ=255)
                                            agtaggctttttccgtaaattcatgtagcgcaaggcgctctgccatatcgctcattaattctgattcctcaacttattcgccagccttcaaaggggcaatattgccggagatactacattatctgccgcgtcgagtcacaggagatcaacgtgaatagcgccgccggacgtgctccgacggcgataaaggcatcaGAATAACCGTATTCTGCAGCCAGGGATTTAATCCTGCCTTGTTTGCCCGGCCATcc                                                                                                                                                                                               <  1:356591‑M2/56‑1 (MQ=255)
                                             gtaggcgttttccgtaaattcatgtagcgcaaggcgctctgccatatcgctcattaattctgattcctcaacttattcgccagccttcaaaggggcaatattgccggagatactacattatctgccgcgtcgagtcacaggagatcaacgtgaatagcgccgccggacgtgctccgacggcgataaaggcatcaGAATAACCGTATTCTGCAGCCAGGGATTTAATCCTGCCTTGTTTGCCCGGCCATCCt                                                                                                                                                                                              >  1:319967‑M2/195‑251 (MQ=255)
                                             gtaggcgttttccgtaaattcatgtagcgcaaggcgctctgccatatcgctcattaattctgattccacaacttattcgccagccttcaaaggggcaatattgccggagatactacattatctgccgcgtcgagtcacaggagatcaacgtgaatagcgccgccggacgtgctccgacggcgataaaggcatcaGAATAACCGTATTCTGCAGCCAGGGATTTAATCCTGCCTTGTTTGCCCGGCCATCCt                                                                                                                                                                                              <  2:225168‑M2/57‑1 (MQ=255)
                                                      ttccgtaaattcatgtagcgcaaggcgctctgccatatcgctcattaattctgattcctcaacttattcgccagccttcaaaggggcaatattgccggagatactacattatctgccgcgtcgagtcacaggagatcaacgtgaatagcgccgccggacgtgctccgacggcgataaaggcatcaGAATAACCGTATTCTGCAGCCAGGGATTTAATCCTGCCTTGTTTGCCCGGCCATCCTGACCGGGc                                                                                                                                                                                      <  2:465722‑M2/65‑1 (MQ=255)
                                                                 catgtagcgcaaggcgctctgccatatcgctcattaattctgattcctcaacttattcgccagccttcaaaggggcaatattgccggagatactacattatctgccgcgtcgagtcacaggagatcaacgtgaatagcgccgccggacgtgctccgacggcgataaaggcatcaGAATAACCGTATTCTGCAGCCAGGGATTTAATCCTGCCTTGTTTGCCCGGCCATCCTGACCGGGCAATGTTCTTTcc                                                                                                                                                                          <  1:637858‑M2/77‑1 (MQ=255)
                                                                 catgtagcgcaaggcgctctgccatatcgctcattaattctgattcctcaacttattcgccagccttcaaaggggcaatagtgccggagatactacattatctgccgcgtcgagtcacaggagatcaacgtgaatagcgccgccggacgtgctccgacggcgataaaggcatcaGAATAACCGTATTCTGCAGCCAGGGATTTAATCCTGCCTTGTTTGCCCGGCCATCCTGACCGGGCAATGTTCTTTc                                                                                                                                                                           <  1:239108‑M2/76‑1 (MQ=255)
                                                                         gcaaggcgctctgccatatcgctcattaattctgattcctcaacttattcgccagccttcaaaggggcaatattgccggagatactacattatctgccgcgtcgagtcacaggagatcaacgtgaatagcgccgccggacgtgctccgacggcgataaaggcatcaGAATAACCGTATTCTGCAGCCAGGGATTTAATCCTGCCTTGTTTGCCCGGCCATCCTGACCGGGCAATGTTCTTTCCTTTAAAcc                                                                                                                                                                  >  2:126201‑M2/167‑251 (MQ=255)
                                                                         gcaaggcgctctgccatatcgcacattaattctgattcctcaacttattcgccagccttcaaaggggcaatattgccggaaatactacattatctgccgcgtcgagtcacaggagatcaacgtgaatagcgccgccggacgtgctccgacggcgataaaggcatcaGAATAACCGTATTCTGCAGCCAGGGATTTAATCCTGCCTTGTTTGCCCGGCCATCCTGACCGGGCAATGTTCTTTCCTTTAAAcc                                                                                                                                                                  >  1:199609‑M2/167‑251 (MQ=255)
                                                                          caaggcgctctgccatatcgctcattaattctgattcctcaacttattcgccagccttcaaaggggcaatattgccggagatactacattatctgccgcgtcgagtcacaggagatcaacgtgaatagcgccgccggacgtgctccgacggcgataaaggcatcaGAATAACCGTATTCTGCAGCCAGGGATTTAATCCTGCCTTGTTTGCCCGGCCATCCTGACCGGGCAATGTTCTTTCCTTTAAACCt                                                                                                                                                                 >  2:370887‑M2/166‑251 (MQ=255)
                                                                          caaggcgctctgccatatcgctcattaattctgattcctcaacttattcgccagccttcaaaggggcaatattgccggagatactacattatctgccgcgtcgagtcacaggagatcaacgtgaatagcgccgccggacgtgctccgacggcgataaaggcatcaGAATAACCGTATTCTGCAGCCAGGGATTTAATCCTGCCTTGTTTGCCCGGCCATCCTGACCGGGCAATGTTCTTTCCTTTAAACCt                                                                                                                                                                 >  2:89996‑M2/166‑251 (MQ=255)
                                                                          caaggcgctctgccatatcgctcattaattctgattcctcaacttattcgccagccttcaaaggggcaatattgccggagatactacattatctgccgcgtcgagtcacaggagatcaacgtgaatagcgccgccggacgtgctccgacggcgataaaggcatcaGAATAACCGTATTCTGCAGCCAGGGATTTAATCCTGCCTTGTTTGCCCGGCCATCCTGACCGGGCAATGTTCTTTCCTTTAAACCt                                                                                                                                                                 >  2:541295‑M2/166‑251 (MQ=255)
                                                                           aaggcgctctgccatatcgctcattaattctgattcctcaacttattcgccagccttcaaaggggcaatattgccggagatactacattatctgccgcgtcgagtcacaggagatcaacgtgaatagcgccgccggacgtgctccgacggcgataaaggcatcaGAATAACCGTATTCTGCAGCCAGGGATTTAATCCTGCCTTGTTTGCCCGGCCATCCTGACCGGGCAATGTTCTTTCCTTTAAACCTc                                                                                                                                                                <  1:547433‑M2/87‑1 (MQ=255)
                                                                               cgctctgccatatcgctcattaattctgattcctcaacttattcgccagccttcaaaggggcaatattgccggagatactacattatctgccgcgtcgagtcacaggagatcaacgtgaatagcgccgccggacgtgctccgacggcgataaaggcatcaGAATAACCGTATTCTGCAGCCAGGGATTTAATCCTGCCTTGTTTGCCCGGCCATCCTGACCGGGCAATGTTCTTTCCTTTAAACCTCAAtc                                                                                                                                                            <  2:605161‑M2/91‑1 (MQ=255)
                                                                                gctctgccatatcgctcattaattctgattcctcaacttattcgccagccttcaaaggggcaatattgccggagatactacattatctgccgcgtcgagtcacaggagatcaacgtgaatagcgccgccggacgtgctccgacggcgataaaggcatcaGAATAACCGTATTCTGCAGCCAGGGATTTAATCCTGCCTTGTTTGCCCGGCCATCCTGACCGGGCAATGTTCTTTCCTTTAAACCTCAAtct                                                                                                                                                           <  1:511354‑M2/92‑1 (MQ=255)
                                                                                     gccatatcgctcattaattctgattcctcaacttattcgccagccttcaaaggggcaatattgccggagatactacattatctgccgcgtcgagtcacaggagatcaacgtgaatagcgccgccggacgtgctccgacggcgataaaggcatcaGAATAACCGTATTCTGCAGCCAGGGATTTAATCCTGCCTTGTTTGCCCGGCCATCCTGACCGGGCAATGTTCTTTCCTTTAAACCTCAATCTcc                                                                                                                                                         <  2:83343‑M2/94‑1 (MQ=255)
                                                                                        atatcgcttattaattctgattcctcaactgattcgccagacttcaaaggggcaatattgccggggatcctacattatctgccgcgtcgagtcacaggagatcaaagtgaatagcgccgccggacgttctcagacggcgataaaggcatcaGAATAACCGTATTCTGCAGCCAGGGATTTAATCCTGCCTTGTTTGCCCGGCCATCCTGACCGGGCAATGGTCTTTCATTAAAACCTCAATCTCCGCCATg                                                                                                                                                   <  2:467955‑M2/100‑1 (MQ=255)
                                                                                           tcgctcattcgttctgattcctcaacttattcgccagccttcaaaggggcaatattgccggagatactacattatctgccgcgtcgagtcacaggagatcaacgtgaatagcgccgccggacgtgctccgacggcgataaaggcatcaGAATAACCGTATTCTGCAGCCAGGGATTTAATCCTGCCTTGTTTGCCCGGCCATCCTGACCGGGCAATGTTCTTTCCTTTAAACCTCAATCTCCGCCATGTCg                                                                                                                                                >  2:185754‑M2/149‑251 (MQ=255)
                                                                                            cgctcattaattctgattcctcaacttattcgccagccttcaaaggggcaatattgccggagatactacattatctgccgcgtcgagtcacaggagatcaacgtgaatagcgccgccggacgtgctccgacggcgataaaggcatcaGAATAACCATATTCTGCAGCCAGGGATTTAATCCTGCCTTGTTTGCCCGGCCATCCTGACCGGGCAATGTTCTTTCCTTTAAACCTCAATCTCCGCCATGTCGc                                                                                                                                               >  1:216640‑M2/148‑251 (MQ=255)
                                                                                                              cctcaacttattcgccagccttcaaaggggcaatattgccggagatactacattatctgccgcgtcgagtcacaggagatcaacgtgaatagcgccgccggacgtgctccgacggcgataaaggcatcaGAATAACCGTATTCTGCAGCCAGGGATTTAATCCTGCCTTGTTTGCCCGGCCATCCTGACCGGGCAATGTTCTTTCCTTTAAACCTCAATCTCCGCCATGTCGCCTTTCTCTTGCAACCAGt                                                                                                                             <  2:40040‑M2/122‑1 (MQ=255)
                                                                                                                tcaacttattcgccagccttcaaaggggcaatattgccggagatactacattatctgccgcgtcgagtcacaggagatcaacgtgaatagcgccgccggacgtgctccgacggcgataaaggcatcaGAATAACCGTATTCTGCAGCCAGGGATTTAATCCTGCCTTGTTTGCCCGGCCATCCTGACCGGGCAATGTTCTTTCCTTTAAACCTCAATCTCCGCCATGTCGCCTTTCTCTTGCAACCAGTTg                                                                                                                           >  1:51815‑M2/128‑251 (MQ=255)
                                                                                                                ccaacttattcgccagccttcaaaggggaaatattgccggagatactacattatctgccgcgtcgagtcacaggagatcaacgtgaatagcgccgccggacgtgctccgacggcgataaaggcatcaGAATAACCGTATTCTGCAGCCAGGGATTTAATCCTGCCTTGTTTGCCCGGCCATCCTGACCGGGCAATGTTCTTTCCTTTAAACCTCAATCTCCGCCATGTCGCCTTTCTCTTGCAACCAGTTg                                                                                                                           <  2:113525‑M2/124‑1 (MQ=255)
                                                                                                                  aacttattcgccagccttcaaaggggcaatattgccggagatactacattatctgccgcgtcgagtcacaggagatcaacgtgaatagcgccgccggacgtgctccgacggcgataaaggcatcaGAATAACCGTATTCTGCAGCCAGGGATTTAATCCTGCCTTGTTTGCCCGGCCATCCTGACCGGGCAATGTTCTTTCCTTTAAACCTCAATCTCCACCATGTCGCCTTTCTCTTGCAACCAGTTgcg                                                                                                                         >  1:418421‑M2/126‑251 (MQ=255)
                                                                                                                    cttattcgccagccttcaaaggggcaatattgccggagatactacattatctgccgcgtcgagtcacaggagatcaccgtgaatagcgccgccggacgtgctccgacggcgataaaggcatcaGAATAACCGTATTCTGCAGCCAGGGATTTAATCCTGCCTTGTTTGCCCGGCCATCCTGACCGGGCAATGTTCTTTCCTTTACACCTCTATCTCCGCCATGTCGCCTTTCTCTTGCAACCAGTTgcggc                                                                                                                       >  2:273165‑M2/124‑251 (MQ=255)
                                                                                                                         tcgccagccttcaaaggggcaatattgccggagatactacattatctgccgcgtcgagtcacaggagatcaacgtgaatagcgccgccggacgtgctccgacggcgataaaggcatcaGAATAACCGTATTCTGCAGCCAGGGATTTAATCCTGCCTTGTTTGCCCGGCCATCCTGACCGGGCAATGTTCTTTCCTTTAAACCTCAATCTCCGCCATGTCGCCTTTCTCTTGCAACCAGTTGCGGCGATc                                                                                                                   <  1:129584‑M2/132‑1 (MQ=255)
                                                                                                                           gccagccttcaaaggggcaatattgccggagatactacattatctgccgcgtcgagtcacaggagatcaacgtgaatagcgccgccggacgtgctccgacggcgataaaggcatcaGAATAACCGTATTCTGCAGCCAGGGATTTAATCCTGCCTTGTTTGCCCGGCCATCCTGACCGGGCAATGTTCTTTCCTTTAAACCTCAATCTCCGCCATGTCGCCTTTCTCTTGCAACCAGTTGCGGCGATCTTc                                                                                                                >  1:396263‑M2/117‑251 (MQ=255)
                                                                                                                            ccagccttcaaaggggcaatattgccggagatactacattatctgccgcgtcgagtcacaggagatcaacgtgaatagcgccgccggacgtgctccgacggcgataaaggcatcaGAATAACCGTATTCTGCAGCCAGGGATTTAATCCTGCCTTGTTTGCCCGGCCATCCTGACCGGGCAATGTTCTTTCCTTTAAACCTCACTCTCCGCCATGTCGCCTTTCTCTTGCAACCAGTTGCGGCGATCTTcc                                                                                                               >  2:388868‑M2/116‑251 (MQ=255)
                                                                                                                                ccttcaaaggggcaatattgccggagatactacattatctgccgcgtcgagtcacaggagatcaacgtgaatagcgccgccggacgtgctccgacggcgataaaggcatcaGAATAACCGTATTCTGCAGCCAGGGATTTAATCCTGCCTTGTTTGCCCGGCCATCCTGACCGGGCAATGTTCTTTCCTTTAAACCTCAATCTCCGCCATGTCGCCTTTCTCTTGCAACCAGTTGCGGCGATCTTCCGAGc                                                                                                           <  1:375629‑M2/140‑1 (MQ=255)
                                                                                                                                     aaaggggcaatattgccggagatactacattatctgccgcgtcgagtcacaggagatcaacgtgaatagcgccgccggacgtgctccgacggcgataaaggcatcaGAATAACCGTATTCTGCAGCCAGGGATTTAATCCTGCCTTGTTTGCCCGGCCATCCTGACCGGGCAATGTTCTTTCCTTTAAACCTCAATCTCCGCCATGTCGCCTTTCTCTTGCAACCAGTTGCGGCGATCTTCCGAGCGTttc                                                                                                      <  1:325879‑M2/145‑1 (MQ=255)
                                                                                                                                          ggcaatattgccggagatactacattatctgccgcgtcgagtcacaggagatcaacgtgaatagcgccgccggacgtgctccgacggcgataaaggcatcaGAATAACCGTATTCTGCAGCCAGGGATTTAATCCTGCCTTGTTTGCCCGGCCATCCTGACCGGGCAATGTTCTTTCCTTTAAACCTCAATCTCCGCCATGTCGCCTTTCTCTTGCAACCAGTTGCGGCGATCTTCCGAGCGTTTCTTCGc                                                                                                 <  2:343443‑M2/150‑1 (MQ=255)
                                                                                                                                              atattgccggagatactacattatctgccgcgtcgagtcacaggagatcaacgtgaatagcgccgccggacgtgctccgacggcgataaaggcatcaGAATAACCGTATTCTGCAGCCAGGGATTTAATCCTGCCTTGTTTGCCCGGCCATCCTGACCGGGCAATGTTCTTTCCTTTAAACCTCAATCTCCGCCATGTCGCCTTTCTCTTGCAACCAGTTGCGGCGATCTTCCGAGCGTTTCTTCGCcagc                                                                                             >  2:384840‑M2/98‑251 (MQ=255)
                                                                                                                                                      ggagatactacattatctgccgcgtcgagtcacaggagatcaacgtgaatagcgccgccggacgtgctccgacggcgataaaggcatcaGAATAACCGTATTCTGCAGCCAGGGATTTAATCCTGCCTTGTTTGCCCGGCCATCCTGACCGGGCAATGTTCTTTCCTTTAAACCTCAATCTCCGCCATGTCGCCTTTCTCTTGCAACCAGTTGCGGCGATCTTCCGAGCGTTTCTTCGCCAGCAGCATATc                                                                                     >  1:230200‑M2/90‑251 (MQ=255)
                                                                                                                                                       gagatactacattatctgccgcgtcgagtcacaggagatcaacgtgaatagcgccgccggacgtgctccgacggcgataaaggcatcaGAATAACCGTATTCTGCAGCCAGGGATTTAATCCTGCCTTGTTTGCCCGGCCATCCTGACCGGGCAATGTTCTTTCCTTTAAACCTCAATCTCCGCCATGTCGCCTTTCTCTTGCAACCAGTTGCGGCGATCTTCCGAGCGTTTCTTCGCCAGCAGCATATcc                                                                                    <  1:286402‑M2/163‑1 (MQ=255)
                                                                                                                                                          atactccattatctgccgcgtcgagtcacaggagatcaacgtgaatagcgccgccggacgtgctccgacggcgataaaggcatcaGAATAACCGTATTCTGCAGCCAGGGATTTAATCCTGCCTTGTTTGCCCGGCCATCCTGACCGGGCAATGTTCTTTCCTTTAAACCTCAATCTCCGCCATGTCGCCTTTCTCTTGCAACCAGTTGCGGCGATCTTCCGAGCGTTTCTTCGCCAGCAGCATATCcatc                                                                                 >  1:316503‑M2/86‑251 (MQ=255)
                                                                                                                                                               gcatcatctgccgcgtcgagtcacaggagatcaacgtgaatagcgccgccggacgtgctccgacggcgataaaggcatcaGAATAACCGTATTCTGCAGCCAGGGATTTAATCCTGCCTTGTTTGCCCGGCCCTCCTGACCGGGCAATGTTCTTTCCTTTAAACCTCAATCTCCGCCATGTCGCCTTTCTCTTGCAACCAGTTGCGGCGATCTTCCGAGCGTTTATTCGCCAGCAGCATCTCCATCATcg                                                                             <  2:342112‑M2/170‑1 (MQ=255)
                                                                                                                                                               acattatctgccgcgtcgagtcacaggagatcaacgtgaatagcgccgccggacgtgctccgacggcgataaaggcatcaGAATAACCGTATTCTGCAGCCAGGGATTTAATCCTGCCTTGTTTGCCCGGCCATCCTGACCGGGCAATGTTCTTTCCTTTAAACCTCAATCTCCGCCATGTCGCCTTTCTCTTGCAACCAGTTGCGGCGATCTTCCGAGCGTTTCTTCGCCAGCAGCATATCCatcatc                                                                              <  2:53640‑M2/169‑1 (MQ=255)
                                                                                                                                                                    atctgccgcgtcgagtcacaggagatcaacgtgaatagcgccgccggacgtgctccgacggcgataaaggcatcaGCATAACCGTATTCTGCAGCCAGGGATTTAATCCTGCCTTGTTTGCCCGGCCATCCTGACCGGGCAATGTTCTTTCCTTTAAACCCCAATCTCCGCCATGTCGCCTTTCTCTTGCAACCAGTTGCGGCGATCTTCCGAGCGTTTCTTCGCCAGCAGCATATCCATCATCGCGTc                                                                         >  1:30263‑M2/76‑249 (MQ=255)
                                                                                                                                                                    atctgccgcgtcgagtcacaggagatcaacgtgaatagcgccgccggacgtgctccgacggcgataaaggcatcaGAATAACCGTATTCTGCAGCCAGGGATTTAATCCTGCCTTGTTTGCCCGGCCATCCTGACCGGGCAATGTTCTTTCCTTTAAACCTCAATCTCCGCCATGTCGCCTTTCTCTTGCAACCAGTTGCGGCGATCTTCCGAGCGTTTCTTCGCCAGCAGCATATCCATCATCGCGTc                                                                         >  1:563764‑M2/76‑249 (MQ=255)
                                                                                                                                                                      ctgccgcgtcgagtcacaggagatcaacgtgaatagcgccgccggacgtgctccgacgacgataaaggcatcaGAATAACCGTATTCTGCAGCCAGGGATTTAATCCTGCCTTGTTTGCCCGGCCATCCTGACCGGGCAATGTTCTTTCCTTTAAACCTCAATCTCCGCCATGTCGCCTTTCTCTTGCAACCAGTTGCGGCGATCTTCCGAGCGTTTCTTCGCCAGCAGCATATCCATCATCGCGTCAGt                                                                      >  1:549687‑M2/74‑250 (MQ=255)
                                                                                                                                                                       ggccgcgtcgagtcacaggagatcaaagtgaatagcgccgccggacgtgctccgacggcgataaaggcatcaGAATAACCGTATTCTGCAGCCAGGGATTTAATCCTGCCTTGTTTGCCCGGCCATCCTGACCGGGCAATGTTCTTTCCTTTAGACCTCAATCTCCGCCATGTCGCCTTACTCGTGCAACCAGTTGCGGCGATCTTCCGGGCGTTTCTTCGCCATCAGCATAGCCATCATCGCGTCAGTAc                                                                    >  2:435076‑M2/73‑251 (MQ=255)
                                                                                                                                                                         ccgcgtcgagtcacaggagatcaacgtgaatagcgccgccggacgtgctccgacggcgataaaggcatcaGAATAACCGTATTCTGCAGCCAGGGATTTAATCCTGCCTTGTTTGCCCGGCCATCCTGACCGGGCAATGTTCTTTCCTTTAAACCTCAATCTCCGCCATGTCGCCTTTCTCTTGCAACCAGTTGCGGCGATCTTCCGAGCGTTTCTTCGCCAGCAGCATATCCATCATCGCGTCAGTACGc                                                                  <  1:311928‑M2/181‑1 (MQ=255)
                                                                                                                                                                           gcgtcgagtcacaggagatcaacgtgaatagcgccgccggacgtgctccgacggcgataaaggcatcaGAATAACCGTATTCTGCAGCCAGGGATTTAATCCTGCCTTGTTTGCCCGGCCATCCTGACCGGGCAATGTTCTTTCCTTTAAACCTCAATCTCCGCCATGTCGCCTTTCTCTTGCAACCAGTTGCGGCGATCTTCCGAGCGTTTCTTCGCCAGCAGCATATCCATCATCGCGTCAGTACGCTg                                                                <  1:367208‑M2/183‑1 (MQ=255)
                                                                                                                                                                              tcgagtcacaggagatcaacgtgaatagcgccgccggaagtgccccgacggcgataaaggcatcaGAATAACCGTATTCTGCAGCCAGGGATTTAATCCTGCCTTGTTTGCCCGGCCATCCTGACCGGGCAATGTTCTTTCCTTTAAACCTAAATCTCCGCCATGTCGCCTTTCTCTTGCAACCAGTTGCGGCGATCTTCCGCGCGTTTCTTCGCCAGCAGCATATCCATCATCGCGTAAGTACGCTGATc                                                             >  2:163962‑M2/66‑251 (MQ=255)
                                                                                                                                                                               cgagtcacaggagatcaacgtgaatagcgcagccggacgtgctccgacggcgataaaggcatcaGAATAACCGTATTCTGCATCCAGTGATTTATGCCTGCCTTGTTTGCCCGGCCATCCTGACCGGGCAATGTTCTTTCCTTTAAACCTCAATCTCCGCCATGTCGCCTTTCTCTTGCTACCAGTTGCGGCGATCTTCCGAGCGTTTCTTCGCCAGCAGCATATCCATCATCGCGTCAGTACGCTGATc                                                             <  2:271746‑M2/186‑1 (MQ=255)
                                                                                                                                                                               cgaggcacaggagatcaacgtgaatagcgccgccggacgtgctccgacggcgataaaggcatcaGAATAACCGTATTCTGCAGCCAGGGATTTAATCCTGCCTTGTTTGCCCTGCCATCCTGACCGGGCAATGTTCTTTCCTTTAAACCTCAATCTCCGCCATGTCTCCTTTCTCTTGCAACCGGTTGCGGCGATCTTCCGAGCGTTTCTTCGCCAGCAGCATATCCATCATCGCGTCAGTACGCTGATc                                                             <  1:468557‑M2/186‑1 (MQ=255)
                                                                                                                                                                                gagtcacaggagatcaacgtgaatagcgccgccggacgtgctccgacggcgataaaggcatcaGAATAACCGTATTCTGCAGCCAGGGATTTAATCCTGCCTTGTTTGCCCGGCCATCCTGACCGGGCAATGTTCTTTCCTTTAAACCTCAATCTCCGCCATGTCGCCTTTCTCTTGCAACCAGTTGCGGCGATCTTCCGAGCGTTTCTTCGCCAGCAGCATATCCATCATCGCGTCAGTACGCTGATCGt                                                           >  1:207402‑M2/64‑251 (MQ=255)
                                                                                                                                                                                 agtcacaggagatcaacgtgaatagcgccgccggacgtgctccgacggcgataaaggcatcaGAATAACCGTATTCTGCAGCCAGGGATTTAATCCTGCCTTGTTTGCCCGGCCATCCTGACCGGGCAATGTTCTTTCCTTTAAACCTCAATCTCCGCCATGTCGCCTTTCTCTTGCAACCAGTTGCGGCGATCTTCCGAGCGTTTCTTCGCCAGCAGCATATCCATCATCGCGTCAGTACGCTGATCGTc                                                          >  2:511415‑M2/63‑251 (MQ=255)
                                                                                                                                                                                      caggagatcaacgtgaatagcgccgccggacgtgctccgacggcgataaaggcatcaGAATAACCGTATTCTGCAGCCAGGGATTTAATCCTGCCTTGTTTGCCCGGCCATCCTGACCGGGCAATGTTCTTTCCTTTAAACCTCAATCTCCGCCATGTCGCCTTTCTCTTGCAACCAGTTGCGGCAATCTTCCGAGCGTTTCTTCGCCAGCAGCATATCCATCATCGCGTCAGTACGCTGATCGTCTtcat                                                     >  2:149092‑M2/58‑251 (MQ=255)
                                                                                                                                                                                        ggagatcaacgtgaatagcgccgccggccgggctccgacggcgataaaggcatctGAATAACCGTATTCTGCAGCCAGGGATTTAATCCTGCCTTGTTTGCCCGGCCATCCTGCCCGGGCAATGTTCTTTCCTTTAAACCTCAATCTCCGCCATGTCGCCTTTCTCTTGCAACCAGTTGCGGCGATCTTCCGAGCGTTTCTTCGCCAGCAGCATATCCATCATCGCGTCAGTACGCTGATCGTCTtcatca                                                   <  2:77106‑M2/196‑1 (MQ=255)
                                                                                                                                                                                                 cgtgaatagcgccgccggacgtgctccgacggcgataaaggcatcaGAATAACCGTATTCTGCAGCCAGGGATTTAATCCTGCCTTGTTTGCCCGGCCATCCTGACCGGGCAATGTTCTTTCCTTTAAACCTCAATCTCCGCCATGTCGCCTTTCTCTTGCAACCAGTTGCGGCGATCTTCCGAGCGTTTCTTCGCCAGCAGCATATCCATCATCGCGTCAGTACGCTGATCTTCTTCATCATCTATAGTc                                          <  1:149119‑M2/205‑1 (MQ=255)
                                                                                                                                                                                                    gaatagcgccgccggacgtgctccgacggcgataaaggcatcaGAATAACCGTATTCTGCAGCCAGGGATTTAATCCTGCCTTGTTTGCCCGGCCATCCTGACCGGGCAATGTTCTTTCCTTTAAACCTCAATCTCCGCCATGTCGCCTTTCTCTTGCAACCAGTTGCGGCGATCTTCCGAGCGTTTCTTCGCCAGCAGCATATCCATCATCGCGTCAGTACGCTGATCGTCTTCATCATCGATAGTCaa                                        <  1:459510‑M2/207‑1 (MQ=255)
                                                                                                                                                                                                    gaatagcgccgccggacgtgctccgacggcgataaaggcatcaGAATAACCGTATTCTGCAGCCAGGGATTTAATCCTGCCTTGTTTGCCCGGCCATCCTGACCGGGCAATGTTCTTTCCTTTAAACCTCAATCTCCGCCATGTCGCCTTTCTCTTGCAACCAGTTGCGGCGATCTTCCGAGCGTTTCTTCGCCAGCAGCATATCCATCATCGCGTCAGTACGCTGATCGTCTTCATCATCGATAGTCAAc                                       <  1:236460‑M2/208‑1 (MQ=255)
                                                                                                                                                                                                     aatagcgccgccggacgtgctccgacggcgataaaggcatcaGAATAACCGTATTCTGCAGCCAGGGATTTAATCCTGCCTTGTTTGCCCGGCCATCCTGACCGGGCAATGTTCTTTCCTTTAAACCTCAATCTCCGCCATGTCGCCTTTCTCTTGCAACCAGTTGCGGCGATCTTCCGAGCGTTTCTTCGCCAGCAGCATATCCATCATCGCGTCAGTACGCTGATCGTCTTCATCATCGATAGTCAACt                                      <  1:126201‑M2/209‑1 (MQ=255)
                                                                                                                                                                                                           gccgccggacgtgctccgacggcgataaaggcatcaGAATAACCGTATTCTGCAGCCAGGGATTTAATCCTGCCTTGTTTGCCCGGCCATCCTGACCGGGCAATGTTCTTTCCTTTAAACCTCAATCTCCGCCATGTCGCCTTTCTCTTGCAACCAGTTGCGGCGATCTTCCGAGCGTTTCTTCGCCAGCAGCATATCCATCATCGCGTCAGTACGCTGATCGTCTTCATCATCGATAGTCAACTGCAcc                                 >  1:397741‑M2/37‑250 (MQ=255)
                                                                                                                                                                                                                cggacgtgctccgacggcgataaaggcatcaGAATAACCGTATTCTGCAGCCAGGGATTTAATCCTGCCTTGTTTGCCCGGCCATCCTGACCGGGCAATGTTCTTTCCTTTAAACCTCAATCTCCGCCATGTCGCCTTTCTCTTGCAACCAGTTGCGGCGATCTTCCGAGCGTTTCTTCGCCAGCAGCATATCCATCATCGCGTCAGTACGCTGATCGTCTTCATCATCGATAGTCAACTGCACCCGACgg                           >  2:431132‑M2/32‑251 (MQ=255)
                                                                                                                                                                                                                 ggacgtgctccgacggcgctaaaggcatcaGAATAACCGTATTCTGCAGCCAGGGATTTAATCCTGCCTTGTTTGCCCGGCAATCCTGACCGGGCAATGTTCTTTCCTTTAAACCTCAATCTCCGCCATGTCGCCTTTCTCTTGCAACCAGTTGCGGCGATCTTCCGAGCGTTTCTTCGCCAGCAGCATATCCATCATCGCGTCAGAACGCAGATCGTCTTCATCATCGATAGTCACCTGCACCAGACGGc                          >  2:322566‑M2/31‑251 (MQ=255)
                                                                                                                                                                                                                   acgtgctccgacggcgataaaggcatcaGAATAACCGTATTCTGCAGCCAGGGATTTAATCCTGCCTTGTTTGCCCGGCCATCCTGACCGGGCAATGTTCTTTCCTTTAAACCTCAATCTCCGCCATGTCGCCTTTCTCTTGCAACCAGTTGCGGCGATCTTCCGAGCGTTTCTTCGCCAGCAGCATATCCATCATCGCGTCAGTACGCTGATCGTCTTCATCATCGATAGTCAACTGCACCAGACGGCGa                        <  1:613911‑M2/223‑1 (MQ=255)
                                                                                                                                                                                                                               ggcgataaaggcatcaGAATAACCGTATTCTGCAGCCAGGGATTTAATCCTGCCTTGTTTGCCCGGCCATCCTGACCGGGCAATGTTCTTTCCTTTAAACCTCAATCTCCGCCATGTCGCCTTTCTCTTGCAACCAGTTGCGGCGATCTTCCGAGCGTTTCTTCGCCAGCAGCATATCCATCATCGCGTCAGTACGCTGATCGTCTTCATCATCGATAGTCAACTGCACCAGACGGCGAGTGTTCGGATc             >  2:604696‑M2/17‑250 (MQ=255)
                                                                                                                                                                                                                               ggcgataaaggcatcaGAATAACCGTATTCTGCAGCCAGGGATTTAATCCTGCCTTGTTTGCCCGGCCATCCTGACCGGGCAATGTTCTTTCCTTTAAACCTCAATCTCCGCCATGTCGCCTTTCTCTTGCAACCAGTTGCGGCGATCTTCCGAGCGTTTCTTCGCCAGCAGCATATCCATCATCGCGTCAGTACGCTGATCGTCTTCATCATCGATAGTCAACTGCACCAGACGGCGAGTGTTCGGATCa            <  1:231005‑M2/235‑1 (MQ=255)
                                                                                                                                                                                                                                gcgataaaggcatcaGAATAACCGTATTCTGCAGCCAGGGATTTAATCCTGCCTTGTTTGCCCGGCCATCCTTGCCGGGCAATGTTCTTTCCTTTAAACCTCAATCTCCGCCATGTCGCCTTTCTCTTGCAACCAGTTGCGGCGATCTTCCGAGCGTTTCTTCGCCAGCAGCATATCCATCATCGCGTCAGTACGCTGATCGTCTTCATCATCGATAGTCAACTGCACCAGACGGCGAGTGTTCGGATCaa           <  2:532082‑M2/236‑1 (MQ=255)
                                                                                                                                                                                                                                 cgataaaggcatcaGAATAACCGTATTCTGCAGCCAGGGATTTAATCCTGCCTTGTTTGCCCGGCCATCCTGACCGGGCAATGTTCTTTCCTTTAAACCTCAATCTCCGCCATGTCGCCTTTCTCTTGCAACCAGTTGCGGCGATCTTCCGAGCGTTTCTTCGCCAGCAGCATATCCATCATCGCGTCAGTACGCTGATCGTCTTCATCATCGATAGTCAACTGCACCAGACGGCGAGTGATCGGATCa            >  1:108204‑M2/15‑249 (MQ=255)
                                                                                                                                                                                                                                  gataaaggcatcaGAATAACCGTATTCTGCAGCCAGGGATTTAATCCTGCCTTGTTTGCCCGGCCATCCTGACCGGGCAATGTTCTTTCCTTTAAACCTCAATCTCCGCCATGTCGCCTTTCTCTTGCAACCAGTTGCGGCGATCTTCCGAGCGTTTCTTCGCCAGCAGCATATCCATCATCGCGTCAGTACGCTGATCGTCTTCATCATCGATAGTCAACTGCACCAGACGGCGAGTGTTCGGATCAAGc         <  1:573781‑M2/238‑1 (MQ=255)
                                                                                                                                                                                                                                        gtcatcaGAATAACCGTATTCTGCAGCCAGGGATTTAATCCTGCCTTGTTTGCCCGGCCATCCTGACCGGGCAATGTTCTTTCCTTTAAACCTCAATCTCGGCCATGTTGCCTTTCTCTTGCAACCAGTTGCGGCGATCTTCCGAGCGTTTCTTCGCCAGCAGCATATCCATCATCGCGTCAGTACGCTGATCGTCTTCATCATCGATAGTCAACTGCACCAGACGGTGAGCGTTCGGATCAAGCGTGGtt   >  1:336997‑M2/8‑251 (MQ=255)
                                                                                                                                                                                                                                        ggcatcaGAATAACCGTATTCTGCAGCCAGGGATTTAATCCTGCCTTGTTTGCCCGGCCATCCTGACCGGGCAATGTTCTTTCCTTTAAACCTCAATCTCCGCCATGTCGCCTTTCTCTTGCAACCAGTTGCGGCGATCTTCCGAGCGTTTCTTCGCCAGCAGCATATCCATCATCGCGTCAGTACGCTGATCGTCTTCATCATCGATAGTCAACTGCACCAGACGGCGAGTGTTCGGATCAAGCGTGGtt   >  1:337001‑M2/8‑251 (MQ=255)
                                                                                                                                                                                                                                         gcatcaGAATAACCGTATTCTGCAGCCAGGGATTTAATCCTGCCTTGTTTGCCCGGCCATCCTGACCGGGCAATGTTCTTTCCTTTAAACCTCAATCTCCGCCATGTCGCCTTTCTCTTGCATCCAGTTGCGGCGATCTTCCGAGCGTTTCTTCGCCAGCAGCATATCCATCATCGCGTCAGTACGCTGATCGTCTTCATCATCGATAGTCAACTGCACCAGACGGCGAGTGTTCGGATCAAGCGTGGttt  <  2:508511‑M2/245‑1 (MQ=255)
                                                                                                                                                                                                                                         gcatcaGAATAACCGTATTCTGCAGCCAGGGATTTAATCCTGCCTTGTTTGCCCGGCCATCCTGACCGGGCAATGTTCTTTCCTTTAAACCTCAATCTCCGCCATGTCGCCTTTCTCTTGCAACCAGTTGCGGCGATCTTCCGAGCGTTTCTTCGCCAGCAGCATATCCATCATCGCGTCAGTACGCTGATCGTCTTCATCATCGATAGTCAACTGCACCAGACGGCGAGTGTTCGGATCAAGCGTGGttt  <  1:33166‑M2/245‑1 (MQ=255)
                                                                                                                                                                                                                                         gcatcaGAATAACCGTATTCTGCAGCCAGGGATTTAATCCTGCCTTGTTTGCACGGCCATCCTGACCGGGCAATGTTCTTTCCTTTAAACCTCAATCTCGGCCATGTCGCCTTTCTCTTGCAACCAGTTGCGGCGATCTTCCGAGCGTTTCTTCGCCAGCAGCATATCCATCATCGCGTCAGTACGCTGATCGTCTTCATCATCGTTAGTCAACTGCACCAGTCGGCGAGTGTTCGGATCAAGCGTGGttt  >  2:113880‑M2/7‑251 (MQ=255)
                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                               
CACCGTCAGGCGGTATTGGATCAGTTTGCTAAAATCGTTCCAACCGTACTGAAAGAAGAAGGTTGCCACGGCTATGCGCCAATGGTGGATTGCGCAGCTGGCGTGAGTTTCCAGTCTATGGCACCGGATTCTATCGTGATGATTGAGCAGTGGGAAAGCATCGCGCATCTTGAAGCGCATCTGCAAACCCCGCACATGAAGGCGTATAGCGAAGCCGTAAAAGGTGACGTGCTGGAGATGAATATCCGTATTCTGCAGCCAGGGATTTAATCCTGCCTTGTTTGCCCGGCCATCCTGACCGGGCAATGTTCTTTCCTTTAAACCTCAATCTCCGCCATGTCGCCTTTCTCTTGCAACCAGTTGCGGCGATCTTCCGAGCGTTTCTTCGCCAGCAGCATATCCATCATCGCGTCAGTACGCTGATCGTCTTCATCATCGATAGTCAACTGCACCAGACGGCGAGTGTTCGGATCAAGCGTGGTTT  >  NC_000913/3173187‑3173670

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