Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	4,439,967	T→A	16.6%	I34N (ATT→AAT)	yfbS →	putative transporter

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	4,439,967	0	T	A	16.6%	57.0 / 4.3	24	I34N (ATT→AAT)	yfbS	putative transporter
Reads supporting (aligned to +/- strand):  ref base T (11/9);  new base A (2/2);  total (13/11)
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

GGAACACCCCCTCCCCATGAACCATGAGTTGCTCTGGGTCCTCGGCCTGCTGGCCATCGTGATCATCCTGTTCGTCATCAACCGCCCACGCATGGACGTGGTCGCGCTGCTGGTGATTCTCGCCCTGCCACTACTCGGCATCCTTACCGTGGAACAGGCCCTGGCCGGGTTCAGCGACCCCAACGTGGTGCTGATCGCCGCGCTGTTCGTGATTGGCGAAGGCCTGGTGCGCACCGGCATTGC  >  NC_002947/4439851‑4440093                                                                                                                     |                                                                                                                               ggAACACCCCCTCCCCATGAACCATGAGTTGCTCTGGGTCCTCGGCCTGCTGGCCATCGTGATCATCCTGTTCGTCATCAACCGCCCACGCATGGACGTGGTCGCGCTGCTGGTGATTCTCGCCCTGCCACTACTc                                                                                                             <  2:223497/136‑1 (MQ=255)        ccccTCCCCATGAACCATGAGTTGCTCTGGGTCCTCGGCCTGCTGGCCATCGTGATCATCCTGTTCGTCATCAACCGCCCACGCATGGACGTGGTCGCGCTGCTGGTGATTCTCGCCCTGCCACTACTCGGCATcc                                                                                                      >  2:201619/1‑136 (MQ=255)        ccccTCCCCATGAACCATGAGTTGCTCTGGGTCCTCGGCCTGCTGGCCATCGTGATCATCCTGTTCGTCATCAACCGCCCACGCATGGACGTGGTCGCGCTGCTGGTGATTCTCGCCCTGCCACTACTCGGCATcc                                                                                                      <  1:35933/136‑1 (MQ=255)                cATGAACCTTGAGTTGCTCTGGGTCCTCGGCCTGCTGGCCATCGTGATCATCCTGTTCGTCATCAACCGCCCACGCATGGACGTGGTCGCGCTGCTGGGGATTCTCGCCCTGCCACTACTCGGCATCCTTACCGTg                                                                                              <  1:353094/136‑1 (MQ=255)                 aTGAACCATGAGTTGCTCTGGGTCCTCGGCCTGCTGGCCATCGTGATCATCCTGTTCGTCATCAACCGCCCACGCATGGACGTGGTCGCGCTGCTGGTGATTCTCGCCCTGCCACTACTCGGCATCCTTACCGTgg                                                                                             <  1:369136/136‑1 (MQ=255)                    aaCCATGAGTTGCTCTGGGTCCTCGGCCTGCTGGCCATCGTGATCATCCTGTTCGTCATCAACCGCCCACGCATGGACGTGGTCGCGCTGCTGGTGATTCTCGCCCTGCCACTACTCGGCATCCTTACCGTGGAAc                                                                                          <  1:309206/136‑1 (MQ=255)                           aGTTGCTCTGGGTCCTCGGCCTGCTGGCCATCGTGATCATCCTGTTCGTCATCAACCGCCCACGCATGGACGTGGTCGCGCTGCTGGTGATTCTCGCCCTGCCACTACTCGGCATCCTTACCGTGGAACAGGCCCt                                                                                   >  2:410189/1‑136 (MQ=255)                           aGTTGCTCTGGGTCCTCGGCCTGCTGGCCATCGTGATCATCCTGTTCGTCATCAACCGCCCACGCATGGACGTGGTCGCGCTGCTGGTGATTCTCGCCCTGCCACTACTCGGCATCCTTACCGTGGAACAGGCCCt                                                                                   <  1:56825/136‑1 (MQ=255)                           aGTTGCTCTGGGTCCTCGGCCTGCTGGCCATCGTGATCATCCTGTTCGTCATCAACCGCCCACGCATGGACGTGGTCGCGCTGCTGGTGAATCTCGCCCTGCCACTACTCGGCATCCTTACCGTGGAACAGGCCCt                                                                                   >  2:154226/1‑136 (MQ=255)                           aGTTGCTCTGGGTCCTCGGCCTGCTGGCCATCGTGATCATCCTGTTCGTCAACAACCGCCCACGCATGGACGAGGTCGCGCTGCTGGTGAATCTCGCCCAGCCACAACTCCGCAACCATACCGTGGAACAGGcaca                                                                                   >  2:154229/1‑133 (MQ=255)                                        ccTCGGCCTGCTGGCCATCGTGATCATCCTGTTCGTCATCAACCGCCCACGCATGGACGTGGTCGCGCTGCTGGTGATTCTCGCCCTGCCACTACTCGGCATCCTTACCGTGGAACAGGCCCTGGCCGGGTTCAGc                                                                      >  1:165412/1‑136 (MQ=255)                                          gcGGCCTGCTGGCCATCGTGATCATCCTGTTCGTCATCAACCGCCCACGCATGGACGTGGTCGCGCTGCTGGTGATTCTCGCCCTGCCACTACTCGGCATCCTTACCGTGGAACAGGCCCTGGCCGGGTTCAGCGa                                                                    <  1:62110/135‑1 (MQ=255)                                             gCCTGCTGGCCATCGTGATCATCCTGTTCGTCATCAACCGCCCACGCATGGACGTGGTCGCGCTGCTGGTGATTCTCGCCCTGCCACTACTCGGCATCCTTACCGTGGAACAGGCCCTGGCCGGGTTCAGCGcccc                                                                 >  2:303641/1‑136 (MQ=255)                                             gCCTGCTGGCCATCGTGATCATCCTGTTCGTCATCAACCGCCCACGCATGGACGAGGTCGCGCTGCTGGTGATACTCGCCCTGCCACAACTCGCCATCCTAACCGTGGAACAGGCCCTGGCCGGGTTCAGCGAccc                                                                 >  2:197651/1‑136 (MQ=255)                                                    ggCCATCGTGATCATCCTGTTCGTCATCAACCGCCCTCGCATGGACGTGGTCGCGCTGCTGGTGAATCTCGCCCTGCCACTACTCGGCATCCTTACCGTGGAACAGGCCCTGGCCGGGTTCAGCGACCCCAACgtg                                                          <  1:154229/136‑1 (MQ=255)                                                    ggCCATCGTGATCATCCTGTTCGTCATCAACCGCCCACGCATGGACGTGGTCGCGCTGCTGGTGAATCTCGCCCTGCCACTACTCGGCATCCTTACCGTGGAACAGGCCCTGGCCGGGTTCAGCGACCCCAACgtg                                                          <  1:154226/136‑1 (MQ=255)                                                      ccATCGTGATCATCCTGTTCGTCATCAACCGCCCACGCATGGACGTGGTCGCGCTGCTGGTGATTCTCGCCCTGCCACTACTCGGCATCCTTACCGTGGAACAGGCCCTGGCCGGGTTCAGCGACCCCAACgtggt                                                        >  2:339872/1‑136 (MQ=255)                                                                          gTCATCAACCGCCCACGCATGGACGTGGTCGCGCTGCTGGTGATTCTCGCCCTGCCACTACTCGGCATCCTTACCGTGGAACAGGCCCTGGCCGGGTTCAGCGACCCCAACGTGGTGCTGATCGCCGCGCTGTTCg                                    >  1:236747/1‑136 (MQ=255)                                                                          gTCATCAACCGCCCACGCATGGACGTGGTCGCGCTGCTGGTGATTCTCGCCCTGCCACTACTCGGCATCCTTACCGTGGAACAGGCCCTGGCCGGGTTCAGCGACCCCAACGTGGTGCTGATCGCCGCGCTGTTCg                                    >  1:212764/1‑136 (MQ=255)                                                                            catcaACCGCCCACGCATGGACGTGGTCGCGCTGCTGGTGATTCTCGCCCTGCCACTACTCGGCATCCTTACCGTGGAACAGGCCCTGGCCGGGTTCAGCGACCCCAACGTGGTGCTGATCGCCGCGCTGTTCGTg                                  >  2:87077/1‑136 (MQ=255)                                                                               caACCGCCCACGCATGGACGTGGTCGCGCTGCTGGTGATTCTCGCCCTGCCACTACTCGGCATCCTTACCGTGGAACAGGCCCTGGCCGGGTTCAGCGACCCCAACGTGGTGCTGATCGCCGCGCTGTTCGTGAtt                               <  1:410189/136‑1 (MQ=255)                                                                                        aCGCATGGACGTGGTCGCGCTGCTGGTGATTCTCGCCCTGCCACTACTCGGCATCCTTACCGTGGAACAGGCCCTGGCCGGGTTCAGCGACCCCAACGTGGTGCTGATCGCCGCGCTGTTCGTGATTGGCGAAGGc                      >  1:228792/1‑136 (MQ=255)                                                                                             tGGACGTGGTCGCGCTGCTGGTGATTCTCGCCCTGCCACTACTCGGCATCCTTACCGTGGAACAGGCCCTGGCCGGGTTCAGCGACCCCAACGTGGTGCTGATCGCCGCGCTGTTCGTGATTGGCGAAGGCCTGGt                 <  1:269106/136‑1 (MQ=255)                                                                                                            tgctgGTGATTCTCGCCCTGCCACTACTCGGCATCCTTACCGTGGAACAGGCCCTGGCCGGGTTCAGCGACCCCAACGTGGTGCTGATCGCCGCGCTGTTCGTGATTGGCGAAGGCCTGGTGCGCACCGGCATTGc  >  1:69513/1‑136 (MQ=255)                                                                                                                     |                                                                                                                               GGAACACCCCCTCCCCATGAACCATGAGTTGCTCTGGGTCCTCGGCCTGCTGGCCATCGTGATCATCCTGTTCGTCATCAACCGCCCACGCATGGACGTGGTCGCGCTGCTGGTGATTCTCGCCCTGCCACTACTCGGCATCCTTACCGTGGAACAGGCCCTGGCCGGGTTCAGCGACCCCAACGTGGTGCTGATCGCCGCGCTGTTCGTGATTGGCGAAGGCCTGGTGCGCACCGGCATTGC  >  NC_002947/4439851‑4440093

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