Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,376,735	0	T	.	60.9%	3.8 / 11.6	24	intergenic (+36/‑140)	ecnB/gdx	bacteriolytic entericidin B lipoprotein/guanidinium exporter
Reads supporting (aligned to +/- strand):  ref base T (5/4);  new base . (12/2);  total (18/6)
Fisher's exact test for biased strand distribution p-value = 1.62e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.67e-01
Rejected as consensus: E-value score below prediction cutoff.
Rejected as consensus: Frequency below/above cutoff threshold.

TTTTTGGCGATGTTGTCTATTATTAATTTGCTATAGGCAAACATAAATAACATTACCT‑AAAAGGAAGACGTTATGGTGAAGAAGACAATTGCAGCGATCTTTTC‑TGTTCTGGTGC‑TTTCAACAGTATTAACTGCCTGCAACACCACGCGTGGCG‑TTGGTGAAGACA‑TTTCTG‑ATGGCGGTAACGCGATTTCTGGCGCAGCAACG‑AAAGCGCAGCAATAAGCAATAACGGTACGACAGCTGTGTCGTGCCG‑TTTG‑TTTTTTCTG‑CGATAGTCACAAAGGTAATAGTTGAAATT‑‑CCCCTGCCACCTGGCAAAATATCCGTTC‑AACCATCA‑GC‑TTTGCAGGACGACCTGCAAACG‑CCTC‑TTTTCA‑CC‑GGGGAC‑GG‑‑CCCCAA‑TT‑‑CTCCGGAGCCTGAT‑ATGTCCTGGATTATCTTAGTTATTGCTGGTCTGCTGGAAGTGGTATGGGCCG‑TTGGCCTGAAATATACCCACGGCTTTAGTCG‑TTTGACGCCGAGTGTTATTACTGTGACGGCGATGATTGTCAGTATGGC  >  NC_000913/4376481‑4377005                                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                                                                  tttttGGCGATGTTGTCTATTATTAATTTGCTATAGGCAAACATAAATAACATTACCTAAAAAGGAAGACGTTATGGTGAAG‑AGACAATTGCAGCGATCTTTTCTTGTTCTGGTGCTTTTCAACAGTATTAACTGCCTGCAACACCACGCGTGGCGTTTGGTGAAGACATTTTCTG‑ATGGCGGTAACGCGATTTCTGGCGCAGCAACGAAAAGCGCAGCAATAAGCAATAACGGTACGACAGCTGTGTCGTGCCG‑TTTG‑TTTTTTCTG‑CGATAGTCACaaa                                                                                                                                                                                                                                                                             >  1:652193/1‑281 (MQ=255)                    ttattaATTTGCTATAGGCAAACATAAATAACATTACCTAAAAAGGAAGACGTTATGGTGAAGAAGACAATTGCAGCGATCTTTTCTTGTTCTGGTGC‑TTTCAACAGTATTAACTGCCTGCAACACCACGCGTGGCG‑TTGGTGAAGACA‑TTTCTG‑ATGGCGGTAACGCGATTTCTGGCGCAGCAACG‑AAAGCGCAGCAATAAGCAATAACGGTACGACAGCTGTGTCGTGCCG‑TTTG‑‑TTTTTCTGTCGATAGTCACAAAGGTAATAGTTGAAATT‑‑CCCCTGCCa                                                                                                                                                                                                                                                  >  1:652197/1‑284 (MQ=255)                                                taaCATTACCT‑AAAAGGAAGACGTTATGGTGAAGAAGACAATTGCAGCGATCTTTTC‑TGTTCTGGTGC‑TTTCAACAGTATTAACTGCCTGCAACACCACGCGTGGCG‑TTGGTGAAGACA‑TTTCTG‑ATGGCGGTAACGCGATTTCTGGCGCAGCAACG‑AAAGCGCAGCAATAAGCAATAACGGTACGACAGCTGTGTCGTGCCG‑TTTG‑tttt                                                                                                                                                                                                                                                                                                >  1:1580304/1‑211 (MQ=255)                                                            aaaaGGAAGACGTTATGGTGAAGAAGACAATTGCAGCGATCTTTTC‑TGTTCTGGTGC‑TTTCAACAGTATTAACTGCCTGCAACACCACGCGTGGCG‑TTGGTGAAGACA‑TTTCTG‑ATGGCGGTAACGCGATTTCTGGCGCAGCAACG‑AAAGCGCAGCAATAAGCAATAACGGTACGACAGCTGTGTCGTGCCG‑TTTG‑TTTTTTCTG‑CGATAGTCACAAAGGTAATAGTTGAAATT‑‑CCCCTGCCACCTGGCAAAATATCCGTTC‑AACCATCA‑GC‑TTTGCAGGACGACCTGCAAACG‑cc                                                                                                                                                                                         >  1:1580307/1‑296 (MQ=255)                                                                         tATGGTGAAGAAGACAATTGCAGCGATCTTTTC‑TGTTCTGGTGC‑TTTCAACAGTATTAACTGCCTGCAACACCACGCGTGGCG‑TTGGTGAAGACA‑TTTCTG‑ATGGCGGTAACGCGATTTCTGGCGCAGCAACG‑AAAGCGCAGCAATAAGCAATAACGGTACGACAGCTGTGTCGTGCCGTTTTG‑TTTTTTCTG‑CGATAGTCACAAAGGTAATAGTTGAAATT‑‑CCCCTGCCACCTGGCAAAATATCCGTTC‑AACCATCA‑GC‑TTTGCAGGACGACCTGCAAACG‑CCTC‑TTTTCA‑CC‑GGGGa                                                                                                                                                                       >  1:1580309/1‑299 (MQ=255)                                                                                    agaCAATTGCAGCGATCTTTTC‑TGTTCTGGTGC‑TTTCAACAGTATTAACTGCCTGCAACACCACGCGTGGCG‑TTGGTGAAGACA‑TTTCTG‑ATGGCGGTAACGCGATTTCTGGCGCAGCAACG‑AAAGCGCAGCAATAAGCAATAACGGTACGACAGCTGTGTCGTGCCG‑TTTG‑‑TTTTTCTG‑CGATAGTCACAAAGGTAATAGTTGAAATT‑‑CCCCTGCCACCTGGC‑AAATATCCGTTC‑AACCATCA‑GC‑TTTGCAGGACGACCTGCAAACG‑CCTC‑TTTTCA‑CC‑GGGGAC‑GG‑‑c                                                                                                                                                                >  1:1580311/1‑289 (MQ=255)                                                                                      aCAATTGCAGCGATCTTTTC‑TGTTCTGGTGC‑TTTCAACAGTATTAACTGCCTGCAACACCACGCGT‑GCG‑TTGGTGAAGACA‑TTTCTG‑ATGGCGGTAACGCGATTTCTGGCGCAGCAACG‑AAAGCGCAGCAATAAGCAATAACGGTACGACAGCTGTGTCGTGCCG‑TTTG‑TTTTTTCTG‑CGATAGTCACAAAGGTAATAGTTGAAATT‑‑CCCCTGCCACCTGGCAAAATATCCGTTC‑AACCATCA‑GC‑TTTg                                                                                                                                                                                                              <  1:1580308/249‑1 (MQ=255)                                                                                        aaTTGCAGCGATCTTTTC‑TGTTCTGGTGC‑TTTCAACAGTATTAACTGCCTGCAACACCACGCGTGGCG‑TTGGTGAAGACA‑TTTCTG‑ATGGCGGTAACGCGATTTCTGGCGCAGCAACG‑AAAGCGCAGCAATAAGCAATAACGGTACGACAGCTGTGTCGTGCCG‑TTTG‑‑TTTTTCTGTCGATAGTCACAAAGGTAATAGTTGAAATT‑‑CCCCTGCCACCTGGCAAAata                                                                                                                                                                                                                                      >  1:652206/1‑227 (MQ=255)                                                                                                    cTTTTC‑TGTTCTGGTGC‑TTTCAACAGTATTAACTGCCTGCAACACCACGCGTGGCG‑TTGGTGAAGACA‑TTTCTG‑ATGGCGGTAACGCGATTTCTGGCGCAGCAACG‑AAAGCGCAGCAATAAGCAATAACGGTACGACAGCTGTGTCGTGCCG‑TTTG‑‑TTTTTCTG‑CGATAGTCACAAAGGTAATAGTTGAAATT‑‑CCCCTGCCACCTGGCa                                                                                                                                                                                                                                           >  1:1580314/1‑209 (MQ=255)                                                                                                            gttctgGTGC‑TTTCAACAGTATTAACTGCCTGCAACACCACGCGT‑GCG‑TTGGTGAAGACA‑TTTCTG‑ATGGCGGTAACGCGATTTCTGGCGCAGCAACG‑AAAGCGCAGCAATAAGCAATAACGGTACGACAGCTGTGTCGTGCCG‑TTTG‑TTTTTTCTG‑CGATAGTCACAAAGGTAATAGTTGAAATT‑‑CCCCTGCCACCTGGCAAAATATCCGTTC‑AACCATCA‑GC‑TTTGCAGgacg                                                                                                                                                                                                       <  1:1580310/235‑1 (MQ=255)                                                                                                                               gTATTAACTGCCTGCAACACCACGCGTGGCG‑TTGGTGAAGACA‑TTTCTG‑ATGGCGGTAACGCGATTTCTGGCGCAGCAACG‑AAAGCGCAGCAATAAGCAATAACGGTACGACAGCTGTGTCGTGCCG‑TTTG‑tttttctgc                                                                                                                                                                                                                                                                                           >  1:652208/1‑136 (MQ=255)                                                                                                                                    aaCTGCCTGCAACACCACGCGTGGCG‑TTGGTGAAGACA‑TTTCTG‑ATGGCGGTAACGCGATTTCTGGCGCAGCAACG‑AAAGCGCAGCAATAAGCAATAACGGTACGACAGCTGTGTCGTGCCG‑TTTG‑tttttct                                                                                                                                                                                                                                                                                             >  1:1580316/1‑131 (MQ=255)                                                                                                                                    aaCTGCCTGCAACACCACGCGTGGCG‑TTGGTGAAGACA‑TTTCTG‑ATGGCGGTAACGCGATTTCTGGCGCAGCAACG‑AAAGCGCAGCAATAAGCAATAACGGTACGACAGCTGTGTCGTGCCG‑TTTG‑‑TTTTTCTGTCGATAGTCACAAAGGTAATAGTTGAAATT‑‑CCCCTGCCACCTGGC‑AAATATCCGTTC‑AACCATCA‑GC‑TTTGCAGGACGACCTGCAAGCG‑CCTC‑TTTTCA‑CC‑GGGGACGGG‑‑CCCCAA‑TT‑‑CTCCGGAGCCTGAT‑ATGTCCTGGATTATCTTAGTTATTGCTGGTCTGCTGGAAGTGGTATGGGCCG‑TTGGcc                                                                            >  1:652209/1‑324 (MQ=255)                                                                                                                                     aCTGCCTGCAACACCACGCGTGGCG‑TTGGTGAAGACA‑TTTCTG‑ATGGCGGTAACGCGATTTCTGGCGCAGCAACG‑AAAGCGCAGCAATAAGCAATAACGGTACGACAGCTGTGTCGTGCCG‑TTTG‑‑TTTTTCTG‑CGATAGTCACAAAGGTAATAGTTGAAATT‑‑CCCCTGCCACCTGGCAAAATATCCGTTC‑AACCATCA‑GC‑TTTGCAGGACGTCCTGCAAACG‑CCTC‑TTTTCACCC‑GGGGACGGG‑‑CCCCAA‑TTC‑CTCCGGAGCCTGAT‑ATGTCCTGGATTATCTTAGTTATTGCTGGTCTGCTGGAAGTGGTATGGGCCG‑TTGGCCTGAAatat                                                                    >  1:652211/1‑333 (MQ=255)                                                                                                                                     aCTGCCTGCAACACCACGCGTGGCG‑TTGGTGAAGACA‑‑TTCTG‑ATGGCGGTAACGCGATTTCTGGCGCAGCAACG‑AAAGCGCAGCAATAAGCAATAACGGTACGACAGCTGTGTCGTGCCG‑TTTG‑‑TTTTTCTG‑CGATAGTCACAAAGGTAATAGTTGAAATT‑‑CCCCTGCCa                                                                                                                                                                                                                                                  >  1:1580317/1‑170 (MQ=255)                                                                                                                                            gcAACACCACGCGT‑GCG‑TTGGTGAAGACA‑TTTCTG‑ATGGCGGTAACGCGATTTCTGGCGCAGCAACG‑AAAGCGCAGCAATAAGCAATAACGGTACGACAGCTGTGTCGTGCCG‑TTTG‑TTTTTTCTG‑CGATAGTCACAAAGGTAATAGTTGAAATT‑‑CCCCTGCCACCTGGCAAAATATCCGTTC‑AACCATCA‑GC‑TTTGCAGgacg                                                                                                                                                                                                       <  1:1580315/204‑1 (MQ=255)                                                                                                                                                        gTGGCG‑TTGGTGAAGACA‑TTTCTG‑ATGGCGGTAACGCGA‑TTCTGGCGCAGCAACG‑AAAGCGCAGCAATAAGCAATAACGGTACGACAGCTGTGTCGTGCCG‑TTTG‑‑TTTTTCTG‑CGATAGTCACAAAGGTAATAGTTGaa                                                                                                                                                                                                                                                                >  1:652212/1‑139 (MQ=255)                                                                                                                                                                tGGTGAAGACA‑TTTCTG‑ATGGCGGTAACGCGATTTCTGGCGCAGCAACG‑AAAGCGCAGCAATAAGCAATAACGGTACGACAGCTGTGTCGTGCCG‑‑TTG‑‑TTTTTCTG‑CGATAGTCACAAAGGTAATAGTTGAAATT‑‑CCCCTGCCACCTGGCAAAATATCCGTTC‑AACCATCA‑GC‑TTTGCAGGACGACCTGCa                                                                                                                                                                                                <  1:1580318/191‑1 (MQ=255)                                                                                                                                                                     aaGACAGTTTCTGAATGGCGGTAACGCGATTTCTGGCGCAGCAACG‑AAAGCGCAGCAATAAGCAATAACGGTACGACAGCTGTGTCGTGCCG‑TTTG‑‑TTTTTCTG‑CGATAGTCAC‑AAGGTAATAGTTGa                                                                                                                                                                                                                                                                 >  1:652213/1‑128 (MQ=255)                                                                                                                                                                                       cGGTAACGCGATTTCTGGCGCAGCAACG‑AAAGCGCAGCAATAAGCAATAACGGTACGACAGCTGTGTCGTGCCG‑TTTG‑‑TTTTTCTGTCGATAGTCACAAAGGTAATAGTTGAAATT‑‑CCCCTGCCACCTGGCAAAATATCCGTTC‑AACCATCA‑GC‑TTTGCAGGACGACCTGCAAACG‑CCTC‑TTTTCA‑CC‑GGGGAC‑GG‑‑CCCCAA‑TT‑‑CTCCGGAGCCTGAT‑ATGTCCTGGATTATCTTAGTTATTGCTGGTCTGCTGGAAGTGGTATggcc                                                                                     >  1:1580321/1‑266 (MQ=255)                                                                                                                                                                                                        cgcagacaACG‑AAAGCGCAGCAATAAGCAATAACGGTACGACAGCTGTGTCGTG‑CG‑‑TTG‑‑TTTTTCTG‑CGATAGTCAC‑AAGGTAATAGTTGAAATT‑‑CCCCTGCCACCTGGCAAAATATCCGTTC‑AACCATCA‑GC‑TTTGCAGGACGACCTGC‑‑ACG‑CCTC‑TTTTCA‑CC‑GGGGAC‑GG‑‑CCCCAA‑TT‑‑CTCCGGAGCCTGAT‑ATGTCCTGGATTATCTTAGTTATTGCTGGTCTGCTGGAAGTGGTATGGGCCG‑TTGGCCTGAAATATACCCACGGCTTTAGTCG‑TTTGACGCCGAGTGTTATTACTGTGACGGc                    <  1:1580320/302‑1 (MQ=255)                                                                                                                                                                                                                aCG‑AAAGCGCAGCAATAAGCAATAACGGTACGACAGCTGTGTCGTGCCG‑TTTG‑‑TTTTTCTG‑CGATAg                                                                                                                                                                                                                                                                                    >  1:1580322/1‑67 (MQ=255)                                                                                                                                                                                                                                         aCGGTACGACAGCTGTGTCGTGCCGTTTTG‑TTTTTTCTG‑CGATAGTCACAAAGGTAATAGTTGaaa                                                                                                                                                                                                                                                               >  1:1580324/1‑66 (MQ=255)                                                                                                                                                                                                                                         aCGGTACGACAGCTGTGTCGTGCCG‑TTTG‑‑TTTTTCTG‑CGATAGTCACAAAGGTAATAGTTGAAATT‑‑CCCCTGCCACCTGGC‑AAATATCCGTTCAAACCATCA‑GC‑TTTGCAGGACGACCTGCAAACG‑CCTC‑TTTTCA‑CC‑‑GGGAC‑GG‑‑CCCCAA‑TT‑‑CTCCGGAGCCTGAT‑ATGTCCTGGATTATCTTAGTTATTGCTGGTCTGCTGGAAGTGGTATGGGCC‑‑TTGGCCTGAAATATACCCACGGCttt                                                        >  1:1580325/1‑244 (MQ=255)                                                                                                                                                                                                                                                   aGCTGTGTCGTG‑CG‑TTTG‑TTTTTTCTG‑CGATAGTCACAAAGGTAATAGTTGAAATT‑‑CCCCTGCCACCTGGCAAAATATCCGTTC‑AACCATCA‑GC‑TTTGCAGGACGACCTGC‑AACG‑CCTC‑TTTTCA‑CC‑GGGGAC‑GG‑‑CCCCAA‑TT‑‑CTCCGGAGCCTGAT‑ATGTCCTGGATTATCTTAGTTATTGCTGGTCTGCTGGAAGTGGTATGGGCCGTTTGGCCTGaaa                                                                       <  1:1580323/221‑1 (MQ=255)                                                                                                                                                                                                                                                          tCGTGCCG‑TTTT‑GTTTTTCTG‑CGATAGTCACAAAGGTAATAGTTGAAATT‑‑CCCCTGCCACCTGGCAAAATATCCGTTC‑AACCATCA‑GC‑TTTGCAGGACGACCTGCAAACG‑CCTC‑TTTTCA‑CC‑GGGGAC‑GG‑‑CCCCAA‑TT‑‑CTCCGGAGCCTGAT‑ATGTCCTGGATTATCTTa                                                                                                                     >  1:1580327/1‑170 (MQ=255)                                                                                                                                                                                                                                                            gTGCCG‑TTTGTTTTTTTCTG‑CGATAGTCACAAAGGTAATAGTTGAAATT‑‑CCCCTGCCACCTGGCAAAATATCCGTTC‑AACC‑TCACGC‑TTTGCAGGACGCCTGGCAAACG‑CCTC‑TTTTCA‑CC‑GGGGAC‑GG‑‑CCCCAATTT‑‑CTCCGGAGCCTGATCATGTCCTGGATTATCTTa                                                                                                                     >  1:652217/1‑171 (MQ=255)                                                                                                                                                                                                                                                                      gtttttttCTG‑CGATAGTCACAAAGGTAATAGTTGAAATTTCCCCCTGCCACCTGGCAAAATATCCGTTC‑AACCATCA‑GCTTTTGCAGGACGACCTGCAAACGCCCTCTTTTTCA‑CCGGGGGAC‑GGGCCCCCAA‑TTTCCTCCGGAGCCTGAT‑ATGTCCTGGATTATCTTAGTTATTGCTGGTCTGCTGGAAGTGGTATGGGCCG‑TTGGCCTGAAATATACCCACGGCTTTAGTCGTTTTGACGCCGAGTGTTATTACTGTGACGGCGATGATTGTCAGTATGgc  >  1:1580329/3‑284 (MQ=255)                                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                                                                  TTTTTGGCGATGTTGTCTATTATTAATTTGCTATAGGCAAACATAAATAACATTACCT‑AAAAGGAAGACGTTATGGTGAAGAAGACAATTGCAGCGATCTTTTC‑TGTTCTGGTGC‑TTTCAACAGTATTAACTGCCTGCAACACCACGCGTGGCG‑TTGGTGAAGACA‑TTTCTG‑ATGGCGGTAACGCGATTTCTGGCGCAGCAACG‑AAAGCGCAGCAATAAGCAATAACGGTACGACAGCTGTGTCGTGCCG‑TTTG‑TTTTTTCTG‑CGATAGTCACAAAGGTAATAGTTGAAATT‑‑CCCCTGCCACCTGGCAAAATATCCGTTC‑AACCATCA‑GC‑TTTGCAGGACGACCTGCAAACG‑CCTC‑TTTTCA‑CC‑GGGGAC‑GG‑‑CCCCAA‑TT‑‑CTCCGGAGCCTGAT‑ATGTCCTGGATTATCTTAGTTATTGCTGGTCTGCTGGAAGTGGTATGGGCCG‑TTGGCCTGAAATATACCCACGGCTTTAGTCG‑TTTGACGCCGAGTGTTATTACTGTGACGGCGATGATTGTCAGTATGGC  >  NC_000913/4376481‑4377005

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 16 ≤ ATCG/ATCG < 25 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑