Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*ACB1432,094,8170TG58.5% ‑3.5 / 21.0 45Q183P (CAG→CCG) CJ019_01995hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (10/7);  new base G (8/16);  total (22/23)
Fisher's exact test for biased strand distribution p-value = 1.25e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.81e-01
Rejected as consensus: E-value score below prediction cutoff.
Rejected as consensus: Frequency below/above cutoff threshold.

CGGTCCAGTACCTGGGCATTGGCCCGGGCCAACAGTTTCTCAAGATACGTCAGCTGCTCTTCATACACCTTCGGCTCCGGCTGCTTGCGCAAATACTGCACGCCGCGCTCGAAGGCCAGGCGCGCCTGGCCGGGCTGGCCTTTTTGCAGGGCCTGTTGGCCGAGGTTGTTGAAGAACTCTATATGCAGCAGCACCAGGATATGGCGGATTTCCTTGACCCAATGCTTGGCCTCGACAGCGGGCAAGAAGCC  >  ACB143/2094691‑2094941
                                                                                                                              |                                                                                                                            
cGGTCCAGTACCTGGGCATTGGCCCGGGCCAACAGTTTCTCAAGATACGTCAGCTGCTCTTCATACACCTTCGGCTCCGGCTGCTTGCGCAAATACTGCACGCCGCGCTCGAAGGCCCCGCGCGCCGGGCCCGGc                                                                                                                      >  3:238826/1‑135 (MQ=255)
  gTCCAGTACCTGGGCATTGGCCCGGGCCAACAGTTTCTCAAGATACGTCAGCTGCTCTTCATACACCTTCGGCTCCGGCTGCGTGCGCAAAAACTGCACGCCGCGCTCGCAGGCCCCGCCCGCCAGGCCgcgcgg                                                                                                                    >  5:174578/1‑133 (MQ=255)
      aGTACCTGGGCATTGGCCCGGGCCAACAGTTTCTCAAGATACGTCAGCTGCGCTTCATACACCTTCGGCTCCGGCTGCTTGCGCAAATACTGCACGCCGCGCTCGCAGGCCCGGCCCGCCAGGCCGGGAGGGccg                                                                                                                >  4:260922/1‑134 (MQ=255)
               gCATTGGCCCGGGCCAACAGTTTCTCAAGATACGTCAGCTGCTCTTCATACACCTTCGGCTCCGGCTGCTTGCGCAAATACTGCACGCCGCGCTCGAAGGCCCGGCCCGCCGGGCCGGGCGGGGCTTTTTGCAgg                                                                                                       >  3:358401/1‑135 (MQ=255)
                cATTGGCCCGGGCCAACAGTTTCTCAAGATACGTCAGCTGCTCTTCATACACCTTCGGCTCCGGCTGCTTGCGCAAATACTGCACGCCGCGCTCGAAGGCCCGGCGCGCCGGGCCGGGCTGGGCGTTTTTGAggg                                                                                                      >  7:13489/1‑135 (MQ=255)
                cATTGGCCCGGGCCAACAGTTTCTCAAGATACGTCAGCTGCTCTTCATACACCTTCGGCTCCGGCTGCTTGCGCAAATACTGCACGCCGCGCTCGAAGGCCAGGCGCGCCGGGCCGGGCTGGGCGTTTTGGAggg                                                                                                      >  1:205692/1‑135 (MQ=255)
                cATTGGCCCGGGCCAACAGTTTCTCAAGATACGTCAGCTGCTCTTCATACACCTTCGGCTCCGGCTGCTTGCGCAAATACTGCACGCCGCGCTCGAAGGCCACGCGCGCCGGGGCGGGGTGGGCGGTTTGGAggg                                                                                                      >  1:151320/1‑135 (MQ=255)
                   tGGCCCGGGCCAACAGTTTCTCAAGATACGTCAGCTGCTCTTCATACACCTTCGGCTCCGGCTGCTTGCGCAAATACTGCACGCCGCGCTCGAAGGCCCAGCGCCCCCGGCCCGGCGGGGCGTTTTTGAggggcg                                                                                                   >  4:125249/1‑132 (MQ=255)
                     gCCCGGGCCAACAGTTTCTCAAGATACGTCAGCTGCTCTTCATACACCTTCGGCTCCGGCTGCTTGCGCAAATACTGCACGCCGCGCTCGAAGGCCAGGCGCGCCGGGCCGGGCGGGCCGTTTTTCAggggcggt                                                                                                 >  1:285392/1‑130 (MQ=255)
                                    ttCTCAATATCCGTCAGCTGCTATCAATACACCTGCGGCTCCGGCTGCTGGCCAAAATAGTCCCCGCCGCGCGCGACGCCCCGGCGCGCCGGGCCGGGCTGGCCTTTTTGCAGGGCCTGTTGGCCGAGgttgtt                                                                                   <  2:31209/134‑1 (MQ=255)
                                      ataaacaTACGTCAGCGGCTCTTCATACACCTTCGGCGCCGGCTGCTTGCGAAAAACCTGCACGCCGCGCTCGAAGGCCAGGCGCGCCTGGCCGGGCTGGCCTTTTTGCAGGGCCTGTTGGCCGAGGTTGTTga                                                                                 <  7:33870/128‑1 (MQ=255)
                                             tACGTCCGCTGCTATACATACACCTGCGGCTCCGGCTGCTGGCGAAAATACTCCCCGCCGCGCGCGAAGGCCCGGCGCGCCGGGCCGGGCTGGCCTTTTTGCAGGGCCTGTTGGCCGAGGTTGTTGAAGAActct                                                                         <  7:176389/135‑1 (MQ=255)
                                             cccGTCCGCGGCGCTACACACCCCTTCGGCGCCGGCTCCTTGCGCAAAACCTCCCCCCCGCGCGCCAACCCCAGGCGCGCGTGGCCGGGCTGGCCTTTTTGCAGGGCCTGTTGGCCGAGGTTGTTGAAGAActct                                                                         <  4:219973/133‑1 (MQ=255)
                                              ccGTCAGCTGCTCTTCATCCACCTGCGGCGCCGGCTGCTTGCTAAAATCCTCCACGCCGCGCTCGACGGCCAGGCGCGCCGGGCCGGGCTGGCCTTTTTGCAGGGCCTGTTGGCCGAGGTTGTTGAAGAACTCta                                                                        <  8:13489/134‑1 (MQ=255)
                                                  ccgctgatCTCCAAACACCTTCGGCTCCGGCGGCTTGCGCAAATACTCCACGCCGCGCGCTAAGGCCAGGCGCGCCGGGCCGGGCTGGCCTTTTTGCAGGGCCTGTTGGCCGAGGTTGTTGAAGAACTCTATATg                                                                    <  4:75554/133‑1 (MQ=255)
                                                  cAGCTCCTCTTAAAACACCTTGGGCGCCGGCGGCTTGCCAAAAACCTCCCCGCCGCGCTCTACGGCCTGGCGCGCGGGGCCGGGCTGGCCTTTTTGCAGGGCCTGTTGGCCGAGGTTGTTGAAGAACTCTATATg                                                                    <  2:76414/135‑1 (MQ=255)
                                                  cAGCGGCTCTACATACACCTTCGGCTCCGGCGGCGTGCGAAAATCCTCCACGCCGCGCTCCACGGCCCGGCGCGCCTGGCCGGGCTGGCCTTTTTGCAGGGCCTGTTGGCCGAGGTTGTTGAAGAACTCTATATg                                                                    <  2:296892/135‑1 (MQ=255)
                                                   ttgtgttctTCATAAACTTGGGGCTCGGGTTGTTTGCTCAAATACTGCACGCCGCGCTCTACGGCCAGGCGCGCGTGGCCGGGCTGGCCTTTTTGCAGGGCCTGTTGGCCGAGGTTGTTGAAGAACTCTATATg                                                                    <  8:347171/128‑1 (MQ=255)
                                                       gctATTCAACCACCTTCGGCTCCGGCTGCTTGCGAAAATCCTGCACGCCGCGCTCGACGGCCAGGCGCGCCGGGCCGGGCTGGCCTTTTTGCAGGGCCTGTTGGCCGAGGTTGTTGAAGAACTCTATATgcagca                                                               <  1:281989/135‑1 (MQ=255)
                                                        ctctTCATACACCTTCGGCTCCGGCTGCTTGCGCAAATACTGCACGCCGCGCTCGAAGGCCAGGCGCGCCTGGCCGGGCGGGGCTTTTTGCAGGGGCTGTTGGCCGTGGGTGTTGTAGTACTCTATATgcagccg                                                              >  5:102576/1‑133 (MQ=255)
                                                        ctctTCATACACCTTCGGCTCCGGCTGCTTGCGCAAATACTGCACGCCGCGCTCGAAGGCCACGCGCGCCAGGCCGGGCGGGCCTGTTTTGAGGGGCGGGTTGCCGCGGGTGGTGAAGAACTAGATATgaagcag                                                              >  7:347171/1‑135 (MQ=255)
                                                             cATACACCTTCGGCTCCGGCTGCTTGCGCAAATACTGCACGCCGCGCTCGAAGGCCCGGCGCGCCTGGCCGGGCTGGCCTTTTTGCAGGGCCTGTTGGGCGAGGTTGTTGAAGAACTCTATATGCAGCAGCACCa                                                         >  5:168241/1‑135 (MQ=255)
                                                                       cGGCTCCGGCTGCTTGCGAAAATACTCCCCGCCGCGCTCGACGGCCTGGCGCGCGGGGCCGGGCTGGCCTTTTTGCAGGGCCTGTTGGCCGAGGTTGTTGAAGAACTCTATATGCAGCAGCACCAGGATATGGc                                                <  4:122747/134‑1 (MQ=255)
                                                                        ggCTCCGGCGGCTGGCGAACATACTCCCCGCCGCGCTGCCGGCCCTGGCGGGCCTGGCCGGGCTGGCCTTTTTGCAGGGCCTGTTGGCCGAGGTTGTTGAAGAACTCTATATGCAGCAGCACCAGGATATGGCg                                               <  6:68333/134‑1 (MQ=255)
                                                                          gtCCGGCTGCTTGCGCAAATACTGCACGCCGCGCTCGAAGGCCCGGCGCGCCTGGCCGGGCTGGCCTTTTTGCAGGGCCGGTTGGCCGAGGTTGTTGAAGAACTCTATATTCAGCAGCCCCAGGATATGGCGGAt                                            >  8:154308/2‑135 (MQ=255)
                                                                          cTCCGGCTGCTTGCGCAAATACTGCACGCCGCGCTCGAAGGCCAGGCGCGCCTGGCCGGGGTGGCGTTGTTGCAGGGCCTGTTGGCCGAGGTTGTTGAAGAACTCTATATGCAGCAGCACCAGGTTATGGCGGAt                                            >  8:262434/1‑135 (MQ=255)
                                                                          cTCCGGCTGCTTGCGCAAATACTGCACGCCGCGCTCGAAGGCCAGGCGCGCCGGGCCGGGCTGGCCTTTTTGCAGGGCCGGTTGGCCGAGGTTGTTGAAGAACTCTATATTCAGCAGCACCAGGATATGGCGGAt                                            >  1:226719/1‑135 (MQ=255)
                                                                            cgggCTGGTTCCCCATATACGGCCCGCCGCGCGCGACGGCCCGGGGCGGCGGGCCGGGCTGGCCTTTTTGCAGGGCCTGTTGGCCGAGGTTGTTGAAGAACTCTATATGCAGCAGCACCAGGATATGGCGGAttt                                          <  1:192160/133‑1 (MQ=255)
                                                                                 tgctTGCGAAAAACCTCCCCGCCGGGCGCGACGGCCTGGCGGGCGGGGCCGGGCTGGCCTTTTTGCAGGGCCTGTTGGCCGAGGTTGTTGAAGAACTCTATATGCAGCAGCACCAGGATATGGCGGATTTCCTTg                                     <  2:285392/135‑1 (MQ=255)
                                                                                 gggttGCGCAAATACTGCCCGCCGCGCGCGACGGCCTGGCGCGCCTGGCCGGGCTGGCCTTTTTGCAGGGCCTGTTGGCCGAGGTTGTTGAAGAACTCTATATGCAGCAGCACCAGGATATGGCGGATTTCCTTg                                     <  3:327097/132‑1 (MQ=255)
                                                                                                   aCGCCGCGCTCGAAGGCCAGGCGCGCCTGGCCGGGCTGGCCTTTTTGCAGGGCCTGTTGGCCGAGGTTGTTGAAGAACTCTATATGCAGCAGCACCAGGATATGGCGGATTTCCTTGACCCAATGCTTGGCCTCg                   >  2:277942/1‑135 (MQ=255)
                                                                                                    cgccgcGCTCGACGCCCCGGCGCGCGGGGCCGGGCTGGCCTTTTTGCAGGGCCTGTTGGCCGAGGTTGTTGAAGAACTCTATATGCAGCAGCACCAGGATATGGCGGATTTCCTTGACCCAATGCTTGGCCTCGa                  <  2:301346/135‑1 (MQ=255)
                                                                                                    cgccgcGCTCGAAGGCCCGGCGCGCCTGGCCGGGCTGGCCTTTTTGCAGGGCCTGTTGGCCGAGGTTGTTGAAGAACTCTATATGCAGCAGCACCAGGATATGGCGGATTTCCTTGACCCAATGCTTGGCCTCg                   >  8:78749/1‑134 (MQ=255)
                                                                                                    cgccgcGCTCGAAGGCCAGGCGCGCCGGGCCGGGCTGGCCTTTTTGCAGGGCCTGTTGGCCGAGGTTGTTTAAGAACTCTATATGAAGCAGCACCAGGATATGGCGGATTTCCTTGACCCAATGCTTGGCCTCg                   >  1:250124/1‑134 (MQ=255)
                                                                                                     gccgcGCTCGAAGGCCAGGCGCGCCTGGCCGGGCTGGCCTTTTTGCAGGGCCTGTTGGCCGAGGTTGTTGAAGAACTCTATATGCAGCAGCACCAGGATATTGCGGATTTCCTTGCCCCAATGCTTGGCCTCg                   >  2:102643/1‑133 (MQ=255)
                                                                                                          gcTCGAAGGCCAGGCGCGCCTGGCCGGGCTGGCCTTTTTGCAGGGCCTGTTGGCCGAGGGTGTTGAAGAACTCTATATGCAGCAGCACCAGGATATGGCGGATTTCCTTGACCCAATGCTTGGCCTCGACAGCgg            >  7:194861/1‑135 (MQ=255)
                                                                                                          gcTCGAAGGCCAGGCGCGCCTGGCCGGGCTGGCCGTTTTGCAGGGCCTGTTGGCCGAGGTTGTTGAAGAACTCTATATGCAGCAGCACCAGGATATGGCGGATTTACTTGACCCAATGCTTGGCCTCGACAGCgg            >  1:301346/1‑135 (MQ=255)
                                                                                                            gcgcccgCCCGGCGCGCCGGGCCGGGCTGGCCTTTTTGCAGGGCCTGTTGGCCGAGGTTGTTGAAGAACTCTATATGCAGCAGCACCAGGATATGGCGGATTTCCTTGACCCAATGCTTGGCCTCGACAGCGGGc          <  5:46943/129‑1 (MQ=255)
                                                                                                               tcggCCAGGCGCGCGGGTCCGTGCTGGCCTTTTTGCAGGGCCTGTTGGCCGAGGTTGTTGAAGAACTCTATATGCAGCAGCACCAGGATATGGCGGATTTCCTTGAGCCAATGCTTGGCCTCGACAGCGGGCaag       <  4:292319/133‑1 (MQ=255)
                                                                                                               ccggCCCGGCGCGCGGGGCCGGGCTGGCCTTTTTGCAGGGCCTGTTGGCCGAGGTTGTTGAAGAACTCTATATGCAGCAGCACCAGGATATGGCGGATTTCCTTGACCCAATGCTTGGCCTCGACAGCGGGCaag       <  3:304452/133‑1 (MQ=255)
                                                                                                               caGGCCAGGCGCGCCGGGCCGGGCTGGCCTTTTTGCAGGGCCTGTTGGCCGAGGTTGTTGAAGAACTCTATATGCAGCAGCACCAGGATATGGCGGATTTCCTTGACCCAATGCTTGGCCTCGACAGCGGGCaag       <  7:32778/134‑1 (MQ=255)
                                                                                                               acggCCGGGCGGGCCTCCCCGGGCTGGCCTTTTTGCAGGGCCTGTTGGCCGAGGTTGTTGAAGAACTCTATATGCAGCAGCACCAGGATATGGCGGATTTCCTTGACCCAATGCTTGGCCTCGACAGCGGGCaag       <  3:125249/133‑1 (MQ=255)
                                                                                                               acggCCCGGCGCGCCGGGCCGGGCTGGCCTTTTTGCAGGGCCTGTTGGCCGAGGTTGTTGAAGAACTCTATATGCAGCAGCACCAGGATATGGCGGATTTCCTTGACCCAATGCTTGGCCTCGACAGCGGGCaag       <  5:299565/133‑1 (MQ=255)
                                                                                                                    ccggCGCGCCGGGCCGGGCTGGCCTTTTTGCAGGGCCTGTTGGCCGAGGTTGTTGAAGAACTCTATATGCAGCAGCACCAGGATATGGCGGATTTCCTTGACCCAATGCTTGGCCTCGACAGCGGGCAAGAAGcc  <  4:339049/133‑1 (MQ=255)
                                                                                                                    cAGGCGCGCCTGGCCGGGCTGGCCTTTTTGCAGGGCCTGTTGGCCGAGGTTGTTGAAGAACTCTATATGCAGCAGCACCAGGATATGGCGGATTTCCTTGACCCAATGCTTGGCCTCGACAGCGGGCAAGAAGc   >  4:220226/1‑134 (MQ=255)
                                                                                                                              |                                                                                                                            
CGGTCCAGTACCTGGGCATTGGCCCGGGCCAACAGTTTCTCAAGATACGTCAGCTGCTCTTCATACACCTTCGGCTCCGGCTGCTTGCGCAAATACTGCACGCCGCGCTCGAAGGCCAGGCGCGCCTGGCCGGGCTGGCCTTTTTGCAGGGCCTGTTGGCCGAGGTTGTTGAAGAACTCTATATGCAGCAGCACCAGGATATGGCGGATTTCCTTGACCCAATGCTTGGCCTCGACAGCGGGCAAGAAGCC  >  ACB143/2094691‑2094941

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