Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,097,3921.G47.4% 7.4 / 38.7 19intergenic (+563/+173)ycdU/serXputative inner membrane protein/tRNA‑Ser
Reads supporting (aligned to +/- strand):  ref base . (6/4);  new base G (5/4);  total (11/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 = 7.13e-01
Rejected as consensus: E-value score below prediction cutoff.
Rejected as consensus: Frequency below/above cutoff threshold.

TCCAACTGGCTGTCGCCAGTTGTCGAACCCCGGTCGGGGCTTCTCATCCCCCCGGTGTGTGCAATATACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGA‑GGGGGGGATTGACTCGCTTCGCTCGCCCTGCGGGCAGCCCGCTCACTGCGTTCACGGTCTGTCCAAC  >  NC_000913/1097273‑1097459
                                                                                                                        |                                                                   
tCCAACTGGCTGTCGCCAGTTGTCGAACCCCGGTCGGGGCTTCTCATCCCCCCGGTGTGTGCAATATACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGA‑GGGGGGGATTGACtcgctt                                                  >  1:568415/1‑139 (MQ=35)
          tGTCGCCAGTTGTCGAACCCCGGTCGGGGCTTCTCATCCCCCCGGTGTGTGCAATATACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGA‑GGGGGGGATTGACTCGCTTCGCTCGCCCt                                        >  1:734020/1‑139 (MQ=35)
                           ccccGGTCGGGGCTTCTCATCCCCCCGGTGTGTGCAATATACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGA‑GGGGGGGATTGACTCGCTTCGCTCGCCCTGCGGGCAGCCCGCTCAc                       <  1:1366997/139‑1 (MQ=32)
                            cccGGTCGGGGCTTCTCATCCCCCCGGTGTGTGCAATATACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGA‑GGGGGGGATTGACTCGCTTCGCTCGCCCTGCGGGCAGCCCGCTCACt                      <  1:1356131/139‑1 (MQ=32)
                             ccGGTCGGGGCTTCTCATCCCCCCGGTGTGTGCAATATACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGGATTGACTCGCTTCGCTCGCCCTGCGGGCAGCCCGCTCACt                      >  2:493543/1‑139 (MQ=17)
                              cGGTCGGGGCTTCTCATCCCCCCGGTGTGTGCAATATACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGGATTGACTCGCTTCGCTCGCCCTGCGGGCAGCCCGCTCACTg                     <  2:1146480/139‑1 (MQ=17)
                              cGGTCGGGGCTTCTCATCCCCCCGGTGTGTGCAATATACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGGATTGACTCGCTTCGCTCGCCCTGCGGGCAGCCCGCTCACTg                     >  2:387076/1‑139 (MQ=17)
                                gTCGGGGCTTCTCATCCCCCCGGTGTGTGCAATATACGAAAAAAAAGCCCGTACTTTCGTACGCGCTCTTCTTTAAATATGGCGGTGA‑GGGGGGGATTGACTCGCTTCGCTCGCCCTGCGGGCAGCCCGCTCACTGCGt                  >  2:1455922/1‑139 (MQ=17)
                                  cGGGGCTTCTCATCCCCCCGGTGTGTGCAATATACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGA‑GGGGGGGATTGACTCGCTTCGCTCGCCCTGCGGGCAGCCCGCTCACTGCGTTc                <  2:734020/139‑1 (MQ=33)
                                        ttCTCATCCCCCCGGTGTGTGCAATATACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGGATTGACTCGCTTCGCTCGCCCTGTGGGCAGCCCGCTCACTGCGTTCACGGt           <  2:1472674/139‑1 (MQ=12)
                                          ctcATCCCCCCGGTGTGTGCAATATACGAAAAAAAAGCCCGTACTTTCGTCCGAGCTCTTCTTTAAATATGGCGGTGA‑GGGGGGGATTGACTCGCTTCGCTCGCCCTGCGGGCAGCCCGCTCACTGCGTTCACGGTCTg        >  2:1457131/1‑139 (MQ=17)
                                          ctcATCCCCCCGGTGTGTGCAATATACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGGATTGACtcgc                                                    >  2:320736/1‑96 (MQ=255)
                                          ctcATCCCCCCGGTGTGTGCAATATACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGGATTGACtcgc                                                    <  1:320736/96‑1 (MQ=16)
                                              tCCCCCCGGTGTGTGCAATATACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGGATTGACTCGCTTCGCTCGCCCTGCGGGCAGCCCGCTCACTGCGTTCACGGTCt         <  1:650395/135‑1 (MQ=17)
                                              tCCCCCCGGTGTGTGCAATATACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGGATTGACTCGCTTCGCTCGCCCTGCGGGCAGCCCGCTCACTGCGTTCACGGTCt         >  2:650395/1‑135 (MQ=17)
                                               ccccccGGTGTGTGCAATATACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGA‑GGGGGGGATTGACTCGCTTCGCTCGCCCTGCGGGCAGCCCGCTCACTGCGTTCACGGTCTGTCCaa   >  1:739102/1‑139 (MQ=32)
                                               ccccccGGTGTGTGCAATATACG‑AAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGA‑‑GGGGGGATTGACTCGCTTCGCTCGCCCTGCGGGCAGc                               >  2:264203/1‑109 (MQ=255)
                                               ccccccGGTGTGTGCAATATACG‑AAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGA‑‑GGGGGGATTGACTCGCTTCGCTCGCCCTGCGGGCAGc                               <  1:264203/109‑1 (MQ=255)
                                                 ccccGGTGTGTGCAATATACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGGATTGACTCGCTTCGCTCGCCCTGCGGGCAGCCCGCTCACTGCGTTCACGGTCTGTCCAAc  >  1:1255600/1‑139 (MQ=11)
                                                                                                                        |                                                                   
TCCAACTGGCTGTCGCCAGTTGTCGAACCCCGGTCGGGGCTTCTCATCCCCCCGGTGTGTGCAATATACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGA‑GGGGGGGATTGACTCGCTTCGCTCGCCCTGCGGGCAGCCCGCTCACTGCGTTCACGGTCTGTCCAAC  >  NC_000913/1097273‑1097459

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