Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*ACB1221,740,3990TG52.0% 4.2 / 13.7 25T56P (ACC→CCC) CJ019_01643hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (5/7);  new base G (0/13);  total (5/20)
Fisher's exact test for biased strand distribution p-value = 1.49e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.50e-01
Rejected as consensus: E-value score below prediction cutoff.
Rejected as consensus: Frequency below/above cutoff threshold.

CGCCAGGAGCGCCTTTAGCGTGTCCACCTTCTCTTGGAAGGCCGGGTGCATTGTTGTCTGCATGGATGATCCTCGGGTGGGTCAGGCGTGGTATTCGAAGGCCTCGGCCTTACGAACGATTCGAACTTGGGCGGTGCGCCGCTCCGGCGCGCGGCGGTCGCGGCGCATGGGGTCGCTGTCGTTGATGGCGGCATGCATTGCCACAAGACAGACCAGGGCGAGGCAGAGCGGAGTGATGATCTTCTGGTATGCCGCCTCTGCTACGA  >  ACB122/1740265‑1740530
                                                                                                                                      |                                                                                                                                   
cgcCAGGAGCGCCTTTAGCGTGTCCACCTTCTCTTGGAAGGCCGGGTGCATTGTTGTCTGCATGGATGAGCCTCGGGTGGGGCAGGCGTGGTATTCGAAGGCCTCGGCCTTACGAACGATTCGAACTTGGGCGGt                                                                                                                                     <  4:405142/135‑1 (MQ=255)
 gcCAGGAGCGCCTTTAGCGTGTCCACCTTCTCTTGGAAGGCCGGGTGCATTGTTGTCTGCATGGATGATCCTCGGGTGGGTCAGGCGTGGTATTCGAAGGCCTCGGCCTTACGAACGATTCGAACTTGGGCGGTg                                                                                                                                    >  3:11173/1‑135 (MQ=255)
 gcCAGGAGCGCCTTTAGCGTGTCCACCTTCTCTTGGAAGGCCGGGTGCATTGTTGTCTGCATGGATGATCCTCGGGTGGGTCAGGCGTGGTATTCGAAGGCCGCGGCCTTACGAACGATTCGAACTTCGGCGGTg                                                                                                                                    >  4:84450/1‑135 (MQ=255)
 gcCAGGAGCGCCTTTAGCGTGTCCACCTTCTCTTGGAAGGCCGGGTGCATTGTTGTCTGCACGGATGATCCTCGGGTGGGTCAGGCGTGGTATTCGAAGGCCTCGGCCTTACGACCGCTTCGAACTAGTGCGGTg                                                                                                                                    >  5:169482/1‑135 (MQ=37)
   cAGGAGCGCCTTTAGCGTGTCCACCTTCTCTTGGAAGGCCGGGTGCATTGTTGTCTGCATGGATGATCCTCGGGTGGGTCAGGCGTGGTATTCGAAGGCCTCGGCCTTACGAACGATTCGAACTTGGGCGGTggg                                                                                                                                  >  2:338875/1‑133 (MQ=255)
   cAGGAGCGCCTTTAGCGTGTCCACCTTCTCTTGGAAGGCCGGGTGCATTGTTGTCTGCATGGATGATCCTCGGGTGGGTCAGGCGTGGTATTCGAAGGCCTCGGCCTTACGAACGATTCGAACTTGGGCGGTggg                                                                                                                                  >  7:449249/1‑133 (MQ=255)
                               tctTGGAAGGCCGGGTGCATTGTTGTCTGCATGGATGAGCCTCGGGGGGGTCAGGCGTGGTATTCGAAGGCCGCGGCCTTACGAACGATTCGAACTTGGGCGGTGCGCCGCGCCGGCGCGCGGCGGTCgcggc                                                                                                        <  7:534803/133‑1 (MQ=255)
                               gctTGGCCGGCCGGGTGCTGTGTTGTCGGCATGGTGCATCCTCGGGGGGGGCTGGCGTGGTATACCCCGCCCGCGGCCTTAGGAACGATTCTACCTTGGGCGGGGCGCCGCGCCGGCGCGCGGCGGTCgcggc                                                                                                        <  2:106756/132‑1 (MQ=255)
                                    gAAGGCCGGGTGCATTGTTGTCTGCATGGATGATCCTCGGGTGGGTCTGGCGTGGTATTCGAAGGCCTCGGCCTTACGAACGATTCGAACTTGGGCGGTGCGCCGCTCCGGCGCGCGGCGGTCGCGGCGCAt                                                                                                    <  4:415019/132‑1 (MQ=255)
                                     aaGGCCGGGTGCATTGTTGTCTGCATGGATGATCCTCGGGTGGGGCAGGCGTGGTATTCGAAGGCCTCGGCCTTACGAACGATTCGAACTTGGGCGGTGCGCCGCTCCGGCGCGCGGCGGTCGCGGCGCATgggg                                                                                                <  3:276796/135‑1 (MQ=255)
                                                     ttgtCTGATGGGATGAACCTGGGTGGGGTCAGGCGTGGTATTCTACGGCCTCGGCATTACTAACGATGCGAAGTTGGGCGGGGCGCCGGTCGGGGGCGGGGGGTTCGCGGCGCATGGGGTCGCTGTCGTTGATg                                                                                 <  8:353472/134‑1 (MQ=9)
                                                            cagggATGAGCCGGGGGGGGGCCAGGCGTGGTATGCGACGGCCGCGGCCTTACGAACGATGCGACCTGGGGCGGGGCGCCGCGCCGGCGCGCGGCGGTCGCGGCGCATGGGGTCGCTGTCGTTGATGGCGgcatg                                                                         <  3:176398/132‑1 (MQ=17)
                                                            cATGGATGAGCCGCGGGTGGGGCCGGCGGGGTATTCGACGGCCGCGGCCTTACGAACGATTCGAACTTGGGCGGGGCGCCGCTCCGGCGCGCGGCGGTCGCGGCGCATGGGGTCGCTGTCGTTGATGGCGgcatg                                                                         <  8:363122/135‑1 (MQ=17)
                                                            cATGGAGGAGCCGGGGGTGGGTCTGGCGTGGTATGCGACGGCCGCGGCCTTACGAACGATGCTAACTTGGGCGGGGCGCCGCTCGGGCGCGCGGCGGGCGCGGCGCATGGGGTCGCTGTCGTTGATGGCGgcatg                                                                         <  7:282454/135‑1 (MQ=17)
                                                                 atgatCCGCGGGTGGGTCAGGCGTGGTATGCGAAGGCCTCGGCCTTACGAACGATTCGAACTTGGGCGGGGCGCCGCTCCGGCGCGCGGCGGTCGCGGCGCATGGGGTCGCTGTCGTTGATGGCGGCATGCATTg                                                                    <  6:244955/135‑1 (MQ=255)
                                                                      ccTCGGGTGGGTCAGGCGTGTTATGCTAAGGCCTCGGCCTTACGAACGATTCGAACTTGGGCGGTGCGCCGCTCCGGCGCGCGGCGGTCGCGGCGCATGGGGTCGCTGTCGTTGATGGCGGCATGCATTGCcaca                                                               <  8:449249/135‑1 (MQ=255)
                                                                       ggggggggggtCTGTGGTGTTATTCTAAGGCCTGGGCATTAGTAACTTTTCGACATTGGGCGGGGCGCGGCTGGGGGGCGCGGAGGTCGCGGCGCATGGGGTCGCTGTCTTTGATGGCGGCATGCATTGCcaca                                                               <  6:443124/131‑1 (MQ=12)
                                                                        gggggtgtgtCTGGCGTGGACTTCGACGGCCGCGGCCTTACGAACGATTCGAACTTGGGCGGTGCGCCGCTCCGGCGCGCGGCGGGCGCGGCGCATGGGGTCGCTGTCGTTGATGGCGGCATGCATTGCcaca                                                               <  7:300537/131‑1 (MQ=18)
                                                                                gtCTGGGGTGGTTTTCGTTGGCCGGGGCTTAACTAACGATTCGAACTTGGGCGGGGCGCGGCTCCGGCGCGCGGCGGTCGCGGCGCATGGGGTCACTGTCGTTGATGGCGGCATGCATTGCCACAAGACAGACCa                                                     <  2:435183/135‑1 (MQ=21)
                                                                                       ggggTATGCGAAGGCCTCGGCCTTACGAACGATTCGAACTTGGGCGGTGCGCCGCGCCGGCGCGCGGCGGTCGCGGCGCATGGGGTCGCTGTCGTTGATGGCGGCATGCATTGCCACAAGACAGACCAGGGCGAg                                              <  4:392416/133‑1 (MQ=255)
                                                                                         gttATGCGACGGCCGCGGCCTACCCAACGATTCGACGTGGGGCGGGGCGCCGCGCGGGCGCGCGGCGGTCGCGGCGCATGGGGTCGCTGTCGTTGATGGCGGCATGCATTGCCACAAGACAGACCAGGGCGAGGc                                            <  4:11173/133‑1 (MQ=21)
                                                                                             ctccacggCCGCGGCATTACTAACCACCCCACCTGGGGCGGGGCGCCGCGCCGGCGCGCGGCGGGCGCGGCGCATGGGGTCGCTGTCGTTGATGGCGGCATGCATTGCCACAAGACAGACCAGGGCGAGGCaga                                         <  1:282735/128‑1 (MQ=25)
                                                                                                  aGGCCGCGGCCTTACTAACGATTCGAACTTGCGCGGGGCGCCGCTCCGGCGGGCGGCGGTCGCGGCGCATGGGGTCGCTGTCGTTGATGGCGGCATGCATTGCCACAAGACAGACCAGGGCGAGGCAGAGCGGAg                                   <  6:276003/135‑1 (MQ=25)
                                                                                                    gCCTCGGCCTTAGGAAGGATGCGAACTTGGGCGGTGCGCCGGTCCGGCGCGCGGCGGTCGCGGCGCATGGGGTCGCTGTCGTTGATGGCGGCATGCATTGCCACAAGACAGACCAGGGCGAGGCAGAGCGGAGtg                                 <  6:532978/135‑1 (MQ=255)
                                                                                                                              ttGGGCGGGGCGCGGCGCCGGCGCGCGGCGGTCGCGGCGCATGGGGTCGCTGTCGTTGATGGCGGCATGCATTGCCACAAGACAGACCAGGGCGAGGCAGAGCGGAGTGATGATCTTCTGGTATGCCGCCTCTGc       <  8:430/135‑1 (MQ=255)
                                                                                                                                   cGGGGCGCGGCTCCGGCGCGCGGCGGTCGCGGCGCTTGGGGTCGCTGTCGTTGATGGAGGCATGCATTGCCACAAGACAGACCAGGGCGAGGCAGAGCGGAGTGATGATCTTCTGGTATGCCGCCTCTGCTACga  <  8:334253/135‑1 (MQ=255)
                                                                                                                                      |                                                                                                                                   
CGCCAGGAGCGCCTTTAGCGTGTCCACCTTCTCTTGGAAGGCCGGGTGCATTGTTGTCTGCATGGATGATCCTCGGGTGGGTCAGGCGTGGTATTCGAAGGCCTCGGCCTTACGAACGATTCGAACTTGGGCGGTGCGCCGCTCCGGCGCGCGGCGGTCGCGGCGCATGGGGTCGCTGTCGTTGATGGCGGCATGCATTGCCACAAGACAGACCAGGGCGAGGCAGAGCGGAGTGATGATCTTCTGGTATGCCGCCTCTGCTACGA  >  ACB122/1740265‑1740530

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