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

TGCGAGTCGAGCTGCGCAGTCGGCTTGCCGCCACGCGCCGTTTCCATGGTCGCAGAAGCACAACCGGCCAGCATCAGCACCGCAGCCAACGGCAGAATCGCACGCATGTTCCATCCTCGGGGTCGCAAAAAAGAAG‑GCGACTCTAGCACTGGCAGGCATCCCGAGTGAATTGCCGCGCACTTGCCGCTGTCGCAGATAAAACATGCCATCACGAGGAACCCATCATGAAACTGCACTATCTGCTACTCGCCGCCCTGCCTTGCCTGCCC  >  ACB122/4424952‑4425220
                                                                                                                                       |                                                                                                                                      
tgCGAGTCGAGCTGCGCAGTCGGCTTGCCGCCACGCGCCGTTTCCATGGTCGCAGAAGCACAACCGGCCAGCATCAGCACCGCAGCCAACGGCAGAATCGCACGCATGTTCCATCCTCGGGGTCGCAAAAAAGAAg‑g                                                                                                                                      >  2:323188/1‑137 (MQ=255)
                 aGTCGGCTTGCCGCCACGCGCCGTTTCCATGGTCGCAGAAGCACAACCGGCCAGCATCAGCACCGCAGCCAACGGCAGAATCGCACGCATGTTCAATCCTCGGGGTCGCAAAAAAGAAG‑GCGACTCTAGCACTggca                                                                                                                     >  5:333423/1‑137 (MQ=255)
                      gCTTGCCGCCACGCGCCGTTTCCATGGTCGCAGAAGCACAACCGGCCAGCATCAGCACCGCAGCCAACGGCAGAATCGCACGCATGTTCCATCCTCGGGGTCGCAAAAAAGAAG‑GCGACTCTAGCACTGGCAGGCAt                                                                                                                >  3:27390/1‑137 (MQ=255)
                      gCTGGCCCCCCCGCGCGTTTGCCCTGGGCGCAGACCCCCACCCGGCCCCCAGCACCCCCCCACCCCACGGATCAATCGCCCGCTTGTTCCATCCGCGGGTACCAAAAAACCACC‑GCGACTCTAGCACTGGCAGGCAt                                                                                                                <  2:342620/137‑1 (MQ=255)
                        ttGCCGCCACGCGCCGTTTCCATGGTCGCAGAAGCACAACCGGCCAGCATCAGCACCGCAGCCAACGGCAGAATCGCACGCATGTTCCATCCTCGGGGTCGCAAAAAAGAAG‑GCGACTCTAGCACTGGCAGGCATcc                                                                                                              >  8:285677/1‑137 (MQ=255)
                             gccACGCGCGGTTTCCTTGGTCGCAGAAGCACACCCGGCCACCATCAGCCCCCCCGACACGGCCAGAATCGCACGCATTTTCCATCCTCGGGGTCCAAAAAAAGACG‑GCGACTCTAGCACTGGCAGGCATCCCGAGt                                                                                                         <  8:32042/137‑1 (MQ=255)
                             cccACGCGCGGTTTCCATGGGCGCAGAACCACAACCGGCCAGCATCACCACCGCAGCCACCGGCAGAATCCCACGCATGTTCCATCCTGGGGGTCCAAAAAAAGAAG‑GCGACTCTAGCACTGGCAGGCATCCCGAGt                                                                                                         <  7:225120/136‑1 (MQ=255)
                                    gcCGTTTCCATGGTCGCAGAAGCACAACCGGCCAGCATCAGCACCGCAGCCAACGGCAGAATCGCACGCATGTTCCATCCTCGGGGTCGCAAAAAAGAAG‑GCGACTCTAGCACTGGAAGGCATCCCGAGTGAATTGc                                                                                                  >  5:166973/1‑137 (MQ=255)
                                           ccATGGTCGCAGAAGCACAACCGGCCAGCATCAGCACCGCAGCCAACGGCAGAATCGCACGCATGTTCCATCCTCGGGGTCGCAAAAAAGAAG‑GCGACTCTAGCACTGGCAGGCATCCCGAGTGAATTGCCGCGCAc                                                                                           >  7:32042/1‑137 (MQ=255)
                                           ccATGGTCGCAGAAGCACAACCGGCCAGCATCAGCACCGCAGCCAACGGCAGAATCGCACGCATGTTCCATCCTCGGGGTCGCAAAAAAGAAG‑GCGACTCTAGCACTGGCAGGCATCCCGAGTGAATTGCCGCGCAc                                                                                           >  6:393669/1‑137 (MQ=255)
                                                   ccAGAACCACAACCGGCCCCCACCACCCCCGCAGCCACGGGCAGACTCCCCCGCATGGTCCAGCCGCGGGGTCGCAAAAACGACC‑GCGACTCTAGCACTGGCAGGCATCCCGAGTGAATTGCCGCGCACTTGCCGCt                                                                                   <  7:113183/136‑1 (MQ=255)
                                                                    cAGCATCAGCACCGCAGCCAACGGCAGAATCGCACGCATGTTCCATCCTCGGGGTCGCAAAAAAGAAG‑GCGACTCTAGCACTGGCAGGCATCCCGAGTGAATTGCCGCGCACTTGCCGCTGTCGCAGATAAAACATg                                                                  >  6:288853/1‑137 (MQ=255)
                                                                     aGCATCAGCACCGCAGCCAACGGCAGAATCGCACGCATGTTCCATCCTCGGGGTCGCAAAAAAGAAG‑GCGACTCTAGCACTGGCAGGCATCCCGAGTGAATTGCCGCGCACTTGCCGCTGTCGCAGATAAAACATGc                                                                 >  1:265739/1‑137 (MQ=255)
                                                                            ccaccccacccAACGCAAAAATCCCTCGAATTTACCAACCTCGTGGAC‑CAAAAACAAAGCGCGACTCTAGCACTGGCAGGCATCCCGAGTGAATTGCCGCGCACTTGCCGCTGTCGCAGATAAAACATGCCATCACg                                                          <  3:85464/136‑1 (MQ=255)
                                                                              aCCGCAGCCAACGGCAGAATCGCACGCATGTTCCATCCTCGGGGTCGCAAAAAAGAAG‑GCGACTCTAGCACTGGCAGGCATCCCGAGTGAATTGCCGCGCACTTGCCGCTGTCGCAGATAAAACATGCCATCACGAg                                                        >  5:352027/1‑137 (MQ=255)
                                                                                cGCAGCCAACGGCATAATCGCACGCATGTCCCAGCCTGGGGGTCCCAAAAAAGACC‑GCGACTCTAGCACTGGCAGGCATCCCGAGTGAATTGCCGCGCACTTGCCGCTGTCGCAGATAAAACATGCCATCACGAGGa                                                      <  8:58786/137‑1 (MQ=255)
                                                                                cGCACCCAACGCCAGAATCGCACGCATGTTCCAGCCTCGGGGTCCAAAAAAAGACC‑GCGACTCTAGCACTGGCAGGCATCCCGAGTGAATTGCCGCGCACTTGCCGCTGTCGCAGATAAAACATGCCATCACGAGGa                                                      <  7:291686/137‑1 (MQ=255)
                                                                                    gCCAACGGCAGAATCGCACGCATGTTCCATCCTCGGGGTCGCAAAAAAGAAG‑GCGACTCTAGCACTGGCAGGCATCCCGAGTGAATTGCCGCGCACTTGCCCCTGTCGCAGATAAAACATGCCATCACGAGGAAccc                                                  >  1:59215/1‑137 (MQ=255)
                                                                                     ccAACGGCAGAATCGCACGCATGTTCCATCCTCGGGGTCGCAAAAAAGAAG‑GCGACTCTAGCACTGGCAGGCATCCCGAGTGAATTGCCGCGCACTTGCCGCTGTCGCAGATAAAACATGCCATCACGCGGAACCca                                                 >  8:376298/1‑137 (MQ=255)
                                                                                     ccAACGGCAGAATCGCACGCATGTTCCATCCTCGGGGTCGCAAAAAAGAAG‑GCGACTCTAGCACTGGCAGGCATCCCGAGAGAATTGCCGCGCACTTGCCGCTGTCGCAGATAAAACATGCAATCACGAGGAACCca                                                 >  4:153434/1‑137 (MQ=255)
                                                                                               aaccgcccgaaTTTTCCACCCGCGGGGCCACAAAAAAGAAC‑GCGACTCTAGCACTGGCAGGCATCCCGAGTGAATTGCCGCGCACTTGCCGCTGTCGCAGATAAAACATGCCATCACGAGGAACCCATCATGAAACt                                       <  5:58812/134‑1 (MQ=255)
                                                                                               aaTCGCCGGCATGTTCCAGCCTCGGGGTAGCAAAAAAGAAC‑GCGACTCTAGCACTGGCAGGCATCCCGAGTGAATTGCCGCGCACTTGCCGCTGTCGCAGATAAAACATGCCATCACGAGGAACCCATCATGAAACt                                       <  7:363363/137‑1 (MQ=255)
                                                                                                      cgcaTGTACCACCCGCGGGGACCAAAAAACCACC‑GCGACTCTAGCACTGGCAGGCATCCCGAGTGAATTGCCGCGCACTTGCCGCTGTCGCAGATAAAACATGCCATCACGAGGAACCCATCATGAAACTGCACTAt                                <  5:238807/137‑1 (MQ=255)
                                                                                                        caTGTTCCATCCTCGGGGCCTCAAAAAAGAAC‑GCGACTCTAGCACTGGCAGGCATCCCGAGTGAATTGCCGCGCACTTGCCGCTGTCGCAGATAAAACATGCCATCACGAGGAACCCATCATGAAACTGCACTATCt                              <  6:97847/137‑1 (MQ=255)
                                                                                                         aTGTTCCATCCTCGGGGTCGCAAAAAAGAAG‑GCGACTCTAGCACTGGCAGTCATCCCGAGTGAATTGCCGCGCACTTGCCGCTGTCGAAGATAAAACATGCCATAACGAGGAACACATCATGAAACTGCACTATCTg                             >  7:48462/1‑137 (MQ=255)
                                                                                                            ttCCATCCTCGGGGTCGCAAAAAAGAAG‑GCGACTCTAGCACTGGCAGGCATCCCGAGTGAATTGCCGCGCACTTGCCGCTGTCGCAGATAAAACATGCCATCACGAGGAACCCATCATGAAACTGCACTATCTGCTa                          >  2:189962/1‑137 (MQ=255)
                                                                                                                    tCGGGGTCGCAAAAAAGAAG‑GCGACTCTAGCACTGGCAGGCATCCCGAGTGAATTGCCGCGCACTTGCCGCTGTCGCAGATAAAACATGCCATCACGAGGAACCCATCATGCAACTGCACTATCTGCTACTcgccgc                  >  2:135297/1‑137 (MQ=255)
                                                                                                                    tCGGGGTCGCAAAAAAGAAG‑GCGACTCTAGCACTGGCAGGCATCCCGAGTGAATTGCCGCGCACTTGCCGCTGTCGCAGATAAAACATGCCATCACGAGGAACCCATCATGAAACTGCACTATCTGCTACTcgccgc                  >  6:90707/1‑137 (MQ=255)
                                                                                                                    tCGGGGTCGCAAAAAAGAAG‑GCGACTCTAGCACTGGCAGGCATCCCGAGTGAATTGCCGCGCACTTGCCGCTGTCGCAGATAAAACATGCCATCACGAGGAACCCATCATGAAACTGCACTATCTGCTACTcgccgc                  >  4:167002/1‑137 (MQ=255)
                                                                                                                    tCGGGGTCGCAAAAAAGAAG‑GCGACTCTAGCACTGGCAGGCATCCCGAGTGAATTGCCGCGCACTTGCCGCTGTCGCAGATAAAACATGCCATCACGAGGAACCCATCATGAAACTGCACTATCTGCTACTcgccgc                  >  7:39696/1‑137 (MQ=255)
                                                                                                                      ggggTCGCAAAAAAGAAG‑GCGACTCTAGCACTGGCAGGCATCCCGAGTGAATTGCCGCGCACTTGCCGCTGTCGCAGATAAAACATGCCATCACGAGGAACCCATCATGAAACTGCACTATCTGCTACTCGCCGccc                >  6:343074/1‑137 (MQ=255)
                                                                                                                             aaaaaaaCCAC‑GCGACTCTAGCACTGGCAGGCATCCCGAGTGAATTGCCGCGCACTTGCCGCTGTCGCAGATAAAACATGCCATCACGAGGAACCCATCATGAAACTGCACTATCTGCTACTCGCCGCCCtgccttg         <  1:338173/136‑1 (MQ=255)
                                                                                                                                    gaag‑GCGACTCTAGCACTGGCAGGCATCCCGAGTGAATTGCCGCGCACTTGCCGCTGTCGCAGATAAAACATGCCATCACGAGGAACCCATCATGAAACTGCACTATCTGCTACTCGCCGCCCTGCCTTGCCTGccc  <  6:333423/137‑1 (MQ=255)
                                                                                                                                       |                                                                                                                                      
TGCGAGTCGAGCTGCGCAGTCGGCTTGCCGCCACGCGCCGTTTCCATGGTCGCAGAAGCACAACCGGCCAGCATCAGCACCGCAGCCAACGGCAGAATCGCACGCATGTTCCATCCTCGGGGTCGCAAAAAAGAAG‑GCGACTCTAGCACTGGCAGGCATCCCGAGTGAATTGCCGCGCACTTGCCGCTGTCGCAGATAAAACATGCCATCACGAGGAACCCATCATGAAACTGCACTATCTGCTACTCGCCGCCCTGCCTTGCCTGCCC  >  ACB122/4424952‑4425220

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