Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,229,177 A→C 14.3% 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.3% 330.6 / 11.9 180intergenic (+112/+261)yohP/dusCuncharacterized protein/tRNA‑dihydrouridine synthase C
Reads supporting (aligned to +/- strand):  ref base A (28/98);  new base C (15/11);  total (49/131)
Fisher's exact test for biased strand distribution p-value = 5.88e-04
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 187 total aligned reads displayed.

TACTGGTGGTGACTGGCGTATTTAAGATGATATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑GCAGGGGGC‑GGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGCCCTCAGCCCACATCCAT  >  NC_000913/2229028‑2229325
                                                                                                                                                      |                                                                                                                                                     
tACTGGTGGTGACTGGCGTATTTAAGATGATATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑GCa                                                                                                                                                       <  1:198192/150‑1 (MQ=255)
   tggtggTGACTGGCGTATTTAAGATGATATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑GCAggg                                                                                                                                                    >  1:784726/1‑150 (MQ=255)
    ggtggtGACTGGCGTATTTAAGATGATATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCCTATAATCACTGGCGCAAGGAGCGC‑GCAgggg                                                                                                                                                   <  2:255333/150‑1 (MQ=255)
     gtggtgACTGGCGTATTTAAGATGATATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑gggggggg                                                                                                                                                  >  2:5399277/1‑150 (MQ=255)
     gtggtgACTGGCGTATTTAAGATGATATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑GGAggggg                                                                                                                                                  >  2:103937/1‑150 (MQ=255)
       ggtgACTGGCGTATTTAAGATGATATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑GCgggggggg                                                                                                                                                >  2:2248838/1‑148 (MQ=255)
       ggtgACTGGCGTATTTAAGATGATATTTTAAAATTAATGAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑GCAGGGggc‑g                                                                                                                                               <  1:204131/150‑1 (MQ=255)
        gtgACTGGCGTATTAAAGATGATATTATAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTAATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑Ggggggggg‑gg                                                                                                                                              >  2:5543329/1‑150 (MQ=255)
          gACTGGCGTATTTAAGATGATATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTGTGTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑Ggggggggg‑gggc                                                                                                                                            >  1:2406266/1‑148 (MQ=255)
          gACTGGCGTATTTAAGATGATATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑Ggggggggg‑gggc                                                                                                                                            >  1:2100091/1‑148 (MQ=255)
             tGGCGTATTTAAGATGATATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAATGAGCGC‑GCAGGGGGC‑GGCCAAt                                                                                                                                         <  1:559336/150‑1 (MQ=255)
             tGGCGTATTTAAGATGATATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGGCGGGGGGG‑GGccca                                                                                                                                          >  2:5363763/1‑148 (MQ=255)
             tGGCGTATTTAAGATGATATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑GCGGGGggc‑gggccaa                                                                                                                                         >  2:124724/1‑145 (MQ=255)
              ggCGTATTTAAGATGATATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑GCAGGGGGC‑GGCCAATc                                                                                                                                        <  1:1444217/150‑1 (MQ=255)
              ggCGTATTTAAGATGATATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑GCAGGGGGC‑GGCCAATc                                                                                                                                        <  2:2224417/150‑1 (MQ=255)
                cGTATTTAAGATGATATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑GCAGGGGGC‑GGCCAATcgc                                                                                                                                      <  2:4576568/150‑1 (MQ=255)
                   aTTTAAGATGATATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑GCAGGGGGC‑GGCCAATcgccgc                                                                                                                                   <  2:3635321/150‑1 (MQ=255)
                   aTTTAAGATGATATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑GCAGGGGGC‑GGCCAATcgccgc                                                                                                                                   <  1:5003472/150‑1 (MQ=255)
                    tttAAGATGATATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑GGCGGGGGC‑GGCCCATccccccc                                                                                                                                  >  1:2558792/1‑150 (MQ=255)
                      tAAGATGATATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑GCAGGGGGC‑GGCCAATcgccgccgc                                                                                                                                <  2:3885796/150‑1 (MQ=255)
                       aatattatatTTTAAAATTAATTAATGTCATCAGGTCCTAAAATAACTATAATTTTTCAGTCTCTATTCCTGTTGCGCTACTGTCATGTGCATTGCTTCATATAATCACTGGCGCTAGGAGCGG‑GCTGGGGGC‑GGCCAATCgccgccgcc                                                                                                                               <  1:2620866/144‑1 (MQ=37)
                       aaGATGATATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAGCGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑GGACGGGGC‑GGCCAAACgcccccccc                                                                                                                               >  1:3322521/1‑150 (MQ=255)
                       aaGATGATATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑GCAGGGGGC‑GGCCAATCgccgccgcc                                                                                                                               >  2:5671277/1‑150 (MQ=255)
                        aGATTATATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑GCAGGGGGC‑GGCCAATCGCCGCCGccc                                                                                                                              <  1:1055827/150‑1 (MQ=255)
                        aGATGATATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑GGAGGGGGG‑GGCCCCACGCCGcccccc                                                                                                                              >  2:2183367/1‑150 (MQ=255)
                          atgatATTTTAAAATTAATTAATGTCCTCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑GCAGGGGGC‑GGCCAATCGCCGCCGccccc                                                                                                                            <  1:3873626/150‑1 (MQ=255)
                           tgatATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑GGCGGGGGC‑GGCCCAAACCCGccccccccc                                                                                                                           >  2:4973436/1‑149 (MQ=255)
                             atatTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑GCAGGGGGC‑GGCCAAACGCCGccccccccacc                                                                                                                         >  2:3549810/1‑147 (MQ=255)
                              tatTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑GCCGGGGGC‑GGCCAAAACCCGCCGCCCCCAGCa                                                                                                                        >  2:1828088/1‑150 (MQ=38)
                                  ttAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑GGCGGGGGCGGGCCAAACccccccccccccccacccc                                                                                                                     >  1:716740/1‑143 (MQ=255)
                                  ttAAAATTAATTAATGTCATCAGGTACGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCGTCATATAATCACTGGCGCAAGGAGCGC‑GGCGGGGGG‑GGGCCCGCGCCGCCGGCGCCCGcccccc                                                                                                                    >  1:3129533/1‑150 (MQ=255)
                                   tAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑GCAGGGGGC‑GGCCAATCGCCGCCGCCCCCTGCAccccc                                                                                                                   <  2:3231063/150‑1 (MQ=255)
                                   tAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑GCAGGGGGC‑GGCCAATCGCCGCCGCCCCCTGCAccccc                                                                                                                   <  2:78921/150‑1 (MQ=255)
                                     aaaTTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑GCAGGGGGG‑GGCCCAAACCCCCCCCCCCCAGCAccccccc                                                                                                                 >  2:4030531/1‑148 (MQ=38)
                                      aattaattAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑GCCGGGGGC‑GGCCCAACGCCGCCCCCCCCCGCCCCCCCcgg                                                                                                                >  2:4499053/1‑150 (MQ=38)
                                            ttactGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGGCATGTGCATTGCTTCATATAATCACTGGCGCACGGAGCGC‑GCTGGGGGC‑GGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTgg                                                                                                          <  2:2489221/150‑1 (MQ=255)
                                            ttaatGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑GCAGGGGGC‑CGCCAAACGCCGCCCCCCCCCGCACCCCCCGGcacagg                                                                                                          >  2:1439416/1‑145 (MQ=37)
                                             taatGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑GCAGGGGGC‑GGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGc                                                                                                         <  2:599459/150‑1 (MQ=255)
                                             taatGTAATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATTTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑GCAGGGGGC‑GGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGc                                                                                                         <  2:1787544/150‑1 (MQ=255)
                                              aatGTCATCAGGTCCGAAAATAACTAGACTATTTCAGTCTCGCATCCTGTTGCGCTCCTGTCCTGTGCATTGGTTCATATAATCACTGGCGAAAGGAGAGC‑GCAGGGGGC‑GGCCAATCGGCGGCGCCCCCTGCACCCCCGGGCTCTGGCg                                                                                                        >  2:3344151/1‑150 (MQ=25)
                                               atGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑GCAGGGGGC‑GGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGa                                                                                                       <  2:4955959/150‑1 (MQ=255)
                                                tGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑GGGGGGGGG‑GGCCCAACGCCCCCCCCCCCCCCACCCCCCGGCGCAGGCgcc                                                                                                      >  1:5313376/1‑148 (MQ=255)
                                                  tcatcaGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑GCAGGGGGC‑GGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACa                                                                                                    <  2:1279481/150‑1 (MQ=255)
                                                    atcaGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGGCGGGGGGG‑GGCCAAAC‑CCGCCCCCCCCAGCCCCCCCCGGCCCGGGCGGCCaaa                                                                                                  >  1:4110856/1‑150 (MQ=255)
                                                       tggTCCTAAAATAACGATATTATTACAGTCTATAAGCGTTTTGCGCTCCTTTCTTGTGAATTGCTTCATATAATCCCTGGCGCACGGCGGGC‑GCCGGGGGG‑GGCAAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATc                                                                                               <  1:4345967/149‑1 (MQ=16)
                                                       aGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑GGCGGGGGG‑GGCCCAACCCCCCCCCCCCCAGCCCCCCCCCGCGCAGGCCCACaaaaaa                                                                                               >  1:452995/1‑148 (MQ=255)
                                                       aGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑GCAGGGGGC‑GGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATc                                                                                               <  2:1672409/150‑1 (MQ=255)
                                                         gTCCGACCATAACGCGAAGATTTCAGTCTCTCAGCCTGGTGCGCTCCTGTCAGGTGCATTGCTTCATCCAAACCCTGGCGCAAGGAGCGC‑GCAGGGGGC‑GGCCACACGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATcgc                                                                                             <  2:3020090/150‑1 (MQ=18)
                                                         gTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑GCAGGGGGC‑GGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATcgc                                                                                             <  2:4169933/150‑1 (MQ=255)
                                                              aaaaTAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑GCAGGGGGC‑GGCCAATCGCCGCCGCCCCCAGCCCCCCCCGGCACAGGCGAACAAAAAcgccgcag                                                                                        >  2:3863860/1‑148 (MQ=25)
                                                                 aTAACGAGATTATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑GCAGGGGGC‑GGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTgcgc                                                                                     <  2:1145046/150‑1 (MQ=255)
                                                                        gaATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACGGGCGCAAGGAGCGC‑GCCGGGGGC‑GGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGccc                                                                              <  2:1667468/150‑1 (MQ=255)
                                                                           tatTTCTGTCTCTCATCCTGTTGCGCTCCTGTCATGTGAATTGCTACATATAATCACTGGCGCACGGAGCGC‑GCAGGGGGC‑GGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCg                                                                           <  1:3253199/150‑1 (MQ=255)
                                                                            atTTCAGTCTCTCTTCGTGTTGCGCTCCTGTCTTGTGAATTGATTAATATAACCACTGGCGCATGGAGCGC‑GCTGGGGGC‑GGCCAACCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCgg                                                                          <  1:1367621/150‑1 (MQ=17)
                                                                              tcccgcctctcAGCGTGTGGCGCCCCTGTCTTGTGATTTTATTCATATAAACCCGGGCGCAGGGCGCGG‑GCGGGGGGC‑GGCCACCCGCCGCCGCCCCCGGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCt                                                                        <  1:519400/144‑1 (MQ=9)
                                                                                  gtatctaATCCTGTTGCGCTCCTGTAATGTGCATTGCTTAATATAATCACTGGCGAATGGAGCGC‑GCTGGGGGG‑GGCCAACCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATc                                                                    <  1:5681562/147‑1 (MQ=18)
                                                                                   tctctcATCCTGTGGCGCTCCTTTCTTGTGCTTTGTTTCATATAATCACTGGCGCATGGCGCGC‑GCTGGGGGC‑GGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATcc                                                                   <  1:1947138/150‑1 (MQ=17)
                                                                                    ctctaATCCTGTTGCGCGCCTGTCTTGTGCTTGGCTTCATATAACCCCGGGCGCATGGCGCGC‑GCCGGGGGC‑GGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATccc                                                                  <  1:2901353/150‑1 (MQ=17)
                                                                                      ctcATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑GCCGGGGGC‑GGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCtt                                                                <  1:3672562/150‑1 (MQ=255)
                                                                                       tcATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑GCAGGGGGC‑GGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTa                                                               <  2:4477963/150‑1 (MQ=255)
                                                                                         aTCCTTTTGCGCTCCTTTCTTGTGCATTGCTTAATATAATCCCGGGCGCACGGAGCGC‑GCCGGGGGC‑GGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACg                                                             <  1:1249408/150‑1 (MQ=17)
                                                                                         aTCCTGTTGCGCTCCTGGCATGTGAATTGCTTAATATAATCCCGGGCGCATGGAGCGC‑GCAGGGGGC‑GGCCAAGCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACg                                                             <  1:5166465/150‑1 (MQ=17)
                                                                                              gTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑GCCGGGGGC‑GGCCAAACCCCCCCCCCCCCCGCACCCCCCGGCGCAGGCGAACAAAAAAGCCCCTGCGCGGTGCCCACGGCTTATACCATAAGGCTAc                                                        >  2:1920451/1‑150 (MQ=17)
                                                                                              gTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑GCAGGGGGC‑GGCCAAAAGCCGCCCCCCCCAGCCCCCCCCCGCGCATGCGAACAAAATACCCGCGGCGCGGTGCCCTCGGCTTATCACTAACCGCAAc                                                        >  2:2368161/1‑150 (MQ=17)
                                                                                               ttGCGCTCCTGTCTTGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑GCAGGGGGC‑GGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAAC‑AAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCg                                                      <  1:2147739/150‑1 (MQ=255)
                                                                                                    cTCCTGTCATGTGCATTGCTTCATTTAATCACTGGCGAACGGAGGGC‑GCAGGGGGC‑GGCAAATCGCGGCCGCCCCGTGCACCCGCTGGCTCTGGCGATGATAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTc                                                  <  2:66155/150‑1 (MQ=17)
                                                                                                    atcatgacATGTGAATTGATTCATATAATCACTGGCGAAAGGAGCGC‑GCTGGGGGC‑GGCAAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTc                                                  <  1:4379943/143‑1 (MQ=14)
                                                                                                         gTCATGTGCATTGCTTCATATAATCACTGGCGCATGGAGCGC‑GCAGGGGGC‑GGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGgcg                                             <  1:4863410/150‑1 (MQ=255)
                                                                                                            aggtgCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑GCAGGGGGC‑GGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGa                                          <  2:5622371/148‑1 (MQ=255)
                                                                                                            aTGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑GCAGGGGGC‑GGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGa                                          <  2:3925117/150‑1 (MQ=255)
                                                                                                             tttgctttgcttAATATAATCACTGGCGCAAGGAGCGC‑GCTGGGGGC‑GGCCAATCGACGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATGCCTTACGGCTACCGGGTCGGGCGCGAg                                         <  2:4404163/148‑1 (MQ=14)
                                                                                                             tgtgCCTTGCTTCATATACTCACGGGCGCACGGAGCGC‑GCCGGGGGC‑GGCCACGCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAg                                         <  1:4143110/150‑1 (MQ=14)
                                                                                                                 cATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑GCAGGGGGC‑GGCCAATCCCCCCCCCCCCCAGCAACCCCCGGCTCAGGCGAACAAAATCGCCGCTGCGCGGTGGCCACGGCGTATCCCTTAAGGCTACCGGGTGGGGCGCGAGGGaa                                     >  2:4840539/1‑150 (MQ=14)
                                                                                                                      cTTCATATAATCACTGGCGCAAGGAGCGC‑GCAGGGGGC‑GGCCAATCCCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATcc                                <  1:2951294/150‑1 (MQ=255)
                                                                                                                        tCATATAATCACTGGCGCAAGGAGCGC‑GCAGGGGGC‑GGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCt                              <  2:2393302/150‑1 (MQ=39)
                                                                                                                        tCATATAATCACTGGCGCAAGGAGCGC‑GCAGGGGGC‑GGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCt                              <  2:2831845/150‑1 (MQ=39)
                                                                                                                        tCATATAATCACTGGCGCAAGGAGCGC‑GCAGGGGGC‑GGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCt                              <  2:2393522/150‑1 (MQ=39)
                                                                                                                             taATCACTGGCGCAAGGAGCGC‑GCAGGGGGC‑GGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTaaa                         <  1:5228425/150‑1 (MQ=35)
                                                                                                                              aaTCACTGGCGCAAGGAGCGC‑GCAGGGGGC‑GGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCACTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTaaaa                        <  2:1536628/150‑1 (MQ=21)
                                                                                                                              aaTCACGGGCGCAAGGAGCGC‑GCAGGGGGC‑GGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTaaaa                        <  2:716740/150‑1 (MQ=21)
                                                                                                                               aTCACTGGCGCAAGGAGCGC‑GCTGGGGGT‑GGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAAc                       <  1:5370502/150‑1 (MQ=18)
                                                                                                                               aTCACTGGCGCAAGGAGCGC‑GCAGGGGTT‑GGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAAc                       <  1:5370483/150‑1 (MQ=18)
                                                                                                                                 cACTGGCGCAAGGAGCGC‑GCAGGGGGC‑GGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAAcgc                     <  1:4870063/150‑1 (MQ=32)
                                                                                                                                   cTGGCGCAAGGAGCGC‑GCAGGGGGC‑GGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACgcgc                   <  2:3607324/150‑1 (MQ=32)
                                                                                                                                   cTGGCGCAAGGAGCGC‑GCAGGGGGC‑GGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACgcgc                   <  1:3603776/150‑1 (MQ=32)
                                                                                                                                   cTGGCGCAAGGAGCGC‑GCAGGGGGC‑GGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACgcgc                   <  1:2226849/150‑1 (MQ=32)
                                                                                                                                   cTGGCGCAAGGAGCGC‑GCAGGGGGC‑GGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACgcgc                   <  2:5215654/150‑1 (MQ=32)
                                                                                                                                    tGGCGCAAGGAGCGC‑GCTGGGGGC‑GGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGcc                  <  1:1005777/150‑1 (MQ=17)
                                                                                                                                    tGGCGCAAGGAGCGC‑GCAGGGGGC‑GGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGcc                  <  1:197342/150‑1 (MQ=32)
                                                                                                                                          aaGGAGCGC‑GCCGGGGGC‑GGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGCCCTCAGc            <  1:873747/150‑1 (MQ=17)
                                                                                                                                           aGGAGCGC‑GCAGGGGGC‑GGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAACATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGATGTAACATCCCTGTAAAACGCGCCCTCAGcc           <  2:4409903/150‑1 (MQ=17)
                                                                                                                                           aGGAGCGC‑GCAGGGGGC‑GGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGCCCTCAGcc           <  2:2902389/150‑1 (MQ=32)
                                                                                                                                             gAGCGC‑GCTGGGGGC‑GGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGCCCTCAGCCca         <  1:5664912/150‑1 (MQ=17)
                                                                                                                                             gAGCGC‑GCAGGGGGC‑GGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGCCCTCAGCCca         <  1:2683878/150‑1 (MQ=32)
                                                                                                                                              aGCGC‑GCAGGGGGC‑GGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGCCCTCAGCCcac        <  2:4318884/150‑1 (MQ=32)
                                                                                                                                                cgc‑gcAGGGGGC‑GGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGCCCTCAGCCCACAt      <  1:149215/150‑1 (MQ=32)
                                                                                                                                                 gc‑gcAGGGGGC‑GGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGTCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGCCCTCAGCCCACATc     >  2:3471563/1‑150 (MQ=17)
                                                                                                                                                 gc‑gcAGGGGGC‑GGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGCCCTCAGCCCACATc     >  2:3055019/1‑150 (MQ=32)
                                                                                                                                                     cAGGGGGC‑GGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGCCCTCAGCCCACATCCAt  >  2:3727423/1‑150 (MQ=32)
                                                                                                                                                     cAGGGGGC‑GGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTCCCCCCCGCGCCCTCAGCCCACATCCAt  <  1:5363763/150‑1 (MQ=11)
                                                                                                                                                      |                                                                                                                                                     
TACTGGTGGTGACTGGCGTATTTAAGATGATATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑GCAGGGGGC‑GGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGCCCTCAGCCCACATCCAT  >  NC_000913/2229028‑2229325

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