Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	1,215,064	G→A	30.2%	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	30.2%	148.5 / 76.9	96	L291L (CTG→CTA)	yebQ	predicted transporter
Reads supporting (aligned to +/- strand):  ref base G (22/45);  new base A (14/15);  total (36/60)
Fisher's exact test for biased strand distribution p-value = 1.73e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TTTTCACTTTCTATTTGCACATCTGTTTGCTCTTTCTGCGCACAAATGCTGGCAATGGTTTCCCTGCCCTTTTACCTGCAAACCGTGCTCGGGCGTAGTGAAGTCGAAACAGGTTTACTTCTGACACCGTG  >  minE/1214999‑1215129                                                                  |                                                                  ttttCACTTTCTATTTGCACATCTGTTTGCTCTTTCTGCGCACAAATGCTGGCAATGGTTTCCCTAccc                                                                <  1:1091727/69‑1 (MQ=255) ttttCACTTTCTATTTGCACATCTGTTTGCTCTTTCTGCGCACAAATGCTGGCAATGGTTTCCCTAccc                                                                <  1:2533444/69‑1 (MQ=255)   ttCACTTTCTATTTGCACATCTGTTTGCTCTTTCTGCGCACAAATGCTGGCAATGGTTTCCCTGCCCttt                                                             >  1:2624449/1‑70 (MQ=255)    tCACTTTCTATTTGCACATCTGTTTGCTCTTTCTGCGCACAAATGCTGGCAATGGTTTCCCTGCCCtttt                                                            <  1:409207/70‑1 (MQ=255)           cTATTTGCACATCTGTTTGCTCTTTCTGCGCACAAATGCTGGCAATGGTTTCCCTGCCCTTTTACCTGCa                                                     <  1:2284615/70‑1 (MQ=255)           cTATTTGCACATCTGTTTGCTCTTTCTGCGCACAAATGCTGGCAATGGTTTCCCTGCCCTTTTACCTGCa                                                     <  1:1501572/70‑1 (MQ=255)           cTATTTGCACATCTGTTTGCTCTTTCTGCGCACAAATGCTGGCAATGGTTTCCCTACCCTTTTAc                                                          >  1:3026458/1‑65 (MQ=255)                  cacaTCTGTTTGCTCTTTCTGCGCACAAATGCTGGCAATGGTTTCCCTGCCCTTTTACCTGCAAACCGTGc                                             <  1:275222/71‑1 (MQ=255)                  cacaTATGTTTGCTCTTTCTGCGCACAAATGCTGGCAATGGTTTCCCTGCCCTTTTACCTGCAAACCGTGc                                             <  1:3647206/71‑1 (MQ=255)                  cacaTATGTTTGCTCTTTCTGCGCACAAATGCTGGCAATGGTTTCCCTGCCCTTTTACCTGCAAACCGTGc                                             <  1:3642387/71‑1 (MQ=255)                   acaTCTGTTTGCTCTTTCTGCGCACAAATGCTGGCAATGGTTTCCCTGCCCTTTTACCTGCAAACCGTGCt                                            <  1:1092628/71‑1 (MQ=255)                   acaTCTGTTTGCTCTTTCTGCGCACAAATGCTGGCAATGGTTTCCCTGCCCTTTTACCTGCAAACCGTGCt                                            <  1:3555706/71‑1 (MQ=255)                   acaTCTGTTTGCTCTTTCTGCGCACAAATGCTGGCAATGGTTTCCCTGCCCTTTTACCTGCAAACCGTGCt                                            <  1:538137/71‑1 (MQ=255)                      tCTGTTTGCTCTTTCTGCGCACAAATGCTGGCAATGGTTTCCCTGCCCtttt                                                            <  1:1389357/52‑1 (MQ=255)                        tGTTTGCTCTTTCTGCGCACAAATGCTGGCAATGGTTTCCCTACCCTTTTACCTGCAAACCGTGCTCGGGc                                       >  1:1453305/1‑71 (MQ=255)                          tttGCTCTTTCTGCGCACAAATGCTGGCAATGGTTTCCCTACCCTTTTAc                                                          <  1:1806792/50‑1 (MQ=255)                          tttGCTCTTTCTGCGCACAAATGCTGGCAATGGTTTCCCTACCCTTTTAc                                                          <  1:3386117/50‑1 (MQ=255)                          tttGCTCTTTCTGCGCACAAATGCTGGCAATGGTTTCCCTACCCTTTTAc                                                          <  1:3988106/50‑1 (MQ=255)                          tttGCTCTTTCTGCGCACAAATGCTGGCAATGGTTTCCCTACCCTTTTAc                                                          <  1:1076953/50‑1 (MQ=255)                            tGCTCTTTCTGCGCACAAATGCTGGCAATGGTTTCCCTGCCCTTTTACCTGCAAACCGTGCt                                            >  1:3896611/1‑62 (MQ=255)                            tGCTCTTTCTGCGCACAAATGCTGGCAATGGTTTCCCTGCCCTTTTACCTGCAAACCGTGCt                                            >  1:239340/1‑62 (MQ=255)                            tGCTCTTTCTGCGCACAAATGCTGGCAATGGTTTCCCTGCCCTTTTACCTGCAAACCGTGCt                                            >  1:2955398/1‑62 (MQ=255)                            tGCTCTTTCTGCGCACAAATGCTGGCAATGGTTTCCCTGCCCTTTTACCTGCAAACCGTGCTCGGGCGTAg                                   >  1:2132383/1‑71 (MQ=255)                            tGCTCTTTCTGCGCACAAATGCTGGCAATGGTTTCCCTGCCCTTTTACCTGCAAACCGTGCTCGGGCGTAg                                   >  1:2829038/1‑71 (MQ=255)                            tGCTCTTTCTGCGCACAAATGCTGGCAATGGTTTCCCTGCCCTTTTACCTGCAAACCGTGCTCGGGCGTAg                                   >  1:1131329/1‑71 (MQ=255)                            tGCTCTTTCTGCGCACAAATGCTGGCAATGGTTTCCCTACCCttt                                                             >  1:1101171/1‑45 (MQ=255)                            tGCTCTTTCTGCGCACAAATGCTGGCAATGGTTTCCCTACCCTTTTAcc                                                         >  1:891498/1‑49 (MQ=255)                            tGCTCTTTCTGCGCACAAATGCTGGCAATGGTTTCCCTACCCTTTTACCt                                                        >  1:2740308/1‑50 (MQ=255)                            tGCTCTTTCTGCGCACAAATGCTGGCAATGGTTTCCCTACCCTTTTACCt                                                        >  1:3583933/1‑50 (MQ=255)                            tGCTCTTTCTGCGCACAAATGCTGGCAATGGTTTCCCTACCCTTTTACCt                                                        >  1:3088368/1‑50 (MQ=255)                            tGCTCTTTCTGCGCACAAATGCTGGCAATGGTTTCCCTACCCTTTTACCt                                                        >  1:1191280/1‑50 (MQ=255)                            tGCTCTTTCTGCGCACAAATGCTGGCAATGGTTTCCCTACCCTTTTACCt                                                        >  1:1377033/1‑50 (MQ=255)                            tGCTCTTTCTGCGCACAAATGCTGGCAATGGTTTCCCTACCCTTTTACCt                                                        >  1:30422/1‑50 (MQ=255)                              ctctTTCTGCGCACAAATGCTGGCAATGGTTTCCCTGCCCTTTTACCTGCAAACCGTGCTCGGGCGTAGt                                  <  1:910063/70‑1 (MQ=255)                              ctctTTCTGCGCACAAATGCTGGCAATGGTTTCCCTGCCCTTTTACCTGCAAACCGTGCTCGGGCGTAGTg                                 <  1:610628/71‑1 (MQ=255)                                   tCTGCGCACAAATGCTGGCAATGGTTTCCCTGCCCTTTTACCTGCAAACCGTGCTCGGGCGTAGTGAAGTc                            >  1:1916563/1‑71 (MQ=255)                                   tCTGCGCACAAATGCTGGCAATGGTTTCCCTGCCCTTTTACCTGCAAACCGTGCTCGGGCGTAGTGAAGTc                            >  1:4007280/1‑71 (MQ=255)                                   tCTGCGCACAAATGCTGGCAATGGTTTCCCTACCCTTTTACCTGCAAACCGTGCTCGGGCGTAGTGAAGTc                            >  1:2854282/1‑71 (MQ=255)                                   tCTGCGCACAAATGCTGGCAATGGTTTCCCTACCCTTTTACCTGCAAACCGTGCTCGGGCGTAGTGAAGTc                            >  1:155500/1‑71 (MQ=255)                                    cTGCGCACAAATGCTGGCAATGGTTTCCCTACCCTTTTAc                                                          <  1:3669654/40‑1 (MQ=255)                                     tGCGCACAAATGCTGGCAATGGTTTCCCTGCCCTTTTACCTGCAAAcc                                                 <  1:246387/48‑1 (MQ=255)                                       cgcACAAATGCTGGCAATGGTTTCCCTGCCCTTTTACCTGCAAACCGTGCTCGGGCGTAGTGAAGTCGa                          <  1:1038499/69‑1 (MQ=255)                                       cgcACAAATGCTGGCAATGGTTTCCCTGCCCTTTTACCTGCAAACCGTGCTCGGGCGTAGTGAAGTCGa                          <  1:3264766/69‑1 (MQ=255)                                        gcACAAATGCTGGCAATGGTTTCCCTGCCCTTTTACCTGCAAACCGTGCTCGGGCGTAGTGAAGTCGAAAc                       <  1:2716910/71‑1 (MQ=255)                                        gcACAAATGCTGGCAATGGTTTCCCTGCCCTTTTACCTGCAAACCGTGCTCGGGCGTAGTGAAGTCGAAAc                       <  1:2522904/71‑1 (MQ=255)                                        gcACAAATGCTGGCAATGGTTTCCCTGCCCTTTTACCTGCAAACCGTGCTCGGGCGTAGTGAAGTCGAAAc                       <  1:3302000/71‑1 (MQ=255)                                        gcACAAATGCTGGCAATGGTTTCCCTGCCCTTTTACCTGCAAACCGTGCTCGGGCGTAGTGAAGTCGAAAc                       <  1:2172319/71‑1 (MQ=255)                                        gcACAAATGCTGGCAATGGTTTCCCTACCCTTTTACCTGCAAACCGTGCTCGGGCGTAGTGAAGTCGAAAc                       <  1:3153007/71‑1 (MQ=255)                                        gcACAAATGCTGGCAATGGTTTCCCTACCCTTTTACCTGCAAACCGTGCTCGGGCGTAGTGAAGTCGAAAc                       <  1:3650793/71‑1 (MQ=255)                                        gcACAAATGCTGGCAATGGTTTCCCTACCCTTTTACCTGCAAACCGTGCTCGGGCGTAGTGAAGTCGAAAc                       <  1:2436314/71‑1 (MQ=255)                                             aaTGCTGGCAATGGTTTCCCTGCCCTTTTACCTGCAAACCGTGCTCGGGCGTAg                                   >  1:2265608/1‑54 (MQ=255)                                             aaTGCTGGCAATGGTTTCCCTGCCCTTTTACCTGCAAACCGTGCTCGGGCGTAGTGAAGTCGAAACAGGtt                  <  1:2892392/71‑1 (MQ=255)                                             aaTGCTGGCAATGGTTTCCCTGCCCTTTTACCTGCAAACCGTGCTCGGGCGTAGTGAAGTCGAAACAGGtt                  <  1:3670776/71‑1 (MQ=255)                                             aaTGCTGGCAATGGTTTCCCTGCCCTTTTACCTGCAAACCGTGCTCGGGCGTAGTGAAGTCGAAACAGGtt                  <  1:3329302/71‑1 (MQ=255)                                             aaTGCTGGCAATGGTTTCCCTGCCCTTTTACCTGCAAACCGTGCTCGGGCGTAGTGAAGTCGAAACAGGtt                  <  1:2287360/71‑1 (MQ=255)                                             aaTGCTGGCAATGGTTTCCCTGCCCTTTTACCTGCAAACCGTGCTCGGGCGTAGTGAAGTCGAAACAGGtt                  <  1:373738/71‑1 (MQ=255)                                             aaTGCTGGCAATGGTTTCCCTGCCCTTTTACCTGCAAACCGTGCTCGGGCGTAGTGAAGTCGAAACAGGtt                  <  1:4083710/71‑1 (MQ=255)                                             aaTGCTGGCAATGGTTTCCCTGCCCTTTTACCTGCAAACCGTGCTCGGGCGTAGTGAAGTCGAAACAGGtt                  <  1:564948/71‑1 (MQ=255)                                             aaTGCTGGCAATGGTTTCCCTGCCCTTTTACCTGCAAACCGTGCTCGGGCGTAGTGAAGTCGAAACAGGtt                  <  1:3032178/71‑1 (MQ=255)                                             aaTGCTGGCAATGGTTTCCCTGCCCTTTTACCTGCAAACCGTGCTCGGGCGTAGTGAAGTCGAAACAGGtt                  <  1:4030937/71‑1 (MQ=255)                                                 cTGGCAATGGTTTCCCTGCCCTTTTACCTGCAAACCGTGCTCGGGCg                                      <  1:256995/47‑1 (MQ=255)                                                 cTGGCAATGGTTTCCCTGCCCTTTTACCTGCAAACCGTGCTCGGGCg                                      <  1:1125482/47‑1 (MQ=255)                                                 cTGGCAATGGTTTCCCTGCCCTTTTACCTGCAAACCGTGCTCGGGCg                                      <  1:185013/47‑1 (MQ=255)                                                  tGGCAATGGTTTCCCTACCCTTTTACCTGCAAACCGTGCTCGGGCGTAGTGAAGTCGa                          <  1:2357632/58‑1 (MQ=255)                                                  tGGCAATGGTTTCCCTACCCTTTTACCTGCAAACCGTGCTCGGGCGTAGTGAAGTCGa                          <  1:3294805/58‑1 (MQ=255)                                                  tGGCAATGGTTTCCCTACCCTTTTACCTGCAAACCGTGCTCGGGCGTAGTGAAGTCGa                          <  1:2704612/58‑1 (MQ=255)                                                    gCAATGGTTTCCCTACCCTTTTACCTGCTAACCGTGCTCGGGCGTAGTGAAGTCGaaa                        >  1:291463/1‑58 (MQ=255)                                                    gCAATGGTTTCCCTACCCTTTTACCTGCAAACCGTGCTCGGGCGTAGTGAAGTCGaaa                        >  1:2572240/1‑58 (MQ=255)                                                       aTGGTTTCCCTGCCCTTTTACCTGCAAACCGTGCTCGGGCGTAGTGAAGTCGAAACAGGTTTACTTCTGa         >  1:602267/1‑70 (MQ=255)                                                       aTGGTTTCCCTGCCCTTTTACCTGCAAACCGTGCTCGGGCGTAGTGAAGTCGAAACAGGTTTACTTCTGa         >  1:39264/1‑70 (MQ=255)                                                       aTGGTTTCCCTACCCTTTTACCTGCAAACCGTGCTCGGGCGTAGTGAAGTCGAAACAGGTTTACTTCTGac        <  1:2009922/71‑1 (MQ=255)                                                        tGGTTTCCCTGCCCTTTTACCTGCAAACCGTGCTCGGGCGt                                     <  1:288992/41‑1 (MQ=255)                                                        tGGTTTCCCTGCCCTTTTACCTGCAAACCGTGCTCGGGCGt                                     <  1:1923147/41‑1 (MQ=255)                                                        tGGTTTCCCTGCCCTTTTACCTGCAAACCGTGCTCGGGCGt                                     <  1:2612475/41‑1 (MQ=255)                                                        tGGTTTCCCTGCCCTTTTACCTGCAAACCGTGCTCGGGCGTAGTGAAGTCGa                          <  1:3581988/52‑1 (MQ=255)                                                        tGGTTTCCCTGCCCTTTTACCTGCAAACCGTGCTCGGGCGTAGTGAAGTCGa                          <  1:3212149/52‑1 (MQ=255)                                                        tGGTTTCCCTGCCCTTTTACCTGCAAACCGTGCTCGGGCGTAGTGAAGTCGa                          <  1:923767/52‑1 (MQ=255)                                                        tGGTTTCCCTGCCCTTTTACCTGCAAACCGTGCTCGGGCGTAGTGAAGTCGa                          <  1:1603261/52‑1 (MQ=255)                                                        tGGTTTCCCTGCCCTTTTACCTGCAAACCGTGCTCGGGCGTAGTGAAGTCGa                          <  1:4020006/52‑1 (MQ=255)                                                        tGGTTTCCCTGCCCTTTTACCTGCAAACCGTGCTCGGGCGTAGTGAAGTCGa                          <  1:2020666/52‑1 (MQ=255)                                                          gTTTCCCTGCCCTTTTACCTGCAAACCGTGCTCGGGCGTAGTg                                 <  1:3783682/43‑1 (MQ=255)                                                          gTTTCCCTGCCCTTTTACCTGCAAACCGTGCTCGGGCGTAGTg                                 <  1:1812770/43‑1 (MQ=255)                                                          gTTTCCCTGCCCTTTTACCTGCAAACCGTGCTCGGGCGTAGTg                                 <  1:2832985/43‑1 (MQ=255)                                                          gTTTCCCTGCCCTTTTACCTGCAAACCGTGCTCGGGCGTAGTg                                 <  1:2795553/43‑1 (MQ=255)                                                          gTTTCCCTGCCCTTTTACCTGCAAACCGTGCTCGGGCGTAGTGaa                               >  1:3081136/1‑45 (MQ=255)                                                          gTTTCCCTACCCTTTTACCTGCAAACCGTGCTCGGGCGt                                     <  1:3915122/39‑1 (MQ=255)                                                            ttCCCTGCCCTTTTACCTGCAAACCGTGCTCGGGCGTAGTGAAGTCGAAACAGGTTTACTTCTGACACCGt   >  1:3952911/1‑71 (MQ=255)                                                             tCCCTGCCCTTTTACCTGCAAACCGTGCTCGGGCGTAGTGAAGTCGAAACAGGTTTACTTCTGACACCGTg  >  1:3353574/1‑71 (MQ=255)                                                             tCCCTGCCCTTTTACCTGCAAACCGTGCTCGGGCGTAGTGAAGTCGAAACAGGTTTACTTCTGACACCGTg  >  1:2408781/1‑71 (MQ=255)                                                              cccTGCCCTTTTACCTGCAAACCGTGCTCGGGCGTAGTGAAGTCGAAACAGGTTTACTTCTGACACCGTg  <  1:2944683/70‑1 (MQ=255)                                                               ccTGCCCTTTTACCTGCAAACCGTGCTCGGGCGTAGTGaa                               >  1:4009245/1‑40 (MQ=255)                                                               ccTGCCCTTTTACCTGCAAACCGTGCTCGGGCGTAGTGaa                               >  1:3942740/1‑40 (MQ=255)                                                               ccTGCCCTTTTACCTGCAAACCGTGCTCGGGCGTAGTGaa                               >  1:1728899/1‑40 (MQ=255)                                                               ccTGCCCTTTTACCTGCAAACCGTGCTCGGGCGTAGTGaa                               >  1:1326317/1‑40 (MQ=255)                                                               ccTGCCCTTTTACCTGCAAACCGTGCTCGGGCGTAGTGaa                               >  1:1259151/1‑40 (MQ=255)                                                               ccTGCCCTTTTACCTGCAAACCGTGCTCGGGCGTAGTGaa                               >  1:3799082/1‑40 (MQ=255)                                                                  |                                                                  TTTTCACTTTCTATTTGCACATCTGTTTGCTCTTTCTGCGCACAAATGCTGGCAATGGTTTCCCTGCCCTTTTACCTGCAAACCGTGCTCGGGCGTAGTGAAGTCGAAACAGGTTTACTTCTGACACCGTG  >  minE/1214999‑1215129

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