Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,438,4660GC61.5% ‑4.4 / 10.9 13intergenic (‑29/‑89)PP_5648/PP_3928hypothetical protein/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (4/1);  new base C (0/8);  total (4/9)
Fisher's exact test for biased strand distribution p-value = 6.99e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.95e-02
Rejected as consensus: E-value score below prediction cutoff.
Rejected as consensus: Frequency below/above cutoff threshold.

TGCGAGTCGAGCTGCGCAGTCGGCTTGCCGCCACGCGCCGTTTCCATGGTCGCAGAAGCACAACCGGCCAGCATCAGCACCGCAGCCAACGGCAGAATCGCACGCATGTTCCATCCTCGGGGTCGCAAAAAAGAAGGCGACTCTAGCACTGGCAGGCATCCCGAGTGAATTGCCGCGCACTTGCCGCTGTCGCAGATAAAACATGCCATCACGAGGAACCCATCATGAAACTGCACTATCTGCTACTCGCC  >  NC_002947/4438331‑4438581
                                                                                                                                       |                                                                                                                   
tgCGAGTCGAGCTGCGCAGTCGGCTTGCCGCCACGCGCCGTTTCCATGGTCGCAGAAGCACAACCGGCCAGCATCAGCACCGCAGCCAACGGCAGAATCGCACGCATGTTCCATCCTCGGGGTCGCAAAAaagaag                                                                                                                     >  2:299456/1‑136 (MQ=255)
                              ccccgcgCCGTTTCCTTGGCCGCACACGCACAACCGGCCCCCACCACCACCGCCCCCAACGGCATAATCGCCCGCATGTACCACCCGCGGGGCCCCAAAAAACACCGCGACTCTAGCACTGGCAGGCATCCCGAGt                                                                                       <  2:214753/133‑1 (MQ=255)
                                                             acccGGCAATCATCACCACCGCAGCAAAGGGAATAATCCCACGCATGTTCCATCCTCGGGGCCTAAAAATAGACCGCGACTCTAGCACTGGCAGGCATCCCGAGTGAATTGCCGCGCACTTGCCGCTGTCGCAGAt                                                        <  2:357375/134‑1 (MQ=255)
                                                              cccGCCCAGCTTCAGCACCGCCCCCAACGGCAGAACCCCACGATTGTTCCATCCTCGGGGTCGAAAAAAACACCGCGACTCTAGCACTGGCAGGCATCCCGAGTGAATTGCCGCGCACTTGCCGCTGTCGCAGATa                                                       <  2:18383/135‑1 (MQ=255)
                                                                     aGCATCAGCACCGCAGCCAACGGCAGAATCGCACGCATGTTCCATCCTCGGGGTCGCAAAAAAGAAGGCGACTCTAGCACTGGCAGGCATCCCGAGTGAATTGCCGCGCACTTGCCGCTGTCGCAGATAAAACATg                                                >  1:89118/1‑136 (MQ=255)
                                                                               ccGCACCCAAGGTAATAATCCCACGCTTGTTCCATCGTGGGGGTCCAAAAAAATAACGCGACTCTAGCACTGGCAGGCATCCCGAGTGAATTGCCGCGCACTTGCCGCTGTCGCAGATAAAACATGCCATCACGAg                                      <  1:176209/136‑1 (MQ=255)
                                                                                     ccAACGGCAGAATCGCACGCATGTTCCATCCTCGGGGTCGCAAAAAAGAAGGCGACTCTAGCACTGGCAGGCATCCCGAGTGAATTGCCGCGCACTTGCCGCTGTCGCAGATAAAACATGCCATCACGAGGAAccc                                >  2:103727/1‑136 (MQ=255)
                                                                                     ccAACGGCAGAATCGCACGCATGTTCCATCCACGGGGTCGCAAAAAAGAAGGCGACTCTAGCACTGGCAGGCATCCCGAGTGAATTGCCGCGCACTTGCCGCTGTCGCAGATAAAACATGCCATCACGAGGAAccc                                >  1:251203/1‑136 (MQ=255)
                                                                                           ccataaTCGCACGCATTCTCCATCCTCGGGGTCCAAAAAAACACCGCGACTCTAGCACTGGCAGGCATCCCGAGTGAATTGCCGCGCACTTGCCGCTGTCGCAGATAAAACATGCCATCACGAGGAACCCATCATg                          <  2:393282/132‑1 (MQ=255)
                                                                                                aTCGCACGCATGTTCCTTCCTGGGGGTCTAAAAAACGACCGCGACTCTAGCACTGGCAGGCATCCCGAGTGAATTGCCGCGCACTTGCCGCTGTCGCAGATAAAACATGCCATCACGAGGAACCCATCATGAAACt                     <  1:299456/136‑1 (MQ=255)
                                                                                                aTCGCACGCATGTCCCATCCTCGGGGGCCCAAAAAAGAACGCGACTCTAGCACTGGCAGGCATCCCGAGTGAATTGCCGCGCACTTGCCGCTGTCGCAGATAAAACATGCCATCACGAGGAACCCATCATGAAACt                     <  2:206081/136‑1 (MQ=255)
                                                                                                  cgcacgcaTGTTCCATCCTCGGGGTCGAAAAAAAGAACGCGACTCTAGCACTGGCAGGCATCCCGAGTGAATTGCCGCGCACTTGCCGCTGTCGCAGATAAAACATGCCATCACGAGGAACCCATCATGAAACTGc                   <  2:151829/136‑1 (MQ=255)
                                                                                                                           cGCAAAAAAGAAGGCGACTCTAGCACTGGCAGGCATCCCGAGTGAATTGCCGCGCACTTGCCGCTGTCGCAGATAAAACATGCCATCACGAGGAACCCATCATGAAACTGCACTATCTGCTACTcgcc  >  1:208391/1‑128 (MQ=255)
                                                                                                                           cGCAAAAAAGAAGGCGACTCTAGCACTGGCAGGCATCCCGAGTGAATTGCCGCGCACTTGCCGCTGTCGCAGATAAAACATGCCATCACGAGGAACCCATCATGAAACTGCACTATCTGCTACTcgcc  <  2:208391/128‑1 (MQ=255)
                                                                                                                                       |                                                                                                                   
TGCGAGTCGAGCTGCGCAGTCGGCTTGCCGCCACGCGCCGTTTCCATGGTCGCAGAAGCACAACCGGCCAGCATCAGCACCGCAGCCAACGGCAGAATCGCACGCATGTTCCATCCTCGGGGTCGCAAAAAAGAAGGCGACTCTAGCACTGGCAGGCATCCCGAGTGAATTGCCGCGCACTTGCCGCTGTCGCAGATAAAACATGCCATCACGAGGAACCCATCATGAAACTGCACTATCTGCTACTCGCC  >  NC_002947/4438331‑4438581

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