Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,229,177 A→C 14.4% intergenic (+112/+261) yohP → / ← dusC uncharacterized protein/tRNA‑dihydrouridine synthase C

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,229,1770AC14.4% 488.5 / 13.2 269intergenic (+112/+261)yohP/dusCuncharacterized protein/tRNA‑dihydrouridine synthase C
Reads supporting (aligned to +/- strand):  ref base A (65/134);  new base C (27/11);  total (106/163)
Fisher's exact test for biased strand distribution p-value = 1.66e-05
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 273 total aligned reads displayed.

TACTGGTGGTGACTGGCGTATTTAAGATGATATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGCAGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGCCCTCAGCCCACATCCATG  >  NC_000913/2229028‑2229326
                                                                                                                                                     |                                                                                                                                                     
tACTGGTGGTGACTGGCGTATTTAAGATGATATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGCa                                                                                                                                                       >  2:2033549/1‑150 (MQ=255)
   tggtggTGACTGGCGTATTTAAGATGATATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGCCggg                                                                                                                                                    >  2:5719780/1‑150 (MQ=255)
   tggtggTGACTGGCGTATTTAAGATGATATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGCCggg                                                                                                                                                    >  1:8817569/1‑150 (MQ=255)
      tggtgACTGGCGTATTTAAGATGATATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCAGCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCggggggggg                                                                                                                                                 >  1:5851438/1‑149 (MQ=255)
      tggtgACTGGCGTATTTAAGATGATATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCAGCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGGCgggggg                                                                                                                                                 >  1:5961148/1‑149 (MQ=255)
        gtgACTGGCGTATTTAAGATGATATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGCAGGGggcgg                                                                                                                                               <  1:7038730/150‑1 (MQ=255)
        gtgACTGGCGTATTTAAGATGATATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGCAGGGggcgg                                                                                                                                               <  2:6766778/150‑1 (MQ=255)
        gtgACTGGCGTATTTAAGATGATATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGCAGGGggcgg                                                                                                                                               <  1:2388329/150‑1 (MQ=255)
        gtgACTGGCGTATTTAAGATGATATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGCAGGGggcgg                                                                                                                                               <  2:9144960/150‑1 (MQ=255)
        gtgACTGGCGTATTTAAGATGATATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGCAGGGggcgg                                                                                                                                               <  1:4265979/150‑1 (MQ=255)
        gtgACTGGCGTATTTAAGATGATATTTTAAAAGTAATTAATGTCATCAGGTCCGACAATAACGAGAATATTTCAGTCTCTCCTGCTGTTGCGCTCCTGTCATGTGCATTGCTTCAGATAAGCACTGGCGCAAGGAGCGCGCAGGGggcgg                                                                                                                                               >  1:5595649/1‑150 (MQ=255)
         tgACTGGCGTATTTAAGATGATATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGCAGGGggcggc                                                                                                                                              <  1:5334719/150‑1 (MQ=255)
          gACTGGCGTATTTAAGATGATATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGCAGGGGGCGGcc                                                                                                                                             <  1:3258536/150‑1 (MQ=255)
          gACTGGCGTATTTAAGATGATATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGCAGGGGGCGGcc                                                                                                                                             >  2:5686468/1‑150 (MQ=255)
          gACTGGCGTATTTAAGATGATATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTAAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGCCGGGGGCGGcc                                                                                                                                             >  1:9474800/1‑150 (MQ=255)
           aCTGGCGTATTTAAGATGATATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGCAGGGGGCGGCCa                                                                                                                                            <  1:3358689/150‑1 (MQ=255)
             tGGCGTATTTAAGATGATATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGCCGGGGGCGGCCCAt                                                                                                                                          >  1:7186292/1‑150 (MQ=255)
              ggCGTATTTAAGATGATATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGCAGGGGGCGGCCAATc                                                                                                                                         <  1:1550845/150‑1 (MQ=255)
                cGTATTTAAGATGATATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGCAGGGGGCGGCCAATcgc                                                                                                                                       >  2:5563677/1‑150 (MQ=255)
                 gTATTTAAGATGATATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGATCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGCAGGGGGCGGCCAAAcgcc                                                                                                                                      >  2:8839475/1‑150 (MQ=255)
                  tATTTAAGATGATATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGCAGGGGGCGGCCCATcgccg                                                                                                                                     >  2:8655857/1‑150 (MQ=255)
                  tATTTAAGATGATATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGCAGGGGGCGGCCAATcgccg                                                                                                                                     <  1:5197103/150‑1 (MQ=255)
                    tttAAGATGATATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGGCGGGGGGGGCCCAAagccgcc                                                                                                                                   >  1:2052247/1‑150 (MQ=255)
                    tttAAGATGATATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGCCGGGGGCGGCCCAAagccgcc                                                                                                                                   >  1:2933997/1‑150 (MQ=255)
                    tttAAGATGATATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGCAGGGGGCGGCCCAAcgccgcc                                                                                                                                   >  1:6242465/1‑150 (MQ=255)
                    tttAAGATGATATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGCAGGCGGCGGCGATTggccgcc                                                                                                                                   <  2:6554476/150‑1 (MQ=255)
                      tAAGATGATATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGCAGGGGGCGGCGATTggccgcccc                                                                                                                                 <  2:8515122/150‑3 (MQ=255)
                        aGATGATATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGCAGGGGGCGGCCAATCGCCGCCGccc                                                                                                                               <  1:8688921/150‑1 (MQ=255)
                        aGATGATATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGCAGGGGGCGGCCAATCGCCGCCGccc                                                                                                                               <  1:5734029/150‑1 (MQ=255)
                        aGATGATATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGCAGGGGGCGGCCAATCGCCGCCGccc                                                                                                                               <  1:8105534/150‑1 (MQ=255)
                           tgatATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGCAGGGGGCGGCCAAGCGCCGCCGCCCCCt                                                                                                                            <  1:993566/150‑1 (MQ=255)
                            gatATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGCAGGGGGCGGCCAATCGCCGCCGCCCCCTg                                                                                                                           <  1:666816/150‑1 (MQ=255)
                            gatATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGCAGGGGGCGGCCAATCGCCGCCGCCCCCTg                                                                                                                           <  2:4210937/150‑1 (MQ=255)
                              tatTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGCAGGGGGTGGCCAATCGCCGCCGCCCCCTGCa                                                                                                                         <  2:1951581/150‑1 (MQ=255)
                              tatTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGCAGGGGGCGGCCAATCGCCGACGCCCCCAGCa                                                                                                                         >  2:3258041/1‑150 (MQ=255)
                               atTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGGCGGGGGCGGCCAAACccccccccccccccccc                                                                                                                        >  1:3838913/1‑145 (MQ=255)
                                 tttAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGCACGGGGCGGCCCAACCCCCCCCCCCCCCGCAccc                                                                                                                      >  2:3749438/1‑150 (MQ=38)
                                 tttAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGCAAGGGGCGGCCAAAAGCCGCCGCCCCCAGCAccc                                                                                                                      >  2:9673960/1‑150 (MQ=255)
                                   tAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGGGGGGGGGGGGCCATAGCCGGCCCCCCCAGCGccccc                                                                                                                    >  1:7920172/1‑150 (MQ=255)
                                    aaaaTTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGCAGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCg                                                                                                                   <  2:7226463/150‑1 (MQ=255)
                                        ttaattaaTGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGGGGGGGGGGGGCCCACGCCCCCCCCCCCCGCccccccccccg                                                                                                               >  1:8998620/1‑145 (MQ=255)
                                        ttaattaaTGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGGCGGGGGCGGGCGATccccccccccccccccccccccccccc                                                                                                               >  1:1092995/1‑145 (MQ=255)
                                          aattaatGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCAGTGCTTCATATAATCACTGGCGCAAGGAGCGCGCAGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGctct                                                                                                             <  2:3009720/150‑1 (MQ=255)
                                            ttaatGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGCAGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTgg                                                                                                           <  2:9286228/150‑1 (MQ=255)
                                            ttaatGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCAGCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGCAGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTgg                                                                                                           <  2:6235043/150‑1 (MQ=255)
                                             taatGTCATCAGGACCGAAAATAACGAGAATATTACAGTCTCACATCCTGTTGCGCTCCTGTCATGTGCATTGCTACATATAATCACTGGCGCAAGGAGCGCGCAGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGc                                                                                                          <  2:6505255/150‑1 (MQ=255)
                                                tGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGCAGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGaa                                                                                                       <  2:6058533/150‑1 (MQ=255)
                                                tGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGCAGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGaa                                                                                                       <  2:8423623/150‑1 (MQ=255)
                                                tGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGCAGGGGGCGGCCAAAAGCCGCCGCCCCCACCACCCCCCGGCTCAGGCgca                                                                                                       >  2:7096106/1‑148 (MQ=37)
                                                     tcaGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGCCCGGGGCGGCCCCTCGCCCCCCCCCCCCGCCCCCCCCGGGGCAGGCGCACaaaa                                                                                                  >  2:5803186/1‑150 (MQ=25)
                                                     tcaGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGCAGGGGGCGGCCAAACCCCCCCCCCCCCAGCACCCCCCGGCGCTGGCGAACaaaa                                                                                                  >  2:7898525/1‑150 (MQ=25)
                                                          tCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGCACGGGGCGGCCCATACCCGCCGCCCCCAGCACCCCCCGGCACAGGCGAACaaaaaaccc                                                                                             >  2:9511532/1‑145 (MQ=25)
                                                            cGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGGCGGGGGGGGCCAAAAGCCCCCCCCCCCAGCCACCCCCCGCGCTGGCGGAAAAAAtaacccc                                                                                           >  2:186308/1‑144 (MQ=255)
                                                             gAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGCAGGGGGCGGCCAAACCCCGCCGCCCCCAGCACCCCCCGGCTCTGGCGAACAAAATccccccg                                                                                          >  2:1464491/1‑147 (MQ=25)
                                                                    aCGAGAATCTTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCATGGAGCGCGCAGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGt                                                                                   <  2:9104700/150‑1 (MQ=255)
                                                                     cGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGCAGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTg                                                                                  <  2:5298626/150‑1 (MQ=255)
                                                                      gagaATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGCAGGGGGCGGCCAATCGCCGCCCCCCCCAGCACCCCCCGGCACAGGCGAACAAAATCGCCCCAGCGCGGTGc                                                                                 >  2:5479520/1‑150 (MQ=21)
                                                                      gagaATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATAGCTTCATATAAACACTGGCGCAAGGAGCGCTCAGGGGGCGGCCAATAGCCGCCCCCCCCAGCACCCCCGGGCTCAGGAGAAAAAAATCGCCCCAGCGCCGTGc                                                                                 >  2:8014112/1‑150 (MQ=255)
                                                                       ataaTATTTCAGTCTCTCTGCCTGTTGCGCTCCTGTCAGGTGAATTGCTCCTTATAATCAGTGGCGCACGGCGCGCGCTGGGGGCGGCCAACCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGcc                                                                                <  1:6773415/148‑1 (MQ=18)
                                                                         aaTATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGCAGGGggcgg                                                                                                                                               >  1:1452263/1‑85 (MQ=255)
                                                                            atTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGCAGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCgg                                                                           <  2:7702616/150‑1 (MQ=255)
                                                                            atTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGCAGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCgg                                                                           <  2:2250481/150‑1 (MQ=255)
                                                                               tCAGTCTCTCATCCTGTTGCGCTCCTGTCTTGTGATTTGCTTAATATAATCACTGGCGCATGGAGCGCGCTGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCtt                                                                        <  1:7235639/150‑1 (MQ=18)
                                                                                  ggctctcAGCCTGTTGCGCTCCTGTCTGGTGCATGGCTTCATATAATCCCGGGCGCACGGAGCGCGCCGGGGGCGGCCACGCGCCGCCGCCCCCTGCACCCCCGGGCTCGGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATc                                                                     <  1:5047663/148‑1 (MQ=17)
                                                                                   tctctcATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGCAGGGGGCGGCCAAAAGCCGCCCCCCCCAGCACCCCCCGGCACAGGCGAAAAAAAACCCCCCAGCGCGGTGCCCACCGCGTATcc                                                                    >  2:3165400/1‑150 (MQ=18)
                                                                                     tctcATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGCAGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCt                                                                  <  2:4909161/150‑1 (MQ=255)
                                                                                       acATCCTGTTGCGCTCCTGACATGTGAATTGCTTCATATAATCACGGCCGCACGGCGCGCGCTGGGGGCGGCCAAGCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTa                                                                <  1:3441023/149‑1 (MQ=18)
                                                                                         aTCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGCACGGGGCGGCCAAAAGCCCCCCCCCCCAGCACCCCCCCGCACGTGCGAACAAAAAAGCCCCAGCCCGGTGCCCCCCGCGTATACCATAcc                                                              >  2:5336612/1‑149 (MQ=12)
                                                                                          tCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGCGGGGGGCGGGCCAAAGCCGCCCCCCCCAGCACCCCCCGGCGCAGGCCAAAAAAAAACCCCCGGCGCGGTGGCCACCGCGTATACCATAAgg                                                             >  2:8616696/1‑150 (MQ=12)
                                                                                                  cgcTCCTGTCTTGTGCATTGCTTCATATAATCCCTGGCGAAAGGAGCGCGCTGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCggg                                                     <  1:8044717/150‑1 (MQ=255)
                                                                                                    cTCCTTTAATGTGCTTTGATTCATATAATCACTGGCGCAAGGAGCGCGCTGGGGGCGGCAAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTc                                                   <  1:7393055/150‑1 (MQ=14)
                                                                                                     tCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGGCGGGGGCGGCCAAACCCCCCCCCCCCCCCCCCCCCCCGGCGCTGGCGAACAAAAACCCCGCGGCGCGGGGCGCCCGGCGTAACCCGTAAGGCGACCcggggg                                                  >  2:1272482/1‑147 (MQ=12)
                                                                                                     tCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGCAGGGGGCGGCCCAAAGCCCCCCCCCCCCGCACCCCCCGGCGCTGGCGAACAAAAAACCCCCGGCGCGGGGCCCGCCGCGTATACCATACGGCAACCCGGGCg                                                  >  2:6904381/1‑150 (MQ=12)
                                                                                                      ccTGTCATGTGCATTGCTTCATATAATCACTGGCGCATGGAGCGCGCAGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCgg                                                 <  1:9178604/150‑1 (MQ=255)
                                                                                                        tGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGCAGGGGGCGGCCAAAACCCGCCGCCCCCAGCACCCCCCGGCGCTGGCGAACAAAATACCCCCTGCGCGGTGCCCACGGCGTATACCAGAAGGCTAACCgggcgggc                                               >  2:7960634/1‑150 (MQ=17)
                                                                                                            aTGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGCAGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGa                                           <  2:660792/150‑1 (MQ=255)
                                                                                                            aTGTGCATTGCTTAATATAATCACTGGCGCAAGGAGCGCGCAGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGa                                           <  2:1357959/150‑1 (MQ=255)
                                                                                                               tgCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGCCGGGGGCCGCCCAAAGCCCCCCCCCCCCGCACCCCCCGGCGCAGGCCAAAAAAAACGCCCCGGCGCGGTGCCCACCGCGTATCCCATACGGCAAACCGGTCGGGCCCGAGGt                                        >  2:4828854/1‑150 (MQ=255)
                                                                                                                 cATTGCTTCATATAATCACTGGCGCAAGGAGCGCGCAGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTaa                                      <  2:6593602/150‑1 (MQ=255)
                                                                                                                 cATTGCTTCATATAATCACTGGCGCAAGGAGCGCGCAGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTaa                                      <  2:3397244/150‑1 (MQ=255)
                                                                                                                      cTTCATATAATCACTGGCGCAAGGAGCGCGCAGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATcc                                 <  2:2379983/150‑1 (MQ=40)
                                                                                                                      cTTCATATAATCACTGGCGCAAGGAGCGCGCAGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATcc                                 <  2:2684291/150‑1 (MQ=40)
                                                                                                                       ttCATATAATCACTGGCGCAAGGAGCGCGCAGGGGGCGGCCAAAAGCCCCCGCCCCCAGCACCCCCCGGCGCAGGCGCAAAAAATCGCCGCCGCGCCGTGCCCACGGCGTATCCCATACGGCGACCCGGGCGGGCGCGCGGTAACaaacc                                >  2:1607978/1‑146 (MQ=9)
                                                                                                                        tCATATACTCCCGGGCGCCCGGAGCGCGCCGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCt                               <  1:2450711/150‑1 (MQ=11)
                                                                                                                        gcATATAATCCCTGCCGCAAGGAGCGCGCAGGGGGCGGCCAACCGCCGCCGCCCCCTGCACCCCCGGGCTCAGGCGAACAAAATCGCCGCTGCGCGCTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGCAACATCCCt                               <  2:8289386/149‑1 (MQ=11)
                                                                                                                              aaTCACTGGCGCAAGGAGCGCGCAGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTaaaa                         <  2:6257832/150‑1 (MQ=34)
                                                                                                                               aTCACTGGCGCAAGGAGCGCGCAGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAATATCGCCGCTGCGCGGTGCCCTCGGATTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTATAAc                        <  2:5351868/150‑1 (MQ=16)
                                                                                                                               aTCACTGGCGCAAGGAGCGCGCAGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAAc                        <  1:1986993/150‑1 (MQ=33)
                                                                                                                               aTCACTGGCGCAAGGAGCGCGCAGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAAc                        <  1:3894520/150‑1 (MQ=33)
                                                                                                                                  aCTGGCGCAAGGAGCGCGCAGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAAcgcg                     <  1:9001859/150‑1 (MQ=32)
                                                                                                                                   cTGGCGCAAGGAGCGCGCTGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACgcgc                    <  1:4838341/150‑1 (MQ=17)
                                                                                                                                   cTGGCGCAAGGAGCGCGCAGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACgcgc                    <  1:264977/150‑1 (MQ=32)
                                                                                                                                   cTGGCGCAAGGAGCGCGCAGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACgcgc                    <  2:5961148/150‑1 (MQ=32)
                                                                                                                                   cTGGCGCAAGGAGCGCGCAGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACgcgc                    <  2:3838913/150‑1 (MQ=32)
                                                                                                                                      gcgcAAGGAGCGCGCAGGGGGCGGCCCAAAGCCCCCCCCCCCAGCACCCCCCGGCACAGGCGAACAAAAAAGCCCCGTCGCCGTGCCCAACGCATAACCCGTACGGCGACCCGGTCGGGCGCGGGGTTACATCCCGGTAAAAACCGccca                 >  2:2417379/1‑149 (MQ=2)
                                                                                                                                           tggAGCGCGCTGGGGGCGGCCAATCGCCGCCGCCCCCGGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGCCCTCAGcc            <  1:3133835/149‑1 (MQ=17)
                                                                                                                                                cgcgcTGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGCCCTCAGCCCACAt       <  1:1848112/150‑1 (MQ=17)
                                                                                                                                                   gcAGGGGGCGGCCAATCGCCGCCGCCCCCAGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCCGCATATCCCTTACGGCTACCCGGTCGGGCGCGAGGTAACATCCCAGTAAAAAGCGCCCACAGCCCACAACCa    >  2:13105/1‑150 (MQ=11)
                                                                                                                                                    cAGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCGCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACCTCCCTGTAAAACGCGCCCTCAGCCCACATCCAt   <  2:7806325/150‑1 (MQ=17)
                                                                                                                                                    cAGGGGGCGGCCAAGCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGCCCTCAGCCCACATCCAt   <  2:4336265/150‑1 (MQ=17)
                                                                                                                                                     |                                                                                                                                                     
TACTGGTGGTGACTGGCGTATTTAAGATGATATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGCAGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGCCCTCAGCCCACATCCATG  >  NC_000913/2229028‑2229326

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 11 ≤ ATCG/ATCG < 25 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: