Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	2,603,354	C→T	9.0%	A439A (GCG→GCA)	yhiP ←	predicted transporter

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	2,603,354	0	C	T	9.0%	262.8 / 11.3	89	A439A (GCG→GCA)	yhiP	predicted transporter
Reads supporting (aligned to +/- strand):  ref base C (29/52);  new base T (6/2);  total (35/54)
Fisher's exact test for biased strand distribution p-value = 5.31e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.26e-01

ACCAATCTTACCAAACACGTTGGTATAGACGGGCAACGTCTCAAGCGGATCGGTAATGTTGTCCGGCACCGCGGTAAATGTTGCCACATAGCCGCCCAGCAAGAACGCGGCAGCCTGCGTCAGGAACCACATCCCGAGAA  >  minE/2603285‑2603424                                                                      |                                                                       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:1648742/68‑1 (MQ=255)          aCCAAACACGTTGGTATAGACGGGCAACGTCTCAAGCGGATCGGTAATGTTGTCCGGCACCGCGGTaa                                                                 <  1:1843319/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:1982399/70‑1 (MQ=255)          aCCAAACACGTTGGTATAGACGGGCAACGTCTCAAGCGGATCGGTAATGTTGTCCGGCACCGCGGTAAAt                                                               <  1:184411/70‑1 (MQ=255)           ccAAACACGTTGGTATAGACGGGCAACGTCTCAAGCGGATCGGTAATGTTGTCCGGCACCGCGGTAAATGt                                                             >  1:4010549/1‑71 (MQ=255)            cAAACACGTTGGTATAGACGGGCAACGTCTCAAGCGGATCGGTAATGTTGTCCGGCACCGCGGTa                                                                  <  1:3939507/65‑1 (MQ=255)              aaCACGTTGGTATAGACGGGCAACGTCTCAAGCGGATCGGTAATGTTGTCCGGCACCGCGGTAAATGTTGc                                                          >  1:3721277/1‑71 (MQ=255)              aaCACGTTGGTATAGACGGGCAACGTCTCAAGCGGATCGGTAATGTTGTCCGGCACCGCGGTAAATGTTGc                                                          >  1:2288308/1‑71 (MQ=255)              aaCACGTTGGTATAGACGGGCAACGTCTCAAGCGGATCGGTAATGTTGTCCGGCACCGCGGTAAATGTTGc                                                          <  1:3940948/71‑1 (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:1694433/71‑1 (MQ=255)                                 gCAACGTCTCAAGCGGATCGGTAATGTTGTCCGGCACCGCGGTAAATGTTGCCACATAGCCGCCCAGCAAg                                       <  1:3267562/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:3062828/1‑67 (MQ=255)                                    aCGTCTCAAGTGGATCGGTAATGTTGTCTGGTACTGCGGTAAATGTTGCCACATAGCCGCCCAGCaa                                        >  1:1377851/1‑67 (MQ=255)                                    aCGTCTCAAGTGGATCGGTAATGTTGTCTGGTACTGCGGTAAATGTTGCCACATAGCCGCCCAGCaa                                        >  1:2974053/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:759902/71‑1 (MQ=255)                                     cGTCTCAAGCGGATCGGTAATGTTGTCCGGCACCGCGGTAAATGTTGCCACATAGCCGCCCAGCAAGAAcg                                   <  1:3884345/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:1414880/1‑70 (MQ=255)                                         tcAAGCGGATCGGTAATGTTGTCCGGCACCGCGGTAAATGTTGCCACATAGCCGCCCAGCAAGAACGCgg                                >  1:1823065/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:853357/70‑1 (MQ=255)                                           aaGCGGATCGGTAATGTTGTCCGGCACCGCGGTAAATGTTGCCACATAGCCGCCCAGCAAGAACGCGGCa                              <  1:305800/70‑1 (MQ=255)                                           aaGCGGATCGGTAATGTTGTCCGGCACCGCGGTAAATGTTGCCACATAGCCGCCCAGCAAGAACGCGGCAg                             >  1:3295283/1‑71 (MQ=255)                                           aaGCGGATCGGTAATGTTGTCCGGCACCGCGGTAAATGTTGCCACATAGCCGCCCAGCAAGAACGCGGCAg                             >  1:414878/1‑71 (MQ=255)                                           aaGCGGATCGGTAATGTTGTCCGGCACCGCGGTAAATGTTGCCACATAGCCGCCCAGCAAGAACGCGGCAg                             <  1:524230/71‑1 (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:2515223/70‑1 (MQ=255)                                                        aTGTTGTCCGGCACCGCGGTAAATGTTGCCACATAGCCGCCCAGCAAGAACGCGGCAGCCTGCGTCAGGa                 <  1:1720201/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)                                                                     ccGCGGTAAATGTTGCCACATAGCCGCCCAGCAAGAACGCGGCAGCCTGCGTCAGGAAc               >  1:742224/1‑59 (MQ=255)                                                                      cgcgGTAAATGTTGCCACATAGCCGCCCAGCAAGAACGCGGCAGCCTGCGTCAGGAACCACATCCCGAGaa  <  1:2897975/71‑1 (MQ=255)                                                                      |                                                                       ACCAATCTTACCAAACACGTTGGTATAGACGGGCAACGTCTCAAGCGGATCGGTAATGTTGTCCGGCACCGCGGTAAATGTTGCCACATAGCCGCCCAGCAAGAACGCGGCAGCCTGCGTCAGGAACCACATCCCGAGAA  >  minE/2603285‑2603424

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