Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*ACB122319,1200CG57.1% 12.4 / 26.8 22P554P (CCG→CCC) CJ019_00291hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (8/1);  new base G (0/12);  total (9/13)
Fisher's exact test for biased strand distribution p-value = 4.42e-05
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.95e-04
Rejected as consensus: Frequency below/above cutoff threshold.

ACCAACAACGCCGTGCAGCAGAACGGTTACAACGACGACGGCAACGAAGTCCGCGTATTCGTCGAGTACCCGATCAGCATCTTCTGATTCTGATCTGATGTATCGAACAGCCCCGGCCTATGCCGGGGCTGTTTCTTTTCGGCGAACGGTTACTGTTTGAGCCGGCGCTGGGCTTACGGTCTCCGGTGCAAAGCTGTCGCTACGCGCCATCCGCCACATCCGCGCATAAAACTCGCTTTCGATCGAGCCACTGAGTAAC  >  ACB122/318997‑319255
                                                                                                                           |                                                                                                                                       
accaccaccGCCGTGCAGCACACCGTTAACAACCACCCCGCCAACCACGGCCGCGTATTCGGCGAGAACCCCATCAGCATCTTCTGATTCTGATCTGATGTATCCAACACCCCCGCCCTATGCCGGGGCTGTTTCt                                                                                                                             <  4:144028/136‑1 (MQ=255)
                agcagAACGGTTACAACGACGACGGCAACGAAGTCCGCGTATTCGTCGAGTACCCGATCAGCATCTTCTGATTCTGATCTGATGTATCGAACAGCCCCGGCCTATGCCGGGGCTGGTTCGTTTCGGCGAACGGtt                                                                                                              >  6:326767/1‑135 (MQ=255)
                          ttACAACGACGACGGCAACGAAGTCCGCGTATTCGTCGAGTACCCGATCAGCATCTTCTGATTCTGATCTGATGTATCGAACAGCACCCGCCAATGCAGGGGGGGTGTGTTTTAGGCGAACCGTTACTGtttttt                                                                                                    >  2:47248/1‑132 (MQ=255)
                                  acgacgGCAACGAAGTCCGCGTATTCGTCGAGTACCCGATCAGCATCTTCTGATTCTGATCTGATGTATCGAACAGCCCCCGCCCATGCCGGGGCGGGTTGTTTTGGGGGCCCCGGTACTGGTTGGggcggccggg                                                                                           >  3:151399/1‑132 (MQ=255)
                                     acgGCAACGAAGTCCGCGTATTCGTCGAGTACCCGATCAGCATCTTCTGATTCTGATCTGATGTATCGAACAGCCCCGGCCTATGCCGGGGCTGTTTCTTTTCGGCGCACGGTTACTGTTTGAGCCGGCGCTGGGc                                                                                        >  5:156581/1‑136 (MQ=255)
                                           aCGAAGTCCGCGTATTCGTCGAGTACCCGATCAGCATCTTCTGATTCTGATCTGATGTATCGAACAGCCCCCGCCGATTGCGGGGGTGTTTGTTTTGGGGGGACGGGTACTGTTTGAGCCGGCGCGGGGCGTACgg                                                                                  >  5:106723/1‑136 (MQ=255)
                                             catgTCCGGTTATTCGTCTAGTCCCCCATCAGCATCTAATGATTATGATCTTATGTAGCGAACAGCCCCGGCCTAGGCGGGGGCTGTTTCTTTTCGGCGAACGGTTACTGTTTGAGCCGGCGCTGGGCTTACGGt                                                                                 <  2:271238/132‑1 (MQ=255)
                                                   cgctTATTCGCCGACTACCCCCTCAGCATATTCTGATTCTGATCTTATGTACCCACCAGCCCCGGCCTCGGCGGGGGCTGTTTCTTTTCGGCGAACGGTTACTGTTTGAGCCGGCGCTGGGCTTACGGTCTCCGGt                                                                          <  7:297199/136‑1 (MQ=255)
                                                   cgcgTATGCGCCCCGTCCCCGATCCCCATCTGGTTATGCTGACCTGTGGCACCCACCAGCCCCGGCCTAGGCGGGGGCTGTTTCTTTTCGGCGAACGGTTACTGTTTGAGCCGGCGCTGGGCTTACGGTCTCCGGt                                                                          <  4:186755/136‑1 (MQ=255)
                                                    gcgTATTCGTCGAGTACCCGATCAGCATCTTCTGATTCTGATCTGATGTATCGAACAGCCCCCGCCAAAGCCGGGGCAGGTTCATTTGCGCGCACCGGTAAAGTTTGAGCCCGCGCTGGGCGTACGGGCTCCCGTg                                                                         >  4:237839/1‑136 (MQ=255)
                                                           cgtcgATTACCCTATCAGCAGCTTCTGATTCTGTTCTGTTGTATCGAACAGCCCCGCCCTTGGGGGGGGCTGTTTCTTTTCGGCGAACGGTTACTGTTTGAGCCGGCGCTGGGCTTACGGTCTCCGGTGCAAAGc                                                                   <  1:109363/135‑1 (MQ=255)
                                                                       gATCAGCATCTTCTGATTCTGAGCTGTTGACCCCACCGGCCCCGCCCTCGGCGGGGGCTGTTTCTTTTCGGCGAACGGTTACTGTTTGAGCCGGCGCTGGGCTTACGGTCTCCGGTGCAAAGCTGTCGCTAcgcg                                                       <  2:108346/135‑1 (MQ=255)
                                                                             cATCTTCTGATTCTGATCTGATGTATCGAACAGCCCCGGCCGATGCCGGGGGGGTTTGTTTTGGGCGCCCGGGTTCTGTTTTGGCCGGCGCGGGGGTTACGGTCTCCGGTGCAAAGCTGTCGCTACGCGCCCTCCg                                                >  4:289107/1‑136 (MQ=255)
                                                                                  tctgattctgatCTGATGTATCGAACAGCCCCGGCCGATTGCGGGGCGGGTTGTTTTCGGGGCACGGGTTGTGTTTTGGCCCGGGCGGGGCGGTCGGGCGCCGGGGCAAAGCTGTGGCTTCGCGCCCTCCGcccc                                            >  6:74462/1‑133 (MQ=255)
                                                                                      attctgatCTGATGTACCGAACCGCCCCGGCCTAGGCGGGGGCTGTTTCTTTTCGGCGAACGGTTACTGTTTGAGCCGGCGCTGGGCTTACGGTCTCCGGTGCAAAGCTGTCGCTACGCGCCATCCGCCACAt                                          <  3:289107/133‑1 (MQ=255)
                                                                                                       cGAACAGCCCCCGCCGATGCCGGGGCGGTTTGTTTTAGGCGAACGGGTTCTGGTTGGGCCGGCGCGGGGGGTACGGTCTCCGGTGCCAAGCTGTCGCTACGCGCCATCCGCCACATACCCGCCTAAAACTCGCgtt                      >  4:226354/1‑136 (MQ=255)
                                                                                                            cgCCCCGGCCTATGCGGGGGCTGTTTCTTTTCGGCGAACGGTTACTGTTTGAGCCGGCGCTGGGCTTACGGTCTCCGGTGCAAAGCTGTCGCTACGCGCCATCCGCCACATCCGCGCATAAAACTCGCTTtcgatc                 <  5:326767/135‑1 (MQ=255)
                                                                                                             gCCCCGGCCTAGGCGGGGGCTGTTTCTTTTCGGCGAACGGTTACTGTTTGAGCCGGCGCTGGGCTTACGGTCTCCGGTGCAAAGCTGTCGCTACGCGCCATCCGCCACATCCGCGCATAAAACTCGCTTtcgatcg                <  3:292766/136‑1 (MQ=255)
                                                                                                             cccccGGCCTATGCGGGGGCTGTTTCTTTTCGGCGAACGGTTACTGTTTGAGCCGGCGCTGGGCTTACGGTCTCCGGTGCAAAGCTGTCGCTACGCGCCATCCGCCACATCCGCGCATAAAACTCGCTTtcgatcg                <  7:310674/135‑1 (MQ=255)
                                                                                                               cggggCTTATGCGGGGGCTGTTTCTTTTCGGCGAACGGTTACTGTTTGAGCCGGCGCTGGGCTTACGGTCTCCGGTGCAAAGCTGTCGCTACGCGCCATCCGCCACATCCGCGCATAAAACTCGCTTTCGATCGAg              <  3:237839/133‑1 (MQ=255)
                                                                                                               cccGGCGTCTGCGGGGGCTGTTTCTTTTCGGCGAACGGTTACTGTTTGAGCCGGCGCTGGGCTTACGGTCTCCGGTGCATAGCTGTCGCTACGCGCCATCCGCCACATCCGCGCATAAAACTCGCTATCGATCGAg              <  3:270240/136‑1 (MQ=255)
                                                                                                                gcGGCTTCTGCGGGGGCTGTTTATTTTCGGCGAACGGTTACTGTTTGAGCCGTCGCTGGGCTTACTGTCTCCGGTGCAAAGCTGTCTCTACGCGCCATCAGCCACATCCGCGCATAAAACTCGCTTTCGATCGAg              <  1:56931/134‑1 (MQ=255)
                                                                                                                           cGGGGCTGTTTCTTTTCGGCGAACGGTTACTGTTTGAGCCGGCGCTGGGCTTACGGTCTCCGGTGCACAGCTGTCGCTACGCGCCATCCGCCACCTCCGCGCATAAAACTCGCTTGCGATCGAGCCCCTGAGTAAc  >  5:225391/1‑136 (MQ=255)
                                                                                                                           cGGGGCTGTTTCTTTTCGGCGAACGGTTACTGTTTGAGCCGGCGCTGGGCTTACGGTCTCCGGTGCAAAGCTGTCGCTACGCGCCATCCGCCACATCCGCGCATAAAACTCGCTTTCGATCGAGCCACTGAGTAAc  >  5:129524/1‑136 (MQ=255)
                                                                                                                           cGGGGCTGTTTCTTTTCGGCGAACGGTTACTGTTTGAGCCGGCGCTGGGCTTACGGTCTCCGGTGCAAAGCTGTCGCTACGCGCCATCCGCCACATCCGCGCATAAAACTCGCTTTCGATCGAGCCACTGAGTAAc  >  3:112356/1‑136 (MQ=255)
                                                                                                                           cGGGGCTGTTTCTTTTCGGCGAACGGTTACTGTTTGAGCCGGCGCTGGGCTTACGGTCTCCGGTGCAAAGCTGTCGCTACGCGCCATCCGCCACATCCGCGCAGAAAACTCGCTTTCGATCGAGCCACTGAGTAAc  >  5:150626/1‑136 (MQ=255)
                                                                                                                           cGGGGCTGTTTCTTTTCGGCGAACGGTTACTGTTTGAGCCGGCGCTGGGCGTACGGTCTCCGGTGCAAAGCTGTCGCTACGCGCCATCCGCCACCTCCGCGCATAAAACTCGCTTTCGATGGAGCCACTGAGTAAc  >  6:4111/1‑136 (MQ=255)
                                                                                                                           cGGGGCTGGTTCTTTTCGGCGAACGGTTACTGTTTGAGCCGGCGCTGGGCTTACGGTCTCCTGTGCAAACCTGTCCCTACGCGCCATCCGCCACATCAGCCCATAAAACTCACTTTCGATCGAGCCACTGAGTAAc  >  5:285000/1‑136 (MQ=255)
                                                                                                                           |                                                                                                                                       
ACCAACAACGCCGTGCAGCAGAACGGTTACAACGACGACGGCAACGAAGTCCGCGTATTCGTCGAGTACCCGATCAGCATCTTCTGATTCTGATCTGATGTATCGAACAGCCCCGGCCTATGCCGGGGCTGTTTCTTTTCGGCGAACGGTTACTGTTTGAGCCGGCGCTGGGCTTACGGTCTCCGGTGCAAAGCTGTCGCTACGCGCCATCCGCCACATCCGCGCATAAAACTCGCTTTCGATCGAGCCACTGAGTAAC  >  ACB122/318997‑319255

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