Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	1,215,064	G→A	20.1%	L291L (CTG→CTA)	yebQ →	predicted transporter

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	1,215,064	0	G	A	20.1%	97.5 / 19.3	45	L291L (CTG→CTA)	yebQ	predicted transporter
Reads supporting (aligned to +/- strand):  ref base G (22/13);  new base A (3/6);  total (25/20)
Fisher's exact test for biased strand distribution p-value = 1.44e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TTCACTTTCTATTTGCACATCTGTTTGCTCTTTCTGCGCACAAATGCTGGCAATGGTTTCCCTGCCCTTTTACCTGCAAACCGTGCTCGGGCGTAGTGAAGTCGAAACAGGTTTACTTCTGACAC  >  minE/1215001‑1215125                                                                |                                                              ttCACTTTCTATTTGCACATCTGTTTGCTCTTTCTGCGCACAAATGCTGGCAATGGTTTCCGTGCCCttt                                                         >  1:1355796/1‑70 (MQ=255) ttCACTTTCTATTTGCACATCTGTTTGCTCTTTCTGCGCACAAATGCTGGCAATGGTTTCCCTGCCCttt                                                         >  1:731684/1‑70 (MQ=255) ttCACTTTCTATTTGCACATCTGTTTGCTCTTTCTGCGCACAAATGCTGGCAATGGTTTCCCTGCCCttt                                                         >  1:382132/1‑70 (MQ=255)  tCACTTTCTATTTGCACATCTGTTTGCTCTTTCTGCGCACAAATGCTGGCAATGGTTTCCCTGCCCtttt                                                        <  1:672496/70‑1 (MQ=255)       ttCTATTTGCACATCTGTTTGCTCTTTCTGCGCACAAATGCTGGCAATTGTTTCCCTGCCCTTTTACCTg                                                   >  1:754266/1‑70 (MQ=255)       ttCTATTTGCACATCTGTTTGCTCTTTCTGCGCACAAATGCTGGCAATGGTTTCCCTACCCTTTTAc                                                      >  1:1778664/1‑67 (MQ=255)        tCTATTTGCACATCTGTTTGCTCTTTCTGCGCACAAATGCTGGCAATGGTTTCCCTACCCTTTTACCTGc                                                  <  1:946147/70‑1 (MQ=255)         cTATTTGCACATCTGTTTGCTCTTTCTGCGCACAAATGCTGGCAATGGTTTCCCTGCCCTTTTACCTGCa                                                 <  1:273887/70‑1 (MQ=255)         cTATTTGCACATCTGTTTGCTCTTTCTGCGCACAAATGCTGGCAATGGTTTCCCTGCCCTTTTACCTGCa                                                 >  1:2082165/1‑70 (MQ=255)            tttGCACATCTGTTTGCTCTTTCTGCGCACAAATGCTGGCAATGGTTTCCCTGCCCTTTTACCTGCAAAcc                                             <  1:1569971/71‑1 (MQ=255)            tttGCACATCTGTTTGCTCTTTCTGCGCACAAATGCTGGCAATGGTTTCCCTACCCTTTTACCTGCaa                                                <  1:1192610/68‑1 (MQ=255)            tttGCACATCTGTTTGCTCTTTCTGCGCACAAATGCTGGCAATGGTTTCCC‑‑CCCTTTTACCTGCAAACCg                                            <  1:1038620/70‑1 (MQ=255)               gCACATCTGTTTGCTCTTTCTGCGCACAAATGCTGGCAATGGTTTCCCTGCCCTTTTACCTGCAAACCGt                                           >  1:1839675/1‑70 (MQ=255)                cacaTCTGTTTGCTCTTTCTGCGCACAAATGCTGGCAATGGTTTCCCTGCCCTTTTACCTGCAAACCGTGc                                         <  1:18192/71‑1 (MQ=255)                 acaTCTGTTTGCTCTTTCTGCGCACAAATGCTGGCAATGGTTTCCCTGCCCTTTTACCTGCAAACCGTGCt                                        <  1:2099862/71‑1 (MQ=255)                     cTGTTTGCTCTTTCTGCGCACAAATGCTGGCAATGGTTTCCCTGCCCTTTTACCTGCAAACCGTGCTCggg                                    >  1:473142/1‑71 (MQ=255)                     cTGTTTGCTCTTTCTGCGCACAAATGCTGGCAATGGTTTCCCTGCCCTTTTACCTGCAAACCGTGCTCg                                      >  1:272987/1‑69 (MQ=255)                      tGTTTGCTCTTTCTGCGCACAAATGCTGGCAATGGTTTCCCTGCCCTTTTACCTGCAAACCGTGCTCGGGc                                   >  1:1516983/1‑71 (MQ=255)                        tttGCTCTTTCTGCGCACAAATGCTGGCAATGGTTTCCCTGc                                                              <  1:306162/42‑1 (MQ=255)                          tGCTCTTTCTGCGCACAAATGCTGGCAATGGTTTCCCTGCCCTTTTACCTGCAAACCGTGCTCGGGCGTa                                >  1:958363/1‑70 (MQ=255)                             tctTTCTGCGCACAAATGCTGGCAATGGTTTCCCTGCCCTTTTACCTGCAAACCGTGCTc                                       <  1:1850657/60‑1 (MQ=255)                              ctTTCTGCGCACAAATGCTGGCAATGGTTTCCCTGCCCTTTTACCTGCAAACCGTGCTCGGGCGTAGTg                             >  1:481792/1‑69 (MQ=255)                               tttCTGCGCACAAATGCTGGCAATGGTTTCCCTGCCCttt                                                         <  1:536175/40‑1 (MQ=255)                               tttCTGCGCACAAATGCTGCCAATGGTTTCCCTGCCCttt                                                         <  1:1479827/40‑1 (MQ=255)                                ttCTGCGCACAAATGCTGGCAATGGTTTCCCTGCCCTTTTACCTGCAAACCGTGCTCGGGCGTAGTGAAGt                         >  1:1056522/1‑71 (MQ=255)                                 tCTGCGCACAAATGCTGGCAATGGTTTCCCTGCCCTTTTAc                                                      >  1:383264/1‑41 (MQ=255)                                 tCTGCGCACAAATGCTGGCAATGGTTTCCCTGCCCTTTTACCTGCaa                                                >  1:169568/1‑47 (MQ=255)                                   tGCGCACAAATGCTGGCAATGGTTTCCCTGCCCTTTTACCTGCAAACCGTGCTCGGGCGTAGTGAAGTCg                       >  1:1070564/1‑70 (MQ=255)                                     cgcACAAATGCTGGCAATGGTTTCCCTGCCCTTTTACCTGCAAACCGTGCTCGGGCGTAGTGAAGTCGa                      <  1:1517559/69‑1 (MQ=255)                                     cgcACAAATGCTGGCAATGGTTTCCCTGCCCTTTTACCTGCAAACCGTGCTCGGGCGTAGTGAAGTCGa                      <  1:26413/69‑1 (MQ=255)                                     cgcACAAATGCTGGCAATGGTTTCCCTACCCTTTTACCTGCAAACCGTGCTCGGGCGTAGTGAAGTCGaa                     >  1:851798/1‑70 (MQ=255)                                     cgcACAAATGCTGGCAATGGTTTCCCTACCCTTTTACCTGCAAACCGTGCTCGGGCGTAGTGAAGTCGaa                     <  1:1278763/70‑1 (MQ=255)                                        acaAATGCTGGCAATGGTTTCCCTGCCCTTTTACCTGCAAACCGTGCTCGGGCGTAGTGAAGTCGAAACAg                 >  1:439276/1‑71 (MQ=255)                                         caAATGCTGGCAATGGTTTCCCTGCCCTTTTACCTGCAAACCGTGc                                         <  1:449546/46‑1 (MQ=255)                                         caAATGCTGGCAATGGTTTCCCTGCCCTTTTACCTGCAAACCGTGCTCGGGCGTAGTGAAGTCGAAACAgg                >  1:987132/1‑71 (MQ=255)                                         caAATGCTGGCAATGGTTTCCCTACCCTTTTACCTGCAAACCGTGCTCGGGCGTAGt                              <  1:1008783/57‑1 (MQ=255)                                           aaTGCTGGCAATGGTTTCCCTACGCTTTTACCTGCAAACCGTGCTCGGGCGTAGTGAAGTCGAAACAGGtt              <  1:1787865/71‑1 (MQ=255)                                                  gCAATGGTTTCCCTACCCTTTTACCTGCAAACCGTGCTCGGGCGTAGTGAAGTCGAAACAGGTTTACTTCt       <  1:528337/71‑1 (MQ=255)                                                    aaTGGTTTCCCTGCCCTTTTACCTGCAAACCGTGCTCGGGCGTAGTGAAGTCGAAACAGGTTTACTTCTGa     >  1:453560/1‑71 (MQ=255)                                                     aTGGTTTCCCTGCCCTTTTACCTGCAAACCGTGCTCGGGCGTAGTGAAGTCGAAACAGGTTTACTTCTGac    <  1:92987/71‑1 (MQ=255)                                                       ggTTTCCCTGCCCTTTTACCTGCAAACCGTGCTCGGGCGTAGTGAAGTCGAAACAGGTTTACTTCTGacac  >  1:620877/1‑71 (MQ=255)                                                       ggTTTCCCTGCCCTTTTACCTGCAAACCGTGCTCGGGCGTAGTGAAGTCGAAACAGGTTTACTTCTGacac  >  1:83164/1‑71 (MQ=255)                                                       ggTTTCCCTACCCTTTTACCTGCAAACCGTGCTCGGGCGTAGTGAAGTCGAAACAGGTTTACTTCTGacac  >  1:1290405/1‑71 (MQ=255)                                                        gTTTCCCTGCCCTTTTACCTGCAAACCGTGCTCGGGCGTAGTGAAGTCGAAACAGGTTTAc           >  1:694297/1‑61 (MQ=255)                                                               tGCCCTTTTACCTGCAAACCGTGCTCGGGCGTAGTGAAGTCGAAAc                   >  1:794580/1‑46 (MQ=255)                                                                |                                                              TTCACTTTCTATTTGCACATCTGTTTGCTCTTTCTGCGCACAAATGCTGGCAATGGTTTCCCTGCCCTTTTACCTGCAAACCGTGCTCGGGCGTAGTGAAGTCGAAACAGGTTTACTTCTGACAC  >  minE/1215001‑1215125

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