Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	2,603,351	C→T	7.5%	V440V (GTG→GTA)	yhiP ←	predicted transporter

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	2,603,351	0	C	T	7.5%	284.7 / 9.4	93	V440V (GTG→GTA)	yhiP	predicted transporter
Reads supporting (aligned to +/- strand):  ref base C (33/53);  new base T (5/2);  total (38/55)
Fisher's exact test for biased strand distribution p-value = 1.17e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.98e-01

CCAGACCAATCTTACCAAACACGTTGGTATAGACGGGCAACGTCTCAAGCGGATCGGTAATGTTGTCCGGCACCGCGGTAAATGTTGCCACATAGCCGCCCAGCAAGAACGCGGCAGCCTGCGTCAGGAACCACATCCC  >  minE/2603281‑2603419                                                                       |                                                                     ccagaccaATCTTACCAAACACGTTGGTATAGACGGGCAACGTCTCAAGCGGATCGGTAATGTTGTCCGGc                                                                      <  1:3702244/71‑1 (MQ=255) caagaccaATCTTACCAAACACGTTGGTATAGACGGGCAACGTCTCAAGCGGATCGGTAATGTTGTCCGGc                                                                      <  1:1186986/69‑1 (MQ=255)  cagaccaATCTTACCAAACACGTTGGTATAGACGGGCAACGTCTCAAGCGGATCGGTAATGTTGTCCGGCa                                                                     >  1:4038084/1‑71 (MQ=255)     accaATCTTACCAAACACGTTGGTATAGACGGGCAACGTCTCAAGCGGATCGGTAATGTTGTCCGGCACcg                                                                  <  1:1263780/71‑1 (MQ=255)      ccaATCTTACCAAACACGTTGGTATAGACGGGCAACGTCTCAAGCGGATCGGTAATGTTGTCCGGCACcgc                                                                 >  1:4469188/1‑71 (MQ=255)       caATCTTACCAAACACGTTGGTATAGACGGGCAACGTCTCAAGCGGATCGGTAATGTTGTCCGGCACcgcg                                                                >  1:4551243/1‑71 (MQ=255)         aTCTTACCAAACACGTTGGTATAGACGGGCAACGTCTCAAGCGGATCGGTAATGTTGTCCGGCACCGCgg                                                               >  1:970098/1‑70 (MQ=255)           cTTACCAAACACGTTGGTATAGACGGGCAACGTCTCAAGCGGATCGGTAATGTTGTCCGGCACCGCGGTa                                                             >  1:818531/1‑70 (MQ=255)              aCCAAACACGTTGGTATAGACGGGCAACGTCTCAAGCGGATCGGTAATGTTGTCCGGCACCGCGGTaaa                                                           <  1:1043519/69‑1 (MQ=255)              aCCAAACACGTTGGTATAGACGGGCAACGTCTCAAGCGGATCGGTAATGTTGTCCGGCACCGCGGTaaa                                                           <  1:2083810/69‑1 (MQ=255)              aCCAAACACGTTGGTATAGACGGGCAACGTCTCAAGCGGATCGGTAATGTTGTCCGGCACCGCGGTaaa                                                           <  1:1490573/69‑1 (MQ=255)              aCCAAACACGTTGGTATAGACGGGCAACGTCTCAAGCGGATCGGTAATGTTGTCCGGCACCGCGGTaa                                                            <  1:1843319/68‑1 (MQ=255)              aCCAAACACGTTGGTATAGACGGGCAACGTCTCAAGCGGATCGGTAATGTTGTCCGGCACCGCGGTaa                                                            <  1:1648742/68‑1 (MQ=255)              aCCAAACACGTTGGTATAGACGGGCAACGTCTCAAGCGGATCGGTAATGTTGTCCGGCACCGCGGTa                                                             <  1:4652112/67‑1 (MQ=255)              aCCAAACACGTTGGTATAGACGGGCAACGTCTCAAGCGGATCGGTAATGTTGTCCGGCACCGCGGTa                                                             <  1:1999554/67‑1 (MQ=255)              aCCAAACACGTTGGTATAGACGGGCAACGTCTCAAGCGGATCGGTAATGTTGTCCGGCACCGCGGTa                                                             <  1:1861925/67‑1 (MQ=255)              aCCAAACACGTTGGTATAGACGGGCAACGTCTCAAGCGGATCGGTAATGTTGTCCGGCACCGCGGTAAAt                                                          <  1:184411/70‑1 (MQ=255)              aCCAAACACGTTGGTATAGACGGGCAACGTCTCAAGCGGATCGGTAATGTTGTCCGGCACCGCGGTAAAt                                                          <  1:1982399/70‑1 (MQ=255)               ccAAACACGTTGGTATAGACGGGCAACGTCTCAAGCGGATCGGTAATGTTGTCCGGCACCGCGGTAAATGt                                                        >  1:4010549/1‑71 (MQ=255)                cAAACACGTTGGTATAGACGGGCAACGTCTCAAGCGGATCGGTAATGTTGTCCGGc                                                                      <  1:15578/56‑1 (MQ=255)                cAAACACGTTGGTATAGACGGGCAACGTCTCAAGCGGATCGGTAATGTTGTCCGGc                                                                      <  1:2068025/56‑1 (MQ=255)                cAAACACGTTGGTATAGACGGGCAACGTCTCAAGCGGATCGGTAATGTTGTCCGGCACCGCGGTa                                                             <  1:3939507/65‑1 (MQ=255)                  aaCACGTTGGTATAGACGGGCAACGTCTCAAGCGGATCGGTAATGTTGTCCGGCACCGCGGTAAATGTTGc                                                     >  1:3721277/1‑71 (MQ=255)                  aaCACGTTGGTATAGACGGGCAACGTCTCAAGCGGATCGGTAATGTTGTCCGGCACCGCGGTAAATGTTGc                                                     <  1:3940948/71‑1 (MQ=255)                  aaCACGTTGGTATAGACGGGCAACGTCTCAAGCGGATCGGTAATGTTGTCCGGCACCGCGGTAAATGTTGc                                                     >  1:2288308/1‑71 (MQ=255)                    cacGTTGGTATAGACGGGCAACGTCTCAAGCGGATCGGTAATGTTGTCCGGCACCGCGGTAAATGTTGCca                                                   <  1:2230495/71‑1 (MQ=255)                     acGTTGGTATAGACGGGCAACGTCTCAAGTGGATCGGTAATGTTGTCTGGTACTGCGGTAAATGTTGCcac                                                  <  1:93216/71‑1 (MQ=255)                     acGTTGGTATAGACGGGCAACGTCTCAAGTGGAGCGGTAATGTTGTCTGGTACTGCGGTAAATGTTGCcac                                                  <  1:3073503/71‑1 (MQ=255)                       gTTGGTATAGACGGGCAACGTCTCAAGCGGATCGGTAATGTTGTCCGGCACCGCGGTAAATGTTGCCACAt                                                <  1:822036/71‑1 (MQ=255)                       gTTGGTATAGACGGGCAACGTCTCAAGCGGATCGGTAATGTTGTCCGGCACCGCGGTAAATGTTGCCACAt                                                <  1:2887222/71‑1 (MQ=255)                       gTTGGTATAGACGGGCAACGTCTCAAGCGGATCGGTAATGTTGTCCGGCACCGCGGTAAATGTTGCCACAt                                                <  1:4420413/71‑1 (MQ=255)                           gTATAGACGGGCAACGTCTCAAGCGGATCGGTAATGTTGTCCGGCACCGCGGTa                                                             <  1:2229406/54‑1 (MQ=255)                           gTATAGACGGGCAACGTCTCAAGCGGATCGGTAATGTTGTCCGGCACCGCGGTAAATGTTGCca                                                   <  1:3162340/64‑1 (MQ=255)                           gTATAGACGGGCAACGTCTCAAGCGGATCGGTAATGTTGTCCGGCACCGCGGTAAATGTTGCCACATAgc                                             <  1:2413015/70‑1 (MQ=255)                            tataGACGGGCAACGTCTCAAGCGGATCGGTAATGTTGTCCGGCACTGCGGTAAATGTTGCcaca                                                 >  1:4573113/1‑65 (MQ=255)                            tataGACGGGCAACGTCTCAAGCGGATCGGTAATGTTGTCCGGCACCGCGGTAAATGTTGCcaca                                                 >  1:303159/1‑65 (MQ=255)                              taGACGGGCAACGTCTCAAGCGGATCGGTAATGTTGTCCGGCACCGCGGTaa                                                            <  1:1111688/52‑1 (MQ=255)                                gACGGGCAACGTCTCAAGCGGATCGGTAATGTTGTCCGGCACCGCGGTAAATGTTGCCACATAgccgcc                                         >  1:3978275/1‑69 (MQ=255)                                 aCGGGCAACGTCTCAAGCGGATCGGTAATGTTGTCCGGCACCGCGGTAAAt                                                          <  1:3141178/51‑1 (MQ=255)                                 aCGGGCAACGTCTCAAGCGGATCGGTAATGTTGTCCGGCACCGCGGTAAATGTTGCCACATAGCCGCCCAg                                      >  1:2044725/1‑71 (MQ=255)                                    ggCAACGTCTCAAGCGGATCGGTAATGTTGTCCGGCACCGCGGTAAATGTTGGCACATAGCCGCCCAGCaa                                   <  1:855733/71‑1 (MQ=255)                                     gCAACGTCTCAAGCGGATCGGTAATGTTGTCCGGCACCGCGGTAAATGTTGCCACATAgccgcc                                         <  1:3007805/64‑1 (MQ=255)                                     gCAACGTCTCAAGCGGATCGGTAATGTTGTCCGGCACCGCGGTAAATGTTGCCACATAGCCGCCCAGCAAg                                  <  1:3267562/71‑1 (MQ=255)                                     gCAACGTCTCAAGCGGATCGGTAATGTTGTCCGGCACCGCGGTAAATGTTGCCACATAGCCGCCCAGCAAg                                  <  1:1694433/71‑1 (MQ=255)                                      cAACGTCTCAAGCGGATCGGTAATGTTGTCCGGCACCGCGGTAAATGTTGCCACATAGCCGCCCAGCAAg                                  >  1:4345756/1‑70 (MQ=255)                                      cAACGTCTCAAGCGGATCGGTAATGTTGTCCGGCACCGCGGTAAATGTTGCCACATAGCCGCCCAGCAAGa                                 <  1:504914/71‑1 (MQ=255)                                       aaCGTCTCAAGCGGATCGGTAATGTTGTCCGGCACCGCGGTAAATGTTGCCACATa                                               >  1:607096/1‑56 (MQ=255)                                        aCGTCTCAAGTGGATCGGTAATGTTGTCTGGTACTGCGGTAAATGTTGCCACATAGCCGCCCAGCaa                                   >  1:2974053/1‑67 (MQ=255)                                        aCGTCTCAAGTGGATCGGTAATGTTGTCTGGTACTGCGGTAAATGTTGCCACATAGCCGCCCAGCaa                                   >  1:3062828/1‑67 (MQ=255)                                        aCGTCTCAAGTGGATCGGTAATGTTGTCTGGTACTGCGGTAAATGTTGCCACATAGCCGCCCAGCaa                                   >  1:1377851/1‑67 (MQ=255)                                        aCGTCTCAAGTGGATCGGTAATGTTGTCTGGTACTGCGGTAAATGTTGCCACATAGCCGCCCAGCaa                                   >  1:1601761/1‑67 (MQ=255)                                        aCGTCTCAAGCGGATCGGTAATGTTGTCCGGCACCGCGGTa                                                             <  1:4439310/41‑1 (MQ=255)                                        aCGTCTCAAGCGGATCGGTAATGTTGTCCGGCACCGCGGTAAATGTTGCCACATAGCCGCCCAGCAAGAAc                               >  1:3017917/1‑71 (MQ=255)                                         cGTCTCAAGCGGATCGGTAATGTTGTCCGGCACCGCGGTAAATGTTGCCACATAGCCGCCCAGCAAGAAcg                              <  1:3884345/71‑1 (MQ=255)                                         cGTCTCAAGCGGATCGGTAATGTTGTCCGGCACCGCGGTAAATGTTGCCACATAGCCGCCCAGCAAGAAcg                              <  1:759902/71‑1 (MQ=255)                                          gTCTCAAGCGGATCGGTAATGTTGTCCGGCACCTCGGTAAATGTTGCCACATAGCCGCCCa                                       >  1:365496/1‑61 (MQ=255)                                          gTCTCAAGCGGATCGGTAATGTTGTCCGGCACCGCGGTAAATGTTGCCACATAGCCGCCCAGCAAGAAcgc                             >  1:2094868/1‑71 (MQ=255)                                           tctcAAGTGGATCGGTAATGTTGTCTGGTACTGCGGTAAATGTTGCCACAt                                                >  1:857072/1‑51 (MQ=255)                                           tctcAAGCGGATCGGTAATGTTGTCCGGCACCGCTGTa                                                             <  1:3964316/38‑1 (MQ=255)                                           tctcAAGCGGATCGGTAATGTTGTCCGGCACCGCGGTa                                                             <  1:2470818/38‑1 (MQ=255)                                           tctcAAGCGGATCGGTAATGTTGTCCGGCACCGCGGTAAATGTTGCCACATAGCCGCCCAGCAAGAAcgcg                            <  1:2851158/71‑1 (MQ=255)                                           tctcAAGCGGATCGGTAATGTTGTCCGGCACCGCGGTAAATGTTGCCACATAGCCGCCCAGCAAGAAcgcg                            >  1:1500796/1‑71 (MQ=255)                                             tcAAGCGGATCGGTAATGTTGTCCGGCACCGCGGTAAAt                                                          >  1:4173192/1‑39 (MQ=255)                                             tcAAGCGGATCGGTAATGTTGTCCGGCACCGCGGTAAATGTTGCCACATAGCCGCCCAGCAAGAACGCgg                           >  1:1823065/1‑70 (MQ=255)                                             tcAAGCGGATCGGTAATGTTGTCCGGCACCGCGGTAAATGTTGCCACATAGCCGCCCAGCAAGAACGCgg                           >  1:1414880/1‑70 (MQ=255)                                               aaGCGGATCGGTAATGTTGTCCGGCACCGCGGTAAATGTTGCca                                                   >  1:272680/1‑44 (MQ=255)                                               aaGCGGATCGGTAATGTTGTCCGGCACCGCGGTAAATGTTGCCACATAGCCGCCCAGCAAg                                  <  1:2573887/61‑1 (MQ=255)                                               aaGCGGATCGGTAATGTTGTCCGGCACCGCGGTAAATGTTGCCACATAGCCGCCCAGCAAGAAcgcg                            >  1:2387303/1‑67 (MQ=255)                                               aaGCGGATCGGTAATGTTGTCCGGCACCGCGGTAAATGTTGCCACATAGCCGCCCAGCAAGAACGCGGCa                         <  1:305800/70‑1 (MQ=255)                                               aaGCGGATCGGTAATGTTGTCCGGCACCGCGGTAAATGTTGCCACATAGCCGCCCAGCAAGAACGCGGCa                         <  1:853357/70‑1 (MQ=255)                                               aaGCGGATCGGTAATGTTGTCCGGCACCGCGGTAAATGTTGCCACATAGCCGCCCAGCAAGAACGCGGCAg                        <  1:524230/71‑1 (MQ=255)                                               aaGCGGATCGGTAATGTTGTCCGGCACCGCGGTAAATGTTGCCACATAGCCGCCCAGCAAGAACGCGGCAg                        >  1:3295283/1‑71 (MQ=255)                                               aaGCGGATCGGTAATGTTGTCCGGCACCGCGGTAAATGTTGCCACATAGCCGCCCAGCAAGAACGCGGCAg                        >  1:414878/1‑71 (MQ=255)                                                aGCGGATCGGTAATGTTGTCCGGCACCGCGGTAAATGTTGCCACATAGCCGCCCAGCAAGAACGCGGCAGc                       <  1:582622/71‑1 (MQ=255)                                                  cGGATCGGTAATGTTGTCCGGCACCGCGGTAAATGTTGCCACATAGCCGCCCAGCAAGAACGCGGCa                         >  1:1918217/1‑67 (MQ=255)                                                  cGGATCGGTAATGTTGTCCGGCACCGCGGTAAATGTTGCCACATAGCCGCCCAGCAAGAACGCGGCAGCCt                     <  1:3155459/71‑1 (MQ=255)                                                  cGGATCGGTAATGTTGTCCGGCACCGCCGTAAATGTTg                                                      >  1:114533/1‑38 (MQ=255)                                                     aTCGGTAATGTTGTCCGGCACCGCGGTAAATGTTGCCACATAGCCGCCCAGCAAGAACGCGGCAGCCTGCg                  <  1:4319031/71‑1 (MQ=255)                                                      tCGGTAATGTTGTCCGGCACCGCGGTAAATGTTGCCACATAGCCGCCCAGCAAGAACGCGGCAGCCTGCGt                 <  1:3396162/71‑1 (MQ=255)                                                      tCGGTAATGTTGTCCGGCACCGCGGTAAATGTTGCCACATAGCCGCCCAGCAAGAACGCGGCAGCCTGCGt                 <  1:186035/71‑1 (MQ=255)                                                       cGGTAATGTTGTCCGGCACCGCGGTAAATGTTGCCACATAGCCGCCCAGCAAGAACGCGGCAGCCTGCGt                 <  1:1039848/70‑1 (MQ=255)                                                       cGGTAATGTTGTCCGGCACCGCGGTAAATGTTGCCACATAGCCGCCCAGCAAGAACGCGGCAGCCTGCGt                 >  1:3364389/1‑70 (MQ=255)                                                         gaaaTGTTGTCCGGCACCGCGGTAAATGTTGCCACATAGCCGCCCAGCAAGAACGCGGCAGCCTGCGTCAg              <  1:3542976/69‑1 (MQ=255)                                                          tAATGTTGTCCGGCACCGCGGTAAATGTTGCCACATAGCCGCCCAg                                      <  1:1682238/46‑1 (MQ=255)                                                            aTGTTGTCCGGCACCGCGGTAAATGTTGCCACATAGCCGCCCAGCAAGaa                                <  1:2300941/50‑1 (MQ=255)                                                            aTGTTGTCCGGCACCGCGGTAAATGTTGCCACATAGCCGCCCAGCAAGAACGCGGCAGCCTGCGTCAGGaa           <  1:1227463/71‑1 (MQ=255)                                                            aTGTTGTCCGGCACCGCGGTAAATGTTGCCACATAGCCGCCCAGCAAGAACGCGGCAGCCTGCGTCAGGaa           <  1:1241793/71‑1 (MQ=255)                                                            aTGTTGTCCGGCACCGCGGTAAATGTTGCCACATAGCCGCCCAGCAAGAACGCGGCAGCCTGCGTCAGGa            <  1:1720201/70‑1 (MQ=255)                                                            aTGTTGTCCGGCACCGCGGTAAATGTTGCCACATAGCCGCCCAGCAAGAACGCGGCAGCCTGCGTCAGGa            <  1:2515223/70‑1 (MQ=255)                                                              gttgtCCGGCACCGCGGTAAATGTTGCCACATAGCCGCCCAGCAAGAACGCGGCAGCCTGCGTCa               >  1:3852678/1‑65 (MQ=255)                                                               ttgtCCGGCACCGCGGTAAATGTTGCCACATAGCCGCCCa                                       >  1:437763/1‑40 (MQ=255)                                                               ttgtCCGGCACCGCGGTAAATGTTGCCACATAGCCGCCCAGCAAGAACGCGGCAGCCTGCGTCAGGAAcc         <  1:4554849/70‑1 (MQ=255)                                                               ttgtCCGGCACCGCGGTAAATGTTGCCACATAGCCGCCCAGCAAGAACGCGGCAGCCTGCGTCAGGAACca        >  1:1957853/1‑71 (MQ=255)                                                               ttgtCCGGCACCGCGGTAAATGTTGCCACATAGCCGCCCAGCAAGAACGCGGCAGCCTGCGTCAGGAACca        >  1:133142/1‑71 (MQ=255)                                                                 gtCCGGCACCGCGGTAAATGTTGCCACATAGCCGCCCAGCAAGAACGCGGCAGCCTGCGTCAGGAACcac       >  1:1441264/1‑70 (MQ=255)                                                                  tCCGGCACCGCGGTAAATGTTGCCACATAGCCGCCCAGCAAGAACGCGGCAGCCTGCGTCAGGAACCACAt     <  1:4562015/71‑1 (MQ=255)                                                                     ggCACCGCGGTAAATGTTGCCACATAGCCGCCCAGCAAGAACGCGGCAGCCTGCGTCAGGAACCACATccc  <  1:3109409/71‑1 (MQ=255)                                                                       |                                                                     CCAGACCAATCTTACCAAACACGTTGGTATAGACGGGCAACGTCTCAAGCGGATCGGTAATGTTGTCCGGCACCGCGGTAAATGTTGCCACATAGCCGCCCAGCAAGAACGCGGCAGCCTGCGTCAGGAACCACATCCC  >  minE/2603281‑2603419

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