Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*ACB1436,087,1260AC53.3% 1.4 / 11.8 15intergenic (+7/+25)spuD_6/gabR_2Putrescine‑binding periplasmic protein SpuD/HTH‑type transcriptional regulatory proteinGabR
Reads supporting (aligned to +/- strand):  ref base A (3/4);  new base C (0/8);  total (3/12)
Fisher's exact test for biased strand distribution p-value = 7.69e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.48e-02
Rejected as consensus: E-value score below prediction cutoff.
Rejected as consensus: Frequency below/above cutoff threshold.

TTCGCAACAACCCCAACCTGTATCCGACGGCAGAGGCGATGGCCAAGCTGTATACCCTCAAGCCTTTGACGCGAAATGCCGAACGGGCCAGGACCCGCGCCTGGACGCGGATAAAGTCCGGTACGTGATGTTGAACCGGGGGCTGCTTTGCAGCCCCCTTCATTTTCTCCAGACTTTGCGTAATTTCATCAAGGCTTCGGCTATCTCGCCCTCCGGCACCGCCGCGAAGCCCAGCACCAGCCCC  >  ACB143/6086992‑6087235
                                                                                                                                      |                                                                                                             
ttCGCAACAACCCCAACCTGTATCCGACGGCAGAGGCGATGGCCAAGCTGTATACCCTCAAGCCTTTGACGCGAAATGCCGAACGGGCCAGGACCCGCGCCTGGACGCGGATAAAGTCCGGTACGTGATGTTGAAcc                                                                                                             >  8:24847/1‑137 (MQ=255)
  cGCAACAACCCCAACCTGTATCCGACGGCAGAGGCGATGGCCAAGCTGTATACCCTCAAGCCTTTGACGCGAAATGCCGAACGGGCCAGGACCCGCGCCTGGACGCGGATAAAGTCCGGTACGTGATGTTGAACCgg                                                                                                           <  1:200615/137‑1 (MQ=255)
  cGCAACAACCCCAACCTGTATCCGACGGCAGAGGCGATGGCCAAGCTGTATACCCTCAAGCCTTTGACGCGAAATGCCGAACGGGCCAGGACCCGCGCCTGGACGCGGATAAAGTCCGGTACGTGATGTTGAACCgg                                                                                                           <  3:54198/137‑1 (MQ=255)
                                 agGCGATGGCCAAGCTGTATACCCTCAAGCCTTTGACGCGAAATGCCGAACGGGCCAGGACCCGCGCCTGGACGCGGATAAAGTCCGGTACGTGATGTTGAACCGGGGGCGGCGTTGGCGCCCCCCTCCttttttgc                                                                            >  8:69334/1‑133 (MQ=255)
                                                                    aCGCTAAATGCTGAAGGGGCCCGCCCCCGCGGCTGGACGCTTTAAACGCCCGGTCCGTGAGGTTTACGCGGGGGCTGCTTTGCAGCCCCCTTCATTTTCTCCTGACTTTGCGTAATTTCATCAAGGCTTCGGCTAtc                                         <  2:181393/137‑1 (MQ=255)
                                                                               ctaccGGGCCAGCCCCCGCGCCGGGACGCGTATAACGACGGGAACGTGATGTTCACCGGGGGGCGGCTTGCCAGCCCCCTTCATTTTCTCCAGACTTTGCGTAATTTCATCAAGGCTTCGGCTATCTCGCCCTCCgg                              <  5:168101/133‑1 (MQ=255)
                                                                                     gcccccgACCCGCGCCGGGCCGCGTAAACAGGCCGTGTCGGCATTTTCACCCGGGGGCGGCTTTCCACCCCCCTTCATTTTCTCCAGACTTTGCGTAATTTCATCAAGGCTTCGGCTATCTCGCCCTCCGGCAccgc                        <  2:174761/131‑1 (MQ=255)
                                                                                      gCCAGCACCCGCGCCGGGTCGCGTATAAAGGCCGTAACGTTATGTTCACCCGGGGGTTGCTTGCCACCCCCCTTCATTTTCTCCAGACTTTGCGTAATTTCATCAAGGCTTCGGCTATCTCGCCCTCCGGCAccgcc                       <  5:112608/137‑1 (MQ=255)
                                                                                      ccaaGGCCCCGCGCCTCGCCTCGTATAACGCCAGGTACTTAATGTGCCACCGGGGGCTGCTTTCACCCCCCCTTCATTTTCTCCAGACTTTGCGTAATTTCATCAAGGCGTCGGCTATCTCGCCCTCCGGCAccgcc                       <  4:202739/134‑1 (MQ=255)
                                                                                           cccccGCGCCGGCACGCGTATAAAGCGCGTTCGGTGATTTTGACCGGGGGGCTGCTTTGAAGCCCCCTTCATTTTCTCCAGACTTTGCGTAATTTCATCAAGGCTTCGGCTATCTCGCCCTCCGGCACCGCCGCGaa                  <  8:160892/135‑1 (MQ=255)
                                                                                            ccccGCGCGCGGCCGCGTAAAAAGGCCGTTAGGTGTTGTTCACCGGGGGGCTGCTTTCCCCCCCCCTTCATTTTCTCCAGACTTTGCGTAATTTCATCAAGGCTTCGGCTATCTCGCCCTCCGGCACCGCCGCGAAg                 <  8:10880/136‑1 (MQ=255)
                                                                                                gcgcCTGGCCGCGTAAAAACTCCGGTACGTCTTGTTGAACCGGGGGCTGCTTTCCAGCCCCCTTCATTTTCTCCAGACTTTGCGTAATTTCATCAAGGCTTCGGCTATCTCGCCCTCCGGCACCGCCGCGAAGCCCa             <  2:265350/137‑1 (MQ=255)
                                                                                                   ccTGGACGCGGATAAAGTCCGGTACGTGATGTTGAACCGGGGGCTGCTTTGCAACCCCCTTCATTTTCTCCAGACTTTGCGTAATTTCATCAAGGCTTCGGCTATCTCGCCCTCCGGCACCGCCGCGACGCCCAGcc          >  5:335293/1‑136 (MQ=255)
                                                                                                           gggATAAAGTCGGGTAGGTGATTTTGACCCGGGGGCTGCTTTGCAGCCCCCTTCATTTTCTCCAGACTTTGCGTAATTTCATCAAGGCTTCGGCTATCTCGCCCTCCGGCACCGCCGCGAAGCCCAGCACCAGcccc  <  8:123595/136‑1 (MQ=255)
                                                                                                           cgcaTAACGGCCGTTAGGTTATGTTCACCCGGGGGCGGGTTGCCCGCCCCCTTCATTTTCTCCAGACTTTGCGTAATTTCATCAAGGCTTCGGCTATCTCGCCCTCCGGCACCGCCGCGAAGCCCAGCACCAGcccc  <  7:24847/134‑1 (MQ=255)
                                                                                                                                      |                                                                                                             
TTCGCAACAACCCCAACCTGTATCCGACGGCAGAGGCGATGGCCAAGCTGTATACCCTCAAGCCTTTGACGCGAAATGCCGAACGGGCCAGGACCCGCGCCTGGACGCGGATAAAGTCCGGTACGTGATGTTGAACCGGGGGCTGCTTTGCAGCCCCCTTCATTTTCTCCAGACTTTGCGTAATTTCATCAAGGCTTCGGCTATCTCGCCCTCCGGCACCGCCGCGAAGCCCAGCACCAGCCCC  >  ACB143/6086992‑6087235

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: