Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*ACB1221,542,1020TG53.3% ‑0.3 / 13.7 17intergenic (+31/‑214)CJ019_01454/groShypothetical protein/10 kDa chaperonin
Reads supporting (aligned to +/- strand):  ref base T (3/4);  new base G (8/0);  total (12/5)
Fisher's exact test for biased strand distribution p-value = 2.56e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.38e-02
Rejected as consensus: E-value score below prediction cutoff.
Rejected as consensus: Frequency below/above cutoff threshold.

GCCAGCGTGCATTCCAGGATGACCCCTTCCAGGCGCAGCCGCGTGATGGCGGCCACCGCCCGAAC‑GTGATCGAAGGCGAGTACGAGCGCCGCGACAAGTAAGAAGGACCGGGGCCGCGCAAGCGG‑CCCC‑GGTTTTTTTTTGGCCTGATGAAAAATTTTCTATCTCGAGCCTTGTAATCGGATGTGGCGGCCTCATGTATGTGGTCACCGCAAGGTTTCTGGTGGCAACACCAGACATTACACAGGTGGTTCGCTCACAGC  >  ACB122/1541971‑1542230
                                                                                                                                      |                                                                                                                                
gcCAGCGTGCATTCCAGGATGACCCCTTCCAGGCGCAGCCGCGTGATGGCGGCCACCGCCCGAAC‑GTGATCGAAGGCGAGTACGAGCGCCGCGACAAGTAAGAAGGACCGGGGCCGCGCCACCGG‑GCCC‑GGggtttt                                                                                                                             >  7:26004/1‑137 (MQ=255)
          aTTCCAGGATGACCCCTTCCAGGCGCAGCCGCGTGATGGCGGCCACCGCCCGAAC‑GTGATCGAAGGCGAGTACGAGCGCCGCGACAAGTAAGAAGGACCGGGGCCGCGCCAGCGGCCCCC‑GGTTTTTTTTTgccggg                                                                                                                    >  2:331911/1‑132 (MQ=255)
          aTTCCAGGATGACCCCTTCCAGGCGCAGCCGCGTGATGGCGGCCACCGCCCGAAC‑GTGATCGAAGGCGAGTACGAGCGCCGCGACAAGTAAGAAGGACCGGGGCCGCGCAAGCGG‑GCGC‑GGGTTTTTTTTGGccggg                                                                                                                   >  6:243580/1‑134 (MQ=255)
            tCCAGGATGACCCCTTCCAGGCGCAGCCGCGTGATGGCGGCCACCGCCCGAAC‑GTGATCGAAGGCGAGTACGAGCGCCGCGACAAGTAAGAAGGACCGGGGCCGCCCCACCCG‑CCCC‑GGGTTTTTTTTGGGCggatt                                                                                                                 >  1:53219/1‑136 (MQ=255)
                      ccccTTCCAGGCGCAGCCGCGTGATGGCGGCCACCGCCCGAAA‑GTGATCGAAGGCGAGTACGAGCGCCGCGACAAGTAAGAAGGACCGGGGCCGCGCAAGCGG‑CCCC‑GGGTGTTTTTTTGCGTGATGAAAAAgttta                                                                                                       >  2:288391/1‑136 (MQ=255)
                         cTTCCAGGCGCAGCCGCGTGATGGCGGCCACCGCCCGAAC‑GTGATCGAAGGCGAGTACGAGCGCCGCGACAAGTAAGAAGGACCGGGGCCGCGCCAGGGG‑CCCCGGGTTTTTTTTTGGCGTGTTTAAAAttttgtat                                                                                                     >  4:35113/1‑133 (MQ=255)
                         cTTCCAGGCGCAGCCGCGTGATGGCGGCCACCGCCCGAAC‑GTGATCGAAGGCGAGTACGAGCGCCGCGACAAGTAAGAAGGACCGGGGCCGCGCAAGCGG‑CCCC‑GGGTTTTTTTTGGGCGGATGAAAAATTTTatat                                                                                                    >  4:39625/1‑137 (MQ=255)
                         cTTCCAGGCGCAGCCGCGTGATGGCGGCCACCGCCCGAAC‑GTGATCGAAGGCGAGTACGAGCGCCGCGACAAGTAAGAAGGACCGGGGCCGCGCAAGCGG‑CCCC‑GGGTTTTTTTTGGGCGGATGAAAAATTTTGTAt                                                                                                    >  7:170373/1‑137 (MQ=255)
                         cTTCCAGGCGCAGCCGCGTGATGGCGGCCACCGCCCGAAC‑GTGATCGAAGGCGAGTACGAGCGCCGCGACAAGTAAGAAGGACCGGGGCCGCGCAAGCGG‑CCCC‑GGGTTTTTTTTGGGCGGATGAAAAATTTTGTAt                                                                                                    >  4:82296/1‑137 (MQ=255)
                            ccAGGCGCAGCCGCGTGATGGCGGCCACCGCCCGAAC‑GTGATCGAAGGCGAGTACGAGCGCCGCGACAAGTAAGAAGGACCGGGGGCCCGCCAGCGGCCCCC‑GGTTTTTTTTTGGCGGGTGTAAAAAtttttttatt                                                                                                  >  3:329453/1‑131 (MQ=255)
                                  gcAGCCGCGTGATGGCGGCCACCGCCCCAAC‑GTGAACGAAGGCGAGAACGAGCGCCGCGACAAGTAAGAAGGACACGGGCCGCCCAACCCG‑CACC‑CGATATTTTTTGGCCATCAGTAAAATTTTGATTATCGAGCCt                                                                                           >  8:244959/1‑137 (MQ=255)
                                                 gggcCCCCCCCCGCACGGTGAGCGATGGCGAGACCGAGCGCCGCCAAAAGAAACAAGCCCCGGGCCCGCCCACGCGG‑CCCC‑GG‑TTTTTTTTGGCCTGATGAAAAATTTTCTATCTCGAGCCTTGTAATCGGATGTgg                                                                            <  1:35644/136‑1 (MQ=255)
                                                       cccccccAAC‑GTCATCCACGGCGCGTACGCGCCCCGCGACAATAAACACGCCCCGGGGCCGCGAACGCGG‑CCCC‑GGTTTTTTTTTGGCCTGATGAAAAATTTTCTATCTCGAGCCTTGTAATCGGATGTGGCGGCCt                                                                      <  5:343609/134‑1 (MQ=255)
                                                             ccacggTGAGCGACGGCTAGTCCGCGCGCCGCAACAACTACCAGCCACGGGGCCCGCGCACGCGG‑CCCC‑GGTTTTTTTTTGGCCTGATGAAAAATTTTCTATCTCGAGCCTTGTAATCGGATGTGGCGGCCTCatgt                                                                 <  8:149650/132‑1 (MQ=255)
                                                                                  aCGAGCGCCGCGACAAGTAAGAAGGACCGGGGCCGCGCAAGCGG‑CCCC‑GGGTTTTTTTTTGGCTGATGTAAAATTTTGTTTGTCGAGCCTTGTTATGGGGTGTGGCGGGCGGATGTTAGTGGTACCCGCAAGGtttt                                            >  7:345011/1‑136 (MQ=255)
                                                                                  aCGAGCGCCGCGACAAGTAAGAAGGACCGGGGCCGCGCAAGCGG‑CCCC‑GGGTTTTTTTTGGGCGGATGTAAAATTTTCTTTGTCGGGCCGTGTAATCGGGTTTTGGGGGCGGATGTTTTTGGTGACCGCAAGGtttt                                            >  5:355314/1‑136 (MQ=255)
                                                                                                 attacgaagGCCGGGGCCCGCGCAAGCGG‑CCCC‑GGTTTTTTTTTGGCCTGATGAAAAATTTTCTATCTCTAGCCTTGTAATCGGATGTGGCGGCCTCATGTATGTGGTCACCGCAAGGTTTCTGGTGGCAACACCAg                             <  2:18838/132‑1 (MQ=255)
                                                                                                                             gccccc‑GGTTTTTTTTTGGCCTGATGAAAAATTTTCTATCTCGAGCCTTGTAATCGGATGTGGCGGCCTCATGTATGTGGTCACCGCAAGGTTTCTGGTGGCAACACCAGACATTACACAGGTGGTTCGCTCACAGc  <  1:331911/135‑1 (MQ=255)
                                                                                                                                      |                                                                                                                                
GCCAGCGTGCATTCCAGGATGACCCCTTCCAGGCGCAGCCGCGTGATGGCGGCCACCGCCCGAAC‑GTGATCGAAGGCGAGTACGAGCGCCGCGACAAGTAAGAAGGACCGGGGCCGCGCAAGCGG‑CCCC‑GGTTTTTTTTTGGCCTGATGAAAAATTTTCTATCTCGAGCCTTGTAATCGGATGTGGCGGCCTCATGTATGTGGTCACCGCAAGGTTTCTGGTGGCAACACCAGACATTACACAGGTGGTTCGCTCACAGC  >  ACB122/1541971‑1542230

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: ‑