Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,482,4800C.26.7% 36.7 / 10.1 15coding (1192/1914 nt)yheSABC‑F family protein predicted regulatory ATPase
Reads supporting (aligned to +/- strand):  ref base C (5/6);  new base . (2/2);  total (7/8)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.53e-01
Rejected as consensus: Frequency below/above cutoff threshold.

ATTAATCAAACTGTTAGCCGGTGAACTTGCGCCAGTCAGCGGTGAAATTGGTCTGGCGAAAGGGATCAAACTCGGCTACTTCGCCCAGCATCAACTTGAATACCTGCGCGCCGACGAATCACCTATTCAACATCTGGCACGTTTAGCGCCGCAGGAGCTGGAACAAAAACTGCGTGACTACCTCGGCGGCTTTGGTTTCCAGGGCGATAAAGTAACCGAAGAAACGCGCCGCTTCTCCGGTGGGGAAAAAGCCCGCCTGGTGCTGGCATTA  >  NC_000913/3482347‑3482617
                                                                                                                                     |                                                                                                                                         
aTTAATCAAACTGTTAGCCGGTGAACTTGCGCCAGTCAGCGGTGAAATTGGTCTGGCGAAAGGGATCAAACTCGGCTACTTCGCCCAGCATCAACTTGAATACCTGCGCGCCGACGAATCACCTATTCAACATCTGGCa                                                                                                                                      >  1:565665/1‑139 (MQ=255)
         aCTGTTAGCCGGTGAACTTGCGCCAGTCAGCGGTGAAATTGGTCTGGCGAAAGGGATCAAACTCGGCTACTTCGCCCAGCATCAACTTGAATACCTGCGCGCCGACGAATCACCTATTCAACATCTGGCACGTTTAgcg                                                                                                                             <  1:607251/139‑1 (MQ=255)
                               ccAGTCAGCGGTGAAATTGGTCTGGCGAAAGGGATCAAACTCGGCTACTTCGCCCAGCATCAACTTGAATACCTGCGCGCCGACGAATCACCTATTCAACATCTGGCACGTTTAGCGCCGCAGGAGCTGGAACAAAAAc                                                                                                       <  2:565665/139‑1 (MQ=255)
                                                                 tCAAACTCGGCTACTTCGCCCAGCATCAACTTGAATACCTGCGCGCCGACGAATCACCTATTCAACATTGGCACGTTTAg                                                                                                                               <  1:396458/80‑1 (MQ=255)
                                                                 tCAAACTCGGCTACTTCGCCCAGCATCAACTTGAATACCTGCGCGCCGACGAATCACCTATTCAACATTGGCACGTTTAg                                                                                                                               <  1:642125/80‑1 (MQ=255)
                                                                 tCAAACTCGGCTACTTCGCCCAGCATCAACTTGAATACCTGCGCGCCGACGAATCACCTATTCAACATTGGCACGTTTAg                                                                                                                               >  2:396458/1‑80 (MQ=255)
                                                                 tCAAACTCGGCTACTTCGCCCAGCATCAACTTGAATACCTGCGCGCCGACGAATCACCTATTCAACATTGGCACGTTTAg                                                                                                                               >  2:642125/1‑80 (MQ=255)
                                                                                            aaCTTGAATACCTGCGCGCCGACGAATCACCTATTCAACATCTGGCACGTTTAGCGCCGCAGGAGCTGGAACAAAAACTGCGTGACTACCTCGGCGGCTTTGGTTTCCAGGGCGATAAAGTAACCGAAGAAACGcgccg                                          >  2:363189/1‑139 (MQ=255)
                                                                                                         gcgcgcCGACGAATCACCTATTCAACATCTGGCACGTTTAGCGCCGCAGGAGCTGGAACAAAAACTGCGTGACTACCTCGGCGGCTTTGGTTTCCAGGGCGATaaa                                                              <  1:192032/106‑1 (MQ=255)
                                                                                                         gcgcgcCGACGAATCACCTATTCAACATCTGGCACGTTTAGCGCCGCAGGAGCTGGAACAAAAACTGCGTGACTACCTCGGCGGCTTTGGTTTCCAGGGCGATaaa                                                              >  2:192032/1‑106 (MQ=255)
                                                                                                                 acgaATCACCTATTCAACATCTGGCACGTTTAGCGCCGCAGGAGCTGGAACAAAAACTGCGTGACTACCTCGGCGGCTTTGGTTTCCAGGGCGATAAAGTAACCGAAGAAACGCGCCGCTTCTCCGGTGGGGAAAAAGc                     <  1:554864/139‑1 (MQ=255)
                                                                                                                               cAACATCTGGCACGTTTAGCGCCGCAGGAGCTGGAACAAAAACTg                                                                                                     >  1:345609/1‑45 (MQ=255)
                                                                                                                               cAACATCTGGCACGTTTAGCGCCGCAGGAGCTGGAACAAAAACTg                                                                                                     <  2:345609/45‑1 (MQ=255)
                                                                                                                                aaCATCTGGCACGTTTAGCGCCGCAGGAGCTGGAACAAAAACTGCGTGACTACCTCGGCGGCTTTGGTTTCCAGGGCGATAAAGTAACCGAAGAAACGCGCCGCTTCTCCGGTGGGGAAAAAGCCCGCCTGGTGCTGGc      <  1:363189/139‑1 (MQ=255)
                                                                                                                                    tCTGGCACGTTTAGCGCCGCAGGAGCTGGAACAAAAACTGCGTGACTACCTCGGCGGCTTTGGTTTCCAGGGCGATAAAGTAACCGAAGAAACGCGCCGCTTCTCCGGTGGGGAAAAAGCCCGCCTGGTGCTGGCATTa  >  1:390582/1‑139 (MQ=255)
                                                                                                                                     |                                                                                                                                         
ATTAATCAAACTGTTAGCCGGTGAACTTGCGCCAGTCAGCGGTGAAATTGGTCTGGCGAAAGGGATCAAACTCGGCTACTTCGCCCAGCATCAACTTGAATACCTGCGCGCCGACGAATCACCTATTCAACATCTGGCACGTTTAGCGCCGCAGGAGCTGGAACAAAAACTGCGTGACTACCTCGGCGGCTTTGGTTTCCAGGGCGATAAAGTAACCGAAGAAACGCGCCGCTTCTCCGGTGGGGAAAAAGCCCGCCTGGTGCTGGCATTA  >  NC_000913/3482347‑3482617

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