Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,438,4660GC45.8% 24.8 / 13.7 24intergenic (‑29/‑89)PP_5648/PP_3928hypothetical protein/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (12/1);  new base C (0/11);  total (12/12)
Fisher's exact test for biased strand distribution p-value = 9.61e-06
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

AGTCGGCTTGCCGCCACGCGCCGTTTCCATGGTCGCAGAAGCACAACCGGCCAGCATCAGCACCGCAGCCAACGGCAGAATCGCACGCATGTTCCATCCTCGGGGTCGCAAAAAAGAAGGCGACTCTAGCACTGGCAGGCATCCCGAGTGAATTGCCGCGCACTTGCCGCTGTCGCAGATAAAACATGCCATCACGAGGAACCCATCATGAAACTGCACTATCTGCTACTCGCCGCCCTGCCTTGCC  >  NC_002947/4438348‑4438594
                                                                                                                      |                                                                                                                                
aGTCGGCTTGCCGCCACGCGCCGTTTCCATGGTCGCAGAAGCACAACCGGCCAGCATCAGCACCGCAGCCAACGGCAGAATCGCACGCATGTTCCATCCTCGGGGTCGCAAAAAAGAAGGCGACTCTAGCACTggc                                                                                                                 >  1:470689/1‑136 (MQ=255)
     gCTTGCCGCCACGCGCCGTTTCCATGGTCGCAGAAGCACAACCGGCCAGCATCAGCACCGCAGCCAACGGCAGAATCGCACGCATGTTCCATCCTCGGGGTCGCAAAAAAGAAGGCGACTCTAGCACTggcaggca                                                                                                            >  2:387713/1‑136 (MQ=255)
                                   cacacccacaCCCGGCCACCAACAGCACCGCACCCACCGGCATAATCGCCCGCTTGTTCCATCCGCGGGGTCCCAAAAACGACCGCGACTCTAGCACTGGCAGGCATCCCGAGTGAATTGCCGCGCACTTGCCGCt                                                                              <  1:141337/130‑1 (MQ=255)
                                    aGAAGCACAACCGGCCAGCATCAGCACCGCAGCCAACGGCAGAATCGCACGCATGTTCCATCCTCGGGGTCGCAAAAAAGAAGGCGACTCTAGCACTGGCAGGCATCCCGAGTGAATTGCCGCGCACTTGCCGCTg                                                                             >  1:38660/1‑136 (MQ=255)
                                        caacaACCGGCCAGCATCCGCACCCCCGCAAACGGCAGAATCTCACGCATGTTCCATCGTCGGGGACCAAAAAATGACCGCGACTCTAGCACTGGCAGGCATCCCGAGTGAATTGCCGCGCACTTGCCGCTGTCGc                                                                         <  2:317081/134‑1 (MQ=255)
                                                 gCCAGCATCAGCACCGCAGCCAACGGCAGAATCGCACGCATGTTCCATCCTCGGGGTCGCAAAAAAGAAGGCGACTCTAGCACTGGCAGGCATCCCGAGTGAATTGCCGCGCACTTGCCGCTGTCGCAGATAAAAc                                                                >  1:19630/1‑136 (MQ=255)
                                                  ccaccaTCACCACCCCACCCAACGGCAGAATCCCACGCATGTACCATCCTCGGGGTCCAAAAAAAGACCGCGACTCTAGCACTGGCAGGCATCCCGAGTGAATTGCCGCGCACTTGCCGCTGTCGCAGATAAAACa                                                               <  2:38660/136‑1 (MQ=255)
                                                     caaTCACCCCCCCAGCAAACGCCATAAGCGCACGCTTGTACCATCCTGGGGGTCCAAAAAACCACCGCGACTCTAGCACTGGCAGGCATCCCGAGTGAATTGCCGCGCACTTGCCGCTGTCGCAGATAAAACATGc                                                            <  1:529681/134‑1 (MQ=255)
                                                                     cAACGTCATAATCGCACGTTTGTGCAAGGGTGGGGGTCGAAAAAAACCCCGCGACTCTAGCACTGGCAGGCATCCCGAGTGAATTGCCGCGCACTTGCCGCTGTCGCAGATAAAACATGCCATCACGAGGAACCca                                            <  1:254208/136‑1 (MQ=255)
                                                                            aGAATCGCACGCATGTTCCATCCTCGGGGTCGCAAAAAAGAAGGCGACTCTAGCACTGGCAGGCATCCCGAGTGAATTGCCGCGCACTTGCCGCTGTCGCAGATAAAACATGCCATCACGAGGAACCCATCATGaa                                     >  1:472403/1‑136 (MQ=255)
                                                                               aTCTCACTCTTGTTCCATCCTGGGGTCCAAAAAAAATACCGCGACTCTAGCACTGGCAGGCATCCCGAGTGAATTGCCGCGCACTTGCCGCTGTCGCAGATAAAACATGCCATCACGAGGAACCCATCATGAAACt                                  <  1:300545/136‑1 (MQ=255)
                                                                                 cgcacgcaTTTTCCATCCTCGGGGTCGAAAAAAAGATCGCGACTCTAGCACTGGCAGGCATCCCGAGTGAATTGCCGCGCACTTGCCGCTGTCGCAGATAAAACATGCCATCACGAGGAACCCATCATGAAACTGc                                <  2:537507/136‑1 (MQ=255)
                                                                                 cgcacgcaTGTTCCATCCTCGGGGTCGCAAAAAAGAAGGCGACTCTAGCACTGGCAGGCATCCCGAGTGAATTGCCGCGCACTTGCCGCTGTCGCAGATAAAACATGCCATCACGAGGAACCCATCATGAAACTGc                                >  2:485662/1‑136 (MQ=255)
                                                                                 cgcacgcaTGTCCCATCCTCGGGGTCCCAAAAAAGAACGCGACTCTAGCACTGGCAGGCATCCCGAGTGAATTGCCGCGCACTTGCCGCTGTCGCAGATAAAACATGCCATCACGAGGAACCCATCATGAAACTGc                                <  2:451425/136‑1 (MQ=255)
                                                                                   cacgcaTGTTCAATCCTGGGGGTCTAAAAAAAGACCGCGACTCTAGCACTGGCAGGCATCCCGAGTGAATTGCCGCGCACTTGCCGCTGTCGCAGATAAAACATGCCATCACGAGGAACCCATCATGAAACTGCAc                              <  2:19630/136‑1 (MQ=255)
                                                                                      gcaTGTTCCATCCTCGGGGTCGCAAAAAAGAAGGCGACTCTAGCACTGGCAGGCATCCCGAGTGAATTGCCGCGCACTTGCCGCTGTCGCAGATAAAACATGCCATCACGAGGAACCCATCATGAAACTGCACTAt                           >  1:483456/1‑136 (MQ=255)
                                                                                                  gtgggggTCTAAAAAAATAACGCGACTCTAGCACTGGCAGGCATCCCGAGTGAATTGCCGCGCACTTGCCGCTGTCGCAGATAAAACATGCCATCACGAGGAACCCATCATGAAACTGCACTATCTGCTACTcgcc               <  2:423490/133‑1 (MQ=255)
                                                                                                  ccgggggTCGCAAAAAACAACGCGACTCTAGCACTGGCAGGCATCCCGAGTGAATTGCCGCGCACTTGCCGCTGTCGCAGATAAAACATGCCATCACGAGGAACCCATCATGAAACTGCACTATCTGCTACTcgcc               <  2:201375/133‑1 (MQ=255)
                                                                                                   tCGGGGTCGCAAAAAAGAAGGCGACTCTAGCACTGGCAGGCATCCCGAGTGAATTGCCGCGCACTTGCCGCTGTCGCAGATAAAACATGCCATCACGAGGAACCCATCATGAAACTGCACTATCTGCTACTcgccg              >  1:193124/1‑136 (MQ=255)
                                                                                                   tCGGGGTCGCAAAAAAGAAGGCGACTCTAGCACTGGCAGGCATCCCGAGTGAATTGCCGCGCACTTGCCGCTGTCGCAGATAAAACATGCCATCACGAGGAACCCATCATGAAACTGCACTATCTGCTACTcgccg              >  2:422532/1‑136 (MQ=255)
                                                                                                   tCGGGGTCGCAAAAAAGAAGGCGACTCTAGCACTGGCAGGCATCCCGAGTGAATTGCCGCGCACTTGCCGCTGTCGCAGATAAAACATGCCATCACGAGGAACCCATCATGAAACTGCACTATCTGCTACTcgccg              >  1:423490/1‑136 (MQ=255)
                                                                                                     ggggTCGCAAAAAAGAAGGCGACTCTAGCACTGGCAGGCATCCCGAGTGAATTGCCGCGCACTTGCCGCTGTCGCAGATAAAACATGCCATCACGAGGAACCCATCATGAAACTGCACTATCTGCTACTCgccgcc            >  1:359728/1‑136 (MQ=255)
                                                                                                              aaaaaGAAGGCGACTCTAGCACTGGCAGGCATCCCGAGTGAATTGCCGCGCACTTGCCGCTGTCGCAGATAAAACATGCCATCACGAGGAACCCATCATGAAACTGCACTATCTGCTACTCGCCGCCCtgccttgc   <  2:359728/136‑1 (MQ=255)
                                                                                                               aaaaGAAGGCGACTCTAGCACTGGCAGGCATCCCGAGTGAATTGCCGCGCACTTGCCGCTGTCGCAGATAAAACATGCCATCACGAGGAACCCATCATGAAACTGCACTATCTGCTACTCGCCGCCCtgccttgcc  >  1:7292/1‑136 (MQ=255)
                                                                                                                      |                                                                                                                                
AGTCGGCTTGCCGCCACGCGCCGTTTCCATGGTCGCAGAAGCACAACCGGCCAGCATCAGCACCGCAGCCAACGGCAGAATCGCACGCATGTTCCATCCTCGGGGTCGCAAAAAAGAAGGCGACTCTAGCACTGGCAGGCATCCCGAGTGAATTGCCGCGCACTTGCCGCTGTCGCAGATAAAACATGCCATCACGAGGAACCCATCATGAAACTGCACTATCTGCTACTCGCCGCCCTGCCTTGCC  >  NC_002947/4438348‑4438594

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