Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,812,327	0	G	A	60.0%	2.0 / 10.1	15	A78V (GCC→GTC)	yicR	RadC‑like JAB domain‑containing protein YicR
Reads supporting (aligned to +/- strand):  ref base G (4/2);  new base A (5/4);  total (9/6)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: E-value score below prediction cutoff.
Rejected as consensus: Frequency below/above cutoff threshold.

GATGGACTTCAACATGGTTTAGCGTGCCGGAAAAAAGACGCCTATGCGTTATAACCCGGTGTTGGGAGTCGAGAAAGATCACCATAAAGATCTCCCGCTCCTCACCCGTGAGCTGGCTTTGTAAAAATTCCCTCGTCATCTCCGGGCTGAGTAAAGGGCTTTCTTCACGCATCCGCACGTTGTAGTAACGCCGCGCCAGTTCAGCAATCCCCTTTAACTGGGCAAATTTCGCCACGCCAATTCCATGAACGCCACTAAATTGTTCATATTCAGAGGTTAACAAGCCATAAAGAGAGCCGAAATTCTCCAGCATCTCTTTTGCCAGGGTTAATACATCTTTACCGCGCGTTCCGGTACGCAGAAATAGCGCCAGCAGCTCGACATCCGTTAAGGCGCTAATACCAAACTTCAGCATTTTTTCGCGCGGCATCAAC  >  NC_000913/3812106‑3812539                                                                                                                                                                                                                              |                                                                                                                                                                                                                     gATGGACTTCAACATGGTTTAGCGTGCCGGAAAAAAGACGCCTATGCGTTATAACCCGGTGTTGGGAGTCGAGAAAGATCACCATAAAGATCTCCCGCTCCTCACCCGTGAGCTGGCTTTGTAAAAATTCCCTCGTCATCTCCGGGCTGAGTAAAGGGCTTTCTTCACGCATCCGCACGTTGTAGTAACGCCGCGCCAGTTCAGCAATCCCCTTTAACTGGACAAATTTCGCCACGCCAATTCCATGAACg                                                                                                                                                                                         <  1:259187/251‑1 (MQ=255)                                                                                                                                                               cgttcttcACGCAGCCGCACGGTGTAGTAACGCCGCGCCAGTTCAGCAATCCCCTTGAACTGGACAAATTTCGCCACGCCAATTCCATGAACGCCACTAAATTGTTCATATTCAGAGGTTAACAAGCCATAAAGAGAGCCGAAATTCTCCAGCATCTCTTTTGc                                                                                                                  <  1:244791/162‑1 (MQ=255)                                                                                                                                                               cTTTCNNCANNNATNNGCACGTTGTAGTAACGCCGCGCCNGTTCAGCAATCCCCTTTAACTGGACAAATNTNGNCACGNNAATNCCATGAACGCNACTAAATTGTTCATATTCAGAGGTTAACAAGCCATACAGAGAGCAGAANTNCTCCAGCATCTCTTTTGc                                                                                                                  >  2:244791/1‑164 (MQ=255)                                                                                                                                                                              cGCACGTTGTAGTAACGCCGCGCCAGTTCAGCAATCCCCTTTAACTGGGCAAATTTCGCCACGCCAATTCCATGAACGCCACTAAAGTGTTCATATTCAGAGGGTAACAAGCCATAAAGAGAGCCGAAATTCTCCAGCATCTCTTTTGCCAGGGTTAATACATCTTTACCGCGCTTTCCGTTACGCAGAAATAGCGCCAGCAGCTCGACATCCGTTAAGGCGCTAATACCAAACTTCAGCATTTTTTcg              <  1:311802/249‑1 (MQ=255)                                                                                                                                                                              cGCACGTTGTAGTAACGCCGCGCCAGTTCAGCAATACCCTTTAACTGGGCAAATTTCGCCACGCCAATTCCATGAACGCCACTAAATTGTTCATATTCAGAGGTTAACAAGCCATAAAGAGAGCCGAAATTCTCCAGCATCTCTTTTGCCAGGGTTAATACATCTTTACCGCGCGTTCCGGTACGCAGAAATAGCGCCAGCAGCTCGACATCCGTTAAGGCGCTAATACCAAACTTCAGCATTTTTTcgcg            >  1:269896/1‑251 (MQ=255)                                                                                                                                                                                     tGTAGTAACGCCGCGCCAGTTCAGCAATCCCCTTTAACTGGACAAATGTCGCCACGCCAATTCCATGAACGCCACTAAATTGTTCATATTCAGAGGTTAACAAGCCATAAAGAGAGCCGAAATTCTCCAGCATCTCTTTTg                                                                                                                   <  1:165265/141‑1 (MQ=255)                                                                                                                                                                                        agtaACGACGCGCCAGTTCAGCAATCCACTTTAACTGGAAAAATTTCGCCACGCCACTTCCATGAACGCCACTAAATTGTTCATATTCAGAGGTTAACAAGCCATAAAGAGAGCCGAAATTCTCCAGCATCTCTTTTGCCAGGGTTAATACATCTTTACCGCGCGTTCCGGTACGCAGAAATAGCGCCAGCAGCTCGACATCCGTTAAGGCGCTAATACCAAACTTCAGCATTTTTTCGCGCGGCATCAAc  >  1:275913/1‑251 (MQ=255)                                                                                                                                                                                         gtaACGCCGCGCCAGTTCAGCAATCCCCTTTAACTGGACACATTTCGCCACGCCAATTCCATGAACGCCACTACATTGTTCATATTCAGAGGTTAACAAGCCATAAAGAGAGCCGCAATTCTCCAGCATCTCTTTTGCCAGGGTTAATACATCTTTACCGCGCGTTCCGGTACGCAGAACTAGCGCCAGCAGCTCGACATCCGTTAAGGCGCTAATACCAAACTTCAGCATTTTTTCGCGCGGCATCAAc  >  1:321614/1‑250 (MQ=255)                                                                                                                                                                                         gtaACGCCGCGCCAGTTCAGCAATCCCCTTTAACTGGACAAATTTCGCCACGCCAATTCCATGAACGCCACTAAATTGTTCATATTCAGAGGTTAACAAGCCATAAAGAGAGCCGAAATTCTCCAGCATCTCTTTTGCCAGGGTTAATTCATCTTTACCGCGCGTTCCGGTCCGCAGAAATAGCGCCAGCAGCTCGACATCCGTTAAGGCGCTAATACCAAACTTCAGCATTTTTTCGCGCGGCATCAAc  >  1:321624/1‑250 (MQ=255)                                                                                                                                                                                                  nnccNNNNNAGCNATCCCCTTTAACGGGGCAAATTTCGCCACGCCAATTCCATGAACGCCACTAAATTGTTCATATTCAGAGGTTAACAAGCCATAAAGAGAGCCGAAATTCTCCAGCATCTCTTTTGCCAGGGTTAATACATCtt                                                                                                 <  2:193092/144‑1 (MQ=255)                                                                                                                                                                                                  ngcCAGTTCAGCAATCCCCTTTAACTGGGCAAATTTCGCCACGCCAATTCCATGAACGCCACTAAATTGTTCATATTCAGAGGTTAACAAGCCATAAAGAGAGCCGAAATTCTCCAGCATCTCTTTTGCCAGGGTTAATACATCtt                                                                                                 >  1:193092/2‑146 (MQ=255)                                                                                                                                                                                                   gcCAGTTCAGCAATCCCCTTTAACTGGGCAAATTTCGCCAAGCCAATTCCATGAACGCCACTAAATTGTTCATCTTCAGAGGTTCACAAGCAATAAAGAGAGCCGAAATTCTCCAGAATCTCTTTTGCCAGGGTTAATACATCTTTACTGCGCGTTCCGGAACGCAGAAATAGCGCCAGCAGCTCGACATCCGTTAAGGCGCAAATACCAAACTTCAGCAtttttt                >  1:329430/1‑226 (MQ=255)                                                                                                                                                                                                            gCAATCCCCTTTAACTGGACAAATTTCGCCACGCCAATTCCATGAACGCCACTAAATTGTTCATATTCAGAGATTAACAAGCCATAAAGAGATCCGAAATTCTCCAGCATCTCTTTTGCCAGGNt                                                                                                            <  1:253368/125‑1 (MQ=255)                                                                                                                                                                                                                 ccccTTAAACTGGACAAATTTCGCCACGCCAATTCCATGAACGCCACTCACATGTACATATTCCGAGGTTCACAAGCCATAAAGAGAGCCGAAATTCTCCAGCCTCTCTTTTCCCAGGGATAATACCTCTTTACCGCGCGTTCCGGTACGCAGAAATAGCGCCAGCAGCTCGACATACGTTAAGgcg                                         >  1:53871/1‑187 (MQ=255)                                                                                                                                                                                                                    cTTTAACTGGGCAAATTTCGCCACGCCAATTCCATGAACGCCACTAAATTGTTCATATTCAGAGGTTAACAAGCCATAACGAGAGCCGAAATTCTCCAGCATCTCTTTTGCCAGGGTTAATACATCTTTACAGCGCGTTCCGGTACGCTGAAATAGCGCCagcagc                                                           >  1:324359/1‑166 (MQ=255)                                                                                                                                                                                                                              |                                                                                                                                                                                                                     GATGGACTTCAACATGGTTTAGCGTGCCGGAAAAAAGACGCCTATGCGTTATAACCCGGTGTTGGGAGTCGAGAAAGATCACCATAAAGATCTCCCGCTCCTCACCCGTGAGCTGGCTTTGTAAAAATTCCCTCGTCATCTCCGGGCTGAGTAAAGGGCTTTCTTCACGCATCCGCACGTTGTAGTAACGCCGCGCCAGTTCAGCAATCCCCTTTAACTGGGCAAATTTCGCCACGCCAATTCCATGAACGCCACTAAATTGTTCATATTCAGAGGTTAACAAGCCATAAAGAGAGCCGAAATTCTCCAGCATCTCTTTTGCCAGGGTTAATACATCTTTACCGCGCGTTCCGGTACGCAGAAATAGCGCCAGCAGCTCGACATCCGTTAAGGCGCTAATACCAAACTTCAGCATTTTTTCGCGCGGCATCAAC  >  NC_000913/3812106‑3812539

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