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

TACTGGTGGTGACTGGCGTATTTAAGATGATATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGCAGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGCCCTCAGCCCACATCCAT  >  NC_000913/2229028‑2229325
                                                                                                                                                      |                                                                                                                                                      
tACTGGTGGTGACTGGCGTATTTAAGATGATATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGCa                                                                                                                                                        <  1:29417/150‑1 (MQ=255)
     gtggtgACTGGCGTATTTAAGATGATATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGCAggggg                                                                                                                                                   >  2:1446967/1‑150 (MQ=255)
     gtggtgACTGGCGTATTTAAGATGATATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGCAggggg                                                                                                                                                   >  2:2463269/1‑150 (MQ=255)
     gtggtgACTGGCGTATTTAAGATGATATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGCAggggg                                                                                                                                                   >  2:2038358/1‑150 (MQ=255)
         tgACTGGCGTATTTAAGATGATATTTTAAAATTAATTACTGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGCAGGGggcggc                                                                                                                                              <  1:2536535/150‑1 (MQ=255)
         tgACTGGCGTATTTAAGATGATATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTTATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGCAGGGggcggc                                                                                                                                              <  1:2536725/150‑1 (MQ=255)
         tgACTGGCGTATTTAAGATGATATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGCAGGGggcggc                                                                                                                                              <  1:2536811/150‑1 (MQ=255)
          gACTGGCGTATTTAAGATGATATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGCAGGGggcggcg                                                                                                                                             >  1:875203/1‑149 (MQ=255)
            cTGGCGTATTTAAGATGATATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGCAGGGGGCGGCCaa                                                                                                                                           <  2:355827/150‑1 (MQ=255)
             tGGCGTATTTAAGATGATATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGCAGGGGGCGGCCCAt                                                                                                                                          >  2:1905057/1‑150 (MQ=255)
              ggCGTATTTAAGATGATATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGCAGGGGGCGGCCAATc                                                                                                                                         <  1:1628671/150‑1 (MQ=255)
               gCGTATTTAAGATGATATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGCAGGGGGCGGCCaaaag                                                                                                                                        >  2:210250/1‑147 (MQ=255)
               gCGTATTTAAGATGATATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGCAGGGGGCGGCCAATcg                                                                                                                                        >  2:1290405/1‑150 (MQ=255)
               gCGTATTTAAGATGATATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGCAGGGGGCGGCCAAAcg                                                                                                                                        >  1:1517994/1‑150 (MQ=255)
                  tATTTAAGATGATATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGCCCCGGGCGGCCCAtacccc                                                                                                                                     >  2:2282108/1‑145 (MQ=255)
                  tATTTAAGATGATATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGCAGGGGGCGGCCAAAcgccg                                                                                                                                     >  2:1694088/1‑150 (MQ=255)
                      tAAGATGATATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATACCGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCACGGAGCGCGCAGGGGGCGGCCAATcgccgccgc                                                                                                                                 <  2:1626092/150‑1 (MQ=255)
                      tAAGATGATATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGCAGGGGGCGGCCCAAccccgcccc                                                                                                                                 >  2:1072015/1‑148 (MQ=255)
                         gattatATTTTAAAATTAATTAATGTCATCAGGTCCTAAAATAACGAGAATTTTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGCAGGGGGCGGCCAATCGCCGCCGcccc                                                                                                                              <  1:1021475/150‑1 (MQ=255)
                         gatgatATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGCAGGGGGCGGCCAATCGCCGCCGcccc                                                                                                                              <  2:1640396/150‑1 (MQ=255)
                           tgatATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCTTTGCTTCATATAATCACGGGCGCAAGGAGCGCGCAGGGGGCGGCCACTCGCCGCCGCCCCCt                                                                                                                            <  1:971812/150‑1 (MQ=255)
                           tgatATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGCCCGGGGGGGCCCCACGCccccccccccc                                                                                                                            >  2:3948/1‑149 (MQ=255)
                            gatATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGGCGGGGGCGGGCCAAACACccccccccccc                                                                                                                            >  1:2217746/1‑149 (MQ=255)
                                 tttAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGCGGGGGGGGGGCCAGAGCCCCCCCCCCCAGCAccc                                                                                                                      >  1:1199127/1‑150 (MQ=38)
                                   tAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGGCGGGGGGGGGCAAACGCCCCCCCCCCCCGCAccccc                                                                                                                    >  2:1314931/1‑150 (MQ=255)
                                      aattaattAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGGACGGGGCGGCCCAAAGCCGCCGccccccccccccccccg                                                                                                                 >  2:340941/1‑147 (MQ=255)
                                      aattaattAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGCCGGGGGCGGCCAAAACCCGCCGCCCCCCGCACCCCCcgg                                                                                                                 >  2:417226/1‑150 (MQ=38)
                                      aattaattAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGCAGGGGGCGGCCAAACGCCGCCCCCCCCAGCACCCCCcgg                                                                                                                 >  2:417081/1‑150 (MQ=255)
                                          aattaatGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGCCGGGGGCGGGCCAAAGCCCCCCCCCCCCCCACCCCCCGGcgct                                                                                                             >  1:203438/1‑150 (MQ=37)
                                             taatGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGCAGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGc                                                                                                          <  2:1517994/150‑1 (MQ=255)
                                              aatGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGCAGGGGGCCGCCCAACGCCGCCCCCCCCAGCACCCCCCGGcagatgcg                                                                                                         >  2:792972/1‑143 (MQ=37)
                                                tgacatcaGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGCAGGGGGCGGCCCAACCCCGCCGCCCCCAGCCCCCCCCGGCGCTGGCGaa                                                                                                       >  2:818878/4‑150 (MQ=37)
                                                  tcatcaGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGCAGGGGGGGGCCAATAGCCGCCGCCCCCAGCACCCCCCGGCGCAGGCGCAcc                                                                                                     >  2:1037554/1‑149 (MQ=37)
                                                        ggTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGGCGGGGGGGGGCCAAACCCCCCGCCCCCACCCCCCCCCGGCGCAGGCCAACaaaaaac                                                                                               >  1:2148003/1‑147 (MQ=255)
                                                        ggTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGGCGGGGGGGGGCCAACGCCCCCCCCCCCCCCCCCCCCCCGCACAGGCGAACaaaaaac                                                                                               >  1:809494/1‑147 (MQ=255)
                                                        ggTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGGGGGGGGGGGCCAAAAGCCCCCCCCCCCCGCCCCCCCCCGCGCTGGCGAAAAAAAAcg                                                                                               >  1:2098730/1‑150 (MQ=255)
                                                        ggTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGGCCGGGGGGGCCCAAAGCCCCCCCCCCCCCCCCCCCCCCCCGCAGGCCAACaaaaaac                                                                                               >  1:2147821/1‑147 (MQ=255)
                                                           ccGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGGAGGGGGGGGGCCAACCCCGCCCCCCCCAGCACCCCCCCGCGCGGGCGGACaaaaaagccc                                                                                            >  1:1624574/1‑144 (MQ=255)
                                                                  tAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGCAGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGcgcg                                                                                     <  2:634837/150‑1 (MQ=255)
                                                                          atatTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGCAGGGGGCGGCCAAAAGCCCCCCCCCCCAGCCACCCCCGGCACAGGCGAACAAAAAAGCCCCAGCCCGGTGcccac                                                                             >  2:1279107/1‑148 (MQ=21)
                                                                            atTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGCAGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCgg                                                                           <  2:1872472/150‑1 (MQ=255)
                                                                            atTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGCAGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCgg                                                                           <  2:2421251/150‑1 (MQ=255)
                                                                            atTTCAGTCTCTAATCCTGTTGCGCTCGTTTCTTGTGATTTGATTCATATAAACACTGGCGCAAGGAGCGCGCCGGGGGCGGCCAAGCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCgg                                                                           <  1:1812050/150‑1 (MQ=17)
                                                                                  gTCTCTCACCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACGGGCGCAAGGCGCGCGCCGGGGGCGGCCAAACGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATc                                                                    <  1:396443/150‑1 (MQ=18)
                                                                                    ctatcATCCTTTTGCGCTCCTTTATTGTGCATTGCTTAATATAATCACTGGCGCATGGTGCGCGCTGGGGGCGGCAAACCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATccc                                                                  <  1:210250/147‑1 (MQ=17)
                                                                                     tctcATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGCAGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCt                                                                 <  2:1008613/150‑1 (MQ=255)
                                                                                      ctcAGCGTGTTGCGCTCCTGTCATGTGATTTGCTACATAAAACCCCTGGCGCACGGCGCGGGCGGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCtt                                                                <  1:1072015/150‑1 (MQ=14)
                                                                                       tcttCCTTTTGGGCTCTTGTATTGTGAATTGTTTAATATAATCACTGGCGCTTGGAGGGCGCTGGGGGGGGCAAATCCCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCTGTGCCCTCGGCTTATCCCTTa                                                               <  1:203623/147‑1 (MQ=14)
                                                                                       tcATCCTGTTGGGCTCCTGTCTTGTGCTTTGCTTAATATAATCACTGGCGAATGGAGCGCGCTGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTa                                                               <  1:194639/150‑1 (MQ=17)
                                                                                       acATCCTGTTGCGCGCCTGTCTTGTGCTTGGCTTCATATAAACACTGGCGCATGGAGCGCGCTGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTa                                                               <  1:203899/149‑1 (MQ=17)
                                                                                                           cATGTGCATTGCTTCCTATAATCCCGGGCGCACGGCGCGCGCAGGGGGCGGCCAAGCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGcgcg                                           <  1:727021/150‑1 (MQ=14)
                                                                                                           cATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGCAGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGcgcg                                           <  1:1667966/150‑1 (MQ=255)
                                                                                                           cATGTGCATTGATAAATATAAACACTGGCGCAAGGAGCGCGCAGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGcgcg                                           <  1:1111447/150‑1 (MQ=255)
                                                                                                                 cATTGCTTCATATAATCACTGGCGCTAGGAGCGCGCAGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTaa                                     <  1:1421959/150‑1 (MQ=255)
                                                                                                                 cATTGCTTCATATAATCACTGGCGCAAGGAGCGCGCAGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTaa                                     <  1:1421952/150‑1 (MQ=255)
                                                                                                                  aTTGCTTCATATAATCACTGGCGCAAGGAGCGCGCAGGGGGGGGCCCAAAGCCCCCCCCCCCCCCCCCCCCCGGCGCAGGCCCAAAAAAAACCCGCGTCGCCGTGCCCACCGCGGTAACCCTAAGGCAAACCGGTCGGGCCCGAGGTAAc                                    >  2:1267481/1‑150 (MQ=255)
                                                                                                                             taATCACTGGCGCAAGGAGCGCGCAGGGGGCCGCCAATACCCCCCGCCCCCAGCACCCCCCGGCTCTGGCGCACAAAAAACCCGCGGCCCGGTGCCCACCGCGTATCCCGTACGGCTACCCGGGCGGGCGCGAGGTAACAACCCCGTaaa                         >  2:2547922/1‑150 (MQ=2)
                                                                                                                              aaTCACTGGCGCAAGGAGCGCGCAGGGGGCGGCCAAACGCCCCCGCCCCCAGCACCCCCCGGCTCTGGCGAAAAAAAAAGCCGCTGCGCGGTGCCCTCGGCGTATACCTTACGGCAACCCGGTCGGGCGCCCGGTAACATCCCGGTaaaa                        >  2:1495109/1‑150 (MQ=11)
                                                                                                                                  aCTGGCGCAAGGAGCGCGCAGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAAcgcg                    <  2:1577176/150‑1 (MQ=32)
                                                                                                                                  aCTGGCGCAAGGAGCGCGCAGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAAcgcg                    <  1:652516/150‑1 (MQ=32)
                                                                                                                                  aCTGGCGCAAGGAGCGCGCAGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCGGGGCGGTGCCCACGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAAcgcg                    <  2:1343462/150‑1 (MQ=12)
                                                                                                                                    tGGCGCAAGGAGCGCGCAGGGGGCGGCCCAACGCCCCCGCCCCCGGCCCCCCCCGGCGCAGGCGAACAAAAAACCCCCTGCGCGGGGCCCCCCGCGTATACCGTAACGCTACCCGGGCGGGGGCGAGGGTACAACCCCGTAAAACCCGcc                  >  2:395072/1‑150 (MQ=2)
                                                                                                                                      gcgcAAGGAGCGCGCATGGGGCGGCCAAAAGCCGCCGCCCCCAGCACCCCCGGGCACTTGCGAACAAAAACGCCGCAGCGCGGTGCCCTAGGCTTATCCCATACGGCTACCGGGTCGGGCGCGAGGTAACATCCCAGTAAAAAGCGccca                >  2:1731551/1‑149 (MQ=11)
                                                                                                                                         cAAGGAGCGCGCAGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGCCCTCAg             <  1:596403/150‑1 (MQ=32)
                                                                                                                                         cAAGGAGCGCGCAGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGCCCTCAg             <  2:692974/150‑1 (MQ=32)
                                                                                                                                                 gcgcAGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGCCCTCAGCCCACATc     <  1:1624535/150‑1 (MQ=32)
                                                                                                                                                     cAGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGCCCTCAGCCCACATCCAt  <  1:1596417/150‑1 (MQ=32)
                                                                                                                                                      |                                                                                                                                                      
TACTGGTGGTGACTGGCGTATTTAAGATGATATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGCAGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGCCCTCAGCCCACATCCAT  >  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: