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

Only 100 of 121 total aligned reads displayed.

TGGTGGTGACTGGCGTATTTAAGATGATATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGCAGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGCCCTCAGCCCACATCCATG  >  NC_000913/2229031‑2229326
                                                                                                                                                   |                                                                                                                                                     
tggtggTGACTGGCGTATTTAAGATGATATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGCAggg                                                                                                                                                    <  1:3487506/150‑1 (MQ=255)
tggtggTGACTGGCGTAATTAAGATGATATTTTAAAATTAAATAATGACATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGCAggg                                                                                                                                                    <  2:2822836/150‑1 (MQ=255)
 ggtggtGACTGGCGTATTTAAGATGATATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGCCgggg                                                                                                                                                   >  1:1165052/1‑150 (MQ=255)
  gtggtgACTGGCGTATTTAAGATGATATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGCAggggg                                                                                                                                                  >  1:3778710/1‑150 (MQ=255)
  gtggtgACTGGCGTATTTAAGATGATATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGCAggggg                                                                                                                                                  >  1:793234/1‑150 (MQ=255)
   tggtgACTGGCGTATTTAAGATGATATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGTATATTTCAGTCTGTCTTCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCggggggggg                                                                                                                                                 >  2:238890/1‑149 (MQ=255)
   tggtgACTGGCGTATTTAAGATGATATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGCAGGGggc                                                                                                                                                 <  1:336247/150‑1 (MQ=255)
   tggtgACTGGCGTATTTAAGATGATATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGCAGGGggc                                                                                                                                                 <  1:3831435/150‑1 (MQ=255)
   tggtgACTGGCGTATTTAAGATGATATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGCAGGGggc                                                                                                                                                 <  1:3831526/150‑1 (MQ=255)
   tggtgACTGGCGTATTTAAGATGATATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATAGTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCggggggggg                                                                                                                                                 >  2:1890189/1‑149 (MQ=255)
    ggtgACTGGCGTATTTAAGATGATATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGCAGGGggcg                                                                                                                                                >  1:2313396/1‑150 (MQ=255)
     gtgACTGGCGTATTTAAGATGATATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGCAGGGggcgg                                                                                                                                               <  2:2451764/150‑1 (MQ=255)
     gtgACTGGCGTATTTAAGATGATATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGCAGGGggcgg                                                                                                                                               <  2:1769875/150‑1 (MQ=255)
     gtgACTGGCGTATTTAAGATGATATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGCAGGGggcgg                                                                                                                                               <  1:2068420/150‑1 (MQ=255)
      tgACTGGCGTATTTAAGATGATATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGCCGGGggcggg                                                                                                                                              >  1:2545903/1‑149 (MQ=255)
      tgACTGGCGTATTTAAGATGATATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGCCGGGggcggg                                                                                                                                              >  1:2545822/1‑149 (MQ=255)
        aCTGGCGTATTTAAGATGATATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGGGGGGggcgggcc                                                                                                                                            >  2:1215888/1‑147 (MQ=255)
         cTGGCGTATTTAAGATGATATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGCAGGGGGCGGCCaa                                                                                                                                           <  2:2460207/150‑1 (MQ=255)
          tGGCGTATTTAAGATGATATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGCAGGGGGCGGCCaaa                                                                                                                                          >  2:1318415/1‑149 (MQ=255)
          tGGCGTATTTAAGATGATATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGCAGGGGGCGGCCAAt                                                                                                                                          <  1:2057185/150‑1 (MQ=255)
          tGGCGTATTTAAGATGATATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGCAGGGGGCGGCCAAt                                                                                                                                          <  1:2056440/150‑1 (MQ=255)
          tGGCGTATTTAAGATGATATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGCAGGGGGCGGCCAAt                                                                                                                                          >  2:3615289/1‑150 (MQ=255)
          tGGCGTATTTAAGATGATATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTATCTCATCCTGTTGCGATCCTGTCATGTGCATAGCTTCATATAATAACTTGCGCAAGGAGCGCGCAAGGGGGGGGCCaa                                                                                                                                           >  2:365962/1‑150 (MQ=255)
           ggCGTATTTAAGATGATATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGCAGGGGGCGGCCAATc                                                                                                                                         <  1:1534156/150‑1 (MQ=255)
            gCGTATTTAAGATGATATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGgggggggggggcccacc                                                                                                                                        >  2:1836009/1‑143 (MQ=255)
            gCGTATTTAAGATGATATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGGCGGGGGGGGcccaaag                                                                                                                                        >  1:2508043/1‑145 (MQ=255)
              gTATTTAAGATGATATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGCAGGGGGCGGCCAATcgcc                                                                                                                                      <  2:682154/150‑1 (MQ=255)
                aTTTAAGATGATATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGCAGGGGGCGGCCAATcgccgc                                                                                                                                    <  1:2034652/150‑1 (MQ=255)
                aTTTAAGATGATATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGCAGGGGGCGGCCAATcgccgc                                                                                                                                    <  1:1129539/150‑1 (MQ=255)
                  ttAAGATGATATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGCAGGGGGCGGCCAAAcgccgccc                                                                                                                                  >  2:2627634/1‑149 (MQ=255)
                      tattatatTTTAAATTTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTAATATAATCACTGGCGCAAGGAGCGCGCTGGGGGCGGCCAATCGCCGCCGcccc                                                                                                                              <  1:3970425/146‑1 (MQ=255)
                      tatgatATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGCAGGGGGCGGCCAATCGCCGCCGcccc                                                                                                                              <  1:3970417/149‑1 (MQ=255)
                      gatgatATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGGCGGGGGCGGGCCAACGcccccccccc                                                                                                                              >  2:1847186/1‑150 (MQ=255)
                      gatgatATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGCCCGGGGGGGCCAAACGcccccccccc                                                                                                                              >  1:1279085/1‑150 (MQ=39)
                         gatATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGGCGGGGGCGGGCCAAAGCcccccccccccg                                                                                                                           >  1:1098246/1‑148 (MQ=255)
                         gatATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGCCGGGGGCGGCCAAAAgcccccgcccccag                                                                                                                           >  2:2083568/1‑148 (MQ=255)
                          atatTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGCAGGGGGCGGCCAAAAGCCGccccccccacc                                                                                                                          >  1:2702901/1‑147 (MQ=255)
                           ttttttAAAATTAATTATTGTCATCTGGTCCGAAAATAACGTTAATTTTTCTGTCTCTCTTCGTGTTGCGCTCCTTTATTGTGCTTTGCTTCATATAATCACTGGCGCATGGTGCGGGCTGGGGGCGGCCAATCGCCGCCGCCCCCTggg                                                                                                                         <  1:2099801/148‑3 (MQ=21)
                           tatTTTAAAATTAATTAATGTAATCAGGTCCCAAAATAACGAGAATATTACAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGCAGGGGGCGGCCAATCGCCGCCGCCCCCTGCa                                                                                                                         <  1:3601677/150‑1 (MQ=255)
                             ttttAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGCCGGGGGGGGCCCAACcccccccccccccccccc                                                                                                                       >  2:3183586/1‑144 (MQ=255)
                                  aaaTTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGCCGGGGGCGGGCCCTCGCCCCCCCCCCCAcccccccccg                                                                                                                  >  1:302138/1‑148 (MQ=38)
                                         ttaatGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGCAGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTgg                                                                                                           <  2:1816496/150‑1 (MQ=255)
                                             tGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGTGCGCGCAGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGaa                                                                                                       <  2:1447020/150‑1 (MQ=255)
                                             tGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGCAGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGaa                                                                                                       <  2:1483714/150‑1 (MQ=255)
                                              gTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGCAGGGGGCGGCCAAACGCCCCCGCCCCCACCACCCCCCGGCCAAGGCGAAc                                                                                                      >  2:2049454/1‑150 (MQ=37)
                                                 atcaGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGCAGGGGGCGGCCAATCGCCGCCGCCCCCAGCACCCCCCGGCGCAGGCGAACaaa                                                                                                   >  2:3775858/1‑150 (MQ=255)
                                                    aGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGCAGGGGGCGGCCAATAGCCCCCGCCCCCAGCACCCCCGGGCTCAGGCGGAAAAAATc                                                                                                >  2:1002669/1‑150 (MQ=25)
                                                                   gagaATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGCCGGGGGCGGCCAATAGCCCCCGCCCCCAGCACCCCCCGGCTCAGGCGAACAAAAACGCCCCGGCGCGGTGc                                                                                 >  2:2079169/1‑150 (MQ=21)
                                                                   gagaATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGCAGGGGGCGGCCAATAGCCCCCCCCCCCAGCACCCCCCGGCAATGGCGAACAAAAACCCCCCAGCGCGGTGc                                                                                 >  2:2079143/1‑150 (MQ=21)
                                                                                  tctcATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGCAGGGGGCGGCCCCAAGCCGCCGCCCCCAGCCCCCCCCCGCTCGGGCGCACAAAAAAGCCGCCGCGCGGGGGCCGCGGCGTATcccc                                                                  >  2:1235576/1‑149 (MQ=18)
                                                                                  tctcATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGCAGGGGGCGGCCAAGCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCt                                                                  <  2:1165052/150‑1 (MQ=255)
                                                                                  tctcATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGCACGGGGCGGCCAAAAGCCCCCCCCCCCAGCCCCCCCCCGCTCAGGCGAACAAAAAAGCCCCTGCGCGGTGCCCACGGCGTATccca                                                                  >  2:1338482/1‑149 (MQ=18)
                                                                                    tcATCGTGTTGCGCTCCTTTCTTGTGCTTTGATTAATATAATCAGTGGCGCATGGAGCGCGCTGGGGGCGGCCATTCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTa                                                                <  1:11117/150‑1 (MQ=14)
                                                                                    tcATCCTGTTGGGCTCCTGTCATGTGCTTTGCTTCATATAATCAGTGGCGCATGGAGCGCGCAGGGGGGGGCCAATCGCCGCGGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTa                                                                <  2:1048217/150‑1 (MQ=17)
                                                                                    tcATCCTGTTGCGCTCCTGTCTGGTGCATTGCTTCATAAAAACACGGGCGCACGGCGCGCGCCGGGGGCGGCCAAGCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTa                                                                <  1:11273/150‑1 (MQ=17)
                                                                                    tcATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGCAGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTa                                                                <  2:3294603/150‑1 (MQ=255)
                                                                                    tcATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGCAGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTa                                                                <  2:1222559/150‑1 (MQ=255)
                                                                                          tGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGCACGGGGCGGCCCAACCCCGCCCCCCCCAGCACCCCCCGGCTCTGGCGGACAAAAAAGCCCCGGCGCGGTGCCCTCGGCGTATCCCATAAGGCTa                                                          >  2:3344427/1‑150 (MQ=17)
                                                                                            ttGCGCTCCTGTCATGTGCATTGCTTCATATAATCAGTGGCGCAAGGAGCGCGCCGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTAcc                                                        <  2:1888559/150‑1 (MQ=255)
                                                                                                gcTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGCTGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGt                                                    <  1:1190052/150‑1 (MQ=255)
                                                                                                gcTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGCAGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGt                                                    <  1:4026018/150‑1 (MQ=255)
                                                                                                    gtGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGCAGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCggg                                                <  2:3865411/149‑1 (MQ=255)
                                                                                                       tCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGCAGGGGGCGGCCCAACGCCGCCGCCCCCAGCAACCCCCGGCTCATGCGGACAAAAACGCCGCGGCCCGGTGCCCACCGCGTATACCATACCGCGACCCGGGCGGgcgc                                             >  2:641159/1‑150 (MQ=12)
                                                                                                        cATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGCAGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGcgcg                                            <  1:2329487/150‑1 (MQ=255)
                                                                                                         aTGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGCAGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGa                                           <  2:3778710/150‑1 (MQ=255)
                                                                                                         aTGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGCAGGGGGCGGCCAAACGCCGCCGCCCCCAGCACCCCCCGGCTCAGGCGAACAAAATAGCCCCTGCGCGGTGCCCTCGGCGTATACCATACGGCTACCCGGTCGGGCGCGa                                           >  2:3423724/1‑150 (MQ=17)
                                                                                                           gtgCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGCAGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAgg                                         <  1:1409930/150‑1 (MQ=255)
                                                                                                           gtgCATTGCTTAATATAATCACTGGCGCAAGGAGCGCGCAGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAgg                                         <  1:3864807/150‑1 (MQ=255)
                                                                                                            tgCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGCAGGGGGCGGCCAAACGCCCCCCCCCCCCCCACCCCCCCGCGCGGGCGAACAAAAACGCCCCGGCGCGGGGCCCGCCGCGGATACCGTACGGCGACCCGGGCGGgcgcccggg                                        >  2:3315888/1‑145 (MQ=9)
                                                                                                             gCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGCAGGGGGCGGCCAATCGCCGCCGCCCCCAGCACCCCCCGGCGCTGGCGAACAAAAACGCCGCGGCGCGGGGCCCACGGCGTAAACCCTACGGCTACCGGGTCGGGCGCGGGGTa                                       >  2:347240/1‑150 (MQ=14)
                                                                                                              cATTGCTTCATATAATCACTGGCGCAAGGAGCGCGCAGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTaa                                      <  2:793234/150‑1 (MQ=255)
                                                                                                              cATTGCTTCATATAATCACTGGCGCAAGGAGCGCGCAGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTaa                                      <  2:2435922/150‑1 (MQ=255)
                                                                                                                 tGCTTCATATAATCACTGGCGCACGGAGCGCGCCGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACAt                                   <  1:1211712/150‑1 (MQ=255)
                                                                                                                   cTTCATATAATCACTGGCGCCAGGAGCGCGCAGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACACAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATcc                                 <  2:3891565/150‑1 (MQ=255)
                                                                                                                   cTTCATATAATCACTGGCGCAAGGAGCGCGCAGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATcc                                 <  2:162083/150‑1 (MQ=40)
                                                                                                                   cTTCATATAATCACTGGCGCAAGGAGCGCGCAGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATcc                                 <  2:3675979/150‑1 (MQ=40)
                                                                                                                     taatatAATCACTGGCGCACGGTGCGCGCTGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCt                               <  1:3823971/148‑1 (MQ=255)
                                                                                                                     tCATATAATCACTGGCGCAAGGAGCGCGCAGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCt                               <  1:1151994/150‑1 (MQ=39)
                                                                                                                     tCATATAATCACTGGCGCAAGGAGCGCGCAGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCt                               <  1:2977541/150‑1 (MQ=39)
                                                                                                                     tCATATAATCACTGGCGCAAGGAGCGCGCAGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCt                               <  1:1836009/150‑1 (MQ=39)
                                                                                                                     tCATATAATCACTGGCGCAAGGAGCGCGCAGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCt                               <  2:3892701/150‑1 (MQ=39)
                                                                                                                         ataATCACTGGCGCATGGAGCGCGCCGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTaa                           <  1:2947827/150‑1 (MQ=25)
                                                                                                                         ataATCACTGGCGCAAGGAGCGCGCTGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTaa                           <  1:3183586/150‑1 (MQ=25)
                                                                                                                              cACTGGCGCAAGGAGCGCGCAGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAAcgc                      <  2:1044021/150‑1 (MQ=32)
                                                                                                                              cACTGGCGCAAGGAGCGCGCAGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAAcgc                      <  2:1391447/150‑1 (MQ=32)
                                                                                                                               ccTGGCGAATGTAGCGCGCTGGGGGCGGCCAATCGCCGCCGCCCCATGCACCCCCGGGCTCTGGCTAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGATACCTGGTCGGGCGCGAGGTAACATCCCTGTAAAAcgcg                     <  1:1479110/149‑1 (MQ=11)
                                                                                                                                  ggCGCAAGGAGCGCGCAGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGccc                  <  1:3081482/150‑1 (MQ=32)
                                                                                                                                   gcgcAAGGAGCGCGCAGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGCCCt                 <  1:3456348/150‑1 (MQ=32)
                                                                                                                                   gcgcAAGGAGCGCGCAGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGCCCt                 <  1:3456336/150‑1 (MQ=32)
                                                                                                                                    cgcAAGGAGCGCGCAGGGGGCGGCCAAAAGCCCCCGCCCCCAGCACCCCCCGGCGCTGGCGGACAAAAACCCCGCGGCGCGGTGCCCACGGCGTATCCCTTAAGGCTACCCGGGCGGGCGCGGGGGAACAAACCAGTAAAAAGCGcccac                >  2:267458/1‑148 (MQ=2)
                                                                                                                                     gcAAGGAGCGCGCAGGGGGCGGCCAAACGCCGCCGCCCCCGGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCGGCGCGGTGCCCGCGGCGTATACCGTAAGGCTACCGGGGCGGGCGCGAGGGAACATCCCAGTAAAACGCGCCCCCa               >  2:3100832/1‑150 (MQ=11)
                                                                                                                                       aaGGAGCGCGCTGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGCCCTCAGc             <  2:810953/150‑1 (MQ=17)
                                                                                                                                        aGGAGCGCGCAGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGCCCTCAGcc            <  1:3969373/150‑1 (MQ=32)
                                                                                                                                        aGGAGCGCGCAGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGCCCTCAGcc            <  1:261961/150‑1 (MQ=32)
                                                                                                                                        aGGAGCGCGCAGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGCCCTCAGcc            <  1:261833/150‑1 (MQ=32)
                                                                                                                                          gAGCGCGCAGGGGGCGGCCAAGCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGCCCTCAGCCca          <  1:1086373/150‑1 (MQ=17)
                                                                                                                                             cgcgcAGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGCCCTCAGCCCACAt       <  1:1918210/150‑1 (MQ=32)
                                                                                                                                                 cAGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGCCCTCAGCCCACATCCAt   <  2:2536364/150‑1 (MQ=32)
                                                                                                                                                 cAGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGCCCTCAGCCCACATCCAt   <  1:1472613/150‑1 (MQ=32)
                                                                                                                                                   aGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGCCCTCAGCCCACCTCCAtg  >  2:3570428/1‑150 (MQ=17)
                                                                                                                                                   |                                                                                                                                                     
TGGTGGTGACTGGCGTATTTAAGATGATATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGCAGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGCCCTCAGCCCACATCCATG  >  NC_000913/2229031‑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: