Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	W3110S.gb	3,780,496	G→A	76.6%	G66R (GGG→AGG)	yidK →	predicted transporter

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	W3110S.gb	3,780,496	0	G	A	76.6%	87.5 / 25.0	47	G66R (GGG→AGG)	yidK	predicted transporter
Reads supporting (aligned to +/- strand):  ref base G (7/4);  new base A (9/27);  total (16/31)
Fisher's exact test for biased strand distribution p-value = 2.92e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TATTGCCGCTTCGTTAATGCTAACCAACCTTTCCACGGAACAACTGGTCGGCCTTTCCGGGCAGGCCTACAAAAGCGGCATGTCGGTGATGGGCTGGGAAGTGACTTCAGCGGTGACGCTG  >  W3110S.gb/3780438‑3780558                                                           |                                                               tATTGCCGCTTCGTTAATGCTAACCAACCTTTCCACGGAACAACTGGTCGGCCTTTCCGGGCAGGCCTAc                                                     >  1:2061196/1‑70 (MQ=255) tATTGCCGCTTCGTTAATGCTAACCAACCTTTCCACGGAACAACTGGTCGGCCTTTCCGGGCAGGCCTAc                                                     >  1:659255/1‑70 (MQ=255) tATTGCCGCTTCGTTAATGCTAACCAACCTTTCCACGGAACAACTGGTCGGCCTTTCCGGGCAGGCCTAc                                                     >  1:1806513/1‑70 (MQ=255) tATTGCCGCTTCGTTAATGCTAACCAACCTTTCCACGGAACAACTGGTCGGCCTTTCCAGGCAGGCCTAc                                                     >  1:3305146/1‑70 (MQ=255)   ttGCCGCTTCGTTAATGCTAACCAACCTTTCCACGGAACAACTGGTCGGCCTTTCCGGGCAGGCCTACaaa                                                  <  1:1735425/71‑1 (MQ=255)   ttGCCGCTTCGTTAATGCTAACCAACCTTTCCACGGAACAACTGGTCGGCCTTTCCAGGCAGGCCTACaa                                                   <  1:3183355/70‑1 (MQ=255)   ttGCCGCTTCGTTAATGCTAACCAACCTTTCCACGGAACAACTGGTCGGCCTTTCCAGGCAGGCCTACaa                                                   <  1:1757815/70‑1 (MQ=255)     gCCGCTTCGTTAATGCTAACCAACCTTTCCACGGAACAACTGGTCGGCCTTTCCAGGCAGGCCTACaa                                                   >  1:2193674/1‑68 (MQ=255)        gCTTCGTTAATGCTAACCAACCTTTCCACGGAACAACTGGTCGGCCTTTCCGGGCAGGCCTACAAAAGCgg                                             <  1:510248/71‑1 (MQ=255)        gCTTCGTTAATGCTAACCAACCTTTCCACGGAACAACTGGTCGGCCTTTCCGGGCAGGCCTACAAAAGCgg                                             <  1:1466383/71‑1 (MQ=255)        gCTTCGTTAATGCTAACCAACCTTTCCACGGAACAACTGGTCGGCCTTTCCAGGCAGGCCTACAAAAGCgg                                             <  1:1286687/71‑1 (MQ=255)        gCTTCGTTAATGCTAACCAACCTTTCCACGGAACAACTGGTCGGCCTTTCCAGGCAGGCCTACAAAAGCgg                                             <  1:227921/71‑1 (MQ=255)        gCTTCGTTAATGCTAACCAACCTTTCCACGGAACAACTGGTCGGCCTTTCCAGGCAGGCCTACAAAAGCgg                                             <  1:2823819/71‑1 (MQ=255)        gCTTCGTTAATGCTAACCAACCTTTCCACGGAACAACTGGTCGGCCTTTCCAGGCAGGCCTACAAAAGCgc                                             <  1:1285518/71‑2 (MQ=255)         cTTCGTTAATGCTAACCAACCTTTCCACGGAACAACTGGTCGGCCTTTCCAGGCAGGCCTACAAAAGCGGc                                            <  1:2286234/71‑1 (MQ=255)         cTTCGTTAATGCTAACCAACCTTTCCACGGAACAACTGGTCGGCCTTTCCAGGCAGGCCTACAAAAGCGGc                                            <  1:1849091/71‑1 (MQ=255)           tCGTTAATGCTAACCAACCTTTCCACGGAACAACTGGTCGGCCTTTCCAGGCAGGCCTACAAAAGCGGCAt                                          >  1:3579067/1‑71 (MQ=255)             gTTAATGCTAACCAACCTTTCCACGGAACAACTGGTCGGCCTTTCCAGGCAGGCCTACAAAAGCGGCATGt                                        >  1:1162179/1‑71 (MQ=255)             gTTAATGCTAACCAACCTTTCCACGGAACAACTGGTCGGCCTTTCCAGGCAGGCCTACAAAAGCGGCATGt                                        >  1:388024/1‑71 (MQ=255)              ttAATGCTAACCAACCTTTCCACGGAACAACTGGTCGGCCTTTCCAGGCAGGCCTACAAAAGCGGCATGTc                                       <  1:549055/71‑1 (MQ=255)                 aTGCTAACCAACCTTTCCACGGAACAACTGGTCGGCCTTTCCGGGCAGGCCTACAAAAGCGGCATGTCGGt                                    <  1:1824458/71‑1 (MQ=255)                 aTGCTAACCAACCTTTCCACGGAACAACTGGTCGGCCTTTCCAGGCAGGCCTACAAAAGCGGCATGTCGGt                                    <  1:243090/71‑1 (MQ=255)                 aTGCTAACCAACCTTTCCACGGAACAACTGGTCGGCCTTTCCAGGCAGGCCTACAAAAGCGGCATGTCGGt                                    <  1:1100580/71‑1 (MQ=255)                  tGCTAACCAACCTTTCCACGGAACAACTGGTCGGCCTTTCCAGGCATGCCTACAAAAGCGGCATGTCGGt                                    >  1:1288599/1‑70 (MQ=255)                  tGCTAACCAACCTTTCCACGGAACAACTGGTCGGCCTTTCCAGGCAGGCCTACAAAAGCGGCATGTCGGt                                    <  1:1212906/70‑1 (MQ=255)                  tGCTAACCAACCTTTCCACGGAACAACTGGTCGGCCTTTCCAGGCAGGCCTACAAAAGCGGCATGTCGGt                                    <  1:2843228/70‑1 (MQ=255)                  tGCTAACCAACCTTTCCACGGAACAACTGGTCGGCCTTTCCAGGCAGGCCTACAAAAGCGGCATGTCGGTg                                   <  1:3172667/71‑1 (MQ=255)                  tGCTAACCAACCTTTCCACGGAACAACTGGGCGGCCTTTCCAGGCAGGCCTACAAAAGCGGCATGTCGGt                                    <  1:918335/70‑1 (MQ=255)                    cTAACCAACCTTTCCACGGAACAACTGGTCGGCCTTTCCAGGCAGGCCTACAAAAGCGGCATGTCGGTGa                                  <  1:2918251/70‑1 (MQ=255)                    cTAACCAACCTTTCCACGGAACAACTGGTCGGCCTTTCCAGGCAGGCCTACAAAAGCGGCATGTCGGTGa                                  <  1:3144386/70‑1 (MQ=255)                       accaaccTTTCCACGGAACAACTGGTCGGCCTTTCCAGGCAGGCCTACAAAAGCGGCATGTCGGTGATg                                <  1:2674432/69‑1 (MQ=255)                          aaccTTTCCACGGAACAACTGGTCGGCCTTTCCAGGCAGGCCTACAAAAGCGGCa                                           >  1:2146985/1‑55 (MQ=255)                          aaccTTTCCACGGAACAACTGGTCGGCCTTTCCAGGCAGGCCTACAAAAGCGGCa                                           >  1:2588035/1‑55 (MQ=255)                              tttCCACGGAACAACTGGTCGGCCTTTCCGGGCAGGCCTACAAAAGCGGCATGTCGGTGATGGGCTggg                         >  1:3184810/1‑69 (MQ=255)                              tttCCACGGAACAACTGGTCGGCCTTTCCGGGCAGGCCTACAAAAGCGGCATGTCGGTGATGGGCTggg                         >  1:1627270/1‑69 (MQ=255)                              tttCCACGGAACAACTGGTCGGCCTTTCCGGGCAGGCCTACAAAAGCGGCATGTCGGTGATGGGCTggg                         >  1:1580377/1‑69 (MQ=255)                              tttCCACGGAACAACTGGTCGGCCTTTCCGGGCAGGCCTACAAAAGCGGCATGTCGGTGATGGGCTgg                          >  1:819698/1‑68 (MQ=255)                              tttCCACGGAACAACTGGTCGGCCTTTCCAGGCAGGCCTACAAAAGCGGCATGTCGGTGATGGGCTggg                         >  1:397632/1‑69 (MQ=255)                                  cACGGAACAACTGGTCGGCCTTTCCAGGCAGGCCTACAAAAGCGGCATGTCGGTGATGGGCTGGGAAGTGa                   <  1:1862068/71‑1 (MQ=255)                                  cACGGAACAACTGGTCGGCCTTTCCAGGCAGGCCTACAAAAGCGGCATGTCGGTGATGGGCTGGGAAGTGa                   <  1:3355301/71‑1 (MQ=255)                                         caacTTGTCGGCCTTTCCAGGCAGGCCTACAAAAGCGGCATGTCGGt                                    <  1:1544652/47‑1 (MQ=39)                                         caacTGGTCGGCCTTTCCAGGCAGGCCTACAAAAGCGGCATGTCGGt                                    <  1:3280739/47‑1 (MQ=255)                                         caacTGGTCGGCCTTTCCAGGCAGGCCTACAAAAGCGGCATGTCGGt                                    <  1:3522192/47‑1 (MQ=255)                                         caacTGGTCGGCCTTTCCAGGCAGGCCTACAAAAGCGGCATGTCGGt                                    <  1:3172588/47‑1 (MQ=255)                                                   gCCTTTCCAGGCAGGCCTACAAAAGCGGCATGTCgg                                     <  1:463530/36‑1 (MQ=255)                                                   gCCTTTCCAGGCAGGCCTACAAAAGCGGCATGTCGGTGATGGGCTGGGAAGTGACTTCAGCGGTGACGCTg  <  1:234125/71‑1 (MQ=255)                                                   gCCTTTCCAGGCAGGCCTACAAAAGCGGCATGTCGGTGATGGGCTGGGAAGTGACTTCAGCGGTGACGCTg  <  1:2323760/71‑1 (MQ=255)                                                           gggCAGGCCTACTAAAGCGGCATGTCGGTGATGGGCTGGGAAGTGACtt                <  1:1399709/49‑1 (MQ=255)                                                           |                                                               TATTGCCGCTTCGTTAATGCTAACCAACCTTTCCACGGAACAACTGGTCGGCCTTTCCGGGCAGGCCTACAAAAGCGGCATGTCGGTGATGGGCTGGGAAGTGACTTCAGCGGTGACGCTG  >  W3110S.gb/3780438‑3780558

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