Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*ACB1223,440,3480TG57.1% ‑4.8 / 14.6 45A19A (GCT→GCGCJ019_03266hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (16/2);  new base G (16/8);  total (34/11)
Fisher's exact test for biased strand distribution p-value = 1.47e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.96e-01
Rejected as consensus: E-value score below prediction cutoff.
Rejected as consensus: Frequency below/above cutoff threshold.

ATCAATGGTGGTAGCAACGGTCTGGAGGATCGCAAGGCCCTCTACAAGCGAGCGCTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCAGCGCTGGCGGCGCCCGGGCAGCCTGGGTGTGGCAGGCCAACGCCTATAAGGCGCAGCTCGCTGATCAGGCCGATGACTTCGGCAAGCAGTTGGCGGAGAAGGATCGCGCGCACGGTCTGGAGCGG  >  ACB122/3440224‑3440468
                                                                                                                            |                                                                                                                        
ctccttgggggTACCCAGGGTCTGGCGCAGCGCCACTCCCTCTACAAGCGAGCGCGCGCGGGGCTGCCGGGCGGGCGTGGCCCCCCCGCGGGGGCGGCGGGGCGCGGGGGCTGGCCGCCAGCGCGGGCGGCGccc                                                                                                                <  4:447420/130‑1 (MQ=255)
    aTGGTGGTAGCAACGGTCTGGAGGATCGCAAGGCCCTCTACAAGCGAGCGCTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGGCGCCAGCGCGGGAGGCGCCCGGGc                                                                                                            >  6:7800/1‑135 (MQ=255)
    aTGGTGGTAGCAACGGTCTGGAGGATCGCAAGGCCCTCTACAAGCGAGCGCTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCCGCGCGGGCGGCGCCCgcgc                                                                                                            >  1:217384/1‑135 (MQ=255)
    aTGGTGGTAGCAACGGTCTGGAGGATCGCAAGGCCCTCTACAAGCGAGCGCTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCCCCGCTGGCGGCGCCCGGGc                                                                                                            >  2:17983/1‑135 (MQ=255)
    aTGGTGGTAGCAACGGTCTGGAGGATCGCAAGGCCCTCTACAAGCGAACGCTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCCGCGCGGGCGGCGCCCGGGc                                                                                                            >  7:424920/1‑135 (MQ=255)
             gCAACGGTCTGGAGGATCGCAAGGCCCTCTACAAGCGAGCGCTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCCGCGCTGGCGGCGCCCGGGCCGCCggggg                                                                                                   >  4:459883/1‑134 (MQ=255)
                 cGGTCTGGAGGATCGCAAGGCCCTCTACAAGCGAGCGCTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCCGCGCTGGCGGCGCCCGGGCCGCCGGGGTGTgg                                                                                               >  3:513905/1‑135 (MQ=255)
                    tCTGGAGGATCGCAAGGCCCTCTACAAGCGAGCGCTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCCGCGCTGGCGGCGCCCGGGCAGCCGGGGTGTGGCAg                                                                                            >  8:440628/1‑135 (MQ=255)
                        gaggaTCGCAAGGCCCTCTACAAGCGAGCGCTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCCCCGCTGGCGGCGCCCGGGCCGCCGGGGTGTGGCAGGccc                                                                                        >  3:177855/1‑134 (MQ=255)
                         aggaTCGCAAGGCCCTCTACAAGCGAGCGCTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCCCCAGCGCTGGCGGCGCCCGGGCAGCCTGGGTGTGGCAGGCCaa                                                                                       >  4:529901/1‑135 (MQ=255)
                                cAAGGCCCTCTACAAGCGAGCGCTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCTCTGCTGCTGGGCGCCCGCGCGGGCGGCGCCCCGGCGCCCGGGGGGTGGGGGGCCCCCGccga                                                                                 >  5:15236/1‑132 (MQ=255)
                                cAAGGCCCTCTACAAGCGAGCGCTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCCCCGCTCGCCGCGCCCGGGCCGCCCGGGGGTGGGAGGGCCACCCCtc                                                                                 >  6:164115/1‑133 (MQ=255)
                                 aaGGCCCTCTACAAGCGAGCGCTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGGCGCCCGCAATGGCGGCGCCCGGGCAGCCGCGGGGTGTGAGGCCACCGCCaata                                                                               >  7:453102/1‑135 (MQ=255)
                                 aaGGCCCTCTACAAGCGAGCGCTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCCGCGCTGGCGGCGCCCGGGCCGCCGGGGTGTGGCAGGCCACCGccgc                                                                                 >  3:193923/1‑131 (MQ=255)
                                 aaGGCCCTCTACAAGCGAGCGCTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCCGCGCTGGCGGCGCCCCGGCAGCCGGGGGGTGGGGGGCCCACGccga                                                                                 >  5:262843/1‑131 (MQ=255)
                                 aaGGCCCTCTACAAGCGAGCGCTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCCGCGCGGGCGGCGCCCCGGCGGCCGGGGTGGGGGAGGCCCACGccga                                                                                 >  7:372088/1‑131 (MQ=255)
                                     cccTCTACAAGCGAGCGCTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCCGCGCTGGCGGCGCCCGGGCAGCCTGGGTGTGGCAGGCCACCGCCTATAAGgc                                                                           >  3:469562/1‑135 (MQ=255)
                                        tctACAAGCGAGCGCTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCCGCGCGGGCGGCGCCCCGGCCGCCTGGGGGGGGCAGGCCACCGCCGATAAGGCGCa                                                                        >  5:239147/1‑135 (MQ=255)
                                          tACAAGCGAGCGCTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCCGCGCGGGCGGCGCCCGGGCCGCCGGGGGGTGGGAGGCCCCCCCCGATAAGGcgccgc                                                                      >  2:184406/1‑135 (MQ=255)
                                           aCAAGCGAGCGCTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCCGCCCAGGCCGCGCCCGGGCCGCCTGGGGGTGGCAGGCCCACGCCGAAAAGGcgccgcg                                                                     >  1:380902/1‑134 (MQ=255)
                                           aCAAGCGAGCGCTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCCCCGCCGGCGGCGCCCCGGCAGCCGGGGGGTGGGCGGCCCACCCCGATAACGCGCAGCt                                                                     >  6:104661/1‑135 (MQ=255)
                                               gcgagcgCGCGCGGGGCTGCCGTGCGCGCGGGGGCCCCCCGGCTGGGCGCCGGGGGGCTGCGGCGGGCCGCCCGCGCGGGCGGCGCCCGGGCAGCCTGGGTGTGGCAGGCCAACGCCTATAAGGCGCAGCTCGCt                                                                 <  3:208603/135‑1 (MQ=255)
                                                  agcgCTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGGCGCCCGCGCGGGCGGCGCCCCGGCCGCCGGGGGGTGGCCGGCCCCCCCCCATAAGTCGCCGCTCGCGGa                                                               >  6:238570/1‑134 (MQ=255)
                                                  agcgCTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCCGCGCGGGCGGCGCCCGGGCCGCCGGGGTGTGGCGGGCCCACGCCCATAAGGCGCACCgcgcggc                                                               >  3:225383/1‑131 (MQ=255)
                                                  agcgCTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCCCCGCGGGCGGCGCCCCGGCCGCCGGGGGGTGGCAGGCCCACCCCGATAAGGCGCCGCGCGCGGa                                                               >  8:103909/1‑134 (MQ=255)
                                                       tCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCCGCGCTGGCGGCGCCCGGGCCGCCTGGGTGTGGCAGGCCCCCGCCGCTAAGGCGCAGCGCGCTGCTCAGGc                                                         >  7:331411/1‑135 (MQ=255)
                                                       tCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCCGCGCTGGCGGCGCCCGGGCAGCCTGGGTGTGGCAGGCCAACGCCTATAAGGCGCAGCTCGCTGCTCAGGc                                                         >  7:362066/1‑135 (MQ=255)
                                                       tCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCCCCGCTGGCGGCGCCCCGGCCGCCTCGGTGGGGGCGGCCAACCCCGATAAGGCGCCGCGCGCGGATccgcg                                                         >  4:153969/1‑131 (MQ=255)
                                                       tCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCACCCCTGGCGGCGCCCCGCCGGCCCGGGTGGGGCCGGCCCACGCCCATAAACCGCACCACCCGGATCacaa                                                         >  6:310377/1‑132 (MQ=255)
                                                       tCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCGGCTGGCCGCCCGCGCGGGCGGCGCCCGGGCCGCCGGGGTGGGGGCGGCCCACGCCGAAAAGGCGCCGCGCGCGGATCCGGc                                                         >  7:289145/1‑135 (MQ=255)
                                                         gAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGGCGCCCGCGCGGGCGGCGGCCGGGCCGCCGGGGTGTGGGGGGCCCACGCCGATAAGGCGCCGCGCGCGGCTCAGGCCg                                                       >  1:444003/1‑135 (MQ=255)
                                                                  cAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCCGCCCGGGCGGCCCCCGGGCCGCCGGGGGGTGGGCGGCCACCCCCGAAAACGCGCACCGCGCTGCTCAGGCCGATGACATCg                                              >  8:313013/1‑135 (MQ=255)
                                                                          gcTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGGCGCCCGCCCGGGCGGCGCCCGGGGCGCCGGGGGGTGGGCGGCCCCCCCCGCTAAAGCGCCGCGCGCGGGTCAGGCCGCGGACTTCGGGcagcag                                      >  6:160476/1‑135 (MQ=255)
                                                                          gcTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCCGCGCGGGCGGCGCCCGGGCCGCCGGGGGGTGGGAGGCCCCCGCCGATAAGGCGCCGCGCGCGGATCAGGCCGCTGACGGCGGCGAGCAg                                      >  7:411501/1‑135 (MQ=255)
                                                                           cgtGGGCCCCCCGCCGGGGCGCCGGGGCGCTGCGGCGGGCCGCCCGCGCTGGCGGCGCCCGGGCAGCCTGGGTGTGGCAGGCCAACGCCTATAAGGCGCAGCTCGCTGATCAGGCCGATGACTTCGGCAAGCAGt                                     <  7:69764/133‑1 (MQ=255)
                                                                                              ggcggtggcggtggcgCCGCCCGCCCGCGCGGGCGGCGCCCGGGCAGCCTGGGTGTGGCAGGCCAACGCCTATAAGGCGCAGCTCGCTGATCAGGCCGATGACTTCGGCAAGCAGTTGGCGGAGAAGGATcgcgc                  <  5:7800/131‑1 (MQ=255)
                                                                                              cgcgggggCGCGGCGGCGGGCCGCGCGCGCGGGCGGCGCCCGGGCAGCCTGGGTGTGGCAGGCCAACGCCTATAAGGCGCAGCTCGCTGATCAGGCCGATGACTTCGGCAAGCAGTTGGCGGAGAAGGATcgcgc                  <  5:200640/129‑1 (MQ=255)
                                                                                               gCCGTGGCGCTGCTGCTGGGCGCCCGCGCGGGGGGCGCCCGGGCCGCCGGGGGGGGGCCGGGCCCCGCCGCTAAGGCGCCGCTCGCGGATCAGGCCGCTGACTTCGGCGCGCCGGTGGGGGCGAAGGCTcgcgcg                 >  8:285610/1‑135 (MQ=255)
                                                                                                gggTGGGGCTGCGGCGGGCCGCCAGCGCGGCCGGCGCCCGGGCAGCCTGGGTGTGGCAGGCCAACGCCTATAAGGCGCAGCTCGCTGATCAGGCCGATGACTTCGGCAAGCAGTTGGCGGAGAAGGATCgcgcgc                <  1:76311/133‑1 (MQ=255)
                                                                                                ccGTGGCGCTGCTGCTGGCCGCCCGCGCTGGCGGCGCCCGGGGCGCCGGGGTGTGGGCGGCCAACGCCGATAAGGCGCCGCTCGCGGGTCAGGCCGATGCCGTGGGCGAGCAGTTGGGGGAGAAGGCTCgcgcgc                >  6:482279/1‑135 (MQ=255)
                                                                                                        ctgctgctgGCCGCCAGGGCTGGCGGCGCCCGGGCAGCCTGGGTGTGGCAGGCCAACGCCTATAAGGCGCAGCTCGCTGATCAGGCCGATGACTTCGGCAAGCAGTTGGCGGAGAAGGATCGCGCGCACGGTCTg        <  8:98568/135‑1 (MQ=255)
                                                                                                             gctgGCCGCCCGCGCCGCCGGCGCCCGGGCAGCCTGGGTGTGGCAGGCCAACGCCTATAAGGCGCAGCTCGCTGATCAGGCCGATGACTTCGGCAAGCAGTTGGCGGAGAAGGATCGCGCGCACGGTCTGGAGCg   <  1:204726/135‑1 (MQ=255)
                                                                                                              gtgGCCGCCCGCGCGGGCCGCGCCCGGGCAGCCTGGGTGTGGCAGGCCAACGCCTATAAGGCGCAGCTCGCTGATCAGGCCGATGACTTCGGCAAGCAGTTGGCGGAGAAGGATCGCGCGCACGGTCTGGAGCgg  <  5:104661/134‑1 (MQ=255)
                                                                                                              ggggCCGCCGGCGCGGGCCGCGCCCGGGCAGCCTGGGTGTGGCAGGCCAACGCCTATAAGGCGCAGCTCGCTGATCAGGCCGATGACTTCGGCAAGCAGTTGGCGGAGAAGGATCGCGCGCACGGTCTGGAGCgg  <  2:89471/133‑1 (MQ=255)
                                                                                                              cgggCCGCCCGGGGGGGCGGCGCCCGGGCAGCCTGGGTGTGGCAGGCCAACGCCTATAAGGCGCAGCTCGCTGATCAGGCCGATGACTTCGGCAAGCAGTTGGCGGAGAAGGATCGCGCGCACGGTCTGGAGCgg  <  7:103909/133‑1 (MQ=255)
                                                                                                                            |                                                                                                                        
ATCAATGGTGGTAGCAACGGTCTGGAGGATCGCAAGGCCCTCTACAAGCGAGCGCTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCAGCGCTGGCGGCGCCCGGGCAGCCTGGGTGTGGCAGGCCAACGCCTATAAGGCGCAGCTCGCTGATCAGGCCGATGACTTCGGCAAGCAGTTGGCGGAGAAGGATCGCGCGCACGGTCTGGAGCGG  >  ACB122/3440224‑3440468

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: