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.

TGGTGGTGACTGGCGTATTTAAGATGATATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGC‑AGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGCCCTCAGCCCACATCCATG  >  NC_000913/2229031‑2229326
                                                                                                                                                   |                                                                                                                                                     
tggtggTGACTGGCGTATTTAAGATGATATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGC‑Aggg                                                                                                                                                    <  1:3487506/150‑1 (MQ=255)
tggtggTGACTGGCGTAATTAAGATGATATTTTAAAATTAAATAATGACATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGC‑Aggg                                                                                                                                                    <  2:2822836/150‑1 (MQ=255)
 ggtggtGACTGGCGTATTTAAGATGATATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGC‑Cgggg                                                                                                                                                   >  1:1165052/1‑150 (MQ=255)
  gtggtgACTGGCGTATTTAAGATGATATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGC‑Aggggg                                                                                                                                                  >  1:3778710/1‑150 (MQ=255)
  gtggtgACTGGCGTATTTAAGATGATATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGC‑Aggggg                                                                                                                                                  >  1:793234/1‑150 (MQ=255)
   tggtgACTGGCGTATTTAAGATGATATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGTATATTTCAGTCTGTCTTCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCgg‑ggggggg                                                                                                                                                 >  2:238890/1‑149 (MQ=255)
   tggtgACTGGCGTATTTAAGATGATATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGC‑AGGGggc                                                                                                                                                 <  1:336247/150‑1 (MQ=255)
   tggtgACTGGCGTATTTAAGATGATATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGC‑AGGGggc                                                                                                                                                 <  1:3831435/150‑1 (MQ=255)
   tggtgACTGGCGTATTTAAGATGATATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGC‑AGGGggc                                                                                                                                                 <  1:3831526/150‑1 (MQ=255)
   tggtgACTGGCGTATTTAAGATGATATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATAGTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCgg‑ggggggg                                                                                                                                                 >  2:1890189/1‑149 (MQ=255)
    ggtgACTGGCGTATTTAAGATGATATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGC‑AGGGggcg                                                                                                                                                >  1:2313396/1‑150 (MQ=255)
     gtgACTGGCGTATTTAAGATGATATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGC‑AGGGggcgg                                                                                                                                               <  2:2451764/150‑1 (MQ=255)
     gtgACTGGCGTATTTAAGATGATATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGC‑AGGGggcgg                                                                                                                                               <  2:1769875/150‑1 (MQ=255)
     gtgACTGGCGTATTTAAGATGATATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGC‑AGGGggcgg                                                                                                                                               <  1:2068420/150‑1 (MQ=255)
      tgACTGGCGTATTTAAGATGATATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGC‑CGGGggcggg                                                                                                                                              >  1:2545903/1‑149 (MQ=255)
      tgACTGGCGTATTTAAGATGATATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGC‑CGGGggcggg                                                                                                                                              >  1:2545822/1‑149 (MQ=255)
        aCTGGCGTATTTAAGATGATATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGG‑GGGGggcgggcc                                                                                                                                            >  2:1215888/1‑147 (MQ=255)
         cTGGCGTATTTAAGATGATATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGC‑AGGGGGCGGCCaa                                                                                                                                           <  2:2460207/150‑1 (MQ=255)
          tGGCGTATTTAAGATGATATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGC‑AGGGGGCGGCCaaa                                                                                                                                          >  2:1318415/1‑149 (MQ=255)
          tGGCGTATTTAAGATGATATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGC‑AGGGGGCGGCCAAt                                                                                                                                          <  1:2057185/150‑1 (MQ=255)
          tGGCGTATTTAAGATGATATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGC‑AGGGGGCGGCCAAt                                                                                                                                          <  1:2056440/150‑1 (MQ=255)
          tGGCGTATTTAAGATGATATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGC‑AGGGGGCGGCCAAt                                                                                                                                          >  2:3615289/1‑150 (MQ=255)
          tGGCGTATTTAAGATGATATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTATCTCATCCTGTTGCGATCCTGTCATGTGCATAGCTTCATATAATAACTTGCGCAAGGAGCGCGCAAGGGGGGGGCCaa                                                                                                                                           >  2:365962/1‑150 (MQ=255)
           ggCGTATTTAAGATGATATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGC‑AGGGGGCGGCCAATc                                                                                                                                         <  1:1534156/150‑1 (MQ=255)
            gCGTATTTAAGATGATATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGg‑ggggggggggcccacc                                                                                                                                        >  2:1836009/1‑143 (MQ=255)
            gCGTATTTAAGATGATATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGG‑CGGGGGGGGcccaaag                                                                                                                                        >  1:2508043/1‑145 (MQ=255)
              gTATTTAAGATGATATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGC‑AGGGGGCGGCCAATcgcc                                                                                                                                      <  2:682154/150‑1 (MQ=255)
                aTTTAAGATGATATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGC‑AGGGGGCGGCCAATcgccgc                                                                                                                                    <  1:2034652/150‑1 (MQ=255)
                aTTTAAGATGATATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGC‑AGGGGGCGGCCAATcgccgc                                                                                                                                    <  1:1129539/150‑1 (MQ=255)
                  ttAAGATGATATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGC‑AGGGGGCGGCCAAAcgccgccc                                                                                                                                  >  2:2627634/1‑149 (MQ=255)
                      tattatatTTTAAATTTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTAATATAATCACTGGCGCAAGGAGCGCGC‑TGGGGGCGGCCAATCGCCGCCGcccc                                                                                                                              <  1:3970425/146‑1 (MQ=255)
                      tatgatATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGC‑AGGGGGCGGCCAATCGCCGCCGcccc                                                                                                                              <  1:3970417/149‑1 (MQ=255)
                      gatgatATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGG‑CGGGGGCGGGCCAACGcccccccccc                                                                                                                              >  2:1847186/1‑150 (MQ=255)
                      gatgatATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGC‑CCGGGGGGGCCAAACGcccccccccc                                                                                                                              >  1:1279085/1‑150 (MQ=39)
                         gatATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGG‑CGGGGGCGGGCCAAAGCcccccccccccg                                                                                                                           >  1:1098246/1‑148 (MQ=255)
                         gatATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGC‑CGGGGGCGGCCAAAAgcccccgcccccag                                                                                                                           >  2:2083568/1‑148 (MQ=255)
                          atatTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGC‑AGGGGGCGGCCAAAAGCCGccccccccacc                                                                                                                          >  1:2702901/1‑147 (MQ=255)
                           ttttttAAAATTAATTATTGTCATCTGGTCCGAAAATAACGTTAATTTTTCTGTCTCTCTTCGTGTTGCGCTCCTTTATTGTGCTTTGCTTCATATAATCACTGGCGCATGGTGCGGGC‑TGGGGGCGGCCAATCGCCGCCGCCCCCTggg                                                                                                                         <  1:2099801/148‑3 (MQ=21)
                           tatTTTAAAATTAATTAATGTAATCAGGTCCCAAAATAACGAGAATATTACAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGC‑AGGGGGCGGCCAATCGCCGCCGCCCCCTGCa                                                                                                                         <  1:3601677/150‑1 (MQ=255)
                             ttttAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGC‑CGGGGGGGGCCCAACcccccccccccccccccc                                                                                                                       >  2:3183586/1‑144 (MQ=255)
                                  aaaTTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGC‑CGGGGGCGGGCCCTCGCCCCCCCCCCCAcccccccccg                                                                                                                  >  1:302138/1‑148 (MQ=38)
                                         ttaatGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGC‑AGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTgg                                                                                                           <  2:1816496/150‑1 (MQ=255)
                                             tGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGTGCGCGC‑AGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGaa                                                                                                       <  2:1447020/150‑1 (MQ=255)
                                             tGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGC‑AGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGaa                                                                                                       <  2:1483714/150‑1 (MQ=255)
                                              gTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGC‑AGGGGGCGGCCAAACGCCCCCGCCCCCACCACCCCCCGGCCAAGGCGAAc                                                                                                      >  2:2049454/1‑150 (MQ=37)
                                                 atcaGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGC‑AGGGGGCGGCCAATCGCCGCCGCCCCCAGCACCCCCCGGCGCAGGCGAACaaa                                                                                                   >  2:3775858/1‑150 (MQ=255)
                                                    aGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGC‑AGGGGGCGGCCAATAGCCCCCGCCCCCAGCACCCCCGGGCTCAGGCGGAAAAAATc                                                                                                >  2:1002669/1‑150 (MQ=25)
                                                                   gagaATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGC‑CGGGGGCGGCCAATAGCCCCCGCCCCCAGCACCCCCCGGCTCAGGCGAACAAAAACGCCCCGGCGCGGTGc                                                                                 >  2:2079169/1‑150 (MQ=21)
                                                                   gagaATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGC‑AGGGGGCGGCCAATAGCCCCCCCCCCCAGCACCCCCCGGCAATGGCGAACAAAAACCCCCCAGCGCGGTGc                                                                                 >  2:2079143/1‑150 (MQ=21)
                                                                                  tctcATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGC‑AGGGGGCGGCCCCAAGCCGCCGCCCCCAGCCCCCCCCCGCTCGGGCGCACAAAAAAGCCGCCGCGCGGGGGCCGCGGCGTATcccc                                                                  >  2:1235576/1‑149 (MQ=18)
                                                                                  tctcATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGC‑AGGGGGCGGCCAAGCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCt                                                                  <  2:1165052/150‑1 (MQ=255)
                                                                                  tctcATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGC‑ACGGGGCGGCCAAAAGCCCCCCCCCCCAGCCCCCCCCCGCTCAGGCGAACAAAAAAGCCCCTGCGCGGTGCCCACGGCGTATccca                                                                  >  2:1338482/1‑149 (MQ=18)
                                                                                    tcATCGTGTTGCGCTCCTTTCTTGTGCTTTGATTAATATAATCAGTGGCGCATGGAGCGCGC‑TGGGGGCGGCCATTCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTa                                                                <  1:11117/150‑1 (MQ=14)
                                                                                    tcATCCTGTTGGGCTCCTGTCATGTGCTTTGCTTCATATAATCAGTGGCGCATGGAGCGCGC‑AGGGGGGGGCCAATCGCCGCGGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTa                                                                <  2:1048217/150‑1 (MQ=17)
                                                                                    tcATCCTGTTGCGCTCCTGTCTGGTGCATTGCTTCATAAAAACACGGGCGCACGGCGCGCGC‑CGGGGGCGGCCAAGCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTa                                                                <  1:11273/150‑1 (MQ=17)
                                                                                    tcATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGC‑AGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTa                                                                <  2:3294603/150‑1 (MQ=255)
                                                                                    tcATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGC‑AGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTa                                                                <  2:1222559/150‑1 (MQ=255)
                                                                                          tGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGC‑ACGGGGCGGCCCAACCCCGCCCCCCCCAGCACCCCCCGGCTCTGGCGGACAAAAAAGCCCCGGCGCGGTGCCCTCGGCGTATCCCATAAGGCTa                                                          >  2:3344427/1‑150 (MQ=17)
                                                                                            ttGCGCTCCTGTCATGTGCATTGCTTCATATAATCAGTGGCGCAAGGAGCGCGC‑CGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTAcc                                                        <  2:1888559/150‑1 (MQ=255)
                                                                                                gcTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGC‑TGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGt                                                    <  1:1190052/150‑1 (MQ=255)
                                                                                                gcTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGC‑AGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGt                                                    <  1:4026018/150‑1 (MQ=255)
                                                                                                    gtGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGC‑AGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCggg                                                <  2:3865411/149‑1 (MQ=255)
                                                                                                       tCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGC‑AGGGGGCGGCCCAACGCCGCCGCCCCCAGCAACCCCCGGCTCATGCGGACAAAAACGCCGCGGCCCGGTGCCCACCGCGTATACCATACCGCGACCCGGGCGGgcgc                                             >  2:641159/1‑150 (MQ=12)
                                                                                                        cATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGC‑AGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGcgcg                                            <  1:2329487/150‑1 (MQ=255)
                                                                                                         aTGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGC‑AGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGa                                           <  2:3778710/150‑1 (MQ=255)
                                                                                                         aTGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGC‑AGGGGGCGGCCAAACGCCGCCGCCCCCAGCACCCCCCGGCTCAGGCGAACAAAATAGCCCCTGCGCGGTGCCCTCGGCGTATACCATACGGCTACCCGGTCGGGCGCGa                                           >  2:3423724/1‑150 (MQ=17)
                                                                                                           gtgCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGC‑AGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAgg                                         <  1:1409930/150‑1 (MQ=255)
                                                                                                           gtgCATTGCTTAATATAATCACTGGCGCAAGGAGCGCGC‑AGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAgg                                         <  1:3864807/150‑1 (MQ=255)
                                                                                                            tgCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGC‑AGGGGGCGGCCAAACGCCCCCCCCCCCCCCACCCCCCCGCGCGGGCGAACAAAAACGCCCCGGCGCGGGGCCCGCCGCGGATACCGTACGGCGACCCGGGCGGgcgcccggg                                        >  2:3315888/1‑145 (MQ=9)
                                                                                                             gCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGC‑AGGGGGCGGCCAATCGCCGCCGCCCCCAGCACCCCCCGGCGCTGGCGAACAAAAACGCCGCGGCGCGGGGCCCACGGCGTAAACCCTACGGCTACCGGGTCGGGCGCGGGGTa                                       >  2:347240/1‑150 (MQ=14)
                                                                                                              cATTGCTTCATATAATCACTGGCGCAAGGAGCGCGC‑AGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTaa                                      <  2:793234/150‑1 (MQ=255)
                                                                                                              cATTGCTTCATATAATCACTGGCGCAAGGAGCGCGC‑AGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTaa                                      <  2:2435922/150‑1 (MQ=255)
                                                                                                                 tGCTTCATATAATCACTGGCGCACGGAGCGCGC‑CGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACAt                                   <  1:1211712/150‑1 (MQ=255)
                                                                                                                   cTTCATATAATCACTGGCGCCAGGAGCGCGC‑AGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACACAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATcc                                 <  2:3891565/150‑1 (MQ=255)
                                                                                                                   cTTCATATAATCACTGGCGCAAGGAGCGCGC‑AGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATcc                                 <  2:162083/150‑1 (MQ=40)
                                                                                                                   cTTCATATAATCACTGGCGCAAGGAGCGCGC‑AGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATcc                                 <  2:3675979/150‑1 (MQ=40)
                                                                                                                     taatatAATCACTGGCGCACGGTGCGCGC‑TGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCt                               <  1:3823971/148‑1 (MQ=255)
                                                                                                                     tCATATAATCACTGGCGCAAGGAGCGCGC‑AGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCt                               <  1:1151994/150‑1 (MQ=39)
                                                                                                                     tCATATAATCACTGGCGCAAGGAGCGCGC‑AGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCt                               <  1:2977541/150‑1 (MQ=39)
                                                                                                                     tCATATAATCACTGGCGCAAGGAGCGCGC‑AGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCt                               <  1:1836009/150‑1 (MQ=39)
                                                                                                                     tCATATAATCACTGGCGCAAGGAGCGCGC‑AGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCt                               <  2:3892701/150‑1 (MQ=39)
                                                                                                                         ataATCACTGGCGCATGGAGCGCGC‑CGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTaa                           <  1:2947827/150‑1 (MQ=25)
                                                                                                                         ataATCACTGGCGCAAGGAGCGCGC‑TGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTaa                           <  1:3183586/150‑1 (MQ=25)
                                                                                                                              cACTGGCGCAAGGAGCGCGC‑AGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAAcgc                      <  2:1044021/150‑1 (MQ=32)
                                                                                                                              cACTGGCGCAAGGAGCGCGC‑AGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAAcgc                      <  2:1391447/150‑1 (MQ=32)
                                                                                                                               ccTGGCGAATGTAGCGCGC‑TGGGGGCGGCCAATCGCCGCCGCCCCATGCACCCCCGGGCTCTGGCTAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGATACCTGGTCGGGCGCGAGGTAACATCCCTGTAAAAcgcg                     <  1:1479110/149‑1 (MQ=11)
                                                                                                                                  ggCGCAAGGAGCGCGC‑AGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGccc                  <  1:3081482/150‑1 (MQ=32)
                                                                                                                                   gcgcAAGGAGCGCGC‑AGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGCCCt                 <  1:3456348/150‑1 (MQ=32)
                                                                                                                                   gcgcAAGGAGCGCGC‑AGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGCCCt                 <  1:3456336/150‑1 (MQ=32)
                                                                                                                                    cgcAAGGAGCGCGC‑AGGGGGCGGCCAAAAGCCCCCGCCCCCAGCACCCCCCGGCGCTGGCGGACAAAAACCCCGCGGCGCGGTGCCCACGGCGTATCCCTTAAGGCTACCCGGGCGGGCGCGGGGGAACAAACCAGTAAAAAGCGcccac                >  2:267458/1‑148 (MQ=2)
                                                                                                                                     gcAAGGAGCGCGC‑AGGGGGCGGCCAAACGCCGCCGCCCCCGGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCGGCGCGGTGCCCGCGGCGTATACCGTAAGGCTACCGGGGCGGGCGCGAGGGAACATCCCAGTAAAACGCGCCCCCa               >  2:3100832/1‑150 (MQ=11)
                                                                                                                                       aaGGAGCGCGC‑TGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGCCCTCAGc             <  2:810953/150‑1 (MQ=17)
                                                                                                                                        aGGAGCGCGC‑AGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGCCCTCAGcc            <  1:3969373/150‑1 (MQ=32)
                                                                                                                                        aGGAGCGCGC‑AGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGCCCTCAGcc            <  1:261961/150‑1 (MQ=32)
                                                                                                                                        aGGAGCGCGC‑AGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGCCCTCAGcc            <  1:261833/150‑1 (MQ=32)
                                                                                                                                          gAGCGCGC‑AGGGGGCGGCCAAGCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGCCCTCAGCCca          <  1:1086373/150‑1 (MQ=17)
                                                                                                                                             cgcgc‑AGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGCCCTCAGCCCACAt       <  1:1918210/150‑1 (MQ=32)
                                                                                                                                                 c‑AGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGCCCTCAGCCCACATCCAt   <  2:2536364/150‑1 (MQ=32)
                                                                                                                                                 c‑AGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGCCCTCAGCCCACATCCAt   <  1:1472613/150‑1 (MQ=32)
                                                                                                                                                   aGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGCCCTCAGCCCACCTCCAtg  >  2:3570428/1‑150 (MQ=17)
                                                                                                                                                   |                                                                                                                                                     
TGGTGGTGACTGGCGTATTTAAGATGATATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGC‑AGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGCCCTCAGCCCACATCCATG  >  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: ‑