Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*ACB1223,491,3370GC56.0% 20.4 / 34.4 25intergenic (‑699/+22)CJ019_03306/CJ019_03307hypothetical protein/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (3/8);  new base C (14/0);  total (17/8)
Fisher's exact test for biased strand distribution p-value = 1.53e-04
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.46e-06
Rejected as consensus: Frequency below/above cutoff threshold.

CTACATCGACCGGTGGTGTCGCCACGTTGTAGGAGCTGGTTTACCCGCGAATAGGCCGGACCTGCCGGCCACTATCCCTGGTCATTCACCGTCACAGCCGCCACCGCATTCGCGGGCAGCCTGCCCGCGAACACCGGCGCAGCCGGTTTCATCCAACTCAACGACCATCGACCTGTACGCAAACGTCTCTGCCATCGATTGCCTTGCATGGCCAATTGTCGTTGATGTAATGATCCGGCTTGAAGCACACCCGTTCATCGTT  >  ACB122/3491202‑3491463
                                                                                                                                       |                                                                                                                              
cTACATCGACCGGTGGTGTCGCCAGGTGGTAGGCGCTGGTTAACCCGCGAATAGGCCGGACCTGCCGGCCAATATCCCTGGGCATTCACCGTCACACCCGCCACCGCATTCGCGGGCAGCCTGCCCGCGAACACCg                                                                                                                                <  8:446404/136‑1 (MQ=21)
 tACATCGACCGGTGGTGTCGCCACGTTGTAGGAGCTGGTTTACCCGCGAATAGGCCGGACCTGCCGGCCACTATCCCTGGTCATTCACCGTCACAGCCGCCACCGCATTCGCCGGCCGCCGGCCCGCGAACACCg                                                                                                                                >  6:81992/1‑135 (MQ=255)
           ggtggtGTCGCCACGTTGTAGGAGCTGGTTTATCCGCGAATAGGCCGGACCTGCCGGCCACTATCCCTGGTCATTCACCGTCACCGCCGCCACCGCATTCCCGGGCAGCCGGCCCGCCCACACCCGCGCAGcccg                                                                                                                      >  6:91523/1‑133 (MQ=21)
           ggtggtGTCGCCACGTTGTAGGAGCTGGTTTACCCGCGAATAGGCCGGACCTGCCGGCCACTATCCCTGGTCATTCACCGTCACAGCCGCCACCGCATTCGCGGGCCGCCGGCCCGCGCACACCCgcgccccccg                                                                                                                      >  2:382716/1‑129 (MQ=21)
           ggtggtGTCGCCACGTTGTAGGAGCTGGTTTACCCGCGAATAGGCCGGACCTGCCGGCCACTATCCCTGGTCATTCACCGTCACAGCCGCCACCGCATTCGCCGGCCGCCCGCCCCCCCACACCCGCgccgcccg                                                                                                                      >  5:423963/1‑133 (MQ=18)
                     ccACGTTGTAGGAGCTGGTTTACCCGCGAATAGGCCGGACCTGCCGGCCACTATCCCTGGTCATTCACCGTCACAGCCGCCACCGCATTCGCGGGCCGCCGGCCCGCGCACCCCCGCGCCCCCGGGTTCATCCaa                                                                                                            >  7:201780/1‑135 (MQ=18)
                                   cgggTTTACCCGCTAATAGGCCGGACCTGCCGGCCACAATCCCTGGTCATTCACCGACACCGCCGCCACCGCATTCGCGGGCAGCCTGCCCGCGAACACCGGCGCAGCCGGTTTCATCCAACTCAACGACCATCGa                                                                                             <  7:323047/134‑1 (MQ=37)
                                          ccccGCTAATTGGCCGCACCTGCCGCCAACTACCCGTGGCCATTCCCCGCCACCGCCGCCCCCCCATGCGGGGGCCGCCGGCCCGCGAACACCGGCGCAGCCGGTTTCATCCAACTCAACGACCATCGACCTGTAc                                                                                      <  7:74246/135‑1 (MQ=255)
                                          aCCCGCGAATAGGCCGGACCTGCCGGCCACTATCCCTGGTCATTCACCGTCACAGCCGCCACCGCATTCCCGGGCACCCTGCCCCCGCACAACCGCGCAGCCGGGTTGAGACAACTCAACGCCCCTCGcccaggcc                                                                                      >  7:121321/1‑131 (MQ=18)
                                              gcgAATAGGCCGGACCTGCCGGCCACTATCCCTGGTCATTCACCGTCACAGCCGCCACCGCATTCGCGGGCCGCCGGCCCGCGAACAACCGCGCCGCCGGGTTCATCCAACTCAACGCCCAACCCCCTGTACGc                                                                                    >  4:186083/1‑134 (MQ=21)
                                                   tAGGCCGGACCTGCCGGCCACTATCCCTGGTCATTCACCGTCACAGCCGCCACCGCATTCCCGGGGACCCGCGCCGCGAACACCCGCCCACCCGCTGTCATCCAACAGCACGACCAAGGACCGGTGCGCAAAccg                                                                              >  4:230603/1‑133 (MQ=12)
                                                                   gCCACTATCCCTGGTCATTCACCGTCACAGCCGCCACCGCATTCGCGGGCCGCCTGCCCGCCAACACCCGCCCCGCCGGGGTCATCCAACTCAACGACCATCGACCAGTAAGCAAAAGTCTCTGCCCTCGaattc                                                              >  6:407556/1‑131 (MQ=25)
                                                                             cTGGTCATTCACCGTCACAGCCGCCACCGCATTCGCCGGCCGCCTGCCCGCGAACACCCGCGCAGCCGGTTTGATCCAACACAACCACCCTCGCCCGGGACGCCAACGTCTCTGCCAACGATTTCCATGCATggg                                                    >  2:82994/1‑134 (MQ=255)
                                                                                   aTTCACCGTCACAGCCGCCACCGCATTCGCGGGCAGCCTGCCCGCGAACACCCGCGCCGCCGGGTTGATCCAACACAACGACCATCGACCTGTACGCAAACGTCTCTGCCATCGAATGCCTTGCATGGGCAAtag                                              >  4:46207/1‑133 (MQ=255)
                                                                                     tCACCGTCACAGCCGCCACCGCATTCGCCGGCACCCTGACCGCGCAAAACCCCGCAGACCGTGTCAGACAACAAAAAGACCACCCACCGGTTCCACAACCTCTCTGCGAAACCTTGCCTTGGAAGGCCCATTGta                                            >  6:151356/1‑134 (MQ=255)
                                                                                       aCCGTCACAGCCGCCACCGCATTCGCGGGCCGCCGGCCCGCGAACACCCGCGCAGCCCGTTTGATACAACTCAACGCCCATCGACCAGTTCGCAAACCTCTCTGCCATCGATTGGCGTGGATGGCCAATTGTCgg                                          >  2:47731/1‑134 (MQ=255)
                                                                                       aCCGTCACAGCCGCCACCGCATTCGCGGGCAGCCTGCCCGCGAACACCCGCGCCGCCGGGTTCATCCAACACAACGCCCCTCGCCCTGTACGCCAACGTCTATGCCATCGATTGGCTTGCATGGCCAATTGTGGt                                          >  4:126171/1‑135 (MQ=255)
                                                                                        ccGTCACAGCCGCCACCGCATTCGCGGGCAGCATGCCCCCGAAAACCGGCGCAGCCCGGTTCATACAACACAACGCCCATAGACCAGTACGCCAACGTAGCTGCCATCGCTTGCCTTGCAGGGACAATTTTCGTTg                                        >  8:96564/1‑136 (MQ=255)
                                                                                                gcccccACGGCATGCGCGGGCAGCCTGCCCGCGAACACCGGCGCAGCCGGTTTCATCCAACTCAACGACCATCGACCTGTACGCAAACGTCTCTGCCATCGATTGCCTTGCATGGCCAATTGTCGTTGATGTAAt                                 <  3:203123/135‑1 (MQ=255)
                                                                                                cccgccCCCGCTTCGGCGGGCCGCCGGCCCGCCAACACCGGCGCAGCCGGTTTCATCCAACTCAACGACCATCGACCTGTACGCAAACGTCTCTGCCATCGATTGCCTTGCATGGCCAATTGTCGTTGATGTAAt                                 <  3:20283/134‑1 (MQ=255)
                                                                                                      aCCGCATTCGCGGGCAGCCTGCCCGCGAACACCGGCGCAGCCGGTTTCATCCAACTCAACGACCATCGACCTGTACGCAAACGTCTCTGCCATCGATTGCCTTGCATGGCCAATTGTCGTTGATGTAATGATCCg                           <  4:72273/135‑1 (MQ=255)
                                                                                                      aCCGCATTCGCGGGCAGCCGGCCCGCGAACACCCGCGCCGCCGGGTTGATCCAACTCAACGCCCCTCGCCCGGTCCGCCCAAGGCTCGGCCCTCGCTTTCCTTGCCTGGGCCATTTTGGTTGATGTTATGATCCg                           >  2:159119/1‑135 (MQ=255)
                                                                                                                     ccccgccccGCAAACACCGGCGCAGCCGGTTTCATCCAACTCAACGACCATCGACCTGTACGCAAACGTCTCTGCCATCGATTGCCTTGCATGGCCAATTGTCGTTGATGTAATGATCCGGCTTGAAGCACAccc            <  5:39350/129‑1 (MQ=255)
                                                                                                                           ccggcgAACACCGGCGCAGCCGGTTTCATCCAACTCAACGACCATCGACCTGTACGCAAACGTCTCTGCCATCGATTGCCTTGCATGGCCAATTGTCGTTGATGTAATGATCCGGCTTGAAGCACACCCGTTCAt      <  1:382716/132‑1 (MQ=255)
                                                                                                                            ccGCGAACACCGGCGCAGCCGGTTTCATCCAACTCAACGACCATCGACCTGTACGCAAACGTCTCTGCCATCGATTGCCTTGCATGGCCAATTGTCGTTGATGTAATGCTCCGGCTTGAAGCCCACCCGTTCATc     >  2:189164/1‑135 (MQ=255)
                                                                                                                            ccGCGAACACCGGCGCAGCCGGTTTCATCCAACTCAACGACCATCGACCTGTACGCAAACGTCTCTGCCATCGATTGCCTTGCATGGCCAATTGTCGTTGATGTAATGATCCGGCTTGAAGCACACCCGTTCa       <  3:126171/133‑1 (MQ=255)
                                                                                                                              gcgAACACCGGCGCAGCCGGTTTCATCCAACTCAACGACCATCGACCTGTACGCAAACGTCTCTGCCATCGATTGCCTTGCATGGCCAATTGTCGTTGATGTAATGATCCGGCTTGAAGCACACCCGTTCATCGtt  >  8:51112/1‑136 (MQ=255)
                                                                                                                                       |                                                                                                                              
CTACATCGACCGGTGGTGTCGCCACGTTGTAGGAGCTGGTTTACCCGCGAATAGGCCGGACCTGCCGGCCACTATCCCTGGTCATTCACCGTCACAGCCGCCACCGCATTCGCGGGCAGCCTGCCCGCGAACACCGGCGCAGCCGGTTTCATCCAACTCAACGACCATCGACCTGTACGCAAACGTCTCTGCCATCGATTGCCTTGCATGGCCAATTGTCGTTGATGTAATGATCCGGCTTGAAGCACACCCGTTCATCGTT  >  ACB122/3491202‑3491463

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