Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*ACB1433,440,3460AC52.9% 11.3 / 28.4 34A24A (GCA→GCCCJ019_03266hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (8/8);  new base C (18/0);  total (26/8)
Fisher's exact test for biased strand distribution p-value = 7.09e-04
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.45e-02
Rejected as consensus: Frequency below/above cutoff threshold.

ATGGTGGTAGCAACGGTCTGGAGGATCGCAAGGCCCTCTACAAGCGAGCGCTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCAGCGCTGGCGGCGCCCGGGCAGCCTGGGTGTGGCAGGCCAACGCCTATAAGGCGCAGCTCGCTGATCAGGCCGATGACTTCGGCAAGCAGTTGGCGGAGAAGGATCGCGCGCACGGTCTGGAGCGGGAGAAGGCTGCAGCAGCAGCGCTCGATCA  >  ACB143/3440211‑3440480
                                                                                                                                       |                                                                                                                                       
aTGGTGGTAGCAACGGTCTGGAGGATCGCAAGGCCCTCTACAAGCGAGCGCTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCAGCGCGGGCGGCGCCCGGGCAg                                                                                                                                        >  8:283833/1‑137 (MQ=255)
      gtAGCAACGGTCTGGAGGATCGCAAGGCCCTCTACAAGCGAGCGCTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCAGCGCTGGCGGCGCCCGGGCCGCCTggg                                                                                                                                  >  5:28165/1‑137 (MQ=255)
       tAGCAACGGTCTGGAGGATCGCAAGGCCCTCTACAAGCGAGCGCTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCAGCGCTGGCGGCGCCCGGGCAGCCTgggg                                                                                                                                 >  6:124000/1‑136 (MQ=255)
       tAGCAACGGTCTGGAGGATCGCAAGGCCCTCTACAAGCGAGCGCTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCAGCGCTGGCGGCGCCCGGGCACCCggggg                                                                                                                                 >  2:51532/1‑136 (MQ=255)
              ggTCTGGAGGATCGCAAGGCCCTCTACAAGCGAGCGCTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCACCGCTGGCGGCGCCCGGGCCGCCGGGGTGTGGccg                                                                                                                          >  1:221203/1‑135 (MQ=255)
                     aggaTCGCAAGGCCCTCTACAAGCGAGCGCTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCAGCGCTGGCGGCGCCCGGGCAGCCTGGGTGTGGCAGGCCAAcc                                                                                                                   >  5:240260/1‑136 (MQ=255)
                        aTCGCAAGGCCCTCTACAAGCGAGCGCTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCCGCGCTGGCGGCGCCCGGGCAGCCTGGGTGTGGCAGGCCAACGcca                                                                                                                >  6:243026/1‑136 (MQ=255)
                            cAAGGCCCTCTACAAGCGAGCGCTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCGGCCCGCGCTGGCGGCGCCCGGGCCGCCGGGGTGTGGGGGGCCCACGccgagaa                                                                                                            >  1:296640/1‑132 (MQ=255)
                                       aCAAGCGAGCGCTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCCGCGCGGGCGGCGCCCGGGCCGCCGGGGGGTGGCAGGCCCACGCCGATAAGGCGCAGCTCg                                                                                                 >  1:128213/1‑137 (MQ=255)
                                         aaGCGAGCGCTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCCCCCCTGGCGGCGCCCCGGCACCCGGGGGGTGGCAGGCCCACCCCCAAAAGGCGCAGCTCGcg                                                                                               >  3:47856/1‑136 (MQ=255)
                                         aaGCGAGCGCTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCCCCAGCGCTGGCGGCGCCCGGGCCGCCGGGGGGGGGGCGGCCACCGCCAATAAGGACCAGCTCGcg                                                                                               >  3:47857/1‑136 (MQ=255)
                                                 gcTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCCGCGCTGGGGGCGCCCCGGCCGCCGGGGTGTGGCGGGCCACCCCCTATAATGAGCCGCTCGCGGCTAAGGc                                                                                      >  3:48136/1‑137 (MQ=255)
                                                 gcTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCCGCGCTGGCGGCGCCCGGGCCGCCGGGGTGTGGCAGGCCAACGCCGATAAGGCGCAGCTCGCTGATCAGGc                                                                                      >  7:293935/1‑137 (MQ=255)
                                                 gcTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCCGCGCGGGCGGCGCCCCGGCCGCCTGGGTGTGGCAGGCCAACGCCTATAAGGCGCAGCTCGCGGCTCAGGc                                                                                      >  2:96065/1‑137 (MQ=255)
                                                 gcTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCACCGCTGGCGGCGCCCCGGCAGCCTGGGTGGGGGAGGCCACCGCCTATAAGGCGCAGATCCCCGATAAGGc                                                                                      >  6:115075/1‑137 (MQ=255)
                                                 gcTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCCCCCGCGCTGGCGGCGCCCGGGCCGCCTGGGTGTGGGAGGCCAACGCCTATAAGGCGCCGCTCGCGGATCAGGc                                                                                      >  4:206946/1‑137 (MQ=255)
                                                   tCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCCGCGCTGGCGGCGCCCGGGCCGCCGGGGTGTGGGCGGCCAACGCCGATAAGGCGCAGCTCGCTGATCAGGCCg                                                                                    >  7:30078/1‑137 (MQ=255)
                                                   tCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCCGCGCGGGCGGCGCCCGGGCCGCGGGGGTGTGGGAGGCCCCCGCCGAGAAGGCGCCGCTCGCGGAGCAGGccc                                                                                    >  6:282849/1‑136 (MQ=255)
                                                   tCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCAGCGCGGGGGGCGCCCGGGCCGCCGGGGTGTGGGCGGCCCACGCCTATAAGGCGCAGCGCGCTGATCAGGCCg                                                                                    >  4:226642/1‑137 (MQ=255)
                                                     gAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCCGCGCGGGCGGGGCCCGGGGCGCCGGGGGGTGGGCGGCCCACGCCGATAAGGCGCCGCTCGCGGCTGAGGCCgcg                                                                                  >  1:244495/1‑135 (MQ=255)
                                                                      gcTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCCGCGCTGGCGGCGCCCGGGCCGCCGGGGGGTGGGCGGCCACCGCCGATAAGGCGCCGCGCGCGGATCAGGCCGATGACTTCGGCAAGCAGtt                                                                 >  5:283085/1‑137 (MQ=255)
                                                                      gcTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCCGCGCTGGCGGCGCCCCGGCCGCCTGGGTGTGGCAGGCCCACGCCTATAAGGCGCCGCTCGCGGATCAGGCCGCTGACGTCGGCCAGCAGtt                                                                 >  8:170963/1‑137 (MQ=255)
                                                                      gcTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCCGCGCGGGCGGCGCCCGGGCAGCCGGGGGGTGGGCCGCCCACCCCTATAAGGCGCCGCACGCGGCTCAGGCCCGATCCGTCGGCAAGCAggt                                                                 >  8:271464/1‑135 (MQ=255)
                                                                      gcTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCCCCGCGGGCGGCGCCCGGGCCCCCGGGGTGTGGGAGGCCAACGCCGATAAAGCGCAGCGCGCTGCTGCGGCCGATGACTTGGGCCAGCAGtt                                                                 >  2:95807/1‑137 (MQ=255)
                                                                           ggCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCCGCGCGGGCGGCGCCCCGGCCGCCGGGGGGTGGGAGGGCCACCCCGATAAGGCGCAGCGCGCGGGTCAGGCCGCTGCCTTGGGCCAGCAGGTGgggg                                                            >  8:278861/1‑137 (MQ=255)
                                                                                cccGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCCGCGCTGGCGGCGCCCGGGCCGCCGGGGTGTGGCAGGCCCACGCCGATAAGGCGCAGCGCGCGGCGCAGGCCGAGGACGTCGTCCAGCAGGTGGCGGAGaaa                                                       >  8:160760/1‑136 (MQ=255)
                                                                                        gggggggTGGCGCTGGGGCGGCCCGCCCGGGGGGGCGGCGCCCGGGCAGCCTGGGTGTGGCAGGCCAACGCCTATAAGGCGCAGCTCGCTGATCAGGCCGATGACTTCGGCAAGCAGTTGGCGGAGAAGGATcgcgc                                               <  5:243026/131‑1 (MQ=255)
                                                                                                          ggggCCGCCTGGGGGGGCGGCGCCCGGGCAGCCTGGGTGTGGCAGGCCAACGCCTATAAGGCGCAGCTCTCTGATCAGGCCGATGACTTCGGCAAGCAGTTGGCGGAGAAGGATCGCGCGCACGGTCTGGAGCGGga                             <  5:207121/135‑1 (MQ=255)
                                                                                                          cgggCCGCCCGCGCGGGCCGCGCCCGGGCAGCCTGGGTGTGGCAGGCCAACGCCTATAAGGCGCAGCTCGCTGATCAGGCCGATGACTTCGGCAAGCAGTTGGCGGAGAAGGATCACGCGCACGGTCTGGAGCGGga                             <  8:30078/135‑1 (MQ=255)
                                                                                                                     ggggggcggcGCCCGGGCAGCCTGGGTGTGGCAGGCCAACGCCTATAAGGCGCAGCTCGCTGATCAGGCCGATGACTTCGGCAAGCAGTTGGCGGAGAAGGATCGCGCGCACGGTCTGGAGCGGGAGAAGGCTgcag                  <  6:283085/133‑1 (MQ=255)
                                                                                                                         ggcgccGCGCGGGCAGGCTGGGTGTGGCAGGCCAACGCCTATAAGGCGCAGCTCGCTGATCAGGCCGATGACTTCGGCAAGCAGTTGGCGGAGAAGGATCGCGCGCACGGTCTGGAGCGGGAGAAGGCTgcagcagc              <  3:38093/137‑1 (MQ=255)
                                                                                                                              cgcCCGGGCAGCCTGGGTGTGGCAGGCCAACGCCTATAAGGCGCAGCTCGCTGATCAGGCCGATGACTTCGGCAAGCAGTTGGCGGAGAAGGATCGCGCGCACGGTCTGGAGCGGGAGAAGGCTGCAGCAGCAgcgc         <  2:304220/137‑1 (MQ=255)
                                                                                                                                     gCAGCCTGGGTGTGGCAGGCCAACGCCTATAAGGCGCAGCTCGCTGATCAGGCCTATGACTTCGGCAAGCAGTTGGCGGAGAAGGATCGCGCGCACGGTCTTGAGCGGGAGAAGGCTTCAGCAGCAGCGCTCGATCa  <  1:67988/137‑1 (MQ=255)
                                                                                                                                     gCAGCCTGGGTGTGGCAGGCCAACGCCTATAAGGCGCAGCTCGCTGATCAGGCCGATGACTTCGGCAAGCAGTTGGCGGAGAAGGATCGCGCGCACGGTCTGGAGCGGGAGAAGGCTGCAGCAGCAGCGCTCGATCa  <  7:170963/137‑1 (MQ=255)
                                                                                                                                       |                                                                                                                                       
ATGGTGGTAGCAACGGTCTGGAGGATCGCAAGGCCCTCTACAAGCGAGCGCTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCAGCGCTGGCGGCGCCCGGGCAGCCTGGGTGTGGCAGGCCAACGCCTATAAGGCGCAGCTCGCTGATCAGGCCGATGACTTCGGCAAGCAGTTGGCGGAGAAGGATCGCGCGCACGGTCTGGAGCGGGAGAAGGCTGCAGCAGCAGCGCTCGATCA  >  ACB143/3440211‑3440480

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: