Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,440,2970CA24.0% 48.5 / 12.7 25intergenic (‑177/+88)yfcJ/fabBputative arabinose efflux transporter/3‑oxoacyl‑[acyl‑carrier‑protein] synthase I
Reads supporting (aligned to +/- strand):  ref base C (10/9);  new base A (3/3);  total (13/12)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.03e-01
Rejected as consensus: Frequency below/above cutoff threshold.

TCCCCGTCTCATCTCTCTGGTTTCAGGGTTACGGTGCGTTGGCAGGATTTAACGCGTACGTCTTTTCAGAAGGAAATCGACAAAGCGGGAAGTTTGCCTGGAACTGGCGGCGATTGTCAATGATGTGAAAAAGGGAACCATCAGGTTCCCTTTTGCGTTAGTGCCGGAGGCCGCATTGGCGCGTAACGTCGGATGCGACGCTGGCGCGTCTACTCCGACCTACTGCGAATTAATCTTTCAGCTTGCGCATTACCAGCGTGGCGTTGGTGCCGCCGAAGCCGAAGCTGT  >  NC_000913/2440156‑2440443
                                                                                                                                             |                                                                                                                                                  
tCCCCGTCTCATCTCTCTGGTTTCAGGGTTACGGTGCGTTGGCAGGATTTAACGCGTACGTCTTTTCAGAAGGAAATCGACAAAGCGGGAAGTTTGCCTGGAACTGGCGGCGATTGTCAATGATGTGAAAAAGGGAACCATCAGGTTcc                                                                                                                                             <  1:27459/149‑1 (MQ=255)
     gTCTCATCTCTCGGGTTTCAGGGGTACGGTGCGTTGTCATGATTTAACGCGTACGTCTTTTCAGAAGGAAATCGACAAAGCGGGAAGTTTGCCTGGATCTGGCGGCGTTTGTCAATGATGTGAAAAAGGGAACCATCAGGTTCCCtttt                                                                                                                                        <  4:131770/149‑1 (MQ=255)
        tcATCTCTCTGGTTTCAGGGTTACGGTGCGTTGGCAGGATTTAACGCGTACGTCTTTTCAGAAGGAAATCGACAAAGCGGGAAGTTTGCCTGGAACTGGCGGCGATTGTCAATTATGTGAAAAAGGGAACCATCAGGTTcc                                                                                                                                             <  8:53889/141‑1 (MQ=255)
        tcATCTCTCTGGTTTCAGGGTTACGGTGCGTTGGCAGGATTTAACGCGTACGTCTTTTCAGAAGGAAATCGACAAAGCGGGAAGTTTGCCTGGAACTGGCGGCGATTGTCAATGATGTGAAAAAGGGAACCATCAGGTTcc                                                                                                                                             >  7:53889/1‑141 (MQ=255)
                                   gCGTTGGCAGGATTTAACGCGTACGTCTTTTCAGAAGGAAATCGACAAAGCGGGAAGTTTGCCTGGAACTGGCGGCGATTGTCAATGATGTGAAAAAGGGAACCATCAGGTTCCCTTTTGCGTTAGTGCCGGAGGCCGCATTGGcgcg                                                                                                           >  8:92826/1‑148 (MQ=255)
                                     gTTGGCAGGATTTAACGCGTACGTCTTTTCAGAAGGAAATCGACAAAGCGGGAAGTTTGCCTGGAACTGGCGGCGATTGTCAATGATGTGAAAAAGGGAACCTGATGGTTCCCttttt                                                                                                                                       <  5:112430/118‑2 (MQ=255)
                                     gTTGGCAGGATTTAACGCGTACGTCTTTTCAGAAGGAAATCGACAAAGCGGGAAGTTTGCCTGGAACTGGCGGCGATTGTCAATGATGTGAAAAAGGGAACCTGATGGTTCCCttttt                                                                                                                                       >  6:112430/1‑117 (MQ=255)
                                                  aaCGCGTACGTCTTTTCAGCAGGAAATCGACAAACCGGGAAGTTGGCCTGGAACTGGCGGCGATTGTCAATGATGTGAAAAAGGGAACCATCAGGTTCCCTTTTGCGTTAGTGCCGGAGGCCGCATTGGCGCGTAACGTCGGATGCGa                                                                                            <  3:122482/148‑1 (MQ=255)
                                                  aaCGCGTACGTCTTTTCAGAAGGAAATCGACAAAGCGGGAAGTTTGCCTGGAACTGGCGGCGATTGTCAATGATGTGAAAAAGGGAACCATCAGGTTCCCTTTTGCGTTAGTGCCGGAGGCCGCATTGGCGCGTAACGTCGGATGCGAc                                                                                           >  3:37191/1‑149 (MQ=255)
                                                              ttttCAGAAGGAAATCGACAAAGCGGGAAGTTTGCCTGGAACTGGCGGCGATTGTCAATGATGTGAAAAAGGGAACCATCAGGTTCCCTTTTGCGTTAGTGCCGGAGGCCGCATTGGCGCGTAACGTCGGATGCGACGCTGGCGCGTc                                                                                >  2:134392/1‑148 (MQ=255)
                                                                ttCAGAAGGAAATCGACAAAGCGGGAAGTTTGCCTGGAACTGGCGGCGATTGTCAATGATGTGAAAAAGGGAACCATCAGGTTCCCTTTTGCGTTAGTGCCGGAGGCCGCATTGGCGCGTAACGTCGGATGCGACGCTGGCGCGTCTAc                                                                             >  3:67260/1‑149 (MQ=255)
                                                                   aGAAGGAAATCGACAAAGCGGGAAGTTTGCCTGGAACTGGCGGCGATTGTCAATGATGTGAAAAAGGGAACCATCAGGTTCCCTTTTGCGTTAGTGCCGGAGGCCGCATTGGCGCGTAACGTCGGATGCGACGCTGGCGCGTCTACTcc                                                                          >  5:61472/1‑149 (MQ=255)
                                                                              gACAAAGCGGGAAGTTTGCCTGGAACTGGCGGCGATTGTCAATGATGTGAAAAAGGGAACCATCAGGTTCCCTTTTGCGTTAGTGCCGGAGGCCGCATTGGCGCGTAACGTCGGATGCGACGCTGGCGCGTCTACTCCGACCTACTGc                                                                <  7:92826/148‑1 (MQ=255)
                                                                                          aGTTTGCCTGGAACTGGCGGCGATTGTCAATGATGTGAAAAAGGGAACCATCAGGTTCCCTTTTGCGTTAGTGCCGGAGGCCGCATTGGCGCGTAACGTCGGATGCGACGCTGGCGCGTCTACTCCGACCTACTGCGAATTAATCTTTc                                                   <  4:37191/149‑1 (MQ=255)
                                                                                                           cggcgATTGTCAATGATGTGAAAAAGGGAACCATCAGGTTCCCTTTTGCGTTAGTGCCGGAGGCCGCATTGGCGCGTAACGTCGGATGCGACGCTGGCGCGTCTACTCCGACCTACTGCGAATTAATCTTTCAGCTTGCGCATTAc                                     >  8:145080/1‑146 (MQ=255)
                                                                                                           cggcgATTGTCAATGATGTGAAAAAGGGAACCATCAGGTTCCCTTTTGCGTTAGTGCCGGAGGCCGCATTGGCGCGTAACGTCGGATGCGACGCTGGCGCGTCTACTCCGACCTACTGCGAATTAATCTTTCAGCTTGCGCATTAc                                     <  7:145080/146‑1 (MQ=255)
                                                                                                                               aaaaaGGGAACCATCAGGTTCCCTTTTGCGTTAGTGCCGGAGGCCGCATTGGCGCGTAACGTCGGATGCGACGCTGGCGCGTCTACTCCGACCTACTGCGAATTAATCTTTCAGCTTGCGCATTACCAGCGTGGCGTTGGTGccaccg               >  1:186996/1‑148 (MQ=255)
                                                                                                                                aaaaGGGAACCTGATGGTTCCCTTTTGCGTTAGTGCCGGAGGCCGCATTGGCGCGTAACGTCGGATGCGACGCTGGCGCGTCTACTCCGACCTACTGCGAATTAATCTTTCAGCTTGCGCATTACCAGCGt                               <  7:2205/131‑1 (MQ=255)
                                                                                                                                aaaaGGGAACCTGATGGTTCCCTTTTGCGTTAGTGCCGGAGGCCGCATTGGCGCGTAACGTCGGATGCGACGCTGGCGCGTCTACTCCGACCTACTGCGAATTAATCTTTCAGCTTGCGCATTACCAGCGt                               >  8:2205/1‑131 (MQ=255)
                                                                                                                                  aaGGGAACCTGATGGTTCCCTTTTGCGTTAGTGCCGGAGGCCGCATTGGCGCGTAACGTCGGATGCGACGCTGGCGCGTCTACTCCGACCTACTGCGaa                                                             >  6:48382/1‑99 (MQ=38)
                                                                                                                                  aaGGGAACCTGATGGTTCCCTTTTGCGTTAGTGCCGGAGGCCGCATTGGCGCGTAACGTCGGATGCGACGCTGGCGCGTCTACTCCGACCTACTGCGaa                                                             <  5:48382/99‑1 (MQ=38)
                                                                                                                                   aGGGAACCATCAGGTTCCCTTTTGCGTTAGTGCCGGAGGCCGCATTGGCGCGTAACGTCGGATGCGACGCTGGCGCGTCTACTCCGACCTACTGCGAATTAATCTTTCAGCTTGCGCATTACCAGCGTGGCGTTGGTGCCgccgaagcc          >  3:164618/1‑149 (MQ=255)
                                                                                                                                         ccATCAGGTTCCCTTTTGCGTTAGTGCCGGAGGCCGCATTGGCGCGTAACGTCGGATGCGACGCTGGCGCGTCTACTCCGACCTACTGCGAATTAATCTTTCAGCTTGCGCATTACCAGCGTGGCGTTGGTGCCGCCGAAGCCGAAGCt    <  4:67260/149‑1 (MQ=255)
                                                                                                                                          cATCAGGTTCCCTTTTGCGTTAGTGCCGGAGGCCGCATTGGCGCGTAACGTCGGATGCGACGCTGGCGCGTCTACTCCGACCTACTGCGAATTAATCTTTCAGCTTGCGCATTACCAGCGTGGCGTTGGTGCCGCCGAAGCCGAAGCTg   <  4:164618/149‑1 (MQ=255)
                                                                                                                                            tCAGGTTCCCTTTTGCGTTAGTGCCGGAGGCCGCATTGGCGCGTAACGTCGGATGCGACGCTGGCGCGTCTACTCCGACCTACTGCGAATTAATCTTTCAGCTTGCGCATTACCAGCGTGGCGTTGGTGCCGCCGAAGCCGAAGCTGt  >  7:132918/1‑148 (MQ=255)
                                                                                                                                             |                                                                                                                                                  
TCCCCGTCTCATCTCTCTGGTTTCAGGGTTACGGTGCGTTGGCAGGATTTAACGCGTACGTCTTTTCAGAAGGAAATCGACAAAGCGGGAAGTTTGCCTGGAACTGGCGGCGATTGTCAATGATGTGAAAAAGGGAACCATCAGGTTCCCTTTTGCGTTAGTGCCGGAGGCCGCATTGGCGCGTAACGTCGGATGCGACGCTGGCGCGTCTACTCCGACCTACTGCGAATTAATCTTTCAGCTTGCGCATTACCAGCGTGGCGTTGGTGCCGCCGAAGCCGAAGCTGT  >  NC_000913/2440156‑2440443

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