Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,148,8000TG50.0% 0.5 / 12.8 26T151P (ACG→CCG) PP_1007transmembrane sensor
Reads supporting (aligned to +/- strand):  ref base T (8/5);  new base G (13/0);  total (21/5)
Fisher's exact test for biased strand distribution p-value = 3.91e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: E-value score below prediction cutoff.
Rejected as consensus: Frequency below/above cutoff threshold.

AGTAAGCCTCGCCCTGTAGCAGGCGGGCGACCTGGCGGCCGTCCTGGCGTTCGCTGGCGAATGCCGAATTGGTGTTGAGCAGCACTTTTGCGCCGTCTTCCAGTTGCAGGCGCTGGCGCTCGCCGACCACCGTAAGGTGGTCGGCCTGCAGGCGCATGGGCAGGTTACCGACGGTGAATAGCCCGACCAGCAGCAGCGCTGCGGTGGCCAGGGGTTTCCAGTGACTGCGCAGGCGCCCGCGCCATGACCTGCGTTGC  >  NC_002947/1148668‑1148924
                                                                                                                                    |                                                                                                                             
aGTAAGCCTCGCCCTGTAGCAGGCGGGCGACCTGGCGGCCGTCCTGGCGTTCGCTGGCGAATGCCGAATTGGTGTTGAGCAGCACTTTTGCGCCGTCTTCCAGTTGCAGGCGCTGGCGCTCGCCGCCCCCCCTAAg                                                                                                                            >  1:107891/1‑136 (MQ=255)
aGTAAGCCTCGCCCTGTAGCAGGCGGGCGACCTGGCGGCCGTCCTGGCGTTCGCTGGCGAATGCCGAATTGGTGTTGAGCAGCACTTTTGCGCCGTCTTCCAGTTGCAGGCGCTGGCGCTCGCCGACCACCGGAAg                                                                                                                            >  2:730902/1‑136 (MQ=255)
aGTAAGCCTCGCCCTGTAGCAGGCGGGCGACCTGGCGGCCGTCCTGGCGTTCGCTGGCGAATGCCGAATTGGTGTTGAGCAGCACTTTTGCGCCGTCTTCCAGTTGCAGGCGCTGGCGCTCGCCCACCACCGTaat                                                                                                                            >  1:511796/1‑135 (MQ=255)
aGTAAGCCTCGCCCTGTAGCAGGCGGGCGACCTGGCGGCCGTCCTGGCGTTCGCTGGCGAATGCCGAATTGGTGTTGAGCAGCACTTTTGCGCCGGCTTCCAGTTGCAGGCGCTGGCGCTCGCCGCCCACCGGAAg                                                                                                                            >  2:1140845/1‑136 (MQ=255)
             cTGTAGCAGGCGGGCGACCTGGCGGCCGTCCTGGCGTTCGCTGGCGAATGCCGAATTGGTGTTGAGCAGCACTTTTGCGCCGTCTTCCAGTTGCAGGCGCTGGCGCTCGCCGACCACCGTAAGGTGGTCGGCCTGc                                                                                                               >  2:33774/1‑136 (MQ=255)
                        gggcgACCTGGCGGCCGTCCTGGCGTTCGCTGGCGAATGCCGAATTGGTGTTGAGCAGCACTTTTGCGCCGTCTTCCAGTTGCAGGCGCTGGCGCTCGCCGACCCCCGTAAGGTGGTCGGCCAGCAGGCGGATggg                                                                                                   >  2:693824/1‑136 (MQ=255)
                         ggcgACCTGGCGGCCGTCCTGGCGTTCGCTGGCGAATGCCGAATTGGTGTTGAGCAGCACTTTTGCGCTGTCTTCCAGTTGCAGGCGCTGGCGCTCGCCGCCCACCCGTAGGGGGTCGGGCGGCAGGGGGATgggg                                                                                                  >  1:592992/1‑135 (MQ=255)
                         ggcgACCTGGCGGCCGTCCTGGCGTTCGCTGGCGAATGCCGAATTGGTGTTGAGCAGCACTTTTGCGCCGTCTTCCAGTTGCAGGCGCTGGCGCTCGCCGCCCCCCCGAAGGTGGTAGGCCGGCAGGCGCCTgggg                                                                                                  >  1:1005176/1‑135 (MQ=255)
                         ggcgACCTGGCGGCCGTCCTGGCGTTCGCTGGCGAATGCCGAATTGGTGTTGAGCAGCACTTTTGCGCCGTCTTCCAGTTGCAGGCGCTGGCGCTCGCCGCCCACCGGAAAGGGGTGGGCCGGGCGgcgccggggg                                                                                                  >  1:1116978/1‑130 (MQ=255)
                         ggcgACCTGGCGGCCGTCCTGGCGTTCGCTGGCGAATGCCGAATTGGTGTTGAGCAGCACTTTTGCGCCGTCTTCCAGTTGCAGGCGCTGGCGCTCGCCGACCACCGTAAGGTGGTCGGCCTGCAGGCGCATGGGc                                                                                                  >  1:536291/1‑136 (MQ=255)
                             cccgggcggcCGGCCGGGCGTCCGCTGGCGAATGCCGAATTGGTGTTCAGCAGCACTTTTGCGCCGTCTTCCAGTGGCAGGCGCTGGCGCTCGCCGACCACCGTAAGGTGGTCGGCCTGCAGGCGCATGGGCAGGt                                                                                              <  2:498900/132‑1 (MQ=255)
                                     gcCGTCCTGGCGTTCGCTGGCGAATGCCGAATTGGTGTTGAGCAGCACTTTTGCGCCGTCTTCCAGTTGCAGGCGCTGGCGCTCGCCGACCACCGTAAGGTGGTCGGCCTGCAGGCGCATGGGCAGGTTACCGACg                                                                                      >  2:766377/1‑136 (MQ=255)
                                      ccGTCCTGGCGTTCGCTGGCGAATGCCGAATTGGTGTTGAGCAGCACTTTTGCGCCGTCTTCCAGTTGCAGGCGCTGGCGCTCGCCGCCCCCCGGAAGGTGGTCGGCCGGCAGGCGCATGGGCCGGTTACcggcgg                                                                                     >  1:878976/1‑136 (MQ=255)
                                           cTGGCGTTCGCTGGCGAATGCCGAATTGGTGTTGAGCAGCACTTTTGCGCCGTCTTCCAGTTGCAGGCGCTGGCGCTCGCCCCCCACCGGAAGGGGGGCGGCCTGCAGGCGCAGGGGCAGGTTACCCACCGGGaaa                                                                                >  1:1276057/1‑135 (MQ=255)
                                                      tGGCGAATGCCGAATTGGTGTTGAGCAGCACTTTTGCGCCGTCTTCCAGTTGCAGGCGCTGGCGCTCGCCGACCCCCGGAAGGTGGTCGGCCGGCAGGCGGATGGGGCGGTTTCCGACGGGGAATAGCCCGACcac                                                                     >  2:1279389/1‑135 (MQ=255)
                                                      tGGCGAATGCCGAATTGGTGTTGAGCAGCACTTTTGCGCCGTCTTCCAGTTGCAGGCGCTGGCGCTCGCCGACCACCGTAAGGTGGTCGGCCTGCAGGCGCATGGGCAGGTTTCCGCCGGTGAATAGCCCGAcccc                                                                     >  1:314593/1‑134 (MQ=255)
                                                      tGGCGAATGCCGAATTGGTGTTGAGCAGCACTTTTGCGCCGTCTTCCAGTTGCAGGCGCTGGCGCTCGCCCACCCCCGGAATGTGGTGGGCCAGCACGCGGCGTGGGAGGTTTACCGCGGTGTAAATCCCGAcccc                                                                     >  2:942805/1‑134 (MQ=255)
                                                                       gtgtGGTGCTGCACTTAGGCGCCGTCTTCATTTTGCTGGCGCGGGCGGTCGCCGACCACCGTACGGGGGTCGGCCTGCAGGCGCATGGGCAGGTTACCGACGGTGAATAGCCCGACCAGCAGCAGCGCTGCGGTgg                                                    <  1:543866/136‑1 (MQ=255)
                                                                                     ttttGCGCCGTCTTCCAGTTGCAGGCGCTGGCGCTCGCCGCCCCCCCGAAGGTGGTGGGGCGGCAGGCGGAAGGGGCGGTTTGCGAAGGTGTATATCCCGACCCACAGCACCGCAGCGGGGGGCAGGGGGTTTCaa                                      >  2:658990/1‑135 (MQ=255)
                                                                                     ttttGCGCCGTCTTCCAGTTGCAGGCGCTGGCGCTCGCCGCCCCCCCGAAGGTGGTGGGCCGGCACGCGCCGGGGGAGGGTTACGCCGGGGGATATCCCCCCCACCCGCACCGCTGCGGTGGGCAGGGGGTTGCac                                      >  2:235456/1‑135 (MQ=255)
                                                                                     ttttGCGCCGTCTTCCAGTTGCAGGCGCTGGCGCTCGCCGCCCACCCGAAGGTGGTGGGCCGGGAGGCGCATGGGGCGGTTAACGACGGTGAATTTCCCGCCCCGCCGCCCCGCGGGCGTGGGCAGGGGGTTACAg                                      >  2:656855/1‑136 (MQ=255)
                                                                                     ttttGCGCCGTCTTCCAGTTGCAGGCGCTGGCGCTCGCCCCCCCCCGGAAGGGGGGAGGGCAGGACGCGCCTGGGGCGGTTTACGCCCGGGAAAAGCCCCACCAGCCGCCGCGCGGCGGGGGCCAGGGGgttgccc                                      >  1:788463/1‑132 (MQ=255)
                                                                                          cgcCGTCTTCCAGTTGCAGGCGCTGGCGCTCGCCGACCACCGTAAGGTGGTCGGCCTGCAGGCGCATGGGCAGGTTACCGACGGTGAATAGCCCGACCAGCCGCAGCGCTGCGGTGGCCAGGGGTTTCCAGTGACt                                 >  1:477072/1‑136 (MQ=255)
                                                                                                           cggcgctggccccCGCCGACCACCGTACGGGGGGCGGCCTGCAGGCGCATGGGCAGGTTACCGACGGTGAATAGCCCGACCAGCAGCAGCGCTGCGGTGGCCAGGGGTTTCCAGTGACTGCGCAGGCGCCCGCGcc                <  2:347969/135‑1 (MQ=255)
                                                                                                                       ccGCCGACCCCGCTACGGGGGGCGGCCGGCAGGCGCATGGGCAGGTTACCGACGGTGAATAGCCCGACCAGCAGCAGCGCTGCGGTGGCCAGGGGGTTCCAGTGACTGCGCAGGCGCCCGCGCCATGACCTGCGtt    <  2:1280377/135‑1 (MQ=255)
                                                                                                                         gCCGACCACCGTAAGGTGGTCGGCCTGCAGGCGCATGGGCAGGTTACCGACGGTGAATAGCCCGACCAGCAGCAGCGCTGCGGTGGCCAGGGGTTTCCAGTGACTGCGCAGGCGCCCGCGCCATGACCTGCGTtgc  <  2:536291/136‑1 (MQ=255)
                                                                                                                                    |                                                                                                                             
AGTAAGCCTCGCCCTGTAGCAGGCGGGCGACCTGGCGGCCGTCCTGGCGTTCGCTGGCGAATGCCGAATTGGTGTTGAGCAGCACTTTTGCGCCGTCTTCCAGTTGCAGGCGCTGGCGCTCGCCGACCACCGTAAGGTGGTCGGCCTGCAGGCGCATGGGCAGGTTACCGACGGTGAATAGCCCGACCAGCAGCAGCGCTGCGGTGGCCAGGGGTTTCCAGTGACTGCGCAGGCGCCCGCGCCATGACCTGCGTTGC  >  NC_002947/1148668‑1148924

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