Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*ACB1224,425,0870GC37.9% 36.3 / 15.1 29intergenic (‑29/‑89)CJ019_04200/CJ019_04201hypothetical protein/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (15/3);  new base C (0/11);  total (15/14)
Fisher's exact test for biased strand distribution p-value = 1.05e-05
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

TGCGAGTCGAGCTGCGCAGTCGGCTTGCCGCCACGCGCCGTTTCCATGGTCGCAGAAGCACAACCGGCCAGCATCAGCACCGCAGCCAACGGCAGAATCGCACGCATGTTCCATCCTCGGGGTCGCAAAAAAGAAGGCGACTCTAGCACTGGCAGGCATCCCGAGTGAATTGCCGCGCACTTGCCGCTGTCGCAGATAAAACATGCCATCACGAGGAACCCATCATGAAACTGCACTATCTGCTACTCGCCGCCCTGCCTTG  >  ACB122/4424952‑4425213
                                                                                                                                        |                                                                                                                              
tgCGAGTCGAGCTGCGCAGTCGGCTTGCCGCCACGCGCCGTTTCCATGGTCGCAGAAGCACAACCGGCCAGCATCAGCACCGCAGCCAACGGCAGAATCGCACGCATGTTCCAACCTCGGGGTCGCAAAAAAGAAgg                                                                                                                               >  2:161795/1‑137 (MQ=255)
   gAGTCGAGCTGCGCAGTCGGCTTGCCGCCACGCGCCGTTTCCATGGTCGCAGAAGCACAACCGGCCAGCATCAGCACCGCAGCCAACGGCAGAATCGCACGCATGTTCCATCCTCGGGGTCGCAAAAAAGAAGGCgc                                                                                                                            >  6:161166/1‑136 (MQ=255)
         aGCTGCGCAGTCGGCTTGCCGCCACGCGCCGTTTCCATGGTCGCAGAAGCACAACCGGCCAGCATCAGCACCGCAGCCAACGGCAGAATCGCACGCACGTTCCATCCTCGGGGTCGCAATAAAGAAGGCGACTCTAg                                                                                                                      >  6:3039/1‑137 (MQ=255)
            ggcgcAGTCGGCTTGCCGCCCCGCGCCGTTTCCATGGGCCCAGAACCACAACCGGCCAGCATCAGCACCCCCCCCAACGGCAGAAGCCCACGCATGTACCACCCTCGGGGGCGAAAAAAAGACCGCGACTCTAGCAc                                                                                                                   <  6:157357/136‑1 (MQ=255)
                                              tGGTCGCAGATCCACAACCGGCCACCATCCGCACCGCAGCCACCGGCAGAAGCGCACTCATGTTCCATCCTCGGGGTCCAAAAAACGAACGCGACTCTAGCACTGGCAGGCATCCCGAGTGAATTGCCGCGCACTTg                                                                                 <  7:229573/137‑1 (MQ=255)
                                                  cGCAGACCCACAACCGCCCAGAATCACCCCCTCCTCCACCCGCATAATCGCCCGCATGTTCCAACCTGGGGTCCCAAAAAAACAACGCGACTCTAGCACTGGCAGGCATCCCGAGTGAATTGCCGCGCACTTGCCGc                                                                             <  7:119158/137‑1 (MQ=255)
                                                     aGAAGCACAACCGGCCAGCATCAGCACCGCAGCCAACGGCAGAATCGCACGCATGTTCCATCCTCGGGGTCGCAAAAAAGAAGGCGACTCTAGCACTGGCAGGCATCCCGAGTGAATTGCCGCGCACTTGCCGCTGt                                                                          >  4:303102/1‑137 (MQ=255)
                                                                     aGCATCAGCACCGCAGCCAACGGCAGAATCGCACGCATGTTCCATCCTCGGGGTCGCAAAAAAGAAGGCGACTCTAGCACTGGCAGGCATCCCGAGTGAATTGCCGCGCACTTGCCGCTGTCGCAGATAAAACATGc                                                          >  6:317509/1‑137 (MQ=255)
                                                                          cAGCACCGCAGCCAACGGCAGAATCGCACGCATGTTCCATCCTCGGGGTCGCAAAAAAGAAGGCGACTCTAGCACTGGCAGGCATCCCGAGTGAATTGCCGCGCACTTGCCGCTGTCGCAGATAAAACATGCCATcc                                                     >  3:245840/1‑136 (MQ=255)
                                                                              cccGCACCCACCGACAGAAGCCCACGCATGTTCCAGCCTCGGGGCCCAAAAAACCACCGCGACTCTAGCACTGGCAGGCATCCCGAGTGAATTGCCGCGCACTTGCCGCTGTCGCAGATAAAACATGCCATCACGAg                                                 <  7:187134/136‑1 (MQ=255)
                                                                                ccaacccACCGGCATAATCGCACGCATGTACCATCCTCTGGGTCAAAAAAACAACCGCGACTCTAGCACTGGCAGGCATCCCGAGTGAATTGCCGCGCACTTGCCGCTGTCGCAGATAAAACATGCCATCACGAGGa                                               <  7:134804/132‑1 (MQ=255)
                                                                                     cccaCGGCAGAATCGCCCGCATGTGCCATCCTCGGGGTCCAAAAAACCAACGCGACTCTAGCACTGGCAGGCATCCCGAGTGAATTGCCGCGCACTTGCCGCTGTCGCAGATAAAACATGCCATCACGAGGAACCca                                          <  6:65320/134‑1 (MQ=255)
                                                                                     ccAACGGCATAATCGAACGCATTTTCAATCATGGGGGTCCAAAAAAACAACGCGACTCTAGCACTGGCAGGCATCCCGAGTGAATTGCCGCGCACTTGCCGCTGTCGCAGATAAAACATGCCATCACGAGGAACCca                                          <  7:321334/137‑1 (MQ=255)
                                                                                     ccAACGGCAGAATCGCACGCATGTTCCATCCTCGGGGTCTCAAAAAATAAGGCGACTCTAGCACTGGCAGGCATCCCGAGTGAATTGCCGCGCACTTGCCGCTGTCGCAGATAAAACATGCCATCACGAGGAACCca                                          <  6:117747/137‑1 (MQ=255)
                                                                                     ccAACGGCAGAATCGCACGCATGTTCCATCCTCGGGGTCGCAAAAAAGAAGGCGACTCTAGCACTGGCAGGCATCCCGAGTGAATTGCCGCGCACTTGCCGCTGTCGCAGATAAAACATGCCATCACGAGGAACCca                                          >  1:98138/1‑137 (MQ=255)
                                                                                     ccAACGGCAGAATCGCACGCATGTTCCATCCTCGGGGTCGCAAAAAAGAAGGCGACTCTAGCACTGGCAGGCATCCCGAGTGAATTGCCGCGCACTTGCCGCTGTCGCAGATAAAACATGCCATCACGAGGAACCca                                          >  7:283681/1‑137 (MQ=255)
                                                                                     ccAACGGCAGAATCGCACGCATGTTCCATCCTCGGGGTCGCAAAAAAGAAGGCGACTCTAGCACTGGCAGGCATCCCGAGTGAATTGCCGCGCACTTGCCGCTGTCGCAGATAAAACATGCCATCACGAGGAACCca                                          >  8:234980/1‑137 (MQ=255)
                                                                                          ggCAGAAGCGCACGCATGTACCAGCCTCGGGGTCGCAAAAAAGACCGCGACTCTAGCACTGGCAGGCATCCCGAGTGAATTGCCGCGCACTTGCCGCTGTCGCAGATAAAACATGCCATCACGAGGAACCCATCATg                                     <  5:161166/137‑1 (MQ=255)
                                                                                               aaTCGCACGCGTGTGCCATCGGGGGGGTCCAAAAAACGACGGCGACTCTAGCACTGGCAGGCATCCCGAGTGAATTGCCGCGCACTTGCCGCTGTCGCAGATAAAACATGCCATCACGAGGAACCCATCATGAAACt                                <  8:283681/137‑1 (MQ=255)
                                                                                                         aTGTTCCATCCTCGGGGTCGCAAAAAAGAAGGCGACTCTAGCACTGGCAGGCATCCCGAGTGAATTGCCGCGCACTTGCCGCTGTCGCAGATAAAACATGCCATCACGAGGCACCCATCATGAAACTGCACTATCTg                      >  7:116026/1‑137 (MQ=255)
                                                                                                                aTCCTCGGGGTCGCAAAAAAGAAGGCGACTCTAGCACTGGCAGGCATCCCGAGTGAATTGCCGCGCACTTGCCGCTGTCGCAGATAAAACATGCCATCACGAGGAACCCATCATGAAACTGCACTATCTGCTACTcg               >  6:332485/1‑137 (MQ=255)
                                                                                                                  ccTCGGGGTCCAAAAAACCACCGCGACTCTAGCACTGGCAGGCATCCCGAGTGAATTGCCGCGCACTTGCCGCTGTCGCAGATAAAACATGCCATCACGAGGAACCCATCATGAAACTGCACTATCTGCTACTcgcc             <  6:272618/137‑1 (MQ=255)
                                                                                                                  ccTCGGGGCCGAAAAAAACAACGCGACTCTAGCACTGGCAGGCATCCCGAGTGAATTGCCGCGCACTTGCCGCTGTCGCAGATAAAACATGCCATCACGAGGAACCCATCATGAAACTGCACTATCTGCTACTcgcc             <  1:28294/137‑1 (MQ=255)
                                                                                                                    tCGGGGTCGCAAAAAAGAAGGCGACTCTAGCACTGGCAGGCATCCCGAGTGAATTGCCGCGCACTTGCCGCTGTCGCAGATAAAACATGCCATCACGAGGAACCCATCATGAAACTGCACTATCTGCTACTcgccgc           >  1:126988/1‑137 (MQ=255)
                                                                                                                    tCGGGGTCGCAAAAAAGAAGGCGACTCTAGCACTGGCAGGCATCCCGAGTGAATTGCCGCGCACTTGCCGCTGTCGCAGATAAAACATGCCATCACGAGGAACCCATCATGAAACTGCACTATCTGCTACTcgccgc           >  8:73805/1‑137 (MQ=255)
                                                                                                                    tCGGGGTCGCAAAAAAGAAGGCGACTCTAGCACTGGCAGGCATCCCGAGTGAATTGCCGCGCACTTGCCGCTGTCGCAGATAAAACATGCCATCACGAGGAACCCAGCATGAAACTGCACTATCTGCTACTcgccgc           >  5:241264/1‑137 (MQ=255)
                                                                                                                      ggggTCGCAAAAAAGAAGGCGACTCTAGCACTGGCAGGCATCCCGAGTGAATTGCCGCGCACTTGCCGCTGTCGCAGATAAAACATGCCATCACGAGGAACCCATCATGAAACTGCACTATCTGCTACTCGCCGccc         >  6:293313/1‑137 (MQ=255)
                                                                                                                         gTCGCAAAAAAGATGGCGACTCTAGCACTGGCAGGCATCCCGAGTGAATTGCCGCGCACTTGCCGCTGTCGCAGATAAAACATGCCATCACGAGGAACCCATCATGAAACTGCACTATCTGCTACTCGCCGCACtgc      >  7:32905/1‑137 (MQ=255)
                                                                                                                             cAAAAAAGAACGCGACTCTAGCACTGGCAGGCATCCCGAGTGAATTGCCGCGCACTTGCCGCTGTCGCAGATAAAACATGCCATCACGAGGAACCCATCATGAAACTGCACTATCTGCTACTCGCCGCCCtgccttg  <  3:303102/137‑1 (MQ=255)
                                                                                                                              aaaaaaaGAAGGCGACTCTAGCACTGGCAGGCATCCCGAGTGAATTGCCGCGCACTTGCCGCTGTCGCAGATAAAACATGCCATCACGAGGAACCCATCATGAAACTGCACTATCTGCTACTCGCCGCCCtgccttg  <  7:73805/136‑1 (MQ=255)
                                                                                                                                        |                                                                                                                              
TGCGAGTCGAGCTGCGCAGTCGGCTTGCCGCCACGCGCCGTTTCCATGGTCGCAGAAGCACAACCGGCCAGCATCAGCACCGCAGCCAACGGCAGAATCGCACGCATGTTCCATCCTCGGGGTCGCAAAAAAGAAGGCGACTCTAGCACTGGCAGGCATCCCGAGTGAATTGCCGCGCACTTGCCGCTGTCGCAGATAAAACATGCCATCACGAGGAACCCATCATGAAACTGCACTATCTGCTACTCGCCGCCCTGCCTTG  >  ACB122/4424952‑4425213

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