Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	1,215,064	G→A	60.4%	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	60.4%	54.6 / 110.7	91	L291L (CTG→CTA)	yebQ	predicted transporter
Reads supporting (aligned to +/- strand):  ref base G (20/16);  new base A (26/29);  total (46/45)
Fisher's exact test for biased strand distribution p-value = 5.22e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.02e-02

TGTTTTCACTTTCTATTTGCACATCTGTTTGCTCTTTCTGCGCACAAATGCTGGCAATGGTTTCCCTGCCCTTTTACCTGCAAACCGTGCTCGGGCGTAGTGAAGTCGAAACAGGTTTACTTCTGACACCGTGGCCGT  >  minE/1214997‑1215134                                                                    |                                                                       tGTTTTCACTTTCTATTTGCACATCTGTTTGCTCTTTCTGCGCACAAATGCTGGCAATGGTTTCCCTGccc                                                                     >  1:3649714/1‑71 (MQ=255)  gTTTTCACTTTCTATTTGCACATCTGTTTGCTCTTTCTGCGCACAAATGCTGGCAATGGTTTCCCTGCCCt                                                                    >  1:1893621/1‑71 (MQ=255)  gTTTTCACTTTCTATTTGCACATCTGTTTGCTCTTTCTGCGCACAAATGCTGGCAATGGTTTCCCTGCCCt                                                                    >  1:3816234/1‑71 (MQ=255)  gTTTTCACTTTCTATTTGCACATCTGTTTGCTCTTTCTGCGCACAAATGCTGGCAATGGTTTCCCTGCCCt                                                                    >  1:615298/1‑71 (MQ=255)  gTTTTCACTTTCTATTTGCACATCTGTTTGCTCTTTCTGCGCACAAATGCTGGCAATGGTTTCCCTGCCCt                                                                    >  1:3283951/1‑71 (MQ=255)     ttCACTTTCTATTTGCACATCTGTTTGCTCTTTCTGCGCACAAATGCTGGCAATGGTTTCCCTACCCttt                                                                  >  1:830104/1‑70 (MQ=255)     ttCACTTTCTATTTGCACATCTGTTTGCTCTTTCTGCGCACAAATGCTGGCAATGGTTTCCCTACCCttt                                                                  >  1:3051003/1‑70 (MQ=255)      tCACTTTCTATTTGCACATCTGTTTGCTCTTTCTGCGCACAAATGCTGGCAATGGTTTCCCTGCCCtttt                                                                 <  1:3851707/70‑1 (MQ=255)      tCACTTTCTATTTGCACATCTGTTTGCTCTTTCTGCGCACAAATGCTGGCAATGGTTTCCCTGCCCtttt                                                                 <  1:1182239/70‑1 (MQ=255)      tCACTTTCTATTTGCACATCTGTTTGCTCTTTCTGCGCACAAATGCTGGCAATGGTTTCCCTGCCCtttt                                                                 <  1:3861016/70‑1 (MQ=255)      tCACTTTCTATTTGCACATCTGTTTGCTCTTTCTGCGCACAAATGCTGGCAATGGTTTCCCTGCCCtttt                                                                 <  1:3894815/70‑1 (MQ=255)      tCACTTTCTATTTGCACATCTGTTTGCTCTTTCTGCGCACAAATGCTGGCAATGGTTTCCCTACCCtttt                                                                 <  1:2883712/70‑1 (MQ=255)      tCACTTTCTATTTGCACATCTGTTTGCTCTTTCTGCGCACAAATGCTGGCAATGGTTTCCCTACCCtttt                                                                 <  1:311457/70‑1 (MQ=255)        aCTTTCTATTTGCACATCTGTTTGCTCTTTCTGCGCACAAATGCTGGCAATGGTTTCCCTGCCCTTTTAcc                                                              <  1:2623716/71‑1 (MQ=255)           ttCTATTTGCACATCTGTTTGCTCTTTCTGCGCACAAATGCTGGCAATGGTTTCCCTACCCTTTTACCTg                                                            >  1:3113463/1‑70 (MQ=255)           ttCTATTTGCACATCTGTTTGCTCTTTCTGCGCACAAATGCTGGCAATGGTTTCCCTACCCTTTTACCTg                                                            >  1:1225811/1‑70 (MQ=255)                tttGCACATCTGTTTGCTCTTTCTGCGCACAAATGCTGGCAATGGTTTCCCTGCCCTTTTACCTGCAAAc                                                       <  1:1507063/70‑1 (MQ=255)                tttGCACATCTGTTTGCTCTTTCTGCGCACAAATGCTGGCAATGGTTTCCCTACCCtttt                                                                 <  1:639278/60‑1 (MQ=255)                tttGCACATCTGTTTGCTCTTGCTGCGCACAAATGCTGGCAATGGGTTCCCTGCCCTTTTACCTGCAAAcc                                                      <  1:3977063/71‑1 (MQ=255)                  tGCACATCTGTTTGCTCTTTCTGCGCACAAATGCTGGCAATGGTTTCCCTACCCTTTTAc                                                               >  1:646578/1‑60 (MQ=255)                  tGCACATCTGTTTGCTCTTTCTGCGCACAAATGCTGGCAATGGTTTCCCTACCCTTTTAc                                                               >  1:1834166/1‑60 (MQ=255)                     acaTCTGTTTGCTCTTTCTGCGCACAAATGCTGGCAATGGTTTCCCTACCCTTTTACCTGCaaa                                                        <  1:2186860/64‑1 (MQ=255)                     acaTCTGTTTGCTCTTTCTGCGCACAAATGCTGGCAATGGTTTCCCTACCCTTTTACCTGCAAACCGTg                                                   >  1:2544729/1‑69 (MQ=255)                     acaTCTGTTTGCTCTTTCTGCGCACAAATGCTGGCAATGGTTTCCCTACCCTTTTACCTGCAAACCGTGCt                                                 >  1:2263762/1‑71 (MQ=255)                     acaTCTGTTTGCTCTTTCTGCGCACAAATGCTGGCAATGGTTTCCCTACCCTTTTACCTGCAAACCGTGCt                                                 <  1:1003965/71‑1 (MQ=255)                      caTCTGTTTGCTCTTTCTGCGCACAAATGCTGGCAATGGTTTCCCTACCCTTTTACCTGCAAACCGTGCTc                                                <  1:2567639/71‑1 (MQ=255)                         cTGTTTGCTCTTTCTGCGCACAAATGCTGGCAATGGTTTCCCTGCCCTTTTACCTGCAAACCGTGCTCgg                                              >  1:2506431/1‑70 (MQ=255)                            tttGCTCTTTCTGCGCACAAATGCTGGCAATGGTTTCCCTACCCTTTTACCTGCAAACCGTGCTCggg                                             <  1:1097739/68‑1 (MQ=255)                            tttGCTCTTTCTGCGCACAAATGCTGGCAATGGTTTCCCTACCCTTTTACCTGCAAACCGTGCTCggg                                             <  1:3998997/68‑1 (MQ=255)                             ttGCTCTTTCTGCGCACAAATGCTGGCAATGGTTTCCCTGCCCTTTTAc                                                               >  1:3648377/1‑49 (MQ=255)                             ttGCTCTTTCTGCGCACAAATGCTGGCAATGGTTTCCCTGCCCTTTTAc                                                               >  1:4207180/1‑49 (MQ=255)                             ttGCTCTTTCTGCGCACAAATGCTGGCAATGGTTTCCCTGCCCTTTTAc                                                               >  1:4241833/1‑49 (MQ=255)                             ttGCTCTTTCTGCGCACAAATGCTGGCAATGGTTTCCCTGCCCTTTTAc                                                               >  1:1502709/1‑49 (MQ=255)                             ttGCTCTTTCTGCGCACAAATGCT‑GCAATGGTTTCCCTGCCCTTTTAc                                                               >  1:2408576/1‑48 (MQ=39)                              tGCTCTTTCTGCGCACAAATGCTGGCAATGGTTTCCCTGCCCTTTTACCTGCAAACCGTGCTCGGGCGTa                                         <  1:511291/70‑1 (MQ=255)                                   tttCTGCGCACAAATGCTGGCAATGGTTTCCCTGCCCTTTTACCTGCAAACCGTGCTCGGGCGTAGTGAAg                                   >  1:1372496/1‑71 (MQ=255)                                   tttCTGCGCACAAATGCTGGCAATGGTTTCCCTACCCTTTTACCTGCAAACCGTGCTCGGGCGTAgagatg                                   >  1:1534571/1‑69 (MQ=255)                                   tttCTGCGCACAAATGCTGGCAATGGTTTCCCTACCCTTTTACCTGCAAACCGTGCTCGGGCGTAgagatg                                   >  1:2571851/1‑69 (MQ=255)                                   ttgcTGCGCACAAATGCTGGCAATGGTTTCCCTGCCCTTTTACCTGCAAACCGTGCTCGGGCGTAGTGAAg                                   >  1:1821537/4‑71 (MQ=255)                                     tCTGCGCACAAATGCTGGCAATGGTTTCCCTGCCCTTTTAc                                                               >  1:1587224/1‑41 (MQ=255)                                     tCTGCGCACAAATGCTGGCAATGGTTTCCCTGCCCTTTTAc                                                               >  1:3828855/1‑41 (MQ=255)                                     tCTGCGCACAAATGCTGGCAATGGTTTCCCTGCCCTTTTACCTGCAAACCGTGCTCg                                               >  1:1084907/1‑57 (MQ=255)                                     tCTGCGCACAAATGCTGGCAATGGTTTCCCTACCCTTTTa                                                                >  1:3162181/1‑40 (MQ=255)                                     tCTGCGCACAAATGCTGGCAATGGTTTCCCTACCCTTTTACCTGCAAACCGTGCTCGGGCGTAGTGAAGTc                                 >  1:1347481/1‑71 (MQ=255)                                     tCTGCGCACAAATGCTGGCAATGGTTTCCCTACCCTTTTACCTGCAAACCGTGCTCGGGCGTAGTGAAGTc                                 >  1:2240883/1‑71 (MQ=255)                                      cagcgcACAAATGCTGGCAATGGTTTCCCTGCCCTTTTACCTGCAAACCGTGCTCGGGCGTAGTGAAGTc                                 <  1:1450829/68‑1 (MQ=255)                                      cTGCGCACAAATGCTGGCAATGGTTTCCCTACCCTTTTACCTGCAAACCGTGCTCGGGCGTAGTGAAGTc                                 >  1:3454447/1‑70 (MQ=255)                                         cgcACAAATGCTGGCAATGGTTTCCCTACCCTTTTACCTGCAAACCGTGCTCGGGCGTAGTGAAGTCGa                               <  1:3365095/69‑1 (MQ=255)                                           cacaAATGCTGGCAATGGTTTCCCTGCCCTTTTACCTGCAAACCGTGc                                                  >  1:589432/1‑48 (MQ=255)                                           cacaAATGCTGGCAATGGTTTCCCTACCCTTTTACCTGCAAACCGTGCTCGGGCGTAGTGAAGTCGAAACa                           <  1:2748674/71‑1 (MQ=255)                                           cacaAATGCTGGCAATGGTTTCCCTACCCTTTTACCTGCAAACCGTGCTCGGGCGTAGTGAAGTCGAAACa                           <  1:4250399/71‑1 (MQ=255)                                           cacaAATGCTGGCAATGGTTTCCCTACCCTTTTACCTGCAAACCGTGCTCGGGCGTAGTGAAGTCGAAACa                           <  1:4012315/71‑1 (MQ=255)                                           cacaAATGCTGGCAATGGTTTCCCTACCCTTTTACCTGCAAACCGTGCTCGGGCGTAGTGAAGTCGAAACa                           <  1:184047/71‑1 (MQ=255)                                           cacaAATGCTGGCAATGGTTTCCCTACCCTTTTACCTGCAAACCGTGCTCGGGCGTAGTGAAGTCGAAACa                           <  1:3491413/71‑1 (MQ=255)                                            acaAATGCTGGCAATGGTTTCCCTGCCCTTTTACCTGCAAACCGTGCTCGGGCGTAGTGAAGTCGAAAc                            <  1:348016/69‑1 (MQ=255)                                                aTGCTGGCAATGGTTTCCCTGCCCTTTTACCTGCAAACCGTGCTCGGGCGTAGTGAAGTCGAAACAGGttt                      <  1:1087480/71‑1 (MQ=255)                                                aTGCTGGCAATGGTTTCCCTGCCATTTTACCTGCAAACCGTGCTCGGGCGTAGTGAAGTCGAAACAGGttt                      <  1:2834750/71‑1 (MQ=255)                                                 tGCTGGCAATGGTTTCCCTACCCTTTTACCTGCaaa                                                        <  1:2931211/36‑1 (MQ=255)                                                  gCTGGCAATGGTTTCCCTACCCTTTTACCTGCAAACCGTGCTCGGGCGTAGTGAAGTCGAAACAGGTTTAc                    >  1:1750347/1‑71 (MQ=255)                                                  gCTGGCAATGGTTTCCCTACCCTTTTACCTGCAAACCGTGCTCGGGCGTAGTGAAGTCGAAACAGGTTTAc                    >  1:3826715/1‑71 (MQ=255)                                                    tGGCAATGGTTTCCCTACCCTTTTACCTGCAAACCGTGCt                                                 <  1:658300/40‑1 (MQ=255)                                                      gCAATGGTTTCCCTGCCCTTTTACCTGCAAACCGTGCTCGGGCGTAGTGAAGTCGAAACAGGTTTACTTCt                <  1:3182437/71‑1 (MQ=255)                                                      gCAATGGTTTCCCTGCCCTTTTACCTGCAAACCGTGCTCGGGCGTAGTGAAGTCGAAACAGGTTTACTTCt                <  1:59961/71‑1 (MQ=255)                                                      gCAATGGTTTCCCTGCCCTTTTACCTGCAAACCGTGCTCGGGCGTAGTGAAGTCGAAACAGGTTTACTTCt                <  1:2041030/71‑1 (MQ=255)                                                      gCAATGGTTTCCCTACCCTTTTACCTGCAAACCGTGCTCGGTCCTAGTGAAGTCGAAACAGGTTTACTTCt                <  1:2022785/71‑1 (MQ=255)                                                      gCAATGGTTTCCCTACCCTTTTACCTGCAAACCGTGCTCGGGCGTAGTGAAGTCGAAAc                            >  1:391228/1‑59 (MQ=255)                                                      gCAATGGTTTCCCTACCCTTTTACCTGCAAACCGTGCTCGGGCGTAGTGAAGTCGAAAc                            >  1:2798545/1‑59 (MQ=255)                                                      gCAATGGTTTCCCTACCCTTTTACCTGCAAACCGTGCTCGGGCGTAGTGAAGTCGAAAc                            >  1:2854910/1‑59 (MQ=255)                                                      gCAATGGTTTCCCTACCCTTTTACCTGCAAACCGTGCTCGGGCGTAGTGAAGTCGAAAc                            >  1:4241649/1‑59 (MQ=255)                                                      gCAATGGTTTCCCTACCCTTTTACCTGCAAACCGTGCTCGGGCGTAGTGAAGTCGAAACAGGTTTACTTCt                <  1:3227424/71‑1 (MQ=255)                                                      gCAATGGTTTCCCTACCCTTTTACCTGCAAACCGTGCTCGGGCGTAGTGAAGTCGAAACAGGTTTACTTCt                <  1:3882358/71‑1 (MQ=255)                                                      gCAATGGTTTCCCTACCCTTTTACCTGCAAACCGTGCTCGGGCGTAGTGAAGTCGAAACAGGTTTACTTCt                <  1:3393126/71‑1 (MQ=255)                                                      gCAATGGTTTCCCTACCCTTTTACCTGCAAACCGTGCTCGGGCGTAGTGAAGTCGAAACAGGTTTACTTCt                <  1:1356416/71‑1 (MQ=255)                                                      gCAATGGTTTCCCTACCCTTTTACCTGCAAACCGTGCTCGGGCGTAGTGAAGTCGAAACAGGTTTACTTCt                <  1:3318174/71‑1 (MQ=255)                                                      gCAATGGTTTCCCTACCCTTTTACCTGCAAACCGTGCTCGGGCGTAGTGAAGTCGAAACAGGTTTACTTCt                <  1:1723723/71‑1 (MQ=255)                                                      gCAATGGTTTCCCTACCCTTTTACCTGCAAACCGTGCTCGGGCGTAGTGAAGTCGAAACAGGTTTACTTCt                <  1:2423616/71‑1 (MQ=255)                                                      gCAATGGTTTCCCTACCCTTTTACCTGCAAACCGTGCTCGGGCGTAGTGAAGTCGAAACAGGTTTACTTCt                <  1:2884786/71‑1 (MQ=255)                                                      gCAATGGTTTCCCTACCCTTTTACCTGCAAACCGTGCTCGGGCGTAGTGAAGTCGAAACAGGTTTACTTCt                <  1:476360/71‑1 (MQ=255)                                                      gCAATGGTTTCCCTACCCTTTTACCTGCAAACCGTGCTCGGGCGTAGTGAAGTCGAAACAGGTTTACTTCt                <  1:2549676/71‑1 (MQ=255)                                                            gTTTCCCTGCCCTTTTACCTGCAAACCGTGCTCGGGCGTAGTGAAGTCGAAACAGGTTTACTTCTGACAcc          >  1:1987984/1‑71 (MQ=255)                                                            gTTTCCCTACCCTTTTACCTGCAAACCGTGCTCGGGCGTAGTGAAGTCGaaa                             >  1:1612145/1‑52 (MQ=255)                                                            gTTTCCCTACCCTTTTACCTGCAAACCGTGCTCGGGCGTAGTGAAGTCGaaa                             >  1:2895278/1‑52 (MQ=255)                                                            gTTTCCCTACCCTTTTACCTGCAAACCGTGCTCGGGCGTAGTGAAGTCGaaa                             >  1:3520742/1‑52 (MQ=255)                                                            gTTTCCCTACCCTTTTACCTGCAAACCGTGCTCGGGCGTAGTGAAGTCGaaa                             >  1:1693638/1‑52 (MQ=255)                                                            gTTTCCCTACCCTTTTACCTGCAAACCGTGCTCGGGCGTAGTGAAGTCGaaa                             >  1:2592087/1‑52 (MQ=255)                                                             tttCCCTACCCTTTTACCTGCAAACCGTGCTCGGGCGTAGTGAAGTCGAAACAGGTTTACTTCTGACACCg         <  1:1660981/71‑1 (MQ=255)                                                             tttCCCTACCCTTTTACCTGCAAACCGTGCTCGGGCGTAGTGAAGTCGAAACAGGTTTACTTCTGACACCg         <  1:582008/71‑1 (MQ=255)                                                              ttCCCTGCCCTTTTACCTGCAAACCGTGCTCGGGCGTAGTg                                      <  1:1218562/41‑1 (MQ=255)                                                               tCCCTGCCCTTTTACCTGCAAACCGTGCTCGGGCGTAGTGAAGTCGAAACAGGTTTACTTCt                >  1:3416970/1‑62 (MQ=255)                                                               tCCCTGCCCTTTTACCTCCAAACCGTGCTCGGGCGTAGTGAAGTCGAAACAGGTTTACTTc                 >  1:1940773/1‑61 (MQ=255)                                                               tCCCTACCCTTTTACCTGCAAACCGTGCTCGGGCGTAGTGAAGTCGAAACAGGTTTACTTCTGACACCGTg       >  1:3019286/1‑71 (MQ=255)                                                                    gCCCTTTTACCTGCAAACCGTGCTCGGGCGTAGTGAAGTCGAAACAGGTTTACTTCTGACACCGTGGCCGt  >  1:3911608/1‑71 (MQ=255)                                                                    gCCCTTTTACCTGCAAACCGTGCTCGGGCGTAGTGAAGTCGAAACAGGTTTACTTCTGACACCGTGGCCGt  >  1:1745044/1‑71 (MQ=255)                                                                    |                                                                       TGTTTTCACTTTCTATTTGCACATCTGTTTGCTCTTTCTGCGCACAAATGCTGGCAATGGTTTCCCTGCCCTTTTACCTGCAAACCGTGCTCGGGCGTAGTGAAGTCGAAACAGGTTTACTTCTGACACCGTGGCCGT  >  minE/1214997‑1215134

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