Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,229,177 A→C 15.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,1770AC15.4% 107.0 / 4.1 66intergenic (+112/+261)yohP/dusCuncharacterized protein/tRNA‑dihydrouridine synthase C
Reads supporting (aligned to +/- strand):  ref base A (22/26);  new base C (7/2);  total (33/33)
Fisher's exact test for biased strand distribution p-value = 1.44e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TACTGGTGGTGACTGGCGTATTTAAGATGATATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑GCAGGGGGC‑GGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGC‑TTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGCCCTCAGCCCACATCCAT  >  NC_000913/2229028‑2229325
                                                                                                                                                      |                                                                                                                                                      
tACTGGTGGTGACTGGCGTATTTAAGATGATATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑GCa                                                                                                                                                        <  1:29417/150‑1 (MQ=255)
     gtggtgACTGGCGTATTTAAGATGATATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑GCAggggg                                                                                                                                                   >  2:1446967/1‑150 (MQ=255)
     gtggtgACTGGCGTATTTAAGATGATATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑GCAggggg                                                                                                                                                   >  2:2463269/1‑150 (MQ=255)
     gtggtgACTGGCGTATTTAAGATGATATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑GCAggggg                                                                                                                                                   >  2:2038358/1‑150 (MQ=255)
         tgACTGGCGTATTTAAGATGATATTTTAAAATTAATTACTGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑GCAGGGggc‑ggc                                                                                                                                              <  1:2536535/150‑1 (MQ=255)
         tgACTGGCGTATTTAAGATGATATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTTATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑GCAGGGggc‑ggc                                                                                                                                              <  1:2536725/150‑1 (MQ=255)
         tgACTGGCGTATTTAAGATGATATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑GCAGGGggc‑ggc                                                                                                                                              <  1:2536811/150‑1 (MQ=255)
          gACTGGCGTATTTAAGATGATATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑GCAGGGggc‑ggcg                                                                                                                                             >  1:875203/1‑149 (MQ=255)
            cTGGCGTATTTAAGATGATATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑GCAGGGGGC‑GGCCaa                                                                                                                                           <  2:355827/150‑1 (MQ=255)
             tGGCGTATTTAAGATGATATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑GCAGGGGGC‑GGCCCAt                                                                                                                                          >  2:1905057/1‑150 (MQ=255)
              ggCGTATTTAAGATGATATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑GCAGGGGGC‑GGCCAATc                                                                                                                                         <  1:1628671/150‑1 (MQ=255)
               gCGTATTTAAGATGATATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑GCAGGGGGC‑GGCCaaaag                                                                                                                                        >  2:210250/1‑147 (MQ=255)
               gCGTATTTAAGATGATATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑GCAGGGGGC‑GGCCAATcg                                                                                                                                        >  2:1290405/1‑150 (MQ=255)
               gCGTATTTAAGATGATATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑GCAGGGGGC‑GGCCAAAcg                                                                                                                                        >  1:1517994/1‑150 (MQ=255)
                  tATTTAAGATGATATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑GCCCCGGGC‑GGCCCAtacccc                                                                                                                                     >  2:2282108/1‑145 (MQ=255)
                  tATTTAAGATGATATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑GCAGGGGGC‑GGCCAAAcgccg                                                                                                                                     >  2:1694088/1‑150 (MQ=255)
                      tAAGATGATATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATACCGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCACGGAGCGC‑GCAGGGGGC‑GGCCAATcgccgccgc                                                                                                                                 <  2:1626092/150‑1 (MQ=255)
                      tAAGATGATATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑GCAGGGGGC‑GGCCCAAccccgcccc                                                                                                                                 >  2:1072015/1‑148 (MQ=255)
                         gattatATTTTAAAATTAATTAATGTCATCAGGTCCTAAAATAACGAGAATTTTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑GCAGGGGGC‑GGCCAATCGCCGCCGcccc                                                                                                                              <  1:1021475/150‑1 (MQ=255)
                         gatgatATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑GCAGGGGGC‑GGCCAATCGCCGCCGcccc                                                                                                                              <  2:1640396/150‑1 (MQ=255)
                           tgatATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCTTTGCTTCATATAATCACGGGCGCAAGGAGCGC‑GCAGGGGGC‑GGCCACTCGCCGCCGCCCCCt                                                                                                                            <  1:971812/150‑1 (MQ=255)
                           tgatATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑GCCCGGGGG‑GGCCCCACGCccccccccccc                                                                                                                            >  2:3948/1‑149 (MQ=255)
                            gatATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑GGCGGGGGCGGGCCAAACACccccccccccc                                                                                                                            >  1:2217746/1‑149 (MQ=255)
                                 tttAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑GCGGGGGGG‑GGGCCAGAGCCCCCCCCCCCAGCAccc                                                                                                                      >  1:1199127/1‑150 (MQ=38)
                                   tAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑GGCGGGGGG‑GGGCAAACGCCCCCCCCCCCCGCAccccc                                                                                                                    >  2:1314931/1‑150 (MQ=255)
                                      aattaattAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑GGACGGGGC‑GGCCCAAAGCCGCCGccccccccccccccccg                                                                                                                 >  2:340941/1‑147 (MQ=255)
                                      aattaattAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑GCCGGGGGC‑GGCCAAAACCCGCCGCCCCCCGCACCCCCcgg                                                                                                                 >  2:417226/1‑150 (MQ=38)
                                      aattaattAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑GCAGGGGGC‑GGCCAAACGCCGCCCCCCCCAGCACCCCCcgg                                                                                                                 >  2:417081/1‑150 (MQ=255)
                                          aattaatGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑GCCGGGGGC‑GGGCCAAAGCCCCCCCCCCCCCCACCCCCCGGcgct                                                                                                             >  1:203438/1‑150 (MQ=37)
                                             taatGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑GCAGGGGGC‑GGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGc                                                                                                          <  2:1517994/150‑1 (MQ=255)
                                              aatGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑GCAGGGGGC‑CGCCCAACGCCGCCCCCCCCAGCACCCCCCGGcagatgcg                                                                                                         >  2:792972/1‑143 (MQ=37)
                                                tgacatcaGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑GCAGGGGGC‑GGCCCAACCCCGCCGCCCCCAGCCCCCCCCGGCGCTGGCGaa                                                                                                       >  2:818878/4‑150 (MQ=37)
                                                  tcatcaGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑GCAGGGGGG‑GGCCAATAGCCGCCGCCCCCAGCACCCCCCGGCGCAGGCGCAcc                                                                                                     >  2:1037554/1‑149 (MQ=37)
                                                        ggTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGGCGGGGGGG‑GGCCAAAC‑CCCCCGCCCCCACCCCCCCCCGGCGCAGGCCAACaaaaaac                                                                                               >  1:2148003/1‑147 (MQ=255)
                                                        ggTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGGCGGGGGGG‑GGCCAA‑CGCCCCCCCCCCCCCCCCCCCCCCGCACAGGCGAACaaaaaac                                                                                               >  1:809494/1‑147 (MQ=255)
                                                        ggTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑GGGGGGGGG‑GGCCAAAAGCCCCCCCCCCCCGCCCCCCCCCGCGCTGGCGAAAAAAAAcg                                                                                               >  1:2098730/1‑150 (MQ=255)
                                                        ggTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑GGCCGGGGG‑GGCCCAAAGCCCCCCCCCCCCCCCCCCCCCCCCGCAGGCCAACaaaaaac                                                                                               >  1:2147821/1‑147 (MQ=255)
                                                           ccGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑GGAGGGGGG‑GGGCCAACCCCGCCCCCCCCAGCACCCCCCCGCGCGGGCGGACaaaaaagccc                                                                                            >  1:1624574/1‑144 (MQ=255)
                                                                  tAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑GCAGGGGGC‑GGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGcgcg                                                                                     <  2:634837/150‑1 (MQ=255)
                                                                          atatTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑GCAGGGGGC‑GGCCAAAAGCCCCCCCCCCCAGCCACCCCCGGCACAGGCGAACAAAAAAGCCCCAGCCCGGTGcccac                                                                             >  2:1279107/1‑148 (MQ=21)
                                                                            atTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑GCAGGGGGC‑GGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCgg                                                                           <  2:1872472/150‑1 (MQ=255)
                                                                            atTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑GCAGGGGGC‑GGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCgg                                                                           <  2:2421251/150‑1 (MQ=255)
                                                                            atTTCAGTCTCTAATCCTGTTGCGCTCGTTTCTTGTGATTTGATTCATATAAACACTGGCGCAAGGAGCGC‑GCCGGGGGC‑GGCCAAGCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCgg                                                                           <  1:1812050/150‑1 (MQ=17)
                                                                                  gTCTCTCACCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACGGGCGCAAGGCGCGC‑GCCGGGGGC‑GGCCAAACGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGC‑TTATc                                                                    <  1:396443/150‑1 (MQ=18)
                                                                                    ctatcATCCTTTTGCGCTCCTTTATTGTGCATTGCTTAATATAATCACTGGCGCATGGTGCGC‑GCTGGGGGC‑GGCAAACCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGC‑TTATccc                                                                  <  1:210250/147‑1 (MQ=17)
                                                                                     tctcATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑GCAGGGGGC‑GGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGC‑TTATCCCt                                                                 <  2:1008613/150‑1 (MQ=255)
                                                                                      ctcAGCGTGTTGCGCTCCTGTCATGTGATTTGCTACATAAAACCCCTGGCGCACGGCGCGG‑GCGGGGGGC‑GGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGC‑TTATCCCtt                                                                <  1:1072015/150‑1 (MQ=14)
                                                                                       tcttCCTTTTGGGCTCTTGTATTGTGAATTGTTTAATATAATCACTGGCGCTTGGAGGGC‑GCTGGGGGG‑GGCAAATCCCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCTGTGCCCTCGGC‑TTATCCCTTa                                                               <  1:203623/147‑1 (MQ=14)
                                                                                       tcATCCTGTTGGGCTCCTGTCTTGTGCTTTGCTTAATATAATCACTGGCGAATGGAGCGC‑GCTGGGGGC‑GGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGC‑TTATCCCTTa                                                               <  1:194639/150‑1 (MQ=17)
                                                                                       acATCCTGTTGCGCGCCTGTCTTGTGCTTGGCTTCATATAAACACTGGCGCATGGAGCGC‑GCTGGGGGC‑GGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGC‑TTATCCCTTa                                                               <  1:203899/149‑1 (MQ=17)
                                                                                                           cATGTGCATTGCTTCCTATAATCCCGGGCGCACGGCGCGC‑GCAGGGGGC‑GGCCAAGCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGC‑TTATCCCTTACGGCTACCGGGTCGGGcgcg                                           <  1:727021/150‑1 (MQ=14)
                                                                                                           cATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑GCAGGGGGC‑GGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGC‑TTATCCCTTACGGCTACCGGGTCGGGcgcg                                           <  1:1667966/150‑1 (MQ=255)
                                                                                                           cATGTGCATTGATAAATATAAACACTGGCGCAAGGAGCGC‑GCAGGGGGC‑GGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGC‑TTATCCCTTACGGCTACCGGGTCGGGcgcg                                           <  1:1111447/150‑1 (MQ=255)
                                                                                                                 cATTGCTTCATATAATCACTGGCGCTAGGAGCGC‑GCAGGGGGC‑GGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGC‑TTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTaa                                     <  1:1421959/150‑1 (MQ=255)
                                                                                                                 cATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑GCAGGGGGC‑GGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGC‑TTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTaa                                     <  1:1421952/150‑1 (MQ=255)
                                                                                                                  aTTGCTTCATATAATCACTGGCGCAAGGAGCGC‑GCAGGGGGG‑GGCCCAAAGCCCCCCCCCCCCCCCCCCCCCGGCGCAGGCCCAAAAAAAACCCGCGTCGCCGTGCCCACCGCGGTAACCC‑TAAGGCAAACCGGTCGGGCCCGAGGTAAc                                    >  2:1267481/1‑150 (MQ=255)
                                                                                                                             taATCACTGGCGCAAGGAGCGC‑GCAGGGGGC‑CGCCAATACCCCCCGCCCCCAGCACCCCCCGGCTCTGGCGCACAAAAAACCCGCGGCCCGGTGCCCACCGC‑GTATCCCGTACGGCTACCCGGGCGGGCGCGAGGTAACAACCCCGTaaa                         >  2:2547922/1‑150 (MQ=2)
                                                                                                                              aaTCACTGGCGCAAGGAGCGC‑GCAGGGGGC‑GGCCAAACGCCCCCGCCCCCAGCACCCCCCGGCTCTGGCGAAAAAAAAAGCCGCTGCGCGGTGCCCTCGGC‑GTATACCTTACGGCAACCCGGTCGGGCGCCCGGTAACATCCCGGTaaaa                        >  2:1495109/1‑150 (MQ=11)
                                                                                                                                  aCTGGCGCAAGGAGCGC‑GCAGGGGGC‑GGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGC‑TTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAAcgcg                    <  2:1577176/150‑1 (MQ=32)
                                                                                                                                  aCTGGCGCAAGGAGCGC‑GCAGGGGGC‑GGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGC‑TTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAAcgcg                    <  1:652516/150‑1 (MQ=32)
                                                                                                                                  aCTGGCGCAAGGAGCGC‑GCAGGGGGC‑GGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCGGGGCGGTGCCCACGGC‑TTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAAcgcg                    <  2:1343462/150‑1 (MQ=12)
                                                                                                                                    tGGCGCAAGGAGCGC‑GCAGGGGGC‑GGCCCAACGCCCCCGCCCCCGGCCCCCCCCGGCGCAGGCGAACAAAAAACCCCCTGCGCGGGGCCCCCCGC‑GTATACCGTAACGCTACCCGGGCGGGGGCGAGGGTACAACCCCGTAAAACCCGcc                  >  2:395072/1‑150 (MQ=2)
                                                                                                                                      gcgcAAGGAGCGC‑GCATGGGGC‑GGCCAAAAGCCGCCGCCCCCAGCACCCCCGGGCACTTGCGAACAAAAACGCCGCAGCGCGGTGCCCTAGGC‑TTATCCCATACGGCTACCGGGTCGGGCGCGAGGTAACATCCCAGTAAAAAGCGccca                >  2:1731551/1‑149 (MQ=11)
                                                                                                                                         cAAGGAGCGC‑GCAGGGGGC‑GGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGC‑TTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGCCCTCAg             <  1:596403/150‑1 (MQ=32)
                                                                                                                                         cAAGGAGCGC‑GCAGGGGGC‑GGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGC‑TTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGCCCTCAg             <  2:692974/150‑1 (MQ=32)
                                                                                                                                                 gc‑gcAGGGGGC‑GGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGC‑TTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGCCCTCAGCCCACATc     <  1:1624535/150‑1 (MQ=32)
                                                                                                                                                     cAGGGGGC‑GGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGC‑TTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGCCCTCAGCCCACATCCAt  <  1:1596417/150‑1 (MQ=32)
                                                                                                                                                      |                                                                                                                                                      
TACTGGTGGTGACTGGCGTATTTAAGATGATATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑GCAGGGGGC‑GGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGC‑TTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGCCCTCAGCCCACATCCAT  >  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: ‑