Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*ACB1432,668,4730AC57.9% ‑0.1 / 11.4 19G234G (GGT→GGGCJ019_02525hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (4/4);  new base C (11/0);  total (15/4)
Fisher's exact test for biased strand distribution p-value = 1.81e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.68e-02
Rejected as consensus: E-value score below prediction cutoff.
Rejected as consensus: Frequency below/above cutoff threshold.

CAGGGCCTTGTCGAACCAGCCGGGTACCACCTTGTCGGTGAAGTGCGGGTACAGCACCAGGTCGTCGCCGCGATAGGGCAGGTCGAGGTGATAGTCGAGCAGGCCGCCGTCGCGGTAGGTGCCGGCGCCGGCACCGGGGATGTCTTTCACGCCCTGCATGACCATCGGTATCGACCCGGAAGCGAGCAGGGCGTGGCGCAGGTTGGTCAAGTCCAGCGGCAGGTAGCGTGAGGGGAAGTCGGTCAGCGTGTCGAGCGGAGGCGCGGCG  >  ACB143/2668341‑2668608
                                                                                                                                    |                                                                                                                                       
cagggcCTTGTCGAACCAGCCGGGTACCACCTTGTCGGTGAAGTGCGGGTACAGCACCAGGTCGTCGCCGCGATAGGGCAGGTCGAGGTGATAGTCGAGCAGGCCGCCGTCGCGGTAGGTGCCGGCGCCGGCAccc                                                                                                                                      >  6:308193/1‑135 (MQ=255)
      cTTGTCGACCCAGCCGGGTACCACCTTGTCGGTGAAGTGCGGGTACAGCACCAGGTCGTCGCCGCGATAGGGCAGGTCGAGGTGATAGTCGAGCAGGCCGCCGTCGCGGTAGGTGCCGGCGCCGGGCCCGGGGGTg                                                                                                                                >  7:36202/1‑136 (MQ=255)
      cTTGTCGAACCAGCCGGGTACCACCTTGTCGGTGAAGTGCGGGTACAGCACCAGGTCGTCGCCGCGATAGGGCAGGTCGAGGTGATAGTCGAGCAGGCCGCCGTCGCGGTAGGTGGCGGCGGCGGCCCCgggggt                                                                                                                                 >  4:218263/1‑133 (MQ=255)
      cTTGTCGAACCAGCCGGGTACCACCTTGTCGGTGAAGTGCGGGTACAGCACCAGGTCGTCGCCGCGATAGGGCAGGTCGAGGTGATAGTCGAGCAGGCCGCCGTCGCGGTAGGTGCCGGCGCCGGCACCGGGGATg                                                                                                                                >  6:267414/1‑136 (MQ=255)
        tGTCGAACCAGCCGGGTACCACCTTGTCGGTGAAGTGCGGGTACAGCACCAGGTCGTCGCCGCGATAGGGCAGGTCGAGGTGATAGTCGAGCAGGCCGCCGTCGCGGTAGGTGCCGGCGCCCGCCCCcggggggt                                                                                                                               >  4:65467/1‑131 (MQ=255)
            gAACCCGCGGGGAACCACCTTGGCGGTTAAGGGCGGGACCTGCACCAGGGCGGCGCCGCGATAGGGCTGGGCGCGGGTATAGTCGAGCAGGCCGCCGGCGCGTTAGGTGCCGGCGCCGGCACCGGGGATGTCTTTc                                                                                                                          <  8:296343/136‑1 (MQ=255)
                cAGCCGGGTACCACCTTGTCGGTGAAGTGCGGGTACAGCACCAGGTCGTCGCCGCGATAGGGCAGGTCGAGGTGATAGTCGAGCAGGCCGCCGTCGCGGTAGGTTGCGGCGCCGGCCCCCGGGCTGTCATTCaacc                                                                                                                      >  6:58772/1‑133 (MQ=255)
                                                        ccAGGTCGTCGCCGCGATAGGGCAGGTCGAGGTGATAGTCGAGCAGGCCGCCGTCGCGGTAGGTGGCCGCGCCCGCCCCGGGGGGGTGGTTGACGCCCAGCCTGACCCTAGGTTTGGGCCCCGGAGCGCGCAGGGc                                                                              >  1:200015/1‑136 (MQ=255)
                                                              cgtcgCCGCGATAGGGCAGGTCGAGGTGATAGTCGAGCAGGCCGCCGTCGCGGTAGGTGCCGGCGCCGGCCCCCGGGATGTCTTTCCCGCCCTGCATGCCCATCGGGATCGACCCGGGAGCGAGCAgggcggggcg                                                                        >  8:69738/1‑136 (MQ=255)
                                                              cgtcgCCGCGATAGGGCAGGTCGAGGTGATAGTCGAGCAGGCCGCCGTCGCGGTAGGTGCCCGCGCCCGCCCCGGGGGTGTGTTTAACCCCCAGCATGACCATCGGTATAGCCCCGGAAGCGCGCCgggcggggcg                                                                        >  1:237895/1‑136 (MQ=255)
                                                              cgtcgCCGCGATAGGGCAGGTCGAGGTGATAGTCGAGCAGGCCGCCGTCGCGGTAGGTGCCCGCGCCCGCCCCGGGGGAGTATTTAACGCCCCGCATGACCAACGGGATCGCCCCCGCAGCGAGCCCGGCGGGgcg                                                                        >  8:48056/1‑136 (MQ=255)
                                                               ggcgccgcGATAGGGCTGGTCGAGGTGATAGTCGAGCAGGCCGCGGGCGCGTTAGGGGCCGGCGCCGGCACCGGGGATGTCTTTCACGCCCTGCATGACCATCGGTATCGACCCGGAAGCGAGCAGGGCGTGgcgc                                                                       <  8:124159/134‑1 (MQ=255)
                                                                   ccgcGATAGGGCAGGTCGAGGTGATAGTCGAGCAGGCCGCCGTCGCGGTAGGTGCCGGCGCCGGCACCGGGGGTGTCTTTGACCCCCGGCATGACCCTCGGGATAGGCCCGGCAGCGCGCAGGGCGGGGGGCacgg                                                                   >  8:87671/1‑133 (MQ=255)
                                                                   ccgcGATAGGGCAGGTCGAGGTGATAGTCGAGCAGGCCGCCGTCGCGGTAGGTGCCCGCGCCCGGCCCCGGGATGTGTTTAAACCCCCGCAATACCAACGGGATAGGCCCCGCAGCGCGCAGGGCGTGGGGGAgg                                                                    >  3:116057/1‑135 (MQ=255)
                                                                                    gAGGTGATAGTCGAGCAGGCCGCCGTCGCGGTAGGTGCCCGCGCCGGCCCCGGGGATGTATTTCAAGCCCAGCAAGACCATCGGTATAGACCCCGAAGCGAGCAGGGCGGGGGGCAGGTTGGTCAAGTCCAGCGGc                                                  >  1:79458/1‑136 (MQ=255)
                                                                                                         gccgTCGCGGTAGGTGCCGGCGCCGGCCCCGGGGGTGTCTTTAACGCCCTGCATGACCATCGGTATCGACCCGGAAGCGCGCAGGGCGGGGCGGCGGGTGGTCAAGTCCCGCGGCAGGGTGCGGGAGGGGAAGTCg                             >  8:38209/1‑136 (MQ=255)
                                                                                                                     ggggcgggcgcgggcACCGGGGATGTCTTTCACGCCCTGCATGACCATCGGTATCGACCCGGAAGCGAGCAGGGCGTGGCGCAGGTTGGTCAAGTCCAGCGGCAGGTAGCGTGAGGGGAAGTCGGTCAGCGTGTc                  <  2:237895/129‑1 (MQ=255)
                                                                                                                          cggccccgccACCGGGGATGTCTTTCACGCCCTGCATGACCATCGGTATCGACCCGGAAGCGAGCAGGGCGTGGCGCAGGTTGGTCAAGTCCAGCGGCAGGTAGCGTGAGGGGAAGTCGGTCAGCGTGTCGAGCgg            <  6:224865/136‑1 (MQ=255)
                                                                                                                             cgccggcACCGGGGATGTCTTTCACGCCCTGCATGACCATCGGTATCGACCCGGAAGCGAGCAGGGCGTGGCGCAGGTTGGTCAAGTCCAGCGGCAGGTAGCGTGAGGGGAAGTCGGTGAGCGTGTCGAGCGGAgg         >  5:157362/1‑136 (MQ=255)
                                                                                                                                    aCCGGGGATGTCTTTCACGCCCTGCATGACCATCGGTATCGACCCGGAAGCGAGCTGGGCGTGGCGCAGGTTGGTCAAGTCCAGCGGCAGGGAGCGTGAGGGGAAGTCGGTCAGCGTGTCGCGCGGAggcgcggcg  >  5:309456/1‑136 (MQ=255)
                                                                                                                                    aCCGGGGATGTCTTTCACGCCCTGCATGACCATCGGTATCGACCCGGAAGCGAGCAGGGCGTGGCGCAGGTTGGTCaa                                                            >  2:86208/1‑78 (MQ=255)
                                                                                                                                    aCCGGGGATGTCTTTCACGCCCTGCATGACCATCGGTATCGACCCGGAAGCGAGCAGGGCGTGGCGCAGGTTGGTCaa                                                            <  1:86208/78‑1 (MQ=255)
                                                                                                                                    aCCGGGGATGTCTTTCACGCCCTGCATGACCATCGGTATCGACCCGGAAGCGAGCAGGGCGTGGCGCAGGTTGGTCAAGTCCAGCGGCAGGTAGCGTGAGGGGAAGTCGGTCAGCGTGTCGAGCGGAggcgcggcg  >  1:207096/1‑136 (MQ=255)
                                                                                                                                    aCCGGGGATGTCTTTCACGCCCTGCATGACCATCGGTATCGACCCGGAAGCGAGCAGGGCGTGGCGCAGGTTGGTCAAGTCCAGCGGCAGGTAGCGTGAGGGGAAGTCGGGCAGCGTGTCGAGCGGAggcgcggcg  >  7:180783/1‑136 (MQ=255)
                                                                                                                                    aCCGGGGATGTCTTTCACGCCCTGCATGACCATCGGTATCGACCCGGAAGCGAGCAGGGCGAGGAGCAGGTTGGTCAAGTCCACCGGCAGGTAGCGGGAGGGGAAGTCGGTGAGCGTGTCCCGCCGcggcgcggcg  >  8:295770/1‑136 (MQ=255)
                                                                                                                                    |                                                                                                                                       
CAGGGCCTTGTCGAACCAGCCGGGTACCACCTTGTCGGTGAAGTGCGGGTACAGCACCAGGTCGTCGCCGCGATAGGGCAGGTCGAGGTGATAGTCGAGCAGGCCGCCGTCGCGGTAGGTGCCGGCGCCGGCACCGGGGATGTCTTTCACGCCCTGCATGACCATCGGTATCGACCCGGAAGCGAGCAGGGCGTGGCGCAGGTTGGTCAAGTCCAGCGGCAGGTAGCGTGAGGGGAAGTCGGTCAGCGTGTCGAGCGGAGGCGCGGCG  >  ACB143/2668341‑2668608

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: