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.

GCCAGCGTGCATTCCAGGATGACCCCTTCCAGGCGCAGCCGCGTGATGGCGGCCACCGCCCGAACGTGATCGAAGGCGAGTACGAGCGCCGCGACAAGTAAGAAGGACCGGGGCCGCGCAAGCGGCCCCGGTTTTTTTTTGGCCTGATGAAAAATTTTCTATCTCGAGCCTTGTAATCGGATGTGGCGGCCTCATGTATGTGGTCACCGCAAGGTTTCTGGTGGCAACACCAGACATTACACAGGTGGTTCGCTCACAGC  >  ACB122/1541971‑1542230
                                                                                                                                      |                                                                                                                                
gcCAGCGTGCATTCCAGGATGACCCCTTCCAGGCGCAGCCGCGTGATGGCGGCCACCGCCCGAACGTGATCGAAGGCGAGTACGAGCGCCGCGACAAGTAAGAAGGACCGGGGCCGCGCCACCGGGCCCGGggtttt                                                                                                                             >  7:26004/1‑137 (MQ=255)
          aTTCCAGGATGACCCCTTCCAGGCGCAGCCGCGTGATGGCGGCCACCGCCCGAACGTGATCGAAGGCGAGTACGAGCGCCGCGACAAGTAAGAAGGACCGGGGCCGCGCCAGCGGCCCCCGGTTTTTTTTTgccggg                                                                                                                    >  2:331911/1‑132 (MQ=255)
          aTTCCAGGATGACCCCTTCCAGGCGCAGCCGCGTGATGGCGGCCACCGCCCGAACGTGATCGAAGGCGAGTACGAGCGCCGCGACAAGTAAGAAGGACCGGGGCCGCGCAAGCGGGCGCGGGTTTTTTTTGGccggg                                                                                                                   >  6:243580/1‑134 (MQ=255)
            tCCAGGATGACCCCTTCCAGGCGCAGCCGCGTGATGGCGGCCACCGCCCGAACGTGATCGAAGGCGAGTACGAGCGCCGCGACAAGTAAGAAGGACCGGGGCCGCCCCACCCGCCCCGGGTTTTTTTTGGGCggatt                                                                                                                 >  1:53219/1‑136 (MQ=255)
                      ccccTTCCAGGCGCAGCCGCGTGATGGCGGCCACCGCCCGAAAGTGATCGAAGGCGAGTACGAGCGCCGCGACAAGTAAGAAGGACCGGGGCCGCGCAAGCGGCCCCGGGTGTTTTTTTGCGTGATGAAAAAgttta                                                                                                       >  2:288391/1‑136 (MQ=255)
                         cTTCCAGGCGCAGCCGCGTGATGGCGGCCACCGCCCGAACGTGATCGAAGGCGAGTACGAGCGCCGCGACAAGTAAGAAGGACCGGGGCCGCGCCAGGGGCCCCGGGTTTTTTTTTGGCGTGTTTAAAAttttgtat                                                                                                     >  4:35113/1‑133 (MQ=255)
                         cTTCCAGGCGCAGCCGCGTGATGGCGGCCACCGCCCGAACGTGATCGAAGGCGAGTACGAGCGCCGCGACAAGTAAGAAGGACCGGGGCCGCGCAAGCGGCCCCGGGTTTTTTTTGGGCGGATGAAAAATTTTatat                                                                                                    >  4:39625/1‑137 (MQ=255)
                         cTTCCAGGCGCAGCCGCGTGATGGCGGCCACCGCCCGAACGTGATCGAAGGCGAGTACGAGCGCCGCGACAAGTAAGAAGGACCGGGGCCGCGCAAGCGGCCCCGGGTTTTTTTTGGGCGGATGAAAAATTTTGTAt                                                                                                    >  7:170373/1‑137 (MQ=255)
                         cTTCCAGGCGCAGCCGCGTGATGGCGGCCACCGCCCGAACGTGATCGAAGGCGAGTACGAGCGCCGCGACAAGTAAGAAGGACCGGGGCCGCGCAAGCGGCCCCGGGTTTTTTTTGGGCGGATGAAAAATTTTGTAt                                                                                                    >  4:82296/1‑137 (MQ=255)
                            ccAGGCGCAGCCGCGTGATGGCGGCCACCGCCCGAACGTGATCGAAGGCGAGTACGAGCGCCGCGACAAGTAAGAAGGACCGGGGGCCCGCCAGCGGCCCCCGGTTTTTTTTTGGCGGGTGTAAAAAtttttttatt                                                                                                  >  3:329453/1‑131 (MQ=255)
                                  gcAGCCGCGTGATGGCGGCCACCGCCCCAACGTGAACGAAGGCGAGAACGAGCGCCGCGACAAGTAAGAAGGACACGGGCCGCCCAACCCGCACCCGATATTTTTTGGCCATCAGTAAAATTTTGATTATCGAGCCt                                                                                           >  8:244959/1‑137 (MQ=255)
                                                 gggcCCCCCCCCGCACGGTGAGCGATGGCGAGACCGAGCGCCGCCAAAAGAAACAAGCCCCGGGCCCGCCCACGCGGCCCCGGTTTTTTTTGGCCTGATGAAAAATTTTCTATCTCGAGCCTTGTAATCGGATGTgg                                                                            <  1:35644/136‑1 (MQ=255)
                                                       cccccccAACGTCATCCACGGCGCGTACGCGCCCCGCGACAATAAACACGCCCCGGGGCCGCGAACGCGGCCCCGGTTTTTTTTTGGCCTGATGAAAAATTTTCTATCTCGAGCCTTGTAATCGGATGTGGCGGCCt                                                                      <  5:343609/134‑1 (MQ=255)
                                                             ccacggTGAGCGACGGCTAGTCCGCGCGCCGCAACAACTACCAGCCACGGGGCCCGCGCACGCGGCCCCGGTTTTTTTTTGGCCTGATGAAAAATTTTCTATCTCGAGCCTTGTAATCGGATGTGGCGGCCTCatgt                                                                 <  8:149650/132‑1 (MQ=255)
                                                                                  aCGAGCGCCGCGACAAGTAAGAAGGACCGGGGCCGCGCAAGCGGCCCCGGGTTTTTTTTTGGCTGATGTAAAATTTTGTTTGTCGAGCCTTGTTATGGGGTGTGGCGGGCGGATGTTAGTGGTACCCGCAAGGtttt                                            >  7:345011/1‑136 (MQ=255)
                                                                                  aCGAGCGCCGCGACAAGTAAGAAGGACCGGGGCCGCGCAAGCGGCCCCGGGTTTTTTTTGGGCGGATGTAAAATTTTCTTTGTCGGGCCGTGTAATCGGGTTTTGGGGGCGGATGTTTTTGGTGACCGCAAGGtttt                                            >  5:355314/1‑136 (MQ=255)
                                                                                                 attacgaagGCCGGGGCCCGCGCAAGCGGCCCCGGTTTTTTTTTGGCCTGATGAAAAATTTTCTATCTCTAGCCTTGTAATCGGATGTGGCGGCCTCATGTATGTGGTCACCGCAAGGTTTCTGGTGGCAACACCAg                             <  2:18838/132‑1 (MQ=255)
                                                                                                                             gcccccGGTTTTTTTTTGGCCTGATGAAAAATTTTCTATCTCGAGCCTTGTAATCGGATGTGGCGGCCTCATGTATGTGGTCACCGCAAGGTTTCTGGTGGCAACACCAGACATTACACAGGTGGTTCGCTCACAGc  <  1:331911/135‑1 (MQ=255)
                                                                                                                                      |                                                                                                                                
GCCAGCGTGCATTCCAGGATGACCCCTTCCAGGCGCAGCCGCGTGATGGCGGCCACCGCCCGAACGTGATCGAAGGCGAGTACGAGCGCCGCGACAAGTAAGAAGGACCGGGGCCGCGCAAGCGGCCCCGGTTTTTTTTTGGCCTGATGAAAAATTTTCTATCTCGAGCCTTGTAATCGGATGTGGCGGCCTCATGTATGTGGTCACCGCAAGGTTTCTGGTGGCAACACCAGACATTACACAGGTGGTTCGCTCACAGC  >  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: