Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,898,467	C→T	A149V (GCC→GTC)	cbrC →	PF03691 family colicin E2 tolerance protein CbrC

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,898,467	0	C	T	94.7%	45.4 / NA	19	A149V (GCC→GTC)	cbrC	PF03691 family colicin E2 tolerance protein CbrC
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base T (10/8);  minor base G (0/1);  total (10/9)
Fisher's exact test for biased strand distribution p-value = 4.74e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GCATAGAAGGTGTTGAATTTGAGTATGATGAAGAGGACGAATTTGCCGGTATTAAGAACACATATCCTGATGAAATGCTGAAAGAGTTGGTTGAACGCACGCCAGGTTATCATGGATGGCAGCAGGAATTCTGGCTCGCGCATTGTGGCGATTTCTGTGTTTTTATCGGCTATGTGGGCTGGAATGATATAAAAGATCGCCTCGATGAATTTGCCAACCTTGAAGAAGATTGTGAGAATTTCGGTATTAGAAATTCTGATCTAGCTAAATGCCTGCAAAAGGGTGGTCATTGTCAGGGTTATCTCTTCCGCTGTCTCCACTGCGGCAAGCTGAGACTGTGGGGTGATTTTTCGTAGTTATTTAAATAATGAGAACAGGCCGGAGCGTAATTCACACATCCGGCCTTATTTCTTAAGCTTAACCGAGTGCTAATTCTGCGGCAGCGGCT  >  NC_000913/3898254‑3898701                                                                                                                                                                                                                      |                                                                                                                                                                                                                                           ncATAGAAGGTGTTGAATTTGAGTATGATGAAGAGGACGAATTTGCCGGTATTAAGAACACATATCCTGATGAAATGCTGAAAGAGTTGGTTGAACGCACGCCAGGTTATCATGGATGGCAGCAGGAATTCTGGCTCGCGCATTGTGGCGATTTCTGTGTTTTTATCGGCTATGTGGGCTGGAATGATATAAAAGATCGCCTCGATGAATTTGTCAACCTTGAAGAAGATTGTGAGAATTTCGGTAtt                                                                                                                                                                                                          >  1:406122/2‑248 (MQ=255)                                   ggACGAATTTGCCGGTATTAAGAACACATATCCTGATGAAATGCTGAAAGAGTTGGTTGAACGCACGCCAGGTTATCATGGATGGCAGCAGGAATTCTGGCTCGCGCATTGTGGCGATTTCTGTGTTTTTATCGGCTATGTGGGCTGGAATGATATAAAAGATCGCCTCGATGAATTTGTCAACCTTGAAGAAGATTGTGAGAATTTCGGTATTAg                                                                                                                                                                                                        <  1:22346/216‑1 (MQ=255)                                                  tATTAAGAACACATATCCTGATGAAATGCTGAAAGAGTTGGTTGAACGCACGCCAGGTTATCATGGATGGCATCAGGAATTCTGGCTCGCGCATGGTGGCGATTTCTGTGTTTTTATCGGCTATGTGGGCTGGAATGATATAAAAGATCGCCTCGATGAACTTGTCAACCGTCAAGAAGATTGTGAGAATTTCGGt                                                                                                                                                                                                             <  1:166032/196‑1 (MQ=255)                                                             caTATCCTGATGAAATGCTGAAAGAGTTGGTTGAACGCACGCCAGGTTATCATGGATGGCAGCAGGAATTCTGGCTCGCGCATTGTGGCGATTTCTGTGTTTTTATCGGCTATGTGGGCTGGAATGATATAAAAGATCGCCTCGATGAATTTGTCAACCTTGAAGAAGATTGTGAGAATTTCGGTATTAGAAATTCTGATCTAGCTAAATGCCTGCAAAAGGGTGGTCATTGTCAGGGTTATCTCTTCCGc                                                                                                                                           >  1:20618/1‑251 (MQ=255)                                                                                               acgcaagcCCGGTTATCATGGCTGGCAGCAGGAATTCTGGCTCGCGCATTGTGGCGATTTCTGTGTTTTTATCGGCTATGTGGGCTGGAATGATATAAAAGATCGCCTCGATGAATTTGTCAACATTGAAGAAGATTGTGAGAATTTCGGTATTAGAAATTCTGATCTAGCTAAATGCCTGCAAAAGGGTGGTCATTGTCAGGGTTATCTCTTCCGCTGTCTCCACTGCGGCAAGCTGAGACTGTGGGGTg                                                                                                         >  1:401574/1‑251 (MQ=255)                                                                                                       cAGGTTATCATGGATGGCAGCAGGAATTCTGGCTCGCGCATTGTGGCGATTTCTGTGTTTTTATCGGCTATGTGGGCTGGAATGATATAAAAGATCGCCTCGATGAATTTGTCAACCTTGAAGAAGATTGTGAGAATTTCGGTATTAGAAATTCTTATCTAGCTAAATGCCTGCAAAAGGGTGGTCATTGTCAGGGTTATCTCTTCCGCTg                                                                                                                                         <  1:439527/211‑1 (MQ=255)                                                                                                                       gcagcaGGAATTCTGGCTCGCGCATTGTGGCGATTTCTGTGTTTTTATCGGCTATGTGGGCTGGAATGATATAAAAGATCGCCTCGATGAATTTGTCAACCTTGAAGAAGATTGTGAGAATTTCGGTATTAGAAATTCTGATCTAGCTAAATGCCTGCAAACGGGTGGTCATTGTCAGGGTTATCTCTTCCGCTGTCTCCACTGCGGCAAGCTGAGACTGTGGGGAGATTTTTCGTAGTTATTTAaataat                                                                                 >  1:482015/1‑251 (MQ=255)                                                                                                                           cagGAATTCTGGCTCGCGCATTGTGGCGATTTCTGTGTTTTTATCGGCTATGAGGGCTATAACGATATAAAAGATCGCCTCGATGAATTTGTCAACATTGAAGTAGATTGTGCGAATTTCGGTATTAGAACTTCTGATATAGCTAACTGCCTGCAAAAGGGTGGCCATTGTCAGGGTTATCTCTTCCGCTGTCTCCACTGCGGCAAGc                                                                                                                        >  1:359638/1‑208 (MQ=255)                                                                                                                                          gcgcATTGTGGCGATTTCTGTGTTTTTATCGGCTATGTGGGCTGGAATGATATAAAAGATCGCCTCGATGAATTTGTCAACCTTGAAGAAGATTGTGAGAATTTCGGTATTAGAAATTCTGATCTAGCTAAATGCCTGCAAAAGGGTGGTCATTGTCAg                                                                                                                                                          <  1:217774/159‑1 (MQ=255)                                                                                                                                                   ggCGATTTCTGTGTTTTTATCGGCTATGTGGGCTGGAATGATATAAAAGATCGCCTCGATGAATTTGTCAACCTTGAAGAAGATTGTGAGAATTTCGGTATTAGAAATTCTGATCTAGCTAAATGCCTGCAAAAGGGTGGTCATTGTCAGGGTGATCTCTTCCGCTGTCTCCACTGCGGCAAGCTGAGACTGTGGGGTGATTTTTCg                                                                                                 <  1:258045/207‑1 (MQ=255)                                                                                                                                                          tCTGTGTTTTTATCGGCTATGTGGGCTGGAATGATATAAAAGATCGCCTCGATGAATTTGTCAACCTTGAAGAAGATTGTGAGAATTTCGGTATTAGAAATTTTGATCTAGCTAAATGCCTGCAAAAGGGTGGTAATTGTTAGGGTTATCTCTTCCGCt                                                                                                                                          >  1:20311/1‑159 (MQ=255)                                                                                                                                                                                 ggCTGGAATGATATAAAAGATCGCCTCGATGAATTTGTCAACCTTGAAGAAGATTGTGAGAATTTCGGTATTAGAAATTCTGAt                                                                                                                                                                                              >  1:211970/1‑84 (MQ=255)                                                                                                                                                                                   cTGGAATGATATAAAAGATCGCCTCGATGAATTTGTCAACCTTGAAGAAGATTGTGAGAATTTCGGTATTAGAAATTCTGATCTAGCTAAATGCCTGCAAAAGGGTGGTCATTGTCAGGGTTATCTCTTCCGCTGTCTCCACTGCGGCAAGCTGAGCCTGTGGGGTGATTTTTCGTAGTTATTTAAATAATGAGAACAGGCCGGAGCGTAATTCACACATCCGGCCTTATTTCTTAAGCTTAACCGAGTGc                     >  1:199010/1‑251 (MQ=255)                                                                                                                                                                                             ataAAAGATCGCCTCGATGAATTTGGCAACCTTGAAGAAGATTGTGAGAATTTCGGTATTAGAAATTCTGATCCAGCTAAATGCCTGCAAAAGGGTGGTCATTGTCAGGGTTATCTCTTCCGCTGTCTCCACTGCGGCAAGCTGAGACTGTGGGGTGATTTTTCGGAGTTATTTAAATAATGAGAACAGTCCGGAGTGTAATTCACACATCCGGCCTTATTTCTTAAGCTTAACCGAGTGCTAATTCTGCg           <  1:472784/251‑1 (MQ=255)                                                                                                                                                                                               aaaaGATCGCCTCGATGAATTTGTCAACCTTGAAGAAGATTGTGAGAATTTCGGTATTAGAAATTCTGATCTAGCTAAATGCCTGCAAAAGGGTGGTCATTGTCAGGGTTATCTCTTCCGCTGTCTCCACTGCGGCAAGCTGAGACTGTGGGGTGATTTTTCGTAGTTATTTAAATAATGAGAACAGGCCGGAGCGTAATTCACACATCCGGCCTTATTTCTTAAGCTTAACCGAGTGCTAATTCTGCggn         <  1:146073/251‑2 (MQ=255)                                                                                                                                                                                                   gATCGCCTCGATGAATTTGTCAACCTTGAAGAAGATTGTGAGAATTTCGGTATTAGAAATTCTGATCTAGCTAAATGCCTGCAAAAGGGTGGTCATTGTCAGGGTTATCTCTTCCGCTGTCTCCACTGCGGCAAGCTGAGACTGTGGGGTGATTTTTCGTAGTTATTTAAATAATGAGAACAGGCCGGAGCGTAATTCACACATCCGGCCTTATTTCTTAAGCTTAACCGAGTGCTAATTCTGCGGCAGCg     <  1:283815/251‑1 (MQ=255)                                                                                                                                                                                                    aTCGCCTCGATGAATTTGTCAACCTTGAAGAAGATTGTGAGAATTGCGGGATTAGAAATTCTGATCTAGCTAAATGCCTGCAAAAGGGTGGTCATTGTCAGGGTTATCTCTTCCGCTGTCTCCACTGCGGCAAGCTGAGACTGTGGGGTGATTTTTCGTAGTTATTTAAATAATGAGAACAGGCCGGAGCGTAATTCACACATCCGGCCTTATTTCTTAAGCTTAACCCATTGTTAATTCTGCGGCAGCgg    >  1:10708/1‑251 (MQ=255)                                                                                                                                                                                                               gAATTTGTCAACCTTGAAGAAGATTGTGAGAATTTCGGTATTAGAAATTCTGATCTAGCTACATGCCTGCAAAAGGGTGGTCATTGTCAGGGTTATCTCTTCCGCTGTCTCCACTGCGGCAAGCTGAGACTGTGGGGTGATTTTTCGTAGTTATTTAAATAATGAGAACAGGCCGGAGCGTAATTCACACATCCGGCCTTATTTCTTAAGCTTAACCGAGTGCTAATTCTGCGGCAGCGGCt  >  1:113732/1‑242 (MQ=255)                                                                                                                                                                                                                 aTTTGTCAACCTTGAAGAAGATTGTGAGAATTTCGGTATTAGAAATTCtn                                                                                                                                                                                                <  1:93346/50‑2 (MQ=255)                                                                                                                                                                                                                      |                                                                                                                                                                                                                                           GCATAGAAGGTGTTGAATTTGAGTATGATGAAGAGGACGAATTTGCCGGTATTAAGAACACATATCCTGATGAAATGCTGAAAGAGTTGGTTGAACGCACGCCAGGTTATCATGGATGGCAGCAGGAATTCTGGCTCGCGCATTGTGGCGATTTCTGTGTTTTTATCGGCTATGTGGGCTGGAATGATATAAAAGATCGCCTCGATGAATTTGCCAACCTTGAAGAAGATTGTGAGAATTTCGGTATTAGAAATTCTGATCTAGCTAAATGCCTGCAAAAGGGTGGTCATTGTCAGGGTTATCTCTTCCGCTGTCTCCACTGCGGCAAGCTGAGACTGTGGGGTGATTTTTCGTAGTTATTTAAATAATGAGAACAGGCCGGAGCGTAATTCACACATCCGGCCTTATTTCTTAAGCTTAACCGAGTGCTAATTCTGCGGCAGCGGCT  >  NC_000913/3898254‑3898701

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