Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*ACB1221,721,8780AC42.9% 11.7 / 14.2 35G95G (GGT→GGG) CJ019_01617hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (19/1);  new base C (5/10);  total (24/11)
Fisher's exact test for biased strand distribution p-value = 1.47e-04
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.27e-01
Rejected as consensus: Frequency below/above cutoff threshold.

GGGTAAACCCATTCCTACAGAAACGAGCAAGCCTGCCGTTTCATGATGTACGCCAGGTGATTTCCTCTTCACCGTCGGCGCTGATGCGGACCCAGCGGTCAGCATCTTCTTCCCCGTCTTCCTCGACCCAGCCACCCGG‑CGCGCAGCGCACTTCCACACCCAGCGCGGCAAATGCGGCGCGGGCGCAGGCAACGTCATCGGCCCACGGCGTCTGGTCGCTTTCCAGATACAGGCTGTTCCATTTCCCTACAGCCTTGGGT  >  ACB122/1721745‑1722004
                                                                                                                                     |                                                                                                                               
gggTAAACCCATTCCTACAGAAACGAGCAAGCCTGCCGTTTCATGATGTACGCCAGGTGATTTCCTCTTCACCGTCGGCGCTGATGCGGACCCAGCGGTCAGCATCTTCTTCCCCGTCTTCCTCGACCCAGCCAc                                                                                                                                <  6:3988/135‑1 (MQ=255)
      aCCCATTCCTACAGAAACGAGCAAGCCTGCCGTTTCATGATGTACGCCAGGTGATTTCCTCTTCACCGTCGGCGCTGATGCGGCCCCAGCGGTCAGCCTCTTCTTCCCCGTCTTCCTCGCCCCAGCCCCCCGg‑cg                                                                                                                         >  3:142683/1‑135 (MQ=255)
               tACAGAAACGAGCAAGCCTGCCGTTTCATGATGTACGCCAGGTGATTTCCTCTTCACCGTCGGCGCTGATGCGGACCCAGCGGTCAGCATCTTCTTCCCCGTCTTCCTCGACCCAGCCACCCGG‑Cgcgcagcgca                                                                                                                >  6:149936/1‑135 (MQ=255)
               tACAGAAACGAGCAAGCCTGCCGTTTCATGATGTACGCCAGGTGATTTCCTCTTCACCGTCGGCGCTGATGCGGACCCAGCGGTCAGCATCGTCTTCCCCGTCTTCCTCGCCCCAGCCCCCCGG‑cgcgccgcgcc                                                                                                                >  2:119090/1‑134 (MQ=255)
               tACAGAAACGAGCAAGCCTGCCGTTTCATGATGTACGCCAGGGGATTTCCTCTTCACCGTCGGCGCTGATGCGGACCCAGCGGTCAGCATCTTCTTCCCCGTCTTCCTCGACCCAGCCACCCGG‑cgcgccgcgcc                                                                                                                >  5:197568/1‑134 (MQ=255)
                   gAAACGAGCAAGCCTGCCGTTTCATGATGTACGCCAGGTGATTTCCTCTTCACCGTCGGCGCTGATGCGGACCCAGCGGTCAGCATCTTCTTCCCCGTCTTCCTCGCCCCCGCCCCCCGG‑CGCGCCGCGCCCTTc                                                                                                            >  4:356971/1‑135 (MQ=255)
                     aaCGAGCAAGCCTGCCGTGTCATGATGTACGCCAGGTGATTTCCTCTTCACCGTCGGCGCTGCTGCGGACCCCGCGGTGAGCATCATCTTCCCCGTCTACCACGACCCCGCCACCCCG‑CCCGCCCCCCCCTTCca                                                                                                          >  1:515637/1‑135 (MQ=255)
                       cGAGCAAGCCTGCCGTTTCATGATGTACGCCAGGTGATTTCCTCTTCACCGTCGGCGCTGATGCGGACCCAGCGGTCAGCATCTTCTTCCCCGGCTTCCTCGACCCAGCCACCCGG‑CGCGCCGCGCACATCcacc                                                                                                        >  3:273400/1‑134 (MQ=255)
                         aGCAAGCCTGCCGTTTCATGATGTACGCCAGGTGATTTCCTCTTCACCGTCGGCGCTGATGCGGACCCAGCGGTCAGCATCTTCTTCCCCGTCTTCCTCGACCCAGCCACCCGG‑CGCGCAGCGCACGTCCACAcc                                                                                                      >  4:413174/1‑135 (MQ=255)
                              gCCTGACGTTTCATGATGTACGCCAGGTGATTTCCTCTTCACCGTCGGCGCTGATGCGGACCCAGCGGTCAGCATCATCTTACCCGTCTTCCTCGACCCAGCCAACCGG‑CGCCCAGCCCACTTCACCACCCAgcg                                                                                                 >  8:126200/1‑135 (MQ=255)
                                                tACGCCAGGTGATTTCCTCTTCACCGTCGGCGCTGATGCGGACCCAGCGGTCAGCATCTTCTTACCCGTCTTCCTCGACCCAGCCACCCCG‑CGCGCAGCGCCCTTCCACACCCCGCGCGGCAAAAGCGGCGCggg                                                                               >  1:199940/1‑135 (MQ=255)
                                                  cGCCAGGTGATTTCCTCTTCACCGTCGGCGCTGATGCGGACCCAGCGGTCAGCATCTTCTTCCCCGTCTTCCTCGACCCAGCCCCCCGG‑CGCGCAGCGCACTTCCACACCCACCGCGGCAAATGCGGCGCGGgcg                                                                             >  1:421158/1‑135 (MQ=255)
                                                          gATTTCCTCTTCACCGTCGGCGCTGATGCGGACCCAGCGGTCAGCATCTTCTTCCCCGTCTTCCTCGACCCAGCCACCCGG‑CGCGCAGCGCACTTCCACACCCCGCGCGGCAAATGCGGCGCGGGCGCAGGCAAc                                                                     >  4:409601/1‑135 (MQ=255)
                                                             tcccTCTCCCCCGGCGGCGCTGATGCGCCCCCCGCGGCCAGAATGTTAACCCCCGTCTTCCGC‑CCCCCGCCCCCCGGCCGCCCCGCGCACTTCCACACCCAGCGCGGCAAATGCGGCGCGGGCGCAGGCAACGt                                                                   <  4:255268/132‑1 (MQ=255)
                                                                ctctTCACCGTCGGCGCTGATGCGGACCCAGCGGTCAGCATCTTCTTCCCCGTCTTCCTCGACCCAGCCACCCGG‑CGCGCAGCGCACTTCCACACCCAGCGCGGCAAATGCGGCGCGGGCGCAGGCAACGTCATc                                                               >  7:25286/1‑135 (MQ=255)
                                                                ctctTCACCGTCGGCGCTGATGCGGACCCAGCGGTCAGCATCTTCTTCCCCGTCGTCCTCGACCCCGCCACCCGG‑CGCGCCCCGCCCGTCCACCCCCCCCGCGGCCAATGCGGCGCCGGCGCAGGCAACggcctc                                                               >  4:308489/1‑130 (MQ=255)
                                                                  ctTCACCGTCGGCGCTGATGCGGACCCAGCGGTCAGCATCTTCTTCCCCGTCTTCCTCGACCCAGCCACCCGG‑CGCGCAGCGCACTTCCACACCCCGCGCGGCAAATGCGGCGCGGGCGCAGGCAACGTCATCgg                                                             >  3:58077/1‑135 (MQ=255)
                                                                           cGGCGCTGATGCGGACCCAGCGGTCAGCATCTTCTTCCCCGTCTTCCTCGACCCAGCCACCCGG‑CGCGCAGCGCACTTCCACACCCAGCGCGGCAAATGCGGCGCGGGCGCAGGCAACGTCATCGGCCCACGGCg                                                    >  6:15018/1‑135 (MQ=255)
                                                                                     gcgcccccGGCGGCCAGCTTGTGCTGCCCCGTCTCCCCCCCCCCCCCCCCCCGG‑CGGCCCCCCCACTTCCACACCCAGCGCGGCAAATGCGGCGCGGGCGCAGGCAACGTCATCGGCCCACGGCGTCTGGTCGCt                                          <  1:78708/130‑1 (MQ=255)
                                                                                     gCGGACCCAGCGGTCAGCATCTTCTTCCCCGTCTTCCTCGACCCAGCCCCCCGG‑CGCGCAGCGCCCATCCACACCCAGCGCGGCAAATGCGGCGCGGGCGCAGGCAACGTCATCGGCCCACGGCGTCTGGTCGcg                                          >  3:324041/1‑134 (MQ=255)
                                                                                           cgcgCGTTCAGAATATTCTCCCCCGTCTCCCCCCCCCCCCCCCCCGGC‑CGCCCCCCCCACTTCCACACCCAGCGCGGCAAATGCGGCGCGGGCGCAGGCAACGTCATCGGCCCACGGCGTCTAGTCGCTTTCCAg                                    <  8:470735/132‑1 (MQ=255)
                                                                                           cctgCGGTCCGATTCTCCTCCCCCGTCACCCGCCCCCCCCCCCCCCGG‑CGCGCCGCCCACTTCCACACCCAGCGCGGCAAATGCGGCGCGGGCGCAGGCAACGTCATCGGCCCACGGCGTCTGGTCGCTTTCCAg                                    <  3:145374/132‑1 (MQ=255)
                                                                                           ccAGCGGTCAGCATCTTCTTCCCCGTCTTCCTCGACCCAGCCACCCGG‑CGCGCAGCGCACTTCCCCACCCCGCGCGGCAAATGCGGCGCGGGCGCAGGCAACGTCAGCGGCCCCCGGCGGCGGGTCGCttt                                        >  1:77241/1‑131 (MQ=255)
                                                                                                ggTCAGCATCTTCTTCCCCGTCTTCCTCGACCCAGCCACCCGG‑CGCGCAGCGCACTTCCACACCCAGCGCGGCAAATGCGGCGCGGGCGCAGGCAACGTCATCGGCCCACGGCGTCTGGTCGCTTTCCAGATACa                               >  4:379658/1‑135 (MQ=255)
                                                                                                  gccgCATCTGCTCCCCGCTCTCCCCCCACCCCCCCCCCCGG‑CGCCCCGCCCACTTCCACACCCAGCGCGGCAAATGCGGCGCGGGCGCAGGCAACGTCATCGGCCCACGGCGTCTGGTCGCTTTCCAGATACAgg                             <  5:174488/132‑1 (MQ=255)
                                                                                                          ggctgccccGTCTTCCCCCCCCCACCCCCCCGG‑CGCGCCGCGCACTTCCACACCCAGCGCGGCAAATGCGGCGCGGGCGCAGGCAACGTCATCGGCCCACGGCGTCTGGTCGCTTTCCAGATAAAGGCTGTTCCa                     <  8:501350/130‑1 (MQ=255)
                                                                                                             tcccccGTCTCCCCCCCCCCCCCCCCCCGG‑CGCCCCCCCCACTTCCACACCCAGCGCGGCAAATGCGGCGCGGGCGCAGGCAACGTCATCGGCCCACGGCGTCTGGTCGCTTTCCAGATACAGGCTGTTcc                      <  6:211834/129‑1 (MQ=255)
                                                                                                              tcCCCGTCTTCCTCGACCCAGCCACCCGG‑CGCGCAGCGCACTTCCACACCCAGCGCGGCAAATGCGGCGCGGGCGCAGGCAACGTCATCGGCCCACGGCGTCTGGTCGCTTTCCAGAAAAAGGCGGTTCCAatta                 >  1:401773/1‑134 (MQ=255)
                                                                                                               ccccGTCTCCCCCGCCCCCCCCCCCCGG‑CGCCCCGCCCACTTCCACACCCAGCGCGGCAAATGCGGCGCGGGCGCAGGCAACGTCATCGGCCCACGGCGTCTGGTCGCTTTCCAGATACAGGCTGTTCCATTTcc                <  6:100124/135‑1 (MQ=255)
                                                                                                                cccGTCTTCGACGCCCCACCCCCCCGG‑CGCCCACCCCACTTCCACACCCAGCGCGGCAAATGCGGCGCGGGCGCAGGCAACGTCATCGGCCCACGGCGTCTGGTCGCTTTCCAGATACAGGCTGTTCCATTTcc                <  4:273400/134‑1 (MQ=255)
                                                                                                                   gTCTTCCTCGACCCAGCCACCCGG‑CGCGCAGCGCACTTCCACACCCAGCGCGGCAAATGCGGCGCGGGCGCAGGCAACGTCATCGGCCCACGGCGTCTGGTCGCTTTCCAGATACAGGCTGTTCCATTTCCCTAc            >  7:375640/1‑135 (MQ=255)
                                                                                                                       tCCTCGACCCAGCCACCCGG‑CGCGCAGCGCACTTCCACACCCAGCGCGGCAAATGCGGCGCGGGCGCAGGCAACGTCATCGGCCCACGGCGTCTGGTCGCTTTCCAGATACAGGCTGTTCCATTTCCCTACAGcc        >  8:126137/1‑135 (MQ=255)
                                                                                                                        ccTCGACCCAGCCACCCGG‑CGCGCAGCGCACTTCCACACCCAGCGCGGCAAATGCGGCGCGGGCGCAGGCAACGTCATCGGCCCACGGCGTCTGGTCGCTTTCCAGATACAGGCTGTTCCATTTCCCTACAGccg       >  6:149954/1‑134 (MQ=255)
                                                                                                                            gACCCAGCCACCCGG‑CGCGCAGCGCACTCCCACACCCAGCGCGGCAAATGCGGCGCGGGCGCAGGCAACGTCATCGGCCCACGGCGTCTGGTCGCTTTCCAGATACAGGCTGTTCCATTTCCCTACAGCCTTgg    >  4:462782/1‑134 (MQ=255)
                                                                                                                             aCCCCGCCCCCCGG‑CGCCCACCCCACTTCCACACCCAGCGCGGCAAATGCGGCGCGGGCGCAGGCAACGTCATCGGCCCACGGCGTCTGGTCGCTTTCCAGATACAGGCTGTTCCATTTCCCTACAGCCTTGGGt  <  6:306059/135‑1 (MQ=255)
                                                                                                                                     |                                                                                                                               
GGGTAAACCCATTCCTACAGAAACGAGCAAGCCTGCCGTTTCATGATGTACGCCAGGTGATTTCCTCTTCACCGTCGGCGCTGATGCGGACCCAGCGGTCAGCATCTTCTTCCCCGTCTTCCTCGACCCAGCCACCCGG‑CGCGCAGCGCACTTCCACACCCAGCGCGGCAAATGCGGCGCGGGCGCAGGCAACGTCATCGGCCCACGGCGTCTGGTCGCTTTCCAGATACAGGCTGTTCCATTTCCCTACAGCCTTGGGT  >  ACB122/1721745‑1722004

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