Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	605,328	G→A	35.5%	P467P (CCG→CCA)	ycaM →	predicted transporter

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	605,328	0	G	A	35.5%	65.1 / 58.9	62	P467P (CCG→CCA)	ycaM	predicted transporter
Reads supporting (aligned to +/- strand):  ref base G (24/16);  new base A (13/9);  total (37/25)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.21e-01

CCTTCCAGTTGGCGCTGAATGTTGCAACGCCGTTTGTGCTGGTAGGATTAGGACTGATATTCCCGCTGCTGGCGCGTAAAGCGAATAGTAAATAATTATTGTGGTCGGTGACGCTGACTCTGGCAGATCTGA  >  minE/605264‑605395                                                                 |                                                                    ccTTCCAGTTGGCGCTGAATGTTGCAACGCCGTTTGTGCTGGTAGGATTAGGACTGATATTCCCActgctg                                                               >  1:1345470/1‑71 (MQ=255)    tCCAGTTGGCGCTGAATGTTGCAACGCCGTTTGTGCTGGTAGGATTAGGACTGATATTCCCgc                                                                    >  1:3130610/1‑63 (MQ=255)    tCCAGTTGGCGCTGAATGTTGCAACGCCGTTTGTGCTGGTAGGATTAGGACTGATATTCCCGCTGCTGgcg                                                            >  1:281247/1‑71 (MQ=255)    tCCAGTTGGCGCTGAATGTTGCAACGCCGTTTGTGCTGGTAGGATTAGGACTGATATTCCCACTGCTGgcg                                                            >  1:2790692/1‑71 (MQ=255)     ccAGTTGGCGCTGAATGTTGCAACGCCGTTTGTGCTGGTAGGATTAGGACTGATATTCCCGCTGCTGgcgc                                                           >  1:3971673/1‑71 (MQ=255)     ccAGTTGGCGCTGAATGTTGCAACGCCGTTTGTGCTGGTAGGATTAGGACTGATATTCCCACTGCTGgcgc                                                           >  1:2329104/1‑71 (MQ=255)     ccAGTTGGCGCTGAATGTTGCAACGCCGTTTGTGCTGGTAGGATTAGGACTGATATTCCCACTGCTGgcgc                                                           >  1:633136/1‑71 (MQ=255)         ttGGCGCTGAATGTTGCAACGCCGTTTGTGCTGGTAGGATTAGGACTGATATTCCCGCTGCTGGCGCGTaa                                                       <  1:1968794/71‑1 (MQ=255)         ttGGCGCTGAATGTTGCAACGCCGTTTGTGCTGGTAGGATTAGGACTGATATTCCCACTGCTGGCGCGTaa                                                       <  1:2183012/71‑1 (MQ=255)           ggCGCTGAATGTTGCAACGCCGTTTGTGCTGGTAGGATTAGGACTGATATTCCCGCTGCTGGCGCGTAAAg                                                     >  1:1739481/1‑71 (MQ=255)           ggCGCTGAATGTTGCAACGCCGTTTGTGCTGGTAGGATTAGGACTGATATTCCCACTGCTGGCGCGTaaa                                                      <  1:820524/70‑1 (MQ=255)            gcgcTGAATGTTGCAACGCCGTTTGTGCTGGTAGGATTAGGACTGATATTCCCGCTGCTGGCGCGTAAAGc                                                    <  1:2581463/71‑1 (MQ=255)               cTGAATGTTGCAACGCCGTTTGTGCTGGTAGGATTAGGACTGATATTCCCGCTGCTGGCGCGTAAAGCGaa                                                 <  1:1515830/71‑1 (MQ=255)                  aaTGTTGCAACGCCGTTTGTGCTGGTAGGATTAGGACTGATATTCCCGCTGCTGGCGCGTAAAGCGAATAg                                              <  1:4054938/71‑1 (MQ=255)                      ttGCAACGCCGTTTGTGCTGGTAGGATTAGGACTGATATTCCCGCTGCTGGCGCGTAAAGCGAATAGTaaa                                          >  1:66614/1‑71 (MQ=255)                      ttGCAACGCCGTTTGTGCTGGTAGGATTAGGACTGATATTCCCGCTGCTGGCGCGTAAAGCGAATAGTaaa                                          >  1:2369434/1‑71 (MQ=255)                        gCAACGCCGTTTGTGCTGGTAGGATTAGGACTGATATTCCCGCTGCTGGCGCGTAAAGCGa                                                  >  1:3684471/1‑61 (MQ=255)                        gCAACGCCGTTTGTGCTGGTAGGATTAGGACTGATATTCCCACTGCTGGCGCGTAAAGCGAATAGTAaat                                         <  1:237158/70‑1 (MQ=255)                         cAACGCCGTTTGTGCTGGTAGGATTAGGACTGATATTCCCACTGCTGGCGCGTAAAGCGAATAGTAaataa                                       <  1:904837/71‑1 (MQ=255)                          aaCGCCGTTTGTGCTGGTAGGATTAGGACTGATATTCCCACTGCTGGCGCGTAAAGCGAATAGTAaataat                                      <  1:1644257/71‑1 (MQ=255)                           aCGCCGTTTGTGCTGGTAGGATTAGGACTGATATTCCCGCTGCTGGCGCGTAAAGCGAATAGTAaataat                                      >  1:2550038/1‑70 (MQ=255)                             gCCGTTTGTGCTGGTAGGATTAGGACTGATATTCCCgctgc                                                                 >  1:215186/1‑41 (MQ=255)                             gCCGTTTGTGCTGGTAGGATTAGGACTGATATTCCCGCTGCTGGCGCGTAAAGCGAATAGTAAATAattat                                   <  1:2972520/71‑1 (MQ=255)                             gCCGTTTGTGCTGGTAGGATTAGGACTGATATTCCCGCTGCTGGCGCGTAAAGCGAATAGTAAATAattat                                   <  1:1744722/71‑1 (MQ=255)                              ccGTTTGTGCTGGTAGGATTAGGACTGATATTCCCACTGCTGGCGCGTAAAGCGAATAGTAaataa                                       >  1:1463982/1‑66 (MQ=255)                              ccGTTTGTGCTGGTAGGATTAGGACTGATATTCCCACTGCTGGCGCGTAAAGCGAATAGTAAATAattatt                                  <  1:527272/71‑1 (MQ=255)                                  ttGTGCTGGTAGGATTAGGACTGATATTCCCGCTGCTGGCGCGTAAAGc                                                    >  1:3152429/1‑49 (MQ=255)                                   tgtgCTGGTAGGATTAGGACTGATATTCCCGCTGCTGGCGCGTAAAGCGAATAGTAAATAATTATTGTGGt                             >  1:215048/1‑71 (MQ=255)                                     tgCTGGTAGGATTAGGACTGATATTCCCACTGCTGGCGCGTAAAGCGAATAGTAAATAATTATTGTGGTCg                           >  1:951926/1‑71 (MQ=255)                                       cTGGTAGGATTAGGACTGATATTCCCGCTGCTGGCGCGTAAAGCGAATAGTAAATAATTATTGTGGTCGGt                         >  1:2635510/1‑71 (MQ=255)                                            aGGATTAGGACTGATATTCCCGCTGCTGGCGCGTAAAGCGAATa                                               >  1:3878016/1‑44 (MQ=255)                                            aGGATTAGGACTGATATTCCCACTGCTGGCGCGTAAAGCGaa                                                 >  1:3457286/1‑42 (MQ=255)                                             ggATTAGGACTGATATTCCCGCTGCTGGCGCGTAAAGc                                                    <  1:3860888/38‑1 (MQ=255)                                             ggATTAGGACTGATATTCCCGCTGCTGGCGCGTAAAGCGAATAGTAAATAATTATTGTGGTCGGTGACGCt                   <  1:3463029/71‑1 (MQ=255)                                              gATTAGGACTGATATTCCCGCTGCTGGCGCGTAAAGCGAATAGTAAATAATTATTGTGGTCGGTg                        <  1:2774070/65‑1 (MQ=255)                                                ttAGGACTGATATTCCCGCTGCTGGCGCGTAAAGCGAATAGTAAATAATTATTGTGGTCGGTGACGCTGAc                <  1:2253940/71‑1 (MQ=255)                                                ttAGGACTGATATTCCCACTGCTGGCGCGTAAAGCGAATa                                               >  1:334216/1‑40 (MQ=255)                                                 tAGGACTGATATTCCCGCTGCTGGCGCGTAAAGCg                                                   <  1:4122571/35‑1 (MQ=255)                                                     aCTGATATTCCCGCTGCTGGCGCGTAAAGCGAATAGTAAATAATTATTGTGGTCGGTGACGCTGACTCTgg           <  1:4135720/71‑1 (MQ=255)                                                     aCTGATATTCCCGCTGCTGGCGCGTAAAGCGAATAGTAAATAATTATTGTGGTCGGTGACGCTGACTCTgg           <  1:4040762/71‑1 (MQ=255)                                                      cTGATATTCCCGCTGCTGGCGCGTAAAGCGAATAGTAAATAATTATTGTGGTCGGTGACGCTGActc              >  1:1811168/1‑67 (MQ=255)                                                       tGATATTCCCGCTGCTGGCGCGTAAAGCGAATAGTAAATAATTATtgtg                               <  1:2841074/49‑1 (MQ=255)                                                        gATATTCCCACTGCTGGCGCGTAAAGCGAATAGTAAATAATTATTGTGGTCGGTGACGCTGACTCTgg           >  1:2321688/1‑68 (MQ=255)                                                        gATATTCCAACTGCTGGCGCGTAAAGCGAATAGTAAATAATTATTGTGGTCGGAGACGCTGACTCTg            >  1:2641820/1‑67 (MQ=255)                                                         atatTCCCGCTGCTGGCGCGTAAAGCGAATAGTAAATAATTATTGTGGTCGGTGCCGCTGACTCTgg           >  1:1753738/1‑67 (MQ=255)                                                         atatTCCCGCTGCTGGCGCGTAAAGCGAATAGTAAATAATTATTGTGGTCGGTGACGCTGACTCTgg           >  1:4084729/1‑67 (MQ=255)                                                         atatTCCCGCTGCTGGCGCGTAAAGCGAATAGTAAATAATTATTGTGGTCGGTGACGCTGACTCTgg           >  1:2139544/1‑67 (MQ=255)                                                         atatTCCCGCTGCTGGCGCGTAAAGCGAATAGTAAATAATTATTGTGGTCGGTGACGCTGACTCTgg           >  1:713554/1‑67 (MQ=255)                                                         atatTCCCACTGCTGGCGCGTAAAGCGAATAGTAAATAATTATTGTGGTCGGTGACGCTGACTCTgg           >  1:910649/1‑67 (MQ=255)                                                         atatTCCCACTGCTGGCGCGTAAAGCGAATAGTAAATAATTATTGTGGTCGGTGACGCTGACTCTgg           >  1:957042/1‑67 (MQ=255)                                                         atatTCCCACTGCTGGCGCGTAAAGCGAATAGTAAATAATTATTGTGGTCGGTGACGCTGACTCTGGCa         <  1:2506452/69‑1 (MQ=255)                                                         atatTCCCACTGCTGGCGCGTAAAGCGAATAGTAAATAATTATTGTGGTCGGTGACGCTGACTCTGGCAg        <  1:1181687/70‑1 (MQ=255)                                                          tatTCCCGCTGCTGGCGCGTAAAGCGAATAGTAAATAatta                                    <  1:2406940/41‑1 (MQ=255)                                                          tatTCCCGCTGCTGGCGCGTAAAGCGAATAGTAAATAatta                                    <  1:788306/41‑1 (MQ=255)                                                          tatTCCCGCTGCTGGCGCGTAAAGCGAATAGTAAATAATTATTGTGGTCGGTGACGCTGACTCTgg           >  1:3224238/1‑66 (MQ=255)                                                           atTCCCGCTGCTGGCGCGTAAAGCGAATAGTAAATAATTATTGTGGTCGGTGACGCTGACTCTGGCAGATc     >  1:3102918/1‑71 (MQ=255)                                                           atTCCCACTGCTGTCGCGTAAAGCGAATAGTAaat                                         <  1:2045342/35‑1 (MQ=37)                                                             tCCCGCTGCTGGCGCGTAAAGCGAATAGTAAATAATTATTGTGGTCGGTGACGCTg                  >  1:18216/1‑56 (MQ=255)                                                              cccGCTGCTGGCGCGTAAAGCGAATAGTAAATAATTATtgt                                >  1:2579745/1‑41 (MQ=255)                                                              cccGCTGCTGGCGCGTAAAGCGAATAGTAAATAATTATTGTGGTCGGTGACGCTGACTCTGGCAGATCTGa  >  1:4041524/1‑71 (MQ=255)                                                              cccGCTGCTGGCGCGTAAAGCGAATAGTAAATAATTATTGTGGTCGGTGACGCTGACTCTGGCAGATCTGa  >  1:2376881/1‑71 (MQ=255)                                                              cccGCTGCTGGCGCGTAAAGCGAATAGTAAATAATTATTGTGGTCGGTGACGCTGACTCTGGCAGATCTGa  >  1:1758238/1‑71 (MQ=255)                                                              cccGCTGCTGGCGCGTAAAGCGAATAGTAAATAATTATTGTGGTCGGTGACGCTGACTCTGGCAGATCTGa  >  1:1723207/1‑71 (MQ=255)                                                              cccACTGCTGGCGCGTAAAGCGAATAGTAAATAATTATTGTGGTCGGTGACGCTGa                 >  1:1385293/1‑56 (MQ=255)                                                                 |                                                                    CCTTCCAGTTGGCGCTGAATGTTGCAACGCCGTTTGTGCTGGTAGGATTAGGACTGATATTCCCGCTGCTGGCGCGTAAAGCGAATAGTAAATAATTATTGTGGTCGGTGACGCTGACTCTGGCAGATCTGA  >  minE/605264‑605395

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: ‑