Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,054,900 A→G 16.8% intergenic (+40/‑61) quiA → / → oprB‑III quinate dehydrogenase/carbohydrate‑selective porin

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,054,9000AG16.8% 441.6 / 77.1 208intergenic (+40/‑61)quiA/oprB‑IIIquinate dehydrogenase/carbohydrate‑selective porin
Reads supporting (aligned to +/- strand):  ref base A (66/92);  new base G (8/26);  total (74/134)
Fisher's exact test for biased strand distribution p-value = 5.34e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.14e-02

Only 100 of 218 total aligned reads displayed.

TGACCTACGTCTCGCCGAAGACCGGCAAGCAGTACGTGGTGATCACTGCTGGCGGCGCCCGCCAGTCGACTGACCGTGGCGACTATGTGATTTCTTACGCTCTGCCGTAAAACCGCGTCGCCTGCTTCGCGGGCAAGCCCGCTCCTACAAGGCCCCTGTAGGAGCGGGCTTGCCCGCGAATGCCCCTCAAAGACAACCCGAGATTCTGTAATGCCATCCGCAATTCGCATCACCCCCACCCTGCTGCTGGCCCTGGCCAGCACCACTGCCCTTGCCAGCGACGATCTCATGACCCGCA  >  NC_002947/4054751‑4055048
                                                                                                                                                     |                                                                                                                                                    
tGACCTACGTCTCGCCGAAGACCGGCAAGCAGTACGTGGTGATCACTGCTGGCGGCGCCCGCCAGTCGACTGACCGTGGCGACTATGTGATTTCTTACGCTCTGCCGTAAAACCGCGTCGCCTGCTTCGCGGGCAAGCCCGCTCCTACaa                                                                                                                                                      <  1:1117899/150‑1 (MQ=255)
   ccTACGTCTCGCCGAAGACCGGCAAGCAGTACGTGGTGATCACTGCTGGCGGCGCCCGCCAGTCGACTGACCGTGGCGACTATGTGATTTCTTACGCTCTGCCGTAAAACCGCGTCGCCTGCTTCGCGGGCAAGCCCGCTCCTACAAGGc                                                                                                                                                   >  2:4033644/1‑150 (MQ=255)
     tACGTCTCGCCGAAGACCGGCAAGCAGTACGTGGTGATCACTGCTGGCGGCGCCCGCCAGTCGATTGACCGTGGCGACTATGTGATTTCTTACGCTCTGCCGTAAAACCGCGTCGCCTGCTTCGCGGGCAAGCCCGCTCCTACAAGGccc                                                                                                                                                 <  2:1523062/150‑1 (MQ=255)
     tACGTCTCGCCGAAGACCGGCAAGCAGTACGTGGTGATCACTGCTGGCGGCGCCCGCCAGTCGACTGACCGTGGCGACTATGTGATTTCTTACGCTCTGCCGTAAAACCGCGTCGCCTGCTTCGCGGGCAAGCCCGCTCCTACAGGGgcc                                                                                                                                                 <  2:2622049/150‑1 (MQ=255)
     tACGTCTCGCCGAAGACCGGCAAGCAGTACGTGGTGATCACTGCTGGCGGCGCCCGCCAGTCGACTGACCGTGGCGACTATGTGATTTCTTACGCTCTGCCGTAAAACCGCGTCGCCTGCTTCGCGGGCAAGCCCGCTACTACAAGGccc                                                                                                                                                 >  1:5040833/1‑150 (MQ=255)
       cGTCTCGCCGAAGACCGGCAAGCAGTACGTGGTGATCACTGCTGGCGGCGCCCGCCAGTCGACTGACCGTGGCGACTATGTGATTTCTTACGCTCTGCCGTAAAACCGCGTCGCCTGCTTCGCGGGCAAGCCCGCTCCTACAGGGgcctt                                                                                                                                               <  2:5188344/150‑3 (MQ=255)
       cGTCTCGCCGAAGACCGGCAAGCAGTACGTGGTGATCACTGCTGGCGGCGCCCGCCAGTCGACTGACCGTGGCGACTATGTGATTTCTTACGCTCTGCCGTAAAACCGCGTCGCCTGCTTCGCGGGCAAGCCCGCTCCTACAGGGgcctt                                                                                                                                               <  2:1753767/150‑3 (MQ=255)
       cGTCTCGCCGAAGACCGGCAAGCAGTACGTGGTGATCACTGCTGGCGGCGCCCGCCAGTCGACTGACCGTGGCGACTATGTGATTTCTTACGCTCTGCCGTAAAACCGCGTCGCCTGCTTCGCGGGCAAGCCCGCTCCTACAAGGCCCCt                                                                                                                                               <  1:2010657/150‑1 (MQ=255)
         tctcGCCGAAGACCGGCAAGCAGTACGTGGTGATCACTGCTGGCGGCGCCCGCCAGTCGACTGACCGTGGCGACTATGTGATTTCTTACGCTCTGCCGTAAAACCGCGTCGCCTGCTTCGCGGGCAAGCCCGCTCCTACAGGGGCCTTGt                                                                                                                                             <  2:38410/150‑1 (MQ=255)
         tctcGCCGAAGACCGGCAAGCAGTACGTGGTGATCACTGCTGGCGGCGCCCGCCAGTCGACTGACCGTGGCGACTATGTGATTTCTTACGCTCTGCCGTAAAACCGCGTCGCCTGCTTCGCGGGCAAGCCCGCTCCTACAGGGGCCTTGt                                                                                                                                             <  1:1184015/150‑1 (MQ=255)
         tctcGCCGAAGACCGGCAAGCAGTACGTGGTGATCACTGCTGGCGGCGCCCGCCAGTCGACTGACCGTGGCGACTATGTGATTTCTTACGCTCTGCCGTAAAACCGCGTCGCCTGCTTCGCGGGCAAGCCCGCTCCTACAAGGCCCCTGt                                                                                                                                             <  2:4805504/150‑1 (MQ=255)
         tctcGCCGAAGACCGGCAAGCAGTACGTGGTGATCACTGCTGGCGGCGCCCGCCAGTCGACTGACCGTGGCGACTATGTGATTTCTTACGCTCTGCCGTAAAACCGCGTCGCCTGCTTCGCGGGCAAGCCCGCTCCTACAAGGCCCCTGt                                                                                                                                             <  2:1239606/150‑1 (MQ=255)
         tctcGCCGAAGACCGGCAAGCAGTACGTGGTGATCACTGCTGGCGGCGCCCGCCAGTCGACTGACCGTGGCGACTATGTGATTTCTTACGCTCTGCCGTAAAACCGCGTCGCCTGCTTCGCGGGCAAGCCCGCTCCTACAAGGCCCCTGt                                                                                                                                             <  2:4087668/150‑1 (MQ=255)
          ctcGCCGAAGACCGGCAAGCAGTACGTGGTGATCACTGCTGGCGGCGCCCGCCAGTCGACTGACCGTGGCGACTATGTGATTTCTTACGCTCTGCCGTAAAACCGCGTCGCCTGCTTCGCGGGCAAGCCCGCTCCTACAAGGCCCCTGTa                                                                                                                                            <  1:2031521/150‑1 (MQ=255)
            cGCCGAAGACCGGCAAGCAGTACGTGGTGATCACTGCTGGCGGCGCCCGCCAGTCGACTGACCGTGGCGACTATGTGATTTCTTACGCTCTGCCGTAAAACCGCGTCGCCTGCTTTGCGGGCAAGCCCGCTCCTACAGGGGCCTTGTAgg                                                                                                                                          <  2:3421052/150‑1 (MQ=255)
             gCCGAAGACCGGCAAGCAGTACGTGGTGATCACTGGTGGCGGCGCCCGCCTGTCGACTGACCGTGGCGACTATGTGATTTCTTACGCTCTGCCGTAAAACCGCGTCGCCTGCTTCGCGGGCAAGCCCGCTCCTACAGGGGCCTTGTAGGa                                                                                                                                         <  2:2414253/150‑1 (MQ=39)
             gCCGAAGACCGGCAAGCAGTACGTGGTGATCACTGCTGGCGGCGCCCGCCAGTCGACTGACCGTGGCGACTATGTGATTTCTTACGCTCTGCCGTAAAACCGCGTCGCCTGCTTCGCGGGGAAGCCCGCGCCGAAAAGGCCCCAGTAggg                                                                                                                                         >  1:2096962/1‑149 (MQ=255)
             gCCGAAGACCGGCAAGCAGTACGTGGTGATCACTGCTGGCGGCGCCCGCCAGTCGACTGACCGTGGCGACTATGTGATTTCTTACGCTCTGCCGTAAAACCGCGTCGCCTGCTTCGCGGGCAAGCCCGCTCCTACAAGGCCCCTGTAGGa                                                                                                                                         <  2:461990/150‑1 (MQ=255)
               cGAAGACCGGCAAGCAGTACGTGGTGATCACTGCTGGCGGCGCCCGCCAGTCGACTGACCGTGGCGACTATGTGATTTCTTACGCTCTGCCGTAAAACCGCGTCGCCTGCTTCGCGGGCAAGCCCGCTCCTACAGGGGCCTTGTAGGAGc                                                                                                                                       <  1:3976620/150‑1 (MQ=255)
               cGAAGACCGGCAAGCAGTACGTGGTGATCACTGCTGGCGGCGCCCGCCAGTCGACTGACCGTGGCGACTATGTGATTTCTTACGCTCTGCCGTAAAACCGCGTCGCCTGCTTCGCGGGCAAGCCCGCTCCTACAGGGGCCTTGTAGGAGc                                                                                                                                       <  1:5562232/150‑1 (MQ=255)
               cGAAGACCGGCAAGCAGTACGTGGTGATCACTGCTGGCGGCGCCCGCCAGTCGACTGACCGTGGCGACTATGTGATTTCTTACGCTCTGCCGTAAAACCGCGTCGCCTGCTTCGCGGGCAAGCCCGCTCCTACAGGGGCCTTGTAGGAGc                                                                                                                                       <  1:4520425/150‑1 (MQ=255)
               cGAAGACCGGCAAGCAGTACGTGGTGATCACTGCTGGCGGCGCCCGCCAGTCGACTGACCGTGGCGACTATGTGATTTCTTACGCTCTGCCGTAAAACCGCGTCGCCTGCTTCGCGGGCAAGCCCGCTCCTACAGGGGCCTTGTAGGAGc                                                                                                                                       <  2:3374003/150‑1 (MQ=255)
               cGAAGACCGGCAAGCAGTACGTGGTGATCACTGCTGGCGGCGCCCGCCAGTCGACTGACCGTGGCGACTATGTGATTTCTTACGCTCTGCCGTAAAACCGCGTCGCCTGCTTCGCGGGCAAGCCCGCTCCTACAAGGCCCCTGTAGGAGc                                                                                                                                       <  2:5078319/150‑1 (MQ=255)
                gAAGACCGGCAAGCAGTACGTGGTGATCACTGCTGGCGGCGCCCGCCAGTCGACTGACCGTGGCGACTATGTGATTTCTTACGCTCTGCCGTAAAACCGCGTCGCCTGCTTCGCGGGCAAGCCCGCTCCTACAAGGCCCCTGTAGGAGCg                                                                                                                                      <  1:1831708/150‑1 (MQ=255)
                gAAGACCGGCAAGCAGTACGTGGTGATCACCGCTGGCGGCGCCCGCCAGTCGACTGACCGTGGCGACTATGTGATTTCTTACGCTCTGCCGTAAAACCGCGTCGCCTGCTTCGCGGGCAAGCCCGCTCCTACAAGGCCCCTGTAGGAGCg                                                                                                                                      <  1:1407054/150‑1 (MQ=255)
                 caGACCGGCAAGCAGTACGTGGTGATCACTGCTGGCGGCGCCCGCCAGTCGACTGCCCGTGGCGACTATGTGATTTCTTCCGCTCTGCCGTAAAACCGCGTCGCCTGCTTCGCGGGCAAGCCCGCTCCTACAAGGCCCCGGTAGGAGCgg                                                                                                                                     <  2:2285150/149‑1 (MQ=255)
                 aaGACCGGCAAGCAGTACGTGGTGATCACTGCTGGCGGCGCCCGCCAGTCGACTGACCGTGGCGACTATGTGATTTCTTACGCTCTGCCGTAAAACCGCGTCGCCTGCTTCGCGGGCAAGCCCGCTCCTACAGGGGCCTTGTAGGAGCgg                                                                                                                                     <  2:4972738/150‑1 (MQ=255)
                 aaGACCGGCAAGCAGTACGTGGTGATCACTGCTGGCGGCGCCCGCCAGTCGACTGACCGTGGCGACTATGTGATTTCTTACGCTCTGCCGTAAAACCGCGTCGCCTGCTTCGCGGGCAAGCCCGCTCCTACAGGGGCCTTGTAGGAGCgg                                                                                                                                     <  2:5365726/150‑1 (MQ=255)
                 aaGACCGGCAAGCAGTACGTGGTGATCACTGCTGGCGGCGCCCGCCAGTCGACTGACCGTGGCGACTATGTGATTTCTTACGCTCTGCCGTAAAACCGCGTCGCCTGCTTCGCGGGCAAGCCCGCTCCTACAAGGCCCCTGTAGGAGCgg                                                                                                                                     <  2:635017/150‑1 (MQ=255)
                 aaGACCGGCAAGCAGTACGTGGTGATCACTGCTGGCGGCGCCCGCCAGTCGACTGACCGTGGCGACTATGTGATTTCTTACGCTCTGCCGTAAAACCGCGTCGCCTGCTTCGCGGGCAAGCCCGCTCCTACAAGGCCCCTGTAGGAGCgg                                                                                                                                     <  2:2161278/150‑1 (MQ=255)
                 aaGACCGGCAAGCAGTACGTGGTGATCACTGCTGGCGGCGCCCGCCAGTCGACTGACCGTGGCGACTATGTGATTTCTTACGCTCTGCCGTAAAACCGCGTCGCCTGCTTCGCGGGCAAGCCCGCTCCTACAAGGCCCCTGTAGGAGCgg                                                                                                                                     <  2:3447662/150‑1 (MQ=255)
                   gACCGGCAAGCAGTACGTGGTGATCACTGCTGGCGGCGCCCGCCAGTCGACTGACCGTGGCGACTATGTGATTTCTTACGCTCTGCCGTAAAACCGCGTCGCCTGCTTCGCGGGCAAGCCCGCTCCTACAAGGCCCCTGTAGGAGCGGGc                                                                                                                                   <  1:1666077/150‑1 (MQ=255)
                   gACCGGCAAGCAGTACGTGGTGATCACTGCTGGCGGCGCCCGCCAGTCGACTGACCGTGGCGACTATGTGATTTCTTACGCTCTGCCGTAAAACCGCGTCGCCTGCTTCGCGGGCAAGCCCGCTCCTACAAGGCCCCTGTAGGAGCGGGc                                                                                                                                   <  2:2039175/150‑1 (MQ=255)
                     ccGGCAAGCAGTACGTGGTGATCACTGCTGGCGGCGCCCGCCAGTCGACTGACCGTGGCGACTATGTGATTTCTTACGCTCTGCCGTAAAACCGCGTCGCCTGCTTCGCGGGCAAGCCCGCTCCTACAAGGCCCCTGTAGGAGCGGGCtt                                                                                                                                 <  1:127603/150‑1 (MQ=255)
                     ccGGCAAGCAGTACGTGGTGATCACTGCTGGCGGCGCCCGCCAGTCGACTGACCGTGGCGACTATGTGATTTCTTACGCTCTGCCGTAAAACCGCGTCGCCTGCTTCGCGGGCAAGCCCCCGCCGACCAGGCCCCAGTAGGAGCGGGcgt                                                                                                                                 >  1:1265932/1‑148 (MQ=255)
                       ggCAAGCAGTACGTGGTGATCACTGCTGGCGGCGCCCGCCAGTCGACTGACCGTGGCGACTATGTGATTTCTCACGCTCTGCCGTAAAACCGCGTCGCCTGCTTCGCGGGCAAGCCCGCTCCTACAAGGCCCCTGTAGGAGCGGGCTTGc                                                                                                                               <  1:2556107/150‑1 (MQ=255)
                              aGTACGTGTTGATCACTGCTGGCGGCGCCCGCCAGTCGACTGACCGTGGCGACTATGTGATTTCTTACGCTCTGCCTTAAACCCGCGTCGCCTGCTGCGCGGGCAAGCCCGCTCCAACAAGGCCCCGGTAGGAGCGGGCTGGCCcgctcc                                                                                                                        <  2:1401794/150‑4 (MQ=255)
                               gTACGTGGTGATCACTGCTGGCGGCGCCCGCCAGTCGACTGACCGTGGCGACTATGTGATTTCTTACGCTCTGCCGTAAAACCGCGTCGCCTGCTTCGCGGGCAAGCCCGCTCCTACAAGGCCCCTGTAGGAGCGGGCTTGCCCGCGAAt                                                                                                                       <  1:3370549/150‑1 (MQ=255)
                                tACGTGGTGATCACTGCTGGCGGCGCCCGCCAGTCGACTGACCGTGGCGACTATGTGATTTCTTACGCTCTGCCGTAAAACCGCGTCGCCTGCTTCGCGGGCAAGCCCGCTCCTACAAGGCCCCTGTAGGAGCGGGCTTGCCCGCGAATg                                                                                                                      >  1:2218570/1‑150 (MQ=255)
                                   gtggtgATCACTGCTGGCGGCGCCCGCCAGTCGACTGCCCGTGGCGACTATGTGATTTCTTACGCTCTGCCGTAAAACCGCGTCGCCTGCTTCGCGGGCAAGCCCGCTCCTACAAGGCCCCTGTAGGAGCGGGCTTGCCCGCGAATGccc                                                                                                                   <  2:3320823/150‑1 (MQ=255)
                                         aTCACTGCTGGCGGCGCCCGCCAGTCGACTGACCGTGGCGACTATGTGATTTCTTACGCTCTGCCGTAAAACCGCGTCGCCTGCTTCGCGGGCAAGCCCGCTCCTACAAGGCCCCTGTAGGAGCGGGCTTGCCCGCGAATGCCCCTCaaa                                                                                                             <  1:851470/150‑1 (MQ=255)
                                              gactgGCGGGCGCCGCCAGTCGACTGACCGTGGCGACTATGTGATTTCTTACGCTCTGCCGTAAAACCGCGTCGCCTGCTTCGCGGGCAAGCCCGCTCCTACAAGGCCCCTGTAGGAGCGGGCTTGCCCGCGAATGCCCCTCAAAGACaa                                                                                                        >  2:4110074/3‑150 (MQ=255)
                                               gctgGCGGCGCCCGCCAGTCGACGGACCGTGGCGACTATGTGATTTCTTACGCTCTGCCGTAAAACCGCGTCGCCTGCTTCGCGGGCAAGCCCGCTCCTACAAGGCCCCTGTAGGAGCGGGCTTGCCCGCGAATGCCCCTCAAAGACAAc                                                                                                       <  2:2771471/150‑1 (MQ=255)
                                                   gcggcgCCCGCCAGTCGACTGACCGTGGCGACTATGTGATTTCTTACGCTCTGCCGTAAAACCGCGTCGCCTGCTTCGCGGGCAAGCCCGCTCCTACAAGGCCCCTGTAGGAGCGGGCTTGCCCGCGAATGCCCCTCAAAGACAACCCga                                                                                                   <  1:5154717/150‑1 (MQ=255)
                                                     ggcgCCCGCCTGTCTACTGACGGTGGCTACTTTGTTTTTTCTAACGCTCTGCCTAAAAACCGCGTGGCCTGCTTCGCGGGCAAGCCCGCTCATACAAGGCCCCTGTAGGAGGGGGCTTGCCCGCGAATGCCCCTCAAAGACAACCCgaga                                                                                                 <  1:2882/150‑1 (MQ=255)
                                                     ggcgCCCGCCAGTCGAGTGACGGTGGCGACTTTGTTTTTTCTTAGGCTCTGCCGTAAAACCGCGTCGCCTGGTTCGGGGGAAAGCCCGCTCCTACAAGGCCCCTGTAGGTGGGGGCTTGCCCGCGAATGCCCCTCAAAGACAACCCgaga                                                                                                 <  1:1427859/150‑1 (MQ=255)
                                                        gcccgccAGTCGACTGACCGTGGCGACTATGTGATTTCTTACGCTCTGCCGTAAAACCGCGTCGCCTGCTTCGCGGGCAAGCCCGCTCCTACAAGGCCCCTGTAGGAGCGGGCTTGCCCGCGAATGCCCCTCAAAGACAACCCGAGATTc                                                                                              <  1:2769596/150‑1 (MQ=255)
                                                         cccgccAGACCACTCACCGTGGCCACTATGTGATTCCAAACGATCTCCCGTAAAACCGCGTCCCCTGCTTCGCGGCAAACCCCCCACATACAAGCCCCCTGTAGCACCGGGCTTGCCCGCGAATGCCCCTCAAAGACAACCCGAGATTCt                                                                                             <  1:989458/150‑1 (MQ=255)
                                                           cgccAGTCGACTGACCGTGGCGACTATGTGATTTCTTACGCTCTGCCGTAAAACCGCGTCGCCTGCTTCGCGGGCAAGCCCGCTCCTACAAGGCCCCTGTAGGAGCGGGCTTGCCCGCGAATGCCCCTCAAAGACAACCCGAGAGTCTGt                                                                                           <  2:1273795/150‑1 (MQ=255)
                                                              cAGTCGACTGACCGTGGCGACTATGTGATTTCTTACGCTCTGCCGTAAAACCGCGTCGCCTGCTTCGCGGGCAAGCCCGCTCCTACAAGGCCCCTGTAGGAGCGGGCTTGCCCGCGAATGCCCCTCAAAGACAACCCGAGATTCTGTAAt                                                                                        <  2:377613/150‑1 (MQ=255)
                                                                 tCGACTGACCGTGGCGACTATGTGATTTCTTACGCTCTGCCGTAAAACCGCGTCGCCTGCTTCGCGGGCAAGCCCGCTCCTACAAGGCCCCTGTCGGAGCGGGCTTGCCCGCGAAATCCCCTCAAAGACCACCCCAGATTCTGTCATGcc                                                                                     >  2:3030045/1‑150 (MQ=38)
                                                                 tCGACTGACCGTGGCGACTATGTGATTTCTTACGCTCTGCCGTAAAACCGCGTCGCCTGCTTCGCGGGCAACCCCGCTCCTACAAGGCCCCTGTAGGAGCGGGCTTGCCCGCGAATGCCCCTCAAAGACAACCCGAGATTCTGTAATTcc                                                                                     >  2:3544135/1‑150 (MQ=255)
                                                                       gACCGTGGCGACTATGTGATTTCTTACGCTCTGCCGTAAAACCGCGTCGCCTGCTTCGCGGGCAAGCCCGCTCCTACAAGGCCCCTGTAGGAGCGGGCTTGCCCGCGAATGCCCCTCAAAGACAACCCGAGATTCTGTAATGCCATCCGc                                                                               >  2:4085305/1‑150 (MQ=255)
                                                                        aCCGTGGCGACTATGTGATTTCTTACGCTCTGCCGTAAAACCGCGTCGCCTGCTTCGCGGGGCAGCCCCCGCCAAAAAAGCCCCAGGAGGAGCGGGCGTTCCCGCGCAATCCCCCAAAAAAAAACCCCAGATTAAGTAATTCCaaacccc                                                                              >  1:5334361/1‑144 (MQ=255)
                                                                           gTGGCGACTATGTGATTTCTTACGCTCTGCCGTAAAACCGCGTCGCCTGCTTCGCGGGCAAGCCCGCTCCTACAAGGCCCCTGTAGGAGCGGGCTTGCCCGCGAATGCCCCTCAAAGACAACCCGAGATTCTGTAATGCCATCCGCAAtt                                                                           >  2:1321574/1‑150 (MQ=255)
                                                                           gTGGCGACTATGTGATTTCTTACGCTCTGCCGTAAAACCGCGTCGCCTGCTTCGCGGGCAAGCCCGCGCCTAAAAGGCCCCTGGAGGAGCGGGGGTGCCCGCCCATGCCCCGCCAAGACAACCCGCGGTTCTGTTATGCCCTCCGccctt                                                                           >  1:3536284/1‑146 (MQ=255)
                                                                              gCGACTATGTGATTTCTTACGCTCTGCCGTAAAACCGCGTCGCCTGCTTCGCGGGCAAGCCCGCTCCTACAAGGCCCCTGTAGGAGCGGGCTTGCCCGCGAATGCCCCTCAAAGACAACCCGAGATTCTGTAATGCCATCCGCAATTCGc                                                                        <  2:142267/150‑1 (MQ=255)
                                                                                gACTTTGTTATTTCTTACGCTCTGCCGTAAAACCGCGTCGCCTGCTTCGCGGGCAAGCCCGCTCCTACACGGCCCCTGTCGGAGCGGGCTTGCCCGCGAATGCCCCTCAAAGACAACCCGAGATTCTGTAATGCCATCCGCAATTCGCAt                                                                      <  1:4707868/150‑1 (MQ=255)
                                                                                   tATGTGATTTCTTACGCTCTGCCGTAAAACCGCGTCGCCTGCTTCGCGGGCAAGCCCGCTACTACAATGCCCCTGTAGGAGCGGGCTTGCCCGCGAATGCCCCTCAAAGACAACCCGAGATTATGTAATGCCATCCGCAATTAGCATCAc                                                                   >  1:3756250/1‑150 (MQ=255)
                                                                                        gATGTCTTACGCTCTGCCGTAAAACCGCGTCGCCTGCTTCGCGGGCCAGCCCGCTCCTACCAGGCCCCTGTAGGAGCGGGCTTGCCCGCGAATGCCCCGCAAAGACAACCCGAGATTCTGTAATGCCCTCCGCAATTCGCATCACCCCCa                                                              >  2:3522534/1‑150 (MQ=39)
                                                                                         aTTTCTTACGCTCTGCCGTAAAACCGCGTCGCCTGCTTCGCGGGCAAGCCCGCTCCTACAGGGGCCTTGTAGGAGc                                                                                                                                       >  1:765229/1‑76 (MQ=21)
                                                                                            tCTTACGCTCTGCCGTAAAACCGCGTCGCCTGCTTCGCGGGCAAGCCCGCTCCTACACGGCCCCTGTAGGAGCGGGCTTGCCCGCGAATGCCCCTCAAAGACAACCCGAGATTCTGTAATGCCATCCGCAATTCGCATCACCCCCACCct                                                          <  1:339542/150‑1 (MQ=255)
                                                                                                   ctctGCCGTAAACCCGCGTCGCCTGCTTCGCGGGCAAGCCCGCTCCTACAAGGCCCCTGTAGGAGCGGGCTTGCCCGCGAATGCCCCTCAAAGACAACCCGAGATTCTGTAATGCCATCCGCAATTCGCATCACCCCCACCctgctgctg                                                   <  2:1015720/150‑1 (MQ=255)
                                                                                                   ctctGCCGTAAAACCGCGTCGCCTGCTTCGCGGGCAAGCCCGCTCCTACACGGCCCCTGTAGGAGCGGGCTTGCCCGCGAATGCCCCTCAAAGACAACCCGAGATTCTGTAATGCCATCCGCAATTCGCATCACCCCCACCcggctgctg                                                   <  1:1835374/150‑1 (MQ=255)
                                                                                                      ggcggTAAACCCGGGTGGCCTGCTTGGCGGGCAAGCCCGCTCCTCCATGCCCCGTGTAGGAGCGGGGTTGCCCGCGAATGCCCCTCAAAGACAACCCGAGATTCTGTAATGCCATCCGCAATTCGCATCACCCCCACCCTGCTGctggcc                                                <  1:90397/146‑1 (MQ=255)
                                                                                                        ccGTAAAACCGCGTCGCCTGCTTCGCGGGCAAGCCCGCTCCTACAATGCCCCTGTAGGAGCGGGCTTGCCCGCGAATGCCCCTCAAAGACAACCCGAGATTCTGTAATGCCATCCGCAATTCGCATCACCCCCACCCTGCTGctggccct                                              >  2:2223286/1‑150 (MQ=255)
                                                                                                         cggaacaCCGCGGCGCCTGCTTCGCGGGCAAGCCCGCGCCTCCCAGGCCCCTGTAGGAGCGGGCGTGCCCGCGAATGCCCCTCAACGACAACCCGAGATTCTGTAATGCCATCCGCAATTCGCATCACCCCCCCCCTGCTGctggccctg                                             <  2:1610577/144‑1 (MQ=255)
                                                                                                               aCCGCGTCGCCTGCTTCGCGGGCAAGCCCGCTCCTACACGGCCCCTGTAGGAGCGGGCTTGCCCGCGAATGCCCCTCAAAGACAACCCGAGATTCTGTAATGCCATCCGCAATTCGCATCACCCCCACCCTGCTGCTGGCCCTGGCCAGc                                       <  2:1417214/150‑1 (MQ=255)
                                                                                                                ccGCGTCGCCTGCTTCGCGGGCAAGCCCGCTCCTACAAGGCCCCTGTAGGAGCGGGCTTGCCCGCGAATGCCCCTCAAAGACAACCCGAGATTCTGTAATGCCATCCGCAATTCGCATCACCCCCACCCTGCTGCTGGCCAGGGCCAGca                                      <  2:5469323/150‑1 (MQ=255)
                                                                                                                ccGCGTCGCCTGCTTCGCGGGCAAGCCCGCTCCTACAAGGCCCCTGTAGGAGCGGGCTTGCCCGCGAATGCCCCTCAAAGACAACCCGAGATTCTGTAATGCCATCCGCAATTCGCATCACCCCCACCCTGCTGCTGGCCAGGGCCAGca                                      <  1:4093276/150‑1 (MQ=255)
                                                                                                                ccGCGTCGCCTGCTTCGCGGGCAAGCCCCCTCCAACAACGCCCCAGTAGGAGCGGGCATGCCCGCGAATTCCCCAAAAAGACAACCCGAGATTCTGTAAAGCCATCCGCCATTAGCTTCACCCCCACCCGGCTGCTGGCCCAGGCCCGcc                                      >  2:2529937/1‑149 (MQ=255)
                                                                                                                 cgcgTCGCCTGCTTCGCGGGCAAGCCCGCTCCTACAAGGCCCCTGTAGGAGCGGGCTTGCCCGCGAATGCCCCTCAAAGACAACCCGAGATTCTGTAATGCCATCCGCAATTCGCATCACCCCCACCCTGCTGCTGGCCCTGGCCAGcac                                     <  2:3705051/150‑1 (MQ=255)
                                                                                                                  gcgTCGCCTGCTTCGCGGGCAAGCCCGCTCCTACAAGGCCCCTGTAGGAGCGGGCTTGCCCGCGAATGCCCCTCAAAGACAACCCGAGATTCTGTAATGCCATCCGCAATTCGCATCACCCCCACCCTGCTGCTGGCCAGGGCCAGcagc                                    <  1:4116774/150‑3 (MQ=255)
                                                                                                                                gcgGGCAAGCCCGCTCCTACAAGGCCCCTGTAGGAGCGGGCTTGCCCGCGAATGCCCCTCAAAGACAACCCGAGATTCTGTAATGCCATCCGCAATTCGCATCACCCCCACCCTGCTGCTGGCCCTGGCCAGCACCACTGCCCTTGCCAg                      <  2:3217600/150‑1 (MQ=255)
                                                                                                                                gcgGGCAAGCCCGCTCCTACAAGGCCCCTGTAGGAGCGGGCTTGCCCGCGAATGCCCCTCAAAGACAACCCGAGATTCTGTAATGCCATCCGCAATTCGCATCACCCCCACCCTGCTGCTGGCCCTGGCCAGCACCACTGCCCTTGCCAg                      >  2:2216276/1‑150 (MQ=255)
                                                                                                                                  gggCAAGCCCGCTCCTACAAGGCCCCTGTTGGAGCGGGCTTGCCCGCGAATGCCCCTCAAAGACAACCCGAGATTCTGTAATGCCATCCGCAATTCGCATCACCCCCCCCCTGCTGCTGGCCCTGGCCAGCACCACTGCCCTTGCCAGcg                    <  1:1217536/150‑1 (MQ=255)
                                                                                                                                    gCAAGCCCGCTCCTACACGGCCCCTGTAGGAGCGGGCTTGCCCGCGAATGCCCCTCAAAGACAACCCGAGATTCTGTAATGCCATCAGCAATTCGCATCACCCCCCCCCTGCTGCTGGCCCTGGACAGCTCCACTGCCCTTGACAGcgac                  <  2:5061410/150‑1 (MQ=39)
                                                                                                                                      aaGCCCGCTCCTACAGGGGCCTTGTAGGAGCGGGCTTGCCCGCGAATGCCCCTCAAAGACAACCCGAGATTCTGTAATGCCATCCGCAATTCGCATCACCCCCACCCTGCTGCTGGCCCTGGCCAGCACCACTGCCCTTGCCAGcgacga                >  2:4363765/1‑150 (MQ=255)
                                                                                                                                      aaGCCCGCTCCTACAAGGCCCCTGTAGGAGCGGGCTTGCCCGCGAATGCCCCTCAAAGACAACCCGAGATTCTGTAATGCCATCCGCAATTCGCATCACCCCCACCCTGCTGCTGGCCCTGGCCAGCACCACTGCCCTTGCCAGcgacga                >  2:2578410/1‑150 (MQ=255)
                                                                                                                                       aGCCCGCTCCTACAAGGCCCCTGTAGGAGCGGGCTTGCCCGCGAATGCCCCTCCAAGACAACCCGAGATTCTGTAATGCCATCCGCAATTCGCATCACCCCCACCCTGCTGCTGGCCCTGGCCAGCACCACTGCCCTTGCCAGCGACGAt               >  2:1489246/1‑150 (MQ=255)
                                                                                                                                         cccGCTCCTACAGGGGCCTTGTAGGAGCGGGCTTGCCCGCGAATGCCCCTCAAAGACAACCCGAGATTCTGTAATGCCATCCGCAAGTCGCATCACCCCCACCCTGCTGCTGGCCCTGGCCCGCACCACTGCCCTTGCCAGCGACGAtct             >  2:2923371/1‑150 (MQ=39)
                                                                                                                                         cccGCTCCTACAAGGCCCCTGTAGGAGCGGGCTTGCCCGCGAATGCCCCTCAAAGACAACCCGAGATTCTGTAATGCCATCCGCAATTCGCATCACCCCCACCCTGCTGCTGGCCCTGGCCAGCACCACTGCCCTTGCCAGCGACGAtct             >  2:4533642/1‑150 (MQ=255)
                                                                                                                                         cccGCTCCTACAAGGCCCCTGTAGGAGCGGGCTTGCCCGCGAATGCCCCTCAAAGACAACCCGAGATTCTGTAATGCCATCCGCAATTCGCATCACCCCCACCCTGCTGCTGGCCCTGGCCAGCACCACTGCCCTTGCCAGCGACGAtct             >  2:1719282/1‑150 (MQ=255)
                                                                                                                                          ccGCTCCTACAGGGGCCTTGTAGGAGCGGGCTTGCCCGCGAATGCCCCTCAAAGACAACCCGAGATTCTGTAATGCCATCCGCAATTCGCATCACCCCCACCCTGCTGCTGGCCCTGGCCAGca                                      >  1:1497799/1‑124 (MQ=255)
                                                                                                                                           cGCACCTACAAGGCCCCTGTCGGAGCGGGCTTGCCCGCGAATGCCCCTCACAGACAACCCGAGATTCTGTAATGCCATCCGCAATTCGCATCACCCCCACCCTGCTGCTGGCCCTGGCCAGCACCACTGCCCTTGCCAGCGACGATCTCa           <  1:346852/150‑1 (MQ=255)
                                                                                                                                            gCTCCTACAAGGCCCCTGTAGGAGCGGGCTTGCCCGCGAATGCCCCTCAAAGACAACCCGAGATTCTGTAATGCCATCCGCAATTCGCATCACCCCCACCCTGCTGCTGGCCCTGGCCAGCACCCCTGCCCTTGCTAGCGCCGATCTCAt          >  1:2987858/1‑150 (MQ=255)
                                                                                                                                            gCTCCTACAAGGCCCCTGTAGGAGCGGGCTTGCCCGCGAATGCCCCTCAAAGACAACCCGAGATTCTGTAATGCCATCCGCAATTCGCATCACCCCCACCCTGCTGCTGGCCCTGGCCAGCACCACTGCCCtt                           <  1:4578391/133‑1 (MQ=255)
                                                                                                                                            gCTCCTACAAGGCCCCTGTAGGAGCGGGCTTGCCCGCGAATGCCCCTCAAAGACAACCCGAGATTCTGTAATGCCATCCGCAATTCGCATCACCCCCACCCTGCTGCTGGCCCTGGCCAGCACCACTGCCCTTGCTAGCGACGATCTCAt          >  1:2987996/1‑150 (MQ=255)
                                                                                                                                              tCCTACAAGGCCCCTGTAGGAGCGGGCTTGCCCGCGAATGCCCCTCAAAGACAACCCGAGATTCTGTAATGCCATCCGCAATTCGCATCACCCCCACCCTGCTGCTGGCCCTGGCCAGCACCACTGCCCTTGCCAGCGACGATCTCATGa        >  2:56006/1‑150 (MQ=255)
                                                                                                                                              tCCTACAAGGCCCCTGTAGGAGCGGGCTTGCCCGCGAATGCCCCTCAAAGACAACCCGAGATTCTGTAATGCCATCCGCAATTCGCATCACCCCCACCCTGCTGCTGGCCCTGGCCAGCACCACTGCCCTTGCCAGCGACGATCTCATGa        >  1:2658426/1‑150 (MQ=255)
                                                                                                                                               ccTCCAAGGCCCCTGTAGGAGCGGGCTTGCCCGCGAATGCCCCTCAAACACAACCCGAGATTCTGTAATGCCATCCGCAATTCGCATCACCCCCACCCTGCTGCTGGCCCTGGCCAGCACCACTGCCCTTGCCAGCGACGATCTCATGAc       <  1:4033644/150‑1 (MQ=255)
                                                                                                                                               ccTACACGGCCCCTGTAGGAGCGGGCTTGCCCGCGAATGCCCCTCAAAGACAACCCGAGATTCTGTAATGCCATCCGCAATTCGCATCACCCCCACCCTGCTGCTGGCCCTGGCCAGCACCACTGCCCTTGCCAGCGACGATCTCATGAc       <  1:2280931/150‑1 (MQ=255)
                                                                                                                                               ccTACAAGGCCCCTGTAGGAGCGGGCTTGCCCGCGAATGCCCCTCAAAGACAACCCGAGATTCTGTAATGCCATCCGCAATTCGCATCACCCCCACCCTGCTGCTGGCCCTGGCCAGCACCACTGCCCTTGCCAGCGACGATCTCATGAc       <  2:4348288/150‑1 (MQ=255)
                                                                                                                                               ccTACAAGGCCCCTGTAGGAGCGGGCTTGCCCGCAAATGCCCCTCAACGACAACCCGAGATTCTGTAATGCCATCCGCAATTCGCATCACCCCCACCCTGCTGCTGGCCCTGGCCAGCACCACTGCCCTTGCCAGCGACGATCTCATGAc       <  2:732835/150‑1 (MQ=255)
                                                                                                                                                cTACAAGGCCCCTGTAGGAGCGGGCTTGCCCGCGAATGCCCCTCAAAGACAACCCGAGATTCTGTAATGCCATCCGCAATTCGCATCACCCCCACCCTGCTGCTGGCCCTGGCCa                                         <  2:5128268/115‑1 (MQ=255)
                                                                                                                                                 tACAAGGCCCCTGTAGGAGCGGGCTTGCCCGCGAATGCCCCTCAAAGACAACCCGAGATTCTGTAATGCCATCCGCAATTCGCATCACCCCCACCCTGCTGCTGGCCCTGGCCAGCACCACTGCCCTTGCCAGCGACGATCTCATGAccc     >  1:4391238/1‑150 (MQ=255)
                                                                                                                                                  aCAAGGCCCCTGTAGGAGCGGGCTTGCCCGCGAATGCCCCTCAAAGACAACCCGAGATTCTGTAATGCCATCCGCAATTCGCATCACCCCCACCCTGCTGCTGGCCCTGGCCAGCACCACTGCCCTTGCCAGCGACGATCTCATGACCCg    >  2:5452575/1‑150 (MQ=255)
                                                                                                                                                    aaGGCCCCTGTAGGAGCGGGCTTGCCCGCGAATGCCCCTCAAAGACAACCCGAGATTCTGTAATGCCATCCGCAATTCGCATCACCCCCACCCTGCTGCTGGCCCTGGCCAGCACCACTGCCCTTGCCAGCGACGATCTCATGACCCgca  <  1:5065493/150‑1 (MQ=255)
                                                                                                                                                    aaGGCCCCTGTAGGAGCGGGCTTGCCCGCGAATGCCCCTCAAAGACAACCCGAGATTCTGTAATGCCATCCGCAATTCGCATCACCCCCACCCTGCTGCTGGCCCTGGCCAGCACCACTGCCCTTGCCAGCGACGATCTCATGACCCgca  >  1:3216519/1‑150 (MQ=255)
                                                                                                                                                    aaGGCCCCTGTAGGAGCGGGCTTGCCCGCGAATGCCCCTCAAAGACAACCCGAGATTCTGTAATGCCATCCGCAATTCGCATCACCCCCACCCTGCTGCTGGCCCTGGCCAGCACCACTGCCCTTGCCAGCGACGATCTCATGACCCgca  >  1:1124571/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
TGACCTACGTCTCGCCGAAGACCGGCAAGCAGTACGTGGTGATCACTGCTGGCGGCGCCCGCCAGTCGACTGACCGTGGCGACTATGTGATTTCTTACGCTCTGCCGTAAAACCGCGTCGCCTGCTTCGCGGGCAAGCCCGCTCCTACAAGGCCCCTGTAGGAGCGGGCTTGCCCGCGAATGCCCCTCAAAGACAACCCGAGATTCTGTAATGCCATCCGCAATTCGCATCACCCCCACCCTGCTGCTGGCCCTGGCCAGCACCACTGCCCTTGCCAGCGACGATCTCATGACCCGCA  >  NC_002947/4054751‑4055048

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: