Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*ACB1223,047,2170CG31.7% 86.5 / 11.1 45L122L (CTG→CTC) CJ019_02862hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (2/26);  new base G (13/0);  total (18/27)
Fisher's exact test for biased strand distribution p-value = 5.96e-09
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

GGCCGTCGCGGTACCAGTACAGGTTGTTATAGCCCAAGGCTTTCGCGCGTTTTACCGCGTTCCAGCTTAACCAGCAATCGGAGCGGCAGTAAAAGACCAGTGGTTGGTCCAGGTGCCCGGCCGT‑GGC‑TTTTTGCAGGTTACGCACGAAGTAGTCCTGCCATGCCGGGCTCAGTTCGCCGTCGCCGGTGTTGGCCAGCCACAGGCTGCCAGGGAGGTTGCTATGCGGCTGGTCTTCTATGAAGCGACCTTGCAGCCATTGACGGCG  >  ACB122/3047084‑3047348
                                                                                                                                       |                                                                                                                                   
ggCCGTCGCGGTACCAGTACAGGTTGTTATAGCCCAAGGCTTTCGCGCGTTTTACCGCGTTCCAGCTTAACCAGCAATCGGAGCGGCAGTAAAAGACCAGTGGTTGGTCCAGGTGCCCGGCCGT‑GGC‑TTTTTGCa                                                                                                                                    <  8:348151/135‑1 (MQ=255)
     tCGCGGTACCAGTACAGGTTGTTATAGCCCAAGGCTTTCGCGCGTTTTACCGCGTTCCAGCTTAACCAGCAATCGGAGCGGCAGTAAAAGACCAGTGGTTGGTCCAGGTGCCCGGCCGT‑GGC‑TTTTTGCAGGTTa                                                                                                                               <  7:271752/135‑1 (MQ=255)
      cgcgGTACCAGTACAGGTTGTTATAGCCCAAGGCTTTCGCGCGTTTTACCGCGTTCCAGCTTAACCAGCAATCGGAGCGGCAGTAAAAGACCAGTGGTTGGTCCAGGTGCCCGGCCGT‑GGC‑TTTTTGCAGGTTAc                                                                                                                              <  8:218674/135‑1 (MQ=255)
       gcgGTACCAGTACAGGTTGTTATAGCCCAAGGCTTTCGCGCGTTTTACCGCGTTCCAGCTTAACCAGCAATCGGAGCGGCAGTAAAAGACCAGTGGTTGGTCCAGGTGCCCGGCCGT‑GGC‑GTTTTGCAGGTTTCg                                                                                                                             >  2:120276/1‑135 (MQ=255)
        cgGTACCAGTACAGGTTGTTATAGCCCAAGGCTTTCGCGCGTTTTACCGCGTTCCAGCTTAACCAGCAATCGGAGCGGCAGTAAAAGACCAGTGGTTGGTCCAGGTGCCCGGCCGT‑GGC‑TTTTTGTAGGTTACGc                                                                                                                            >  2:76116/1‑135 (MQ=255)
          gTACCAGTACAGGTTGTTATAGCCCAAGGCTTTCGCGCGTTTTACCGCGTTCCAGCTTAACCAGCAATCGGAGCGGCAGTAAAAGACCAGTGGTTGGTCCAGGTGCCCGGCCGT‑GGC‑TTTTTGCAGGTTACGCAc                                                                                                                          <  4:115325/135‑1 (MQ=255)
           tACCAGTACAGGTTGTTATAGCCCAAGGCTTTCGCGCGTTTTACCGCGTTCCAGCTTAACCAGCAATCGGAGCGGCAGTAAAAGACCAGTGGTTGGTCCAGGTGCCCGGCCGT‑GGC‑TTTTTGCAGGTTACGCACg                                                                                                                         <  8:289331/135‑1 (MQ=255)
           tACCAGTACAGGTTGTTATAGCCCAAGGCTTTCGCGCGCTTTACCGCGTTCCAGCTTAACCAGCAATCGGAGCGGCAGTAAAAGACCAGTGGTTGGTCCAGGTGCCCGGCCGT‑GGC‑TTTTTGCAGGTTACGCACg                                                                                                                         <  6:62821/135‑1 (MQ=255)
              cAGTACAGGTTGTTATAGCCCAAGGCTTTCGCGCGTTTTACCGCGTTCCAGCTTAACCAGCAATCGGAGCGGCAGTAAAAGACCAGTGGTTGGTCCAGGTGCCCGGCCGT‑GGC‑TTTTTGCAGGTTACGCACGAag                                                                                                                      <  8:63365/135‑1 (MQ=255)
                gTACAGGTTGTTATAGCCCAAGGCTTTCGCGCGTTTTACCGCGTTCCAGCTTAACCAGCAATCGGAGCGGCAGTAAAAGACCAGTGGTTGGTCCAGGTGCCCGGCCGT‑GGC‑TTTTTGTAGGTTACGCACGAagta                                                                                                                    <  1:76116/135‑1 (MQ=255)
                gTACAGGTTGTTATAGCCCAAGGCTTTCGCGCGTTTTACCGCGTTCCAGCTTAACCAGCAATCGGAGCGGCAGTAAAAGACCAGTGGTTGGTCCAGGTGCCCGGCCGT‑GGC‑TTTTTGCAGGTTACGCACGAagta                                                                                                                    <  6:143308/135‑1 (MQ=255)
                     ggTTGTTATAGCCCAAGGCTTTCGCGCTTTTTACCGCGTTCCAGCTTAACCAGCAATCGGAGCGGCAGTAAATGAGCAGTGGTTGGTCCAGGTGCCCGGCCGT‑GGC‑TTTTTGCAGGTTAAGCACGAAGTAGTCCt                                                                                                               <  7:163563/135‑1 (MQ=255)
                        tgttATAGCCCAAGGCTTTCGCGCGTTTTACCGCGTTCCAGCTTAACCAGCAATCGGAGCGGCAGTATAAGACCAGTGGTTGGTCCAGGTGCCCGGCCGT‑GGG‑GTTTTGGTGGTTTCGACCGAACTTGTACTGc                                                                                                             >  2:383459/1‑134 (MQ=255)
                              aGCCCAAGGCTTTCGCGCGTTTTACCGCGTTCCAGCTTAACCAGCAATCGGAGCGGCAGTAAAAGACCAGTGGTTGGTCCAGGTGCCCGGCCGT‑GGG‑GTTTTGGAGGTTTCGCACGCAGTAGTCCTGCCCTGCCg                                                                                                      >  2:298710/1‑135 (MQ=255)
                              aGCCCAAGGCTTTCGCGCGTTTTACCGCGTTCCAGCTTAACCAGCAATCGGAGCGGCAGTAAAAGACCAGTGGTTGGTCCAGGTGCCCGGCCGT‑GGG‑GTGTTTGAGGTTAAGCCGCAAGTACTCCGGCCATGccc                                                                                                      >  8:361739/1‑134 (MQ=255)
                              aGCCCAAGGCTTTCGCGCGTTTTACCGCGTTCCAGCTTAACCAGCAATCGGAGCGGCAGTAAAAGACCAGTGGTTGGTCCAGGTGCCCGGCCGT‑GGC‑TTTTTGCAGGTTACGCACGAAGTAGTCCTGCCATGCCg                                                                                                      >  4:164709/1‑135 (MQ=255)
                              aGCCCAAGGCTTTCGCGCGTTTTACCGCGTTCCAGCTTAACCAGCAATCGGAGCGGCAGTAAAAGACCAGTGGTTGGTCCAGGTGCCCGGCCGT‑GGC‑GTTTTTGTGGGTAGCCACGCACTAGTGCTGCCATGCCg                                                                                                      >  6:328472/1‑135 (MQ=255)
                              aGCCCAAGGCTTTCGCGCGTTTTACCGCGTTCCAGCTTAACCAGCAATCGGAGCGGCAGTAAAAGACCAGTGGTTGGTCCAGGTGCCCGGCCGT‑GGC‑GTTTTGGATGTTAAGCCCGAAGTAGTACTGCCAAGCCg                                                                                                      >  4:264749/1‑135 (MQ=255)
                              aGCCCAAGGCTTTCGCGCGTTTTACCGCGTTCCAGCTTAACCAGCAATCGGAGCGGCAGTAAAAGACCAGTGGTTGGTCCAGGTCCCCGGCCGT‑GGC‑TTTTTTGAGGGTACGCCCGAAGTAGTGCTGCCATGccc                                                                                                      >  4:337979/1‑134 (MQ=255)
                                   aaGGCTTTCGCGCGTTTTACCGCGTTCCAGCTTAACCAGCAATCGGAGCGGCAGTAAAAGACCAGTGGTTGGTCCAGGTGCCCGGCCGTGGGG‑TTTTTGACGGTT‑CGGCCCACGTAGTCCTGCCAGGCCGGGcgc                                                                                                 >  4:373769/1‑133 (MQ=255)
                                   aaGGCTTTCGCGCGTTTTACCGCGTTCCAGCTTAACCAGCAATCGGAGCGGCAGTAAAAGACCAGTGGTTGGTCCAGGTGCCCGGCCGT‑GGG‑GTTTTTGTGGTTTAGCTCGAAGTATGCCTGGCCTGCCgggggc                                                                                                 >  1:206108/1‑132 (MQ=255)
                                       cTTTCGCGCGTTTTACCGCGTTCCAGCTTAACCAGCAATCGGAGCGGCAGTAAAAGACCAGTGGTTGGTCCAGGTGCCCGGCCGT‑GGG‑GGTTTGGAGGTTACGCACGCAGTAGTCCTGCCATGCCGGGCTCAggt                                                                                             >  3:98620/1‑133 (MQ=255)
                                             cgcgTTTTACCGCGTTCCAGCTTAACCAGCAATCGGAGCGGCAGTAAAAGACCAGTGGTTGGTCCAGGTGCCCGGCCGT‑GGG‑TTTTTGGAGGGTTCGCACGCAGTAGTCCTGCCATTCCGGGCTACGTtcggcgt                                                                                       >  7:219270/1‑135 (MQ=255)
                                             cgcgTTTTACCGCGTTCCAGCTTAACCAGCAATCGGAGCGGCAGTAAAAGACCAGTGGTTGGTCCAGGTGCCCGGCCGT‑GGC‑TTTTTGCAGGTTACGCACGAAGTAGTCCTGCCATGCCGGGCTCAGTtcgccgt                                                                                       <  3:373769/135‑1 (MQ=255)
                                             cgcgTTTTACCGCGTTCCAGCTTAACCAGCAATCGGAGCGGCAGTAAAAGACCAGTGGTTGGTCCAGGTGCCCGGCCGT‑GGC‑TTTTTGCAGGTTACGCACGAAGTAGTCCTGCCATGCCGGGCTCAGTtcgccgt                                                                                       <  4:200015/135‑1 (MQ=255)
                                             cgcgTTTTACCGCGTTCCAGCTTAACCAGCAATCGGAGCGGCAGTAAAAGACCAGTGGTTGGTCCAGGTGCCCGGCCGT‑GGC‑TTTTTGCAGGTTACGCACGAAGTAGTCCTGCCATGCCGGGCTCAGTtcgccgt                                                                                       <  7:289360/135‑1 (MQ=255)
                                                  tttACCGCGTTCCAGCTTAACCAGCAATCGGAGCGGCAGTAAAAGACCAGTGGTTGGTCCAGGTGCCCGGCCGT‑GGG‑GTTTTTGAGGTTTCGCCCGAAGTAGTCCGGGCCAGCCCGGCGGCGTtggccgtcgccg                                                                                  >  3:145773/1‑135 (MQ=255)
                                                  tttACCGCGTTCCAGCTTAACCAGCAATCGGAGCGGCAGTAAAAGACCAGTGGTTGGTCCAGGTGCCCGGCCGT‑GGCGGTTTTGAGGGTTA‑GCCCGCAGTAGTCCTGCCATGCCGGGCGCCGTtcgccggcggcg                                                                                  >  5:131117/1‑135 (MQ=255)
                                                  tttACCGCGTTCCAGCTTAACCAGCAATCGGAGCGGCAGTAAAAGACCAGTGGTTGGTCCAGGTGCCCGGCCGT‑GGC‑GGTTTGGAGGGTTCGCCCGAAGTTGTCCTGCCAAGCCGGGGGCAGTtcgccggcgccg                                                                                  >  5:138059/1‑135 (MQ=255)
                                                       cgcgTTCCAGCTTAACCAGCAATCGGAGCGGCAGTAAAAGACCAGTGGTTGGTCCAGGTGCCCGGCCGT‑GGC‑TTTTTGCAGGTTACGCACGAAGTAGTCCTGCCATGCCGGGCTCAGTTCGCCGTCGCCGGTGtt                                                                             <  8:264888/135‑1 (MQ=255)
                                                       cgcgTTCCAGCTTAACCAGCAATCGGAGCGGCAGTAAAAGACCAGTGGGTGGTCCAGGTGCCCGGCCGT‑GGC‑TTTTTGCAGGTTACGCACGAAGTAGTCCTGCCATGCCGGGCTCAGTTCGCCGTCGCCGGTGtt                                                                             <  8:204679/135‑1 (MQ=255)
                                                               aGCTTAACCAGCAATCGGAGCGGCAGTAAAAGACCAGTGGTTGGTCCAGGTGCCCGGCCGT‑GGC‑GTTTTGGAGGTTACGCACGCAGTAGTCCAGCCATGCCGGGCGCAGTTAGCCCTGGCTGGGGTTGgcccccc                                                                     >  7:125119/1‑131 (MQ=255)
                                                               aGCTTAACCAGCAATCGGAGCGGCAGTAAAAGACCAGAGGTTGGTCCAGGTGCCCGGCCGT‑GGC‑TTTTTGCAGGTTACGCACGAAGTAGTCCTGCCATGCCGGGCTCAGTTCGCCGTCGCCGGTGTTGgccagcc                                                                     <  5:170138/135‑1 (MQ=255)
                                                                 cTTAACCAGCAATCGGAGCGGCAGTAAAAGACCAGTGGTTGGTCCAGGTGCCCGGCCGT‑GGC‑TTTTTGCAGGTTACGCACGAAGTAGTCCTGCCATGCCGGGCTCAGTTCGCCGTCGCCGGTGTTGGCCAGCcac                                                                   <  5:31379/135‑1 (MQ=255)
                                                                  ttAACCAGCAATCGGAGCGGCAGTAAAAGACCAGTGGTTGGTCCAGGTGCCCGGCCGT‑GGC‑TTTTTGCAGGTTACGCACGAAGTAGTCCTGCCATGCCGGGCTCAGTTCGCCGTCGCCGGTGTTGGCCAGCcaca                                                                  <  1:86529/135‑1 (MQ=255)
                                                                             tCGGAGCGGCAGTAAAAGACCAGTGGTTGGTCCAGGTGCCCGGCCGT‑GGC‑TTTTTGCAGGTTACGCACGAAGTAGTCCTGCCATGCCGGGCTCAGTTCGCCGTCGCCGATGTTGGCCAGCCACAGGCTGCCAggg                                                       <  1:298710/135‑1 (MQ=255)
                                                                                        gTAACAGACCAGTGGGTGGTCCAGGTGCCCGGCCGT‑GGC‑TTTTTGCAGGTTACGCACGAAGTAGTCCTGCCATGCCGGGCTCAGTTCGCCGTCGCCGGTGTTGGCCAGCCACAGGCTGCCAGGGAGGTTGCTATg                                            <  2:309658/135‑1 (MQ=255)
                                                                                         gaaaaGACCAGTGGTTGGTCCAGGTGCCCGGCCGT‑GGC‑TTTTTGCAGGTTACGCACGAAGTAGTCCTGCCATGCCGGGCTCAGTGCGCCGTCGCCGGTGTTGGCCAGCCACAGGCTGCCAGGGAGGGTGCTATGc                                           <  8:198244/134‑1 (MQ=255)
                                                                                             aGACCAGTGGTTGGTCCAGGTGCGCGGCCGT‑GGC‑TTTTTGCAGGTTACGCACGAAGTAGTCCTGCCATGCCGGGCTCAGTTCGCCGTCGCCGGTGTTGGCCAGCCACAGGCTGCCAGGGAGGTTGCTATGCGGCt                                       <  3:264749/135‑1 (MQ=255)
                                                                                             aGACCAGTGGTTGGTCCAGGTGCCCGGCCGT‑GGC‑TTTTTGCAGGTTACGCACGAAGTAGTCCTGCCATGCCGGGCTCAGTTCGCCGTCGCCGGTGTTGGCCAGCCACAGGCTGCCAGGGAGGTTGCTATGCGGCt                                       <  3:299191/135‑1 (MQ=255)
                                                                                                   gTGGTTGGTCCAGGTGCCCGGCCGT‑GGC‑TTTTTGCAGGTTACGCACGAAGTAGTCCTGCCATGCCGGGCTCAGTTCGCCGTCGCCGGTGTTGGCCAGCCACAGGCTGCCAGGGAGGTTGCTATGCGGCTGGtctt                                 <  2:77516/135‑1 (MQ=255)
                                                                                                        tggtCCAGGTGCCCGGCCGT‑GGC‑TTTTTGCAGGTTACGCACGAAGTAGTCCTGCCATGCCGGGCTCAGTTCGCCGTCGCCGGTGTTGGCCAGCCACAGGCTGCCAGGGAGGTTGCTATGCGGCTGGTCTTCTATg                            <  3:234246/135‑1 (MQ=255)
                                                                                                        tggtCCAGGTGCCCGGCCGT‑GGC‑TTTTTGCAGGTTACGCACGAAGTAGTCCTGCCATGCCGGGCTCAGTTCGCCGTCGCCGGTGTTGGCCAGCCACAGGCTGCCAGGGAGGTTGCTATGCGGCTGGTCTTCTATg                            <  7:361739/135‑1 (MQ=255)
                                                                                                             cAGGTGCCCGGCCGT‑GGC‑TTTTTGCAGGTTACGCACGAAGTAGTCCTGCCATGCCGGGCTCAGTTCGCCGTCGCCGGTGTTGGCCAGCCACAGGCTG‑CAGGGAGGTTGCTATGCGGCTGGTCTTCTATGAAGCGa                      <  1:78489/135‑1 (MQ=255)
                                                                                                                         cGT‑GGC‑GTTTTGGAGGTTACGCACGAAGTAGTCCTGCCATGACGGGCTCAGTTAGCCGGCGCCCGGGTTTGGCAGCCCCACCCTGCCAGGGAGGTTGCTATTAGGCTGGTCGTGGATGAAGAGCCCTTGAAGCCa           >  2:166144/1‑135 (MQ=255)
                                                                                                                                       cAGGTTACGCACGAAGTAGTCCTGCCATGCCGGGCTCAGTTCGCCGTCGCCGGTGTTGGCCAGCCACAGGCTGCCAGGGAGGTTGCTATGCGGCTGGTCTTCTATGAAGCGACCTTGCAGCCATTGACGGCg  >  2:143655/1‑132 (MQ=255)
                                                                                                                                       |                                                                                                                                   
GGCCGTCGCGGTACCAGTACAGGTTGTTATAGCCCAAGGCTTTCGCGCGTTTTACCGCGTTCCAGCTTAACCAGCAATCGGAGCGGCAGTAAAAGACCAGTGGTTGGTCCAGGTGCCCGGCCGT‑GGC‑TTTTTGCAGGTTACGCACGAAGTAGTCCTGCCATGCCGGGCTCAGTTCGCCGTCGCCGGTGTTGGCCAGCCACAGGCTGCCAGGGAGGTTGCTATGCGGCTGGTCTTCTATGAAGCGACCTTGCAGCCATTGACGGCG  >  ACB122/3047084‑3047348

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