Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*ACB1431,986,6610AC53.8% 1.9 / 16.6 26intergenic (+42/+56)mshA_1/CJ019_01891D‑inositol‑3‑phosphate glycosyltransferase/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (12/0);  new base C (9/5);  total (21/5)
Fisher's exact test for biased strand distribution p-value = 4.25e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.31e-02
Rejected as consensus: E-value score below prediction cutoff.
Rejected as consensus: Frequency below/above cutoff threshold.

TGAAAGTGCAGCTGAACTAAAACGCCTGGTGTTGGATGTTACGGAAAATACTCAGCTGAAGCATAAGGCGCCTGCTCAAGCGCACTGATTGCGTACGCTTGGCCATCCATGAACCCTGCAGTGCGGCGCAGCCG‑CACTGCAGGGCTTTTGCCATATGACTTGTATGAGGCACGCAGTAAACGCGCTTATCCGGCTACGGGCACAGCCTGGCGACGCTGGCGCAGATAGCGTTGAGCCGGGCGTTCCAGCCAGCGGTGGGTCAACA  >  ACB143/1986532‑1986796
                                                                                                                                 |                                                                                                                                        
tGAAAGTGCAGCTGAACTAAAACGCCTGGTGTTGGATGTTACGGAAAATACTCAGCTGAAGCATAAGGCGCCTGCTCAAGCGCACTGATTGCGTACGCTTGGCCATCCATGAACCCTGCAGTGCGGCGCAGCCG‑c                                                                                                                                    >  4:242206/1‑135 (MQ=255)
 gAAAGTGCAGCTGAACTAAAACGCCTGGTGTTGGATGTTACGGAAAATACTCAGCTGAAGCATAAGGCGCCTGCTCAAGCGCACTGATTGCGTACGCTTGGCCATCCATGAACCCTGCAGTGCGGCGCAGCCG‑c                                                                                                                                    >  4:254763/1‑134 (MQ=255)
 gAAAGTGCAGCTGAACTAAAACGCCTGGTGTTGGATGTTACGGAAAATACTCAGCTGAAGCATAAGGCGCCTGCTCAAGCGCACTGATTGCGTACGCTTGGCCATCCATGAACCCTGCAGTGCCGCGCCGACG‑CAc                                                                                                                                  >  8:242138/1‑136 (MQ=255)
                 aaaacCGCGTGGTGTTGTTTGTTACGCAAAATACTCAGCTCACCAAAATGGCGCCGGCACACCCCAACTTTTTGCTTCGGCTGGCCCACCCATCAACCCTCCCGGGCGGCGCCGCCG‑CACTGCAGGGCTTTTGcc                                                                                                                   <  4:302662/134‑1 (MQ=255)
                   aaaCGCCTGGTGTTGGATGTTACGGAAAATACTCAGCTGAAGCATAAGGCGCCTGCTCAAGCGCACTGATAGCGTACGCTTCGCCATCCATGAACCCTGCAGTGCGGCCCAGCCG‑AACTGCAGGGAATTTGACata                                                                                                                >  2:113359/1‑136 (MQ=255)
                                                        tGAAGCATAAGGCGCCTGCTCAAGCGCACTGATTGCGTACGCTTGGCCATCCATGAACCCTGCCGTGCGGCCCAGCCG‑CCCTGCAGGGCGTTTGCCATTTGACATGGAATAGGCCCCCCGTAAACGCGCGTATCCg                                                                           >  8:337414/1‑136 (MQ=255)
                                                         gAAGCATAAGGCGCCTGCTCAAGCGCACTGATTGCGTACGCTTGGCCATCCATGAACCCTGCAGTGCGGCGCAGCCG‑CACTGCAGGGCTTTTGCCATATGACTTGTATGCGGCACGCAGTAAACGCGCTTATCCgg                                                                          >  5:359352/1‑136 (MQ=255)
                                                                    cgcCTGCTCAAGCGCACTGATTGCGTACGCTTGGCCATCCATGAACCCTGCAGTGCGGCGCAGCCG‑CACTGCAGGGCTTTTGCCATATGACTTGTATGAGGCACGCAGTAAACGCGCTTATCCGGCTACGGGcaca                                                               >  7:184038/1‑136 (MQ=255)
                                                                          cTCAAGCGCACTGATTGCGTACGCTTGGCCATCCATGAACCCTGCCGTGCGGGGCCGCCG‑CCCGGCGGGGGGTTTGGCATTTGACTTGTTTGAGGGCCGCAGGAAACGCGCGTATACGGGGACGGGCACCGCCTg                                                          >  4:345260/1‑135 (MQ=255)
                                                                          cTCAAGCGCACTGATTGCGTACGCTTGGCCATCCATGAACCCTGCAGTTGGGCGCCGCCG‑CACTGGCGGGGGTTTGCCCTATGACTTGTATGAGGCCCGCAGTAAACCCGCTTATACGGCGACGGGCCCAGccag                                                          >  4:73215/1‑133 (MQ=255)
                                                                               gcgcACTGATTGCGTACGCTTGGCCATCCATGAACCCTGCAGTGCGGCGCCGCCG‑CCCTGCAGGGCTTTTGCCATATGACTTGTATTAGGCACGCAGTAAACGCGCTTATCCGGCTACGGGCACCGCCTGGCGACg                                                    >  6:332737/1‑136 (MQ=255)
                                                                                   agttaTTGCGTACGCTTGGCCACCCATCACCCCTGCAGGGCGGCGCCGCCG‑CACTGCAGGGCTTTTGCCATATGACTTGTATGAGGCACGCAGTAAACGCGCTTATCCGGCTACGGGCACAGCCTGGCGACGCTgg                                                <  7:62011/132‑1 (MQ=255)
                                                                                          gCGTACGCTTGGCCATCCATGAACCCTGCAGTGCGGCGCAGCCG‑CACTGCAGGGCTTTTGCCATATGACTTGTATGAGGCACGCAGTAAACGCGCTTATCCGGCTACGGGCACAGCCTGGCGACGCTGGCGCAGa                                          >  7:259229/1‑135 (MQ=255)
                                                                                            gTACGCTTGGCCATCCATGAACCCTGCAGTGGGGGGCCGCCG‑CACGGCAGGGGGTTTTGCATATTTCTTGTTTGGGGGCCGCAGTTAAAGCGGGTATACCGCGACCGGGCCAGCCGGGCGACGCGGGCGCAggaag                                       >  2:156729/1‑132 (MQ=255)
                                                                                            gTACGCTTGGCCATCCATGAACCCTGCAGTGCGGCGCAGCCG‑CACTGCAGGGCTTTTGCCATATGACTTGTATGAGGCACGCAGTAAACGCGCTTATCCGGCTACGGGCACAGCCTGGCGACGCTGGCGCAGATAg                                       >  1:277848/1‑136 (MQ=255)
                                                                                                  ttGGCCATCCATGAACCCTGCAGTGCGGCGCAGCCG‑CACTGCAGGGCTTTTGCCATATGACTTGTATGAGGCACGCAGTAAACGCGCTTATCCGGCTACGGGCACAGCCTGGCGACGCTGGCGCAGATAGCGtt                                   >  6:361773/1‑134 (MQ=255)
                                                                                                    ggCCATCCAGCACCCCTGCCGGGGGGCCCCGCCG‑CACTGCAGGGCTTTTGCCATATGACTTGTATGAGGCACGCAGTAAACGCGCTTATCCGGCTACGGGCACAGCCTGGCGACGCTGGCGCAGATAGCGTTGAGc                               <  7:352946/136‑1 (MQ=255)
                                                                                                         tccatCAACCCTCCAGGGCGGCGCCCCCG‑CACTGCAGGGCTTTTGCCATATGACTTGTATGAGGCACGCAGTAAACGCGCTTATCCGGCTACGGGCACAGCCTGGCGACGCTGGCGCAGATAGCGTTGAGCCGGGc                          <  6:331288/136‑1 (MQ=255)
                                                                                                               aaCCCGGCAGTGCGGCGCCGCCG‑CACTGCAGGGCTTTTGCCATATGACTTGTATGAGGCACGCAGTAAACGCGCTTATCCGGCTACGGGCACAGCCTGGCGACGCTGGCGCAGATAGCGTTGAGCCGGGCGTTcc                     <  5:321961/135‑1 (MQ=255)
                                                                                                                       aGTGCGGCGCCGCCG‑CACTGCAGGGCTTTTGCCATATGACGTGTATGAGGCACGCAGTAAACGCGCTTATCCGGCTACGGGCACAGCCTGGCGACGCTGGCGCAGATAGCGTTGAGCCGGGCGTTCCAGCCCGCgg            >  7:264715/1‑136 (MQ=255)
                                                                                                                       aGTGCGGCGCCGCCG‑CAATGCAGGGCTTTTGCCATATGACTTGTATGAGGCACGCAGTAAACGCGCGTATACGGCTACGGGCACAGCCTGGCGACCCTGGCGCAGATAGCGTTGAGCCGGGCGTTGCAGCCAGCgg            >  7:92259/1‑136 (MQ=255)
                                                                                                                       aGTGCGGCGCCGCCC‑CACGGCAGGGCTTTTGCCATATGACGTGTATGAGGCACGCAGTAAACGCGCTTATCCGGCTACGGGCACAGCCTGGCGACGCGGGCGCAGCTAGCGTTGAGCCGGGCGTTCCAGCCAGCgg            >  2:323055/1‑136 (MQ=255)
                                                                                                                       aGTGCGGCGCAGCCG‑CACTGCCGGGCTTTTGCCATATGTCTTGTATGAGGCACGCAGTAAACGCGCGTATCCGGCTACGGGCCCAGCCTGGCGACGCGGGCGCAGATAGCGTTGAGCCGGGCGTTCCAGCCCGCgg            >  5:166455/1‑136 (MQ=255)
                                                                                                                       aGTGCGGCGCAGCCG‑CACTGCAGGGCTTTTTCCATATTACGTGTATGAGGGACGCAGTAAACGCGCTTATACGGCTACGGGCACAGCCTGGCGACGCTGGCGCAGATAGCGTTGAGCCGGGCGTTCCAACCAGCgg            >  6:310931/1‑136 (MQ=255)
                                                                                                                       aGTGCGGCGCAGCCG‑CACTGCAGGGCTTTTGCCATATGACTTGTATGAGGCACGCAGTAAACGCGCTTATACGGCTACGGGCACAGCCTGGCGACGCTGGCGCAGATAGCGTTGAGCCGGGCGTTCCAGCCAGCgg            >  7:259607/1‑136 (MQ=255)
                                                                                                                         tGCGGCGCCGCCG‑CACTGCAGGGCGTTTGCCATATGACTTGTATGAGGCACGCAGTAAACGCGCGTATCCGGCTACGGGCACCGCCTGGCGACGCTGGCGCAGCTAGCGTTGAGCCGGGCGTTCCAGCCAGCGGTg          >  8:74345/1‑136 (MQ=255)
                                                                                                                              cgccgccgcCACTGCAGGGCTTTTGCCATATGACTTGTATGAGGCACGCAGTAAACGCGCTTATCCGGCTACGGGCACAGCCTGGCGACGCTGGCGCAGATAGCGTTGAGCCGGGCGTTCCAGCCAGCGGTGGGt       <  3:242206/135‑1 (MQ=255)
                                                                                                                                 aGCCG‑CACTGCAGGGCTTTTGCCATATGACTTGTATGAGGCACGCAGTAAACGCGCTTATCCGGCTACGGGCACAGCCTGGCGACGCTGGCGCAGATAGCGTTGAGCCGGGCGTTCCAGCCAGCGGTGGGTCAACa  <  1:228244/136‑1 (MQ=255)
                                                                                                                                 |                                                                                                                                        
TGAAAGTGCAGCTGAACTAAAACGCCTGGTGTTGGATGTTACGGAAAATACTCAGCTGAAGCATAAGGCGCCTGCTCAAGCGCACTGATTGCGTACGCTTGGCCATCCATGAACCCTGCAGTGCGGCGCAGCCG‑CACTGCAGGGCTTTTGCCATATGACTTGTATGAGGCACGCAGTAAACGCGCTTATCCGGCTACGGGCACAGCCTGGCGACGCTGGCGCAGATAGCGTTGAGCCGGGCGTTCCAGCCAGCGGTGGGTCAACA  >  ACB143/1986532‑1986796

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