Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,148,0230GC26.1% 32.2 / 14.3 24P68P (CCG→CCCarcD‑IIarginine/ornithine antiporter
Reads supporting (aligned to +/- strand):  ref base G (10/7);  new base C (2/4);  total (12/12)
Fisher's exact test for biased strand distribution p-value = 3.71e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.87e-01
Rejected as consensus: Frequency below/above cutoff threshold.

GCGGGATCTTCTCGCTGCCGCAGAACATGGCCGCGCGTGCTGATGTGGGCGCAGTGCTGATCGGCTGGGGGATCACCGCGGTCGGCATGTTGGCCCTGGCCTTCGTGTTCCAGACCCTGGCCAACCGCAAGCCGGAACTGGACTCGGGGGTATATGCCTACGCCAAGGCCGGCTTTGGCGAGTACATGGGCTTTTCCTCGGCCTGGGGCTACTGGATCAGCGCCTGGCTGGGCAACGTCGGCTATTTCGTGCTGCTGTTCAGCA  >  NC_002947/1147890‑1148153
                                                                                                                                     |                                                                                                                                  
gcggGATCTTCTCGCTGCCGCAGAACATGGCCGCGCGTGCTGATGTGGGCGCAGTGCTGATCGGCTGGGGGATCACCGCGGTCGGCATGTTGGCCCTGGCCTTCGTGTTCCAGACCCTGGCCAACCGCAAGCCCGa                                                                                                                                  <  1:208183/136‑1 (MQ=255)
gcggGATCTTCTCGCTGCCGCAGAACATGGCCGCGCGTGCTGATGTGGGCGCAGTGCTGATCGGCTGGGGGATCACCGCGGTCGGCATGTTGGCCCTGGCCTTCGTGTTCCAGACCCTGGCCAACCGCAAGCCCGa                                                                                                                                  <  2:389129/136‑1 (MQ=255)
              cTGCCGCAGAACATGGCCGCGCGTGCTGATGTGGGCGCAGTGCTGATCGGCTGGGGGATCACCGCGGTCGGCATGTTGGCCCTGGCCTTCGTGTTCCAGACCCTGGCCAACCGCAAGCCGGAACTGGACTCggggg                                                                                                                    <  1:321970/136‑1 (MQ=255)
                              ccGCGCGTGCTGATGTGGGCGCAGTGCTGATCGGCTGGGGGATCACCGCGGTCGGCATGTTGGCCCTGGCCTTCGTGTTCCAGACCCTGGCCAACCGCAAGCCGGAACTGGACTCGGGGGTATATGCCTACGCCaa                                                                                                    >  2:84328/1‑136 (MQ=255)
                              ccGCGCGTGCTGATGTGGGCGCAGTGCTGATCGGCTGGGGGATCACCGCGGTCGGCATGTTGGCCCTGGCCTTCGTGTTCCAGACCCTGGCCAACCGCAAGCCGGAACTGGACTCGGGGGTATATGCCTACGCCaa                                                                                                    >  2:294562/1‑136 (MQ=255)
                              ccGCGCGTGCTGATGTGGGCGCAGTGCTGATCGGCTGGGGGATCACCGCGGTCGGCATGTTGGCCCTGGCCTTCGTGTTCCAGACCCTGGCCAACCGCAAGCCGGAACTGGACTCGGGGGTATATGCCTACGCCaa                                                                                                    >  1:320797/1‑136 (MQ=255)
                              ccGCGCGTGCTGATGTGGGCGCAGTGCTGATCGGCTGGGGGATCACCGCGGTCGGCATGTTGGCCCTGGCCTTCGTGTTCCAGACCCTGGCCAACCGCAAGCCCGa                                                                                                                                  >  1:150261/1‑106 (MQ=255)
                              ccGCGCGTGCTGATGTGGGCGCAGTGCTGATCGGCTGGGGGATCACCGCGGTCGGCATGTTGGCCCTGGCCTTCGTGTTCCAGACCCTGGCCAACCGCAAGCCCGa                                                                                                                                  <  2:307246/106‑1 (MQ=255)
                              ccGCGCGTGCTGATGTGGGCGCAGTGCTGATCGGCTGGGGGATCACCGCGGTCGGCATGTTGGCCCTGGCCTTCGTGTTCCAGACCCTGGCCAACCGCAAGCCCGa                                                                                                                                  >  1:307246/1‑106 (MQ=255)
                              ccGCGCGTGCTGATGTGGGCGCAGTGCTGATCGGCTGGGGGATCACCGCGGTCGGCATGTTGGCCCTGGCCTTCGTGTTCCAGACCCTGGCCAACCGCAAGCCCGa                                                                                                                                  <  2:150261/106‑1 (MQ=255)
                                                              ggCTGGGGGATCACCGCGGTCGGCATGTTGGCCCTGGCCTTCGTGTTCCAGACCCTGGCCAACCGCAAGCCg                                                                                                                                    <  2:254682/72‑1 (MQ=255)
                                                              ggCTGGGGGATCACCGCGGTCGGCATGTTGGCCCTGGCCTTCGTGTTCCAGACCCTGGCCAACCGCAAGCCg                                                                                                                                    >  1:254682/1‑72 (MQ=255)
                                                                 tGGGGGATCACCGCGGTCGGCATGTTGGCCCTGGCCTTCGTGTTCCAGACCCTGGCCAACCGCAAGCCGGAACTGGACTCGGGGGTATATGCCTACGCCAAGGCCGGCTTTGGCGAGTACATGGGCTTTTCCTCgg                                                                 >  2:83294/1‑136 (MQ=255)
                                                                             gcgGTCGGCATGTTGGCCCTGGCCTTCGTGTTCCAGACCCTGGCCAACCGCAAGCCGGAACTGGACTCGGGGGTATATGCCTACGCCAAGGCCGGCTTTGGCGAGTACATGGGCTTTTCCTCGGCCTGGGGCTACt                                                     >  2:196145/1‑136 (MQ=255)
                                                                                 tCGGCATGTTGGCCCTGGCCTTCGTGTTCCAGACCCTGGCCAACCGCAAGCCGGAACTGGACTCGCGGGTATATGCCTACGCCAAGGCCGGCTTTGGCGAGTACATGGGCATTTCCTCGGCCTGGGGCTACTGGAt                                                 >  2:140830/1‑136 (MQ=255)
                                                                                                    cTTCGTGTTCCAGACCCTGGCCAACCGCAAGCCGGAACTGGACTCGGGGGTATATGCCTACGCCAAGGCCgg                                                                                              <  2:30380/72‑1 (MQ=255)
                                                                                                    cTTCGTGTTCCAGACCCTGGCCAACCGCAAGCCGGAACTGGACTCGGGGGTATATGCCTACGCCAAGGCCgg                                                                                              >  1:30380/1‑72 (MQ=255)
                                                                                                    cTTCGTGTTCCAGACCCTGGCCAACCGCAAGCCGGAACTGGACTCGGGGGTATATGCCTACGCCAAGGCCGGCTTTGGCGAGTACATGGGCTTTTCCTCGGCCTGGGGCTACTGGATCAGCGCCTGGCTGGGCAAc                              >  2:16366/1‑136 (MQ=255)
                                                                                                    cTTCGTGTTCCAGACCCTGGCCAACCGCAAGCCGGAACTGGACTCGGGGGTATATGCCTACGCCAAGGCCGGCTTTGGCGAGTACATGGGCTTTTCCTCGGCCTGGGGCTACTGGATCAGCGCCTGGCTGGGCAAc                              >  1:28661/1‑136 (MQ=255)
                                                                                                                         cAACCGCAAGCCGGAACTGGACTCGGGGGTATATGCCTACGCCAAGGCCGGCTTTGGCGAGTTCATGGGCTTTTCCTCGGCCTGGGGCTACTGGATCAGCGCCTGGCTGGGCAACGTCGGCTATTTCGTGctgctg         >  1:330700/1‑136 (MQ=255)
                                                                                                                                aaGCCTGAACTGGACTCGGGGGTATATGCCTACGCCAAGGCCGGCTTTGGCGAGTACATGGGCTTTTCCTCGGCCTGGGGCTACTGGATCAGCGCCTGGCTGGGCAACGTCGGCTACTTCGTGCTGCTGTTCAGCa  <  1:83294/136‑1 (MQ=255)
                                                                                                                                aaGCCGGAACTGGACTCGGGGGTATATGCCTACGCCAAGGCCGGCTTTGGCGAGTTCATGGGCTTTTCCTCGGCCTGGGGCTACTGGATCAGCGCCTGGCTGGGCAACGTCGGCTATTTCGTGCTGCTGTTCAGCa  <  2:330700/136‑1 (MQ=255)
                                                                                                                                aaGCCGGAACTGGACTCGGGGGTATATGCCTACGCCAAGGCCGGCTTTGGCGAGTACATGGGCTTTTCCTCGGCCTGGGGCTACTGGATCAGCGCCTGGCTGGGCAACGTCGGCTATTTCGTGCTGCTGTTCAGCa  <  1:84328/136‑1 (MQ=255)
                                                                                                                                aaGCCGGAACTGGACTCGGGGGTATATGCCTACGCCAAGGCCGGCTTTGGCGAGTACATGGGCTTTTCCTCGGCCTGGGGCTACTGGATCAGCGCCTGGCTGGGCAACGTCGGCTATTTCGTGCTGCTGTTCAGCa  <  2:28661/136‑1 (MQ=255)
                                                                                                                                aaGCCGGAACTGGACTCGGGGGTATATGCCTACGCCAAGGCCGGCTTTGGCGAGTACATGGGCTTTTCCTCGGCCTGGGGCTACTGGATCAGCGCCTGGCTGGGCAACGTCGGCTATTTCGTGCTGCTGTTCAGCa  <  2:320797/136‑1 (MQ=255)
                                                                                                                                aaGCCGGAACTGGACTCGGGGGTATATGCCTACGCCAAGGCCGGCTTTGGCGAGTACATGGGCTTTTCCTCGGCCTGGGGCTACTGGATCAGCGCCTGGCTGGGCAACGTCGGCTATTTCGTGCTGCTGTTCAGCa  <  1:294562/136‑1 (MQ=255)
                                                                                                                                     |                                                                                                                                  
GCGGGATCTTCTCGCTGCCGCAGAACATGGCCGCGCGTGCTGATGTGGGCGCAGTGCTGATCGGCTGGGGGATCACCGCGGTCGGCATGTTGGCCCTGGCCTTCGTGTTCCAGACCCTGGCCAACCGCAAGCCGGAACTGGACTCGGGGGTATATGCCTACGCCAAGGCCGGCTTTGGCGAGTACATGGGCTTTTCCTCGGCCTGGGGCTACTGGATCAGCGCCTGGCTGGGCAACGTCGGCTATTTCGTGCTGCTGTTCAGCA  >  NC_002947/1147890‑1148153

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