Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,859,191	T→C	D309G (GAT→GGT)	yidK ←	putative transporter

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,859,191	0	T	C	92.1%	115.4 / 0.1	38	D309G (GAT→GGT)	yidK	putative transporter
Reads supporting (aligned to +/- strand):  ref base T (1/2);  new base C (18/17);  total (19/19)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.95e-01

CACCGCACCAAATAACACCGCGCCGAAGAAACCCACCATTGGCACTGGCAGAACGTTATTGACCAGCGTCGGGTAGGCCATGTCGGCTTTCGGTAAATCCTGATACAGATGAAACGCAATCAACCCTGGCAGCACCAGTACCAGCGGGTCGAGCATTTTCAGCACCGCCGTTAACAGCGCCCCCTTTTGCCCTT  >  NC_000913/3859094‑3859287                                                                                                  |                                                                                                 cACCGCACCAAATAACACCGCGCCGAAGAAACCCACCATTGGCACTGGCAGAACGTTATTGACCAGCGTCGGGTAGGCCATGTCGGCTTTCGGTAAAccc                                                                                                >  2:1135055/1‑100 (MQ=255)  aCCGCACCAAATAACACCGCGCCGAAGAAACCCACCATTGGCACTGGCAGAACGTTATTGACCAGCGTCGGGTAGGCCATGTCGGCTTTCGGTAAATCCt                                                                                               >  2:1207960/1‑100 (MQ=255)   ccGCACCAAATAACACCGCGCCGAAGAAACCCACCATTGGCACTGGCAGAACGTTATTGACCAGCGTCGGGTAGGCCATGTCGGCTTTCGGTAAACCCTg                                                                                              >  2:816255/1‑100 (MQ=255)     gCACCAAATAACACCGCGCCGAAGAAACCCACCATTGGCACTGGCAGAACGTTATTGACCAGCGTCGGGTAGGCCATGTCGGCTTTCGGTAAACCCTGAt                                                                                            >  2:771741/1‑100 (MQ=255)     gCACCAAATAACACCGCGCCGAAGAAACCCACCATTGGCACTGGCAGAACGTTATTGACCAGCGTCGGGTAGGCCATGTCGGCTTTCGGTAAACCCTGAt                                                                                            <  1:675648/100‑1 (MQ=255)      cACCAAATAACACCGCGCCGAAGAAACCCACCATTGGCACTGGCAGAACGTTATTGACCAGCGTCGGGTAGGCCATGTCGGCTTTCGGTAAACCCTGATa                                                                                           >  2:299376/1‑100 (MQ=255)       aCCAAATAACACCGCGCCGAAGAAACCCACCATTGGCACTGGCAGAACGTTATTGACCAGCGTCGGGTAGGCCATGTCGGCTTTCGGTAAACCCTGATAc                                                                                          <  1:699434/100‑1 (MQ=255)            aTAACACCGCGCCGAAGAAACCCACCATTGGCACTGGCAGAACGTTATTGACCAGCGTCGGGTAGGCCATGTCGGCTTTCGGTAAACCCTGATACAGATg                                                                                     <  2:215720/100‑1 (MQ=255)                    gcgcCGAAGAAACCCACCATTGGCACTGGCAGAACGTTATTGACCAGCGTCGGGTAGGCCATGTCGGCTTTCGGTAAATCCTGATACAGATGAAACGCaa                                                                             <  1:497392/100‑1 (MQ=255)                     cgcCGAAGAAACCCACCATTGGCACTGGCAGAACGTTATTGACCAGCGTCGGGTAGGCCATGTCGGCTTTCGGTAAACCCTGATACAGATGAAACGCAAt                                                                            <  2:862942/100‑1 (MQ=255)                          aagaaACCCACCATTGGCACTGGCAGAACGTTATTGACCAGCGTCGGGTAGGCCATGTCGGCTTTCGGTAAATCCTGATACAGATGAAACGCAATCAAcc                                                                       <  1:177864/100‑1 (MQ=255)                                cccACCATTGGCACTGGCAGAACGTTATTGACCAGCGTCGGGTAGGCCATGTCGGCTTTCGGTAAACCCTGATACAGATGAAACGCAATCAACCCTGgca                                                                 >  1:484407/1‑100 (MQ=255)                                     caTTGGCACTGGCAGAACGTTATTGACCAGCGTCGGGTAGGCCATGTCGGCTTTCGGTAAACCCTGATACAGATGAAACGCAATCAACCCTGGCAGCAcc                                                            <  1:1195115/100‑1 (MQ=255)                                     caTTGGCACTGGCAGAACGTTATTGACCAGCGTCGGGTAGGCCATGTCGGCTTTCGGTAAACCCTGATACAGATGAAACGCAATCAACCCTGGCAGCAcc                                                            <  2:557981/100‑1 (MQ=255)                                     caTTGGCACTGGCAGAACGTTATTGACCAGCGTCGGGTAGGCCATGTCGGCTTTCGGTAAACCCTGATACAGATGAAACGCAATCAACCCTGGCAGCAcc                                                            <  1:944702/100‑1 (MQ=255)                                         ggCACTGGCAGAACGTTATTGACCAGCGTCGGGTAGGCCATGTCGGCTTTCGGTAAACCCTGATACAGATGAAACGCAATCAACCCTGGCAGCACCAGTa                                                        <  1:1082199/100‑1 (MQ=255)                                         ggCACTGGCAGAACGTTATTGACCAGCGTCGGGTAGGCCATGTCGGCTTTCGGTAAACCCTGATACAGATGAAACGCAATCAACCCTGGCAGCACCAGTa                                                        <  1:1295233/100‑1 (MQ=255)                                         gacACTGGCAGAACGTTATTGACCAGCGTCGGGTAGGCCATGTCGGCTTTCGGTAAACCCTGATACAGATGAAACGCAATCAACCCTGGCAGCACCAGTa                                                        >  1:878323/3‑100 (MQ=255)                                          gCACTGGCAGAACGTTATTGACCAGCGTCGGGTAGGCCATGTCGGCTTTCGGTAAACCCTGATACAGATGAAACGCAATCAACCCTGGCAGCACCAGTAc                                                       >  1:1264394/1‑100 (MQ=255)                                             cTGGCAGAACGTTATTGACCAGCGTCGGGTAGGCCATGTCGGCTTTCGGTAAACCCTGATACAGATGAAACGCAATCAACCCTGGCAGCACCAGTACCAg                                                    >  2:271683/1‑100 (MQ=255)                                                 cAGAACGTTATTGACCAGCGTCGGGTAGGCCATGTCGGCTTTCGGTAAACCCTGATACAGATGAAACGCAATCAACCCTGGCAGCACCAGTACCAGCggg                                                >  1:205558/1‑100 (MQ=255)                                                  aGAACGTTATTGACCAGCGTCGGGTAGGCGATGTCGGCTTTCGGTAAACCCTGATACAGATGAAACGCAATCAACCCTGGCAGCACCAGTACCAGCggg                                                <  2:1192407/99‑1 (MQ=255)                                                  aGAACGTTATTGACCAGCGTCGGGTAGGCCATGTCGGCTTTCGGTAAACCCTGATACAGATGAAACGCAATCAACCCTGGCAGCACCAGTACCAGCGGGt                                               <  2:1183288/100‑1 (MQ=255)                                                             gACCAGCGTCGGGTAGGCCATGTCGGCTTTCGGTAAACCCTGATACAGATGAAACGCAATCAACCCTGGCAGCACCAGTACCAGCGGGTCGAGCATTTTc                                    >  2:1138428/1‑100 (MQ=255)                                                             gACCAGCGTCGGGTAGGCCATGTCGGCTTTCGGTAAACCCTGATACAGATGAAACGCAATCAACCCTGGCAGCACCAGTACCAGCGGGTCGAGCATTTTc                                    >  2:237244/1‑100 (MQ=255)                                                                         gTAGGCCATGTCGGCTTTCGGTAAACCCTGATACAGATGAAACGCAATCAACCCTGGCAGCACCAGTACCAGCGGGTCGAGCATTTTCAGCACCGCCGtt                        >  1:279963/1‑100 (MQ=255)                                                                          tAGGCCATGTCGGCTTTCGGTAAACCCTGATACAGATGAAACGCAATCAACCCTGGCAGCACCAGTACCAGCGGGTCGAGCATTTTCAGCACCGCCGTTa                       >  2:1297852/1‑100 (MQ=255)                                                                           aGGCCATGTCGGCTTTCGGTAAACCCTGATACAGATGAAACGCAATCAACCCTGGCAGCACCAGTACCAGCGGGTCGAGCATTTTCAGCACCGCCGTTaa                      >  2:350021/1‑100 (MQ=255)                                                                             gCCATGTCGGCTTTCGGTAAACCCTGATACAGATGAAACGCAATCAACCCTGGCAGCACCAGTACCAGCGGGTCGAGCATTTTCAGCACCGCCGTTAACa                    >  1:89859/1‑100 (MQ=255)                                                                                aTGTCGGCTTTCGGTAAACCCTGATACAGATGAAACGCAATCAACCCTGGCAGCACCAGTACCAGCGGGTCGAGCATTTTCAGCACCGCCGTTAACAgcg                 <  1:335114/100‑1 (MQ=255)                                                                                aTGTCGGCTTTCGGTAAACCCTGATACAGATGAAACGCAATCAACCCTGGCAGCACCAGTACCAGCGGGTCGAGCATTTTCAGCACCGCCGTTAACAgcg                 >  1:100164/1‑100 (MQ=255)                                                                                 tGTCGGCTTTCGGTAAACCCTGATACAGATGAAACGCAATCAACCCTGGCAGCACCAGTACCAGCGGGTCGAGCATTTTCAGCACCGCCGTTAACAgcgc                >  1:609745/1‑100 (MQ=255)                                                                                  gTCGGCTTTCGGTAAACCCTGATACAGATGAAACGCAATCAACCCTGGCAGCACCAGTACCAGCGGGTCGAGCATTTTCAGCACCGCCGTTAACAGCGcc               <  1:229578/100‑1 (MQ=255)                                                                                    cGGCTTTCGGTAAACCCTGATACAGATGAAACGCAATCAACCCTGGCAGCACCAGTACCAGCGGGTCGAGCATTTTCAGCACCGCCGTTAACAGCGcccc             >  2:336553/1‑100 (MQ=255)                                                                                          tCGGTAAACCCTGATACAGATGAAACGCAATCAACCCTGGCAGCACCAGTACCAGCGGGTCGAGCATTTTCAGCACCGCCGTTAACAGCGCCCCCTTTTg       <  1:816261/100‑1 (MQ=255)                                                                                          tCGGTAAACCCTGATACAGATGAAACGCAATCAACCCTGGCAGCACCAGTACCAGCGGGTCGAGCATTTTCAGCACCGCCGTTAACAGCGCCCCCTTTTg       <  2:878316/100‑1 (MQ=255)                                                                                           cGGTAAACCCTGATACAGATGAAACGCAATCAACCCTGGCAGCACCAGTACCAGCGGGTCGAGCATTTTCAGCACCGCCGTTAACAGCGCCCCCTTTTGc      <  1:62658/100‑1 (MQ=255)                                                                                            ggTAAACCCTGATACAGATGAAACGCAATCAACCCTGGCAGCACCAGTACAAGCGGGTCGAGCATTTTCAGCACCGCCGTTCACAGCGCCCCCTTTTGcc     >  2:1258928/1‑100 (MQ=255)                                                                                               aaaCCCTGATACAGATGAAACGCAATCAACCCTGGCAGCACCAGTACCAGCGGGTCGAGCATTTTCAGCACCGCCGTTAACAGCGCCCCCTTTTGCCCtt  <  1:979770/100‑1 (MQ=255)                                                                                                  |                                                                                                 CACCGCACCAAATAACACCGCGCCGAAGAAACCCACCATTGGCACTGGCAGAACGTTATTGACCAGCGTCGGGTAGGCCATGTCGGCTTTCGGTAAATCCTGATACAGATGAAACGCAATCAACCCTGGCAGCACCAGTACCAGCGGGTCGAGCATTTTCAGCACCGCCGTTAACAGCGCCCCCTTTTGCCCTT  >  NC_000913/3859094‑3859287

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