Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,854,595	C→T	A225V (GCC→GTC)	emrD →	multidrug efflux pump EmrD

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,854,595	0	C	T	91.7%	24.4 / NA	12	A225V (GCC→GTC)	emrD	multidrug efflux pump EmrD
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base T (5/6);  minor base G (1/0);  total (6/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

GATACGATGTGGAACTGGCGCGCCTGTTATCTCTTTTTGTTGGTTCTTTGTGCTGGTGTGACCTTCAGTATGGCCCGCTGGATGCCGGAAACGCGTCCGGTCGATGCACCGCGCACGCGCCTGCTTACCAGTTATAAAACGCTTTTCGGTAACAGCGGTTTTAACTGTTATTTGCTGATGCTGATTGGCGGTCTGGCCGGGATTGCCGCCTTTGAAGCCTGCTCCGGCGTGCTGATGGGCGCGGTGTTAGGGCTGAGCAGTATGACGGTCAGTATTTTGTTTATTCTGCCGATTCCGGCAGCGTTTTTTGGCGCATGGTTTGCCGGACGTCCCAATAAACGCTTCTCCACGTTAATGTGGCAGTCGGTTATCTGCTGCCTGCTGGCTGGCTTGCTGATGTGGATCCCCGACTGGTTTGGC  >  NC_000913/3854387‑3854806                                                                                                                                                                                                                 |                                                                                                                                                                                                                    gATACGATGTGGAACTGGCGCGCCTGGTATCTCTTTTTGTTGGTTCTTTGTGCTGGTGTGACCTTCAGTATGGCCCGCTGGATGCCGGAAACGCGTCCGGTCGATGCACCGCGCACGCGCCTGCTTACCAGTTATAAAACTCTTTTCGGTAACAGCGGTTGTAACTGTTATTTGCTGATGCTGATTGGCGGTCTGGCCGGGATTGCCGTCTTTGAAGCCTGCTCCGGCGTGCTGATGGGCGCGGTGTTAgg                                                                                                                                                                           <  1:148329/251‑1 (MQ=255)                   ngcgcCTGTTATCTCTTTTTGTTGGTTCTTTGTGCTGGTGTGACCTTCAGTATGGCCCGCTGGATGCCGGAAACGCGTCCGGTCGATGCACCGCGCACGCGCCTGCTTACCAGTTATAAAACGCTTTTCGGTAACAGCGGTTTTAACTGTTATTTGCTGATGCTGATTGGCGGTCTGGCCGGGATTGCCGTCTTTGAAGCCTGCTCCGGCGTGCTGATGGGCGCGGTGTTAgg                                                                                                                                                                           >  1:54186/2‑233 (MQ=255)                                                           tgACCTTCAGTATGGCCCGCTGGATGCCGGAAACGCGTCCGGTCGATGCACCGCGCACGAGCCTGCTTACCAGTTCTAACACGCTTTTCGGTAACAGCGGTTTTAACTGTTATTTGCTGATGCTGATTGGCGGTCTGTCCGGGATTGCCGTCTTTGAAGCCTGCTCCGGCGTGCTGATGGGCGCGGTGTTAGGGCTGAGCAGTATGACGGTCAGTATTTTGTTTATTCTGCCGATTc                                                                                                                               >  1:88022/1‑237 (MQ=255)                                                                 tCAGTATGGCCCGCTGGATGCCGGAAACGCGTCCGGTCGATGCACCGCGCACGCGCCTGCTTCCCAGTTTTAAAACGCTTTTCGGTAACAGCGGTTTTAACTGTTATTTGCTGATGATGATTGGCGGTCTGGCCGGGATTGCCGGCTTTGAAGCCTGCTCCGGCGTGCTGATGGGCGAGGTGTTAGGGCTGAGCAGTATGACGGTCAGTATTTTGTTTATTCTGCCGATTCCGGCGGCGTTTTTTGgcgcc                                                                                                           >  1:106013/1‑250 (MQ=255)                                                                      aTGGCNNGCNNNATNNCGGAAACGCGTCCGGTCGATGCANCGCGCACGCGCCTGCTTACCAGTTATCAANCNCNTTTCNNTAANAGCGGTTTTANCTGTTATTTGCTGATGCTGATTGGCGGTCTGGCCGGGATTGCCGTCTTNGNAGCCTGCTCCGGCGTGCt                                                                                                                                                                                             >  2:160756/1‑164 (MQ=255)                                                                      aTGGCCCGCTGGATGCCGGAAACGTGTCCGGTCGATGCACCGCGCACGCGCCTGCTTACCAGTTATAAAGCGCTTTTCTGTAACAGTGGGTTTAACTGTTATTTTCGGATTCTTATTGGCGGTCTGGCCGGGATTGCCGTCTTTGAAGCCTGCTCCGGCGTGCt                                                                                                                                                                                             <  1:160756/164‑1 (MQ=255)                                                                                ggATGCCGGAAACGCGTCCGGTCGATGCACCGCGCACGCGCCTGCTTACCAGTTATAAAACGCTTTTCGGTAACAGCGGTTTTAACTGTTATTTGCTGATGCTGATTGGCGGTCTGGCCGGGATTGCCGTCTTTGAAGCCTGCTCCGGCGTGCTGATGGGCGCGGTGTTAGGGCTGAGCAGTATGACGGTCAGTATTTTGTTTATTCTGCCGATTCCGGCAGCGTTTTTTGGCGCATGGTTTGCCGGACNt                                                                                            <  1:171487/251‑1 (MQ=255)                                                                                                                                              cTTTTCGGTAACAGCGGTTTTAACTGTTATTTGCTGATGCTGATTGGCGGGCTGGCCGGGATTGCCGTCTTTGAAGCCGGCTCCGGCGTGCTGATGGGCGCGGTGTTAGGGCTGAGCAGTAGTACGGGCAGTATGTTGTTTAGTCTGCCGAGTCCGGCAGCGTTTTTTGGCGCATGGTTTGCCGGACGTCCCa                                                                                        <  1:193703/193‑1 (MQ=255)                                                                                                                                                                    aCTGTTATTTGCTGATGCTGATTGGCGGTCTGGCCGGGATTGCCGTCTTTGAAGCCTGCTCCGGCGTGCTGATGGGCGCGGTGTTAGGGCTGAGCAGTATGACGGTCAGTATTTTGTTTATTCTGCCGATTCCGGCAGCGTTTTTTGGCGCATGGTTTGCCGGACGTCCCAATAAACGCTTCTCCACGTTAATGTGGCAGTCGGTTATCTGCTGCCTGCTGGCTGGCTTGCTGATGTGGATCCCCGACTgg        <  1:119719/251‑1 (MQ=255)                                                                                                                                                                          aTTTGCTTAGGCTGATTGGCGGTCTGGCCGGGATTGCCGTCTTTGAAGCCTGCTCCGGCGTGTTGATGGGCGCGGTGTTAGGGCTGAGCAGTATTACGGTCAGTATTTTGGTTATTCTGGCGATTCCGGCAGCGTTTTTTGGCGCAGGGGGTGTCGGTCGTCCCAATAAACGCTTCTCCACTTTAAGGTGTCAGTCGGTTATCTGCTGCCTGCTGGCTGGCTTGCTGATGTGGATCCCCGACTGGTTTGGc  <  1:143232/251‑1 (MQ=255)                                                                                                                                                                                                   gnCCGGGATTGACGTCTTTGAAGCCTGCTCCGGCGTGCTGATGGGCGCGGTGTTAGGGCTGAGCAGTATGACGGTCAGTATTTTGTTTATTCTGCCGATTCCGGCAGCGTTTTTTGgc                                                                                                              >  1:171407/1‑118 (MQ=255)                                                                                                                                                                                                    gCCGGGATTGCCGTCTTTGAAGCCTGCTCCGGCGTGCTGATGGGCGCGGTGGTAGGGCTGAGCCGCATGACGGTCAGTATTTTGTTTATTCTGCGGATTCCGGCAGCGTTTTTTGGCGCCTGGTTTGCCGGAc                                                                                              >  1:86630/1‑133 (MQ=255)                                                                                                                                                                                                                 |                                                                                                                                                                                                                    GATACGATGTGGAACTGGCGCGCCTGTTATCTCTTTTTGTTGGTTCTTTGTGCTGGTGTGACCTTCAGTATGGCCCGCTGGATGCCGGAAACGCGTCCGGTCGATGCACCGCGCACGCGCCTGCTTACCAGTTATAAAACGCTTTTCGGTAACAGCGGTTTTAACTGTTATTTGCTGATGCTGATTGGCGGTCTGGCCGGGATTGCCGCCTTTGAAGCCTGCTCCGGCGTGCTGATGGGCGCGGTGTTAGGGCTGAGCAGTATGACGGTCAGTATTTTGTTTATTCTGCCGATTCCGGCAGCGTTTTTTGGCGCATGGTTTGCCGGACGTCCCAATAAACGCTTCTCCACGTTAATGTGGCAGTCGGTTATCTGCTGCCTGCTGGCTGGCTTGCTGATGTGGATCCCCGACTGGTTTGGC  >  NC_000913/3854387‑3854806

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