Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0129671,341,9221.G21.1% 101.9 / 30.1 38intergenic (‑223/+46)osmB/yciTosmotically‑inducible lipoprotein OsmB/DNA‑binding transcriptional regulator YciT
Reads supporting (aligned to +/- strand):  ref base . (18/12);  new base G (8/0);  total (26/12)
Fisher's exact test for biased strand distribution p-value = 3.85e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.99e-01
Rejected as consensus: Frequency below/above cutoff threshold.

TAATTTTATATCTTGAGAGTGTTAATAACAGGTAAATAGTCTTAATTATCAACCAGGAATCATCTTAGAGCGGATGATTTGCCAAACTGCAAATCATCCGTAGAGAAGGGAAATGGTTAAATGTCAACGACGTGAACGGTGA  >  NC_012967/1341852‑1341993
                                                                       |                                                                          
taatTTTATATCTTGAGAGTGTTAATAACAGGTAAATAGTCTTAATTATCAACCAGGAATCATCTTAGAGCg                                                                           >  1:9620286/1‑72 (MQ=255)
taatTTTATATCTTGAGAGTGTTAATAACAGGTAAATAGTCTTAATTATCAACCAGGAATCATATTAGAGCg                                                                           <  1:1880426/72‑1 (MQ=255)
      tataTCTTGAGAGTGTTAATAACAGGTAAATAGTCTTAATTATCAACCAGGAATCATCTTAGAGCGgatgat                                                                    >  2:7153373/1‑72 (MQ=255)
           cTTGAGAGTGTTAATAACAGGTAAATAGTCTTAATTATCAACCAGGAATCATCTTAGAGCGGTTGATTTGcc                                                               <  1:238874/72‑1 (MQ=255)
           cTTGAGAGTGTTAATAACAGGTAAATAGTCTTAATTATCAACCAGGAATCATCTTAGAGCGGGAGAttttgc                                                               >  2:2826884/1‑69 (MQ=255)
            ttGAGAGTGTTAATAACAGGTAAATAGTCTTAATTATCAACCAGGAATCATCTTAGAGCGGATGATTTGCCa                                                              <  1:2787723/72‑1 (MQ=255)
                gagTGTTAATAACAGGTAAATAGTCTTAATTATCAACCAGGAATCATCTTAGAGCGGATGATTTGCCAAACt                                                          <  1:9812850/72‑1 (MQ=255)
                gagTGTTAATAACAGGTAAATAGTCTTAATTATCAACCAGGAATCATCTTAGAGCGGATGATTTGCCAAACt                                                          >  2:7779868/1‑72 (MQ=255)
                      taataaCAGGTAAATAGTCTTAATTATCAACCAGGAATCATCTTAGAGCGGATGATTTGCCAAACTGCAAat                                                    <  1:8989328/72‑1 (MQ=255)
                      taataaCAGGTAAATAGTCTTAATTATCAACCAGGAATCATCTTAGAGCGGATGATTTGCCAAACTGCAAat                                                    <  2:3005734/72‑1 (MQ=255)
                        ataaCAGGTAAATAGTCTTAATTATCAACCAGGAATCATCTTAGAGCGGGAGGATTTTCCAAATGGCaaaca                                                   >  1:5249824/1‑70 (MQ=255)
                        ataaCAGGTAAATAGTCTTAATTATCAACCAGGAATCATCTTAGAGCGGATGATTTGCCAAACTGCAAatca                                                  <  1:7109209/72‑1 (MQ=255)
                                  aaTAGTCTTAATTATCAACCAGGAATCATCTTAGAGCGGAAGGTTTGCCAAACTGCAAATCATCCGTAggga                                      >  1:6584582/1‑72 (MQ=255)
                                    tAGTCTTAATTATCAACCAGGAATCATCTTAGAGCGGGAGGATTTGCCAAAATCTCAAACATCCGCAGGGaa                                     >  2:5116524/1‑72 (MQ=255)
                                     actCTTAATTATCAACCAGGAATCATCTTAGAGCGGGAGGATTTGCCAAATTGCAAATAATCGGTGGGGAAg                                    >  2:8084484/3‑72 (MQ=255)
                                     aGTCTTAATTATCAACCAGGAATCATCTTAGAGCGGGAGAATTTCCCAACTTCCAAACCTCCCGAGAGAAgg                                   >  2:8228305/1‑72 (MQ=255)
                                         ttAATTATCAACCAGGAATCATCTTAGAGCGGATGATTTGCCAAACTGCAAATCATCCGTAGAGAAGGGaaa                               <  2:7793754/72‑1 (MQ=255)
                                           atttATCACCCAGGAATCATCTTAGAGCGGATGATTTGCAAAACTGAAAATCATCCGTAGAGAAGGGAAATg                             <  2:8912504/70‑1 (MQ=255)
                                           aaTTATCAACCAGGAATCATCTTAGAGCGGGATGTTTTGCCAAACTGCAAATAATCCGTGAGAAAGGGaaag                              >  2:8686331/1‑71 (MQ=255)
                                            aTTATCAACCAGGAATCATCTTAGAGCGGGATGATTTGCCCAAAGCCAAAACACCGGTGAAGAAGGGgaagg                             >  2:4538573/1‑70 (MQ=255)
                                            aTTATCAACCAGGAATCATCTTAGAGCGGATGATTTTCCAAACTGCAAATCATCCGTAGAGAAGGGGAATgg                            >  2:9665762/1‑72 (MQ=255)
                                            aTTATCAACCAGGAATCATCTTAGAGCGGAAGATTTTCCAAAACTCAAAACATCCCTAGGGAAGGGGAAAgg                            >  1:7333242/1‑72 (MQ=255)
                                             ttATCAACCAGGAATCATCTTAGAGCGGGAGGATTTGCCCAAATCCCAAACACCCGCAggaaaggggaaagg                            >  2:3081780/1‑72 (MQ=255)
                                             ttATCAACCAGGAATCATCTTAGAGCGGGTGATTTTTCCAAACTGCAATAATCCCTAGAGGAGGGGAAAggg                           >  1:8876837/1‑71 (MQ=255)
                                              tATCAACCAGGAATCATCTTAGAGCGGGATAATTTTCCAAATGGAAAATATCCCTTGGGGAAGGGAAAGGGt                           >  1:5055148/1‑72 (MQ=255)
                                                tCAACCAGGAATCATCTTAGAGCGTATGATTTGCCAAACTGCAAATCATCCGTAGAGAAGGGAAATGGTTaa                        >  1:6902873/1‑72 (MQ=255)
                                                 cAACCAGGAATCATCTTAGAGCGGGGGGATTTTTCCAAATCCAAAACCTCCCTTGAGAAGGGAAAAGGGTaa                        >  1:3014557/1‑72 (MQ=255)
                                                 cAACCAGGAATCATCTTAGAGCGGAGGATTTTCCAAAATGCAAATCATCCGTAGAGAAGGGAAAGGGGTaaa                       >  1:7772050/1‑72 (MQ=255)
                                                  aaCCAGGAATCATCTTAGAGCGGGAGATTTTTCCAAAAGCCAAAAAATCCTCGAGAAAGGGAAAAGGTTaaa                       >  1:5336888/1‑72 (MQ=255)
                                                  aaCCAGGAATCATCTTAGAGCGGATGATTTGCCAAACTGCCAATCATCCGTAGAGAAGGGGAATGGTTAAAt                      >  2:9657997/1‑72 (MQ=255)
                                                    ccAGGAATCATCTTAGAGCGGGATGATTTTCCCAAATCAAAATATCCCCAGGAGAAGGGAAACGGTTaaagt                     >  2:4412266/1‑70 (MQ=255)
                                                     cAGGAATCATCTTAGAGCGGATGATTTGCCAAACTGCAAATCATCCGTAGAGAAGGGGAATGGTTAAATGTc                   >  1:7811685/1‑72 (MQ=255)
                                                      ggaaaTCATTTTAGAGGGGATGATTTGCAAAACGCAAAATCATCCGTAGAGAAGGGAAATGGTTAAATGTCa                  <  2:3927141/69‑1 (MQ=255)
                                                      aGGAATCATCTTAGAGCGGGTGATTTTCCCAACTGCCAATCCTCCGTAGGGAAGGGGAATGGTTTAATGTCa                  >  2:9548573/1‑72 (MQ=255)
                                                      aGGAATCATCTTAGAGCGGATGATTTGCCAAACTGCAAATCATCCGTAGAGAAGGGAAATGGTTAAATGTCa                  <  1:6953903/72‑1 (MQ=255)
                                                      aGGAATCATCTTAGAGCGGATGATTTGCCAAAATGCAAATCATCCGTAGAGAAGGGGAATGGTTAAATGTCa                  >  2:9729135/1‑72 (MQ=255)
                                                       ggAATCATCTTAGAGCGTATGATTTGCCAAACTGCAAATCATCCGTAGAGAAGGGAAATGGTTAAATGTCaa                 <  2:4565418/72‑1 (MQ=255)
                                                           tcatcTTAGAGCGGATGATTTGCAAAACTGCAAATCACCCGTAGAGAAGGGAAATGGTTAAATGTCAacgac             <  1:9313313/72‑1 (MQ=255)
                                                              tcTTATAGCGGATGATTTGCCAAACTGCAAATCATCCGTAGAGAAGGGAAATTGTTAAATGTCAACGACGTg          >  1:2269844/1‑72 (MQ=255)
                                                              tcTTAGAGCGGATGATTTGCCAAACTGCAAATCATCCGTAGAGAAGGGAAATGGTTAAATGTCAACGACGTg          <  2:7172268/72‑1 (MQ=255)
                                                                ttAGAGCGGATGATTTGCAAAACTGCAAATCATCCGTAGAGAAGGGAAATGGTTAAATGTCAACGACGTGaa        <  2:7146255/72‑1 (MQ=255)
                                                                 tAGGGCGGTTGATTGGCAAAACTGCAAATCATCNGTAGAGNAGGGAAATGGTTAAATGTCAACGACGTGAAc       <  1:581572/72‑1 (MQ=255)
                                                                   gagCGGATGATTTGCCAAACTGAAAATCATCCGTAGAGAAGGGAAATGGTTAAATGTCAACGACGTGAACgg     <  1:6643848/72‑1 (MQ=255)
                                                                    agCGGATGATTTCCAAAACTGAAATTAATCCGTAGAGAAGGGAAATGGTTAAATGTCAACGACGTGAACGGt    <  1:2654277/72‑1 (MQ=255)
                                                                     gCGGATGATTTGCCAAACTGCAAATCATCCGTAGAGAAGGGAAATGGTTAAATGTCAACGACGTGAACGGTg   >  1:5159351/1‑72 (MQ=255)
                                                                     gCGGAAGACTTGCCAAACTGCAAATCATCCGTAGAGAAGGGAAATGGTTAAATGTCAACGACGTGAACgggg   >  2:2916989/1‑70 (MQ=255)
                                                                      cGGATGATTTGCCAAACTGCAAATCATCCGTAGAGAAGGGAAATGGTTAAATGTCAACGACGTGAACGGGGa  >  1:4999819/1‑72 (MQ=255)
                                                                       |                                                                          
TAATTTTATATCTTGAGAGTGTTAATAACAGGTAAATAGTCTTAATTATCAACCAGGAATCATCTTAGAGCGGATGATTTGCCAAACTGCAAATCATCCGTAGAGAAGGGAAATGGTTAAATGTCAACGACGTGAACGGTGA  >  NC_012967/1341852‑1341993

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 23 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: