Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	2,080,563	G→A	6.6%	Q86Q (CAG→CAA)	exuT →	hexuronate transporter

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	2,080,563	0	G	A	6.6%	282.4 / 6.6	91	Q86Q (CAG→CAA)	exuT	hexuronate transporter
Reads supporting (aligned to +/- strand):  ref base G (50/35);  new base A (2/4);  total (52/39)
Fisher's exact test for biased strand distribution p-value = 3.96e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.75e-01

ACCTGACGCGTAACACTGTGGCGGCAGCTGCGCCAACTCTGATGGAAGAGTTAAACATCTCCACCCAACAGTATTCCTATATCATCGCAGCCTATTCTGCTGCTTATACGGTCATGCAACCGGTAGCAGGTTATGTACTG  >  minE/2080493‑2080632                                                                       |                                                                      aCCTGACGCGTAACACTGTGGCGGCAGCTGCGCCAACTCTGATGGAAGAGTTAAACATCTCCACCCAACAg                                                                       >  1:3759162/1‑71 (MQ=255)  ccTGACGCGTAACACTGTGGCGGCAGCTGCGCCAACTCTGATGGAAGAGTTAAACATCTCCACCCAACAg                                                                       <  1:248227/70‑1 (MQ=255)   cTGACGCGTAACACTGTGGCGGCAGCTGCGCCAACTCTGATGGAAGAGTTAAACATCTCCACCCAACAATa                                                                     >  1:3752379/1‑71 (MQ=255)    tGACGCGTAACACTGTGGCGGCAGCTGCGCCAACTCTGATGGAAGAGTTAAACATCTCCACCCAACAGTAt                                                                    >  1:2527212/1‑71 (MQ=255)       caggTAACACTGTGGCGGCAGCTGCGCCAACTCTGATGGAAGAGTTAAACATCTCCACCCAACAGTATTcc                                                                 >  1:899343/4‑71 (MQ=255)       caggTAACACTGTGGCGGCAGCTGCGCCAACTCTGATGGAAGAGTTAAACATCTCCACCCAACAGTATTc                                                                  >  1:2669722/4‑70 (MQ=255)            aaCACTGTGGCGGCAGCTGCGCCAACTCTGATGGAAGAGTTAAACATCTCCACCCAACAGTATTCCTatat                                                            >  1:3031764/1‑71 (MQ=255)            aaCACTGTGGCCGCAGCTGCGCCAACTCTGATGGAAGAGTTAAACATCTCCACCCAACAGTATTCCtat                                                              >  1:617001/1‑69 (MQ=255)             acacTGTGGCGGCAGCTGCGCCAACTCTGATGGAAGAGTTAAACATCTCCACCCAACAGTATTCCTatat                                                            >  1:1082871/1‑70 (MQ=255)              cacTGTGGCGGCAGCTGCGCCAACTCTGATGGAAGAGTTAAACATCTCCACCCAACAGTATTCCTATatca                                                          <  1:2283836/71‑1 (MQ=255)              cacTGTGGCGGCAGCTGCGCCAACTCTGATGGAAGAGTTAAACATCTCCACCCAACAGTATTCCTATatca                                                          <  1:3466624/71‑1 (MQ=255)                 tgtgGCGGCAGCTGCGCCAACTCTGATGGAAGAGTTAAACATCTCCACCCAACAGTATTCCTATATCATCg                                                       >  1:1918661/1‑71 (MQ=255)                   tgGCGGCAGCTGCGCCAACTCTGATGGAAGAGTTAAACATCTCCACCCAACAGTATTCCTATATCATCg                                                       <  1:614161/69‑1 (MQ=255)                   tgGCGGCAGCTGCGCCAACTCTGATGGAAGAGTTAAACATCTCCACCCAACAGTATTCCTATATCATCGc                                                      <  1:3919586/70‑1 (MQ=255)                   tgGCGGCAGCTGCGCCAACTCTGATGGAAGAGTTAAACATCTCCACCCAACAATATTCCTATATCATCGc                                                      <  1:1163768/70‑1 (MQ=255)                    ggcggcAGCTGCGCCAACTCTGATGGAAGAGTTAAACATCTCCACCCAACAGTATTCCTATatc                                                           <  1:959086/64‑1 (MQ=255)                     gcggcAGCTGCGCCAACTCTGATGGAAGAGTTAAACATCTCCACCCAACAGTATTCCTATatcatc                                                        >  1:1480852/1‑66 (MQ=255)                     gcggcAGCTGCGCCAACTCTGATGGAAGAGTTAAACATCTCCACCCAACAGTATTCCTATatcatc                                                        >  1:2173637/1‑66 (MQ=255)                     gcggcAGCTGCGCCAACTCTGATGGAAGAGTTAAACATCTCCACCCAACAGTATTCCTATATCATCGCAGc                                                   >  1:3328367/1‑71 (MQ=255)                        gcAGCTGCGCCAACTCTGATGGAAGAGTTAAACATCTCCACCCAACAGTATTCCt                                                                <  1:2780061/55‑1 (MQ=255)                        gcAGCTGCGCCAACTCTGATGGAAGAGTTAAACATCTCCACCCAACAGTATTCCTATATCATCGc                                                      >  1:3469609/1‑65 (MQ=255)                        gcAGCTGCGCCAACTCTGATGGAAGAGTTAAACATCTCCACCCAACAGTATTCCTATATCATCGCa                                                     >  1:642391/1‑66 (MQ=255)                        gcAGCTGCGCCAACTCTGATGGAAGAGTTAAACATCTCCACCCAACAGTATTCCTATATCATCGCAGCCTa                                                >  1:2527910/1‑71 (MQ=255)                        gcAGCTGCGCCAACTCTGATGGAAGAGTTAAACATCTCCACCCAACAGTATTCCTATATCATCGCAGCCTa                                                <  1:2174572/71‑1 (MQ=255)                        gcAGCTGCGCCAACTCTGATGGAAGAGTTAAACATCTCCACCCAACAGTATTCCTATATCATCGCAGCCTa                                                <  1:3418279/71‑1 (MQ=255)                         cAGCTGCGCCAACTCTGATGGAAGAGTTAAACATCTCCACCCAACAGTATTCCTATATCATCGCAGCCt                                                 >  1:1338041/1‑69 (MQ=255)                         cAGCTGCGCCAACTCTGATGGAAGAGTTAAACATCTCCACCCAACAATATTCCTATATCATCGCAGCCTAt                                               >  1:2644982/1‑71 (MQ=255)                          aGCTGCGCCAACTCTGATGGAAGAGTTAAACATCTCCACCCAACAGTATTCCTATATCATCGCAGCCTAt                                               <  1:2579501/70‑1 (MQ=255)                              gcgcCAACTCTGATGGAAGAGTTAAACATCTCCACCCAACAGTATTCCTATATCATCGCAGCCTATTct                                            >  1:17007/1‑69 (MQ=255)                               cgcCAACTCTGATGGAAGAGTTAAACATCTCCACCCAACAGTATTCCTATATCATCGCAGCCTATTctg                                           <  1:3144240/69‑1 (MQ=255)                                  cAACTCTGATGGAAGAGTTAAACATCTCCACCCAACAGTATTCCTATATCATCGCAGCCTATTCTgctgct                                      >  1:1044615/1‑71 (MQ=255)                                   aaCTCTGATGGAAGAGTTAAACATCTCCACCCAACAGTATTCCTATATCATCGCAGCCTATTCTgctgct                                      >  1:3741577/1‑70 (MQ=255)                                   aaCTCTGATGGAAGAGTTAAACATCTCCACCCAACAGTATTCCTATATCATCGCAGCCTATTCTgctgct                                      >  1:2518872/1‑70 (MQ=255)                                   aaCTCTGATGGAAGAGTTAAACATCTCCACCCAACAGTATTCCTATATCATCGCAGCCTATTCTGCTGCtt                                     <  1:2083412/71‑1 (MQ=255)                                    aCTCTGATGGAAGAGTTAAACATCTCCACCCAACAGt                                                                      >  1:1799607/1‑37 (MQ=255)                                     ctctGATGGAAGAGTTAAACATCTCCACCCAACAGTATTCCTATATCATCGCAGCCTATTCTGCTGCTta                                    >  1:2818532/1‑70 (MQ=255)                                      tctGATGGAAGAGTTAAACATCTCCACCCAACAGTATTCCTATATCATCGCAGCCTATTCTGCTGCTtata                                  >  1:1299460/1‑71 (MQ=255)                                        tGATGGAAGAGTTAAACATCTCCACCCAACAGTATTCCTATATCATCGCAGCCTATTCTGCTGCTtata                                  <  1:462992/69‑1 (MQ=255)                                        tGATGGAAGAGTTAAACATCTCCACCCAACAGTATTCCTATATCATCGCAGCCTATTCTGCTGCTTATACg                                >  1:412964/1‑71 (MQ=255)                                         gATGGAAGAGTTAAACATCTCCACCCAACAGTATTCCTATATCATCGCa                                                     <  1:2442253/49‑1 (MQ=255)                                         gATGGAAGAGTTAAACATCTCCACCCAACAGTATTCCTATATCATCGCa                                                     <  1:409863/49‑1 (MQ=255)                                          aTGGAAGAGTTAAACATCTCCACCCAACAGTATTCCTATATCATCGCAGCCTATTCTGCTGCTTATAc                                 <  1:3088413/68‑1 (MQ=255)                                          aTGGAAGAGTTAAACATCTCCACCCAACAGTATTCCTATATCATCGCAGCCTATTCTGCTGCTTATACGGt                              <  1:486051/71‑1 (MQ=255)                                          aTGGAAGAGTTAAACATCTCCACCCAACAATATTCCTATATCATCGCAGCCTATTCTGCTGCTTATACgg                               <  1:1582314/70‑1 (MQ=255)                                           tGGAAGAGTTAAACATCTCCACCCAACAGTATTCCTATATCATCGCAGCCTATTCTGCTGCTTATACGGt                              >  1:1984846/1‑70 (MQ=255)                                           tGGAAGAGTTAAACATCTCCACCCAACAGTATTCCTATATCATCGCAGCCTATTCTGCTGCTTATACGGTc                             <  1:2420341/71‑1 (MQ=255)                                           tGGAAGAGTTAAACATCTCCACCCAACAGTATTCCTATATCATCGCAGCCTATTCTGCTGCTTATACGGTc                             <  1:663734/71‑1 (MQ=255)                                           tGGAAGAGTTAAACATCTCCACCCAACAGTATTCCTATATCATCGCAGCCTATTCTGCTGCTTATACGGTc                             <  1:980240/71‑1 (MQ=255)                                           tGGAAGAGTTAAACATCTCCACCCAACAGTATTCCTATATCATCGCAGCCTATTCTGCTGCTTATACGGTc                             <  1:1823688/71‑1 (MQ=255)                                             gAAGAGTTAAACATCTCCACCCAACAGTATTCCTATatc                                                           >  1:273511/1‑39 (MQ=255)                                               agagTTAAACATCTCCACCCAACAGTATTCCTATATCATCGCAGCCTATTctgctg                                        >  1:193968/1‑56 (MQ=255)                                                gagTTAAACATCTCCACCCAACAGTATTCCTATATCATCg                                                       >  1:3240321/1‑40 (MQ=255)                                                gagTTAAACATCTCCACCCAACAGTATTCCTATATCATCGCAGc                                                   >  1:822254/1‑44 (MQ=255)                                                gagTTAAACATCTCCACCCAACAGTATTCCTATATCATCGCAGc                                                   >  1:880648/1‑44 (MQ=255)                                                gagTTAAACATCTCCACCCAACAGTATTCCTATATCATCGCAGc                                                   >  1:3750805/1‑44 (MQ=255)                                                gagTTAAACATCTCCACCCAACAGTATTCCTATATCATCGCAGCCTATTCTGCTGCTTATACGGTCATGc                         >  1:2946167/1‑70 (MQ=255)                                                gagTTAAACATCTCCACCCAACAGTATTCCTATATCATCGCAGCCTATTCTGCTGCTTATACGGTCATGc                         >  1:3442726/1‑70 (MQ=255)                                                gagTTAAACATCTCCACCCAACAGTATTCCTATATCATCGCAGCCTATTCTGCTGCTTATACGGTCATGCa                        <  1:3252216/71‑1 (MQ=255)                                                gagTTAAACATCTCCACCCAACAGTATTCCTATATCATCGCAGCCTATTCTGCTGCTTATACGGTCATGCa                        <  1:2246171/71‑1 (MQ=255)                                                gagTTAAACATCTCCACCCAACAGTATTCCTATATCATCGCAGCCTATTCTGCTGCTTATACGGTCATGCa                        <  1:2745120/71‑1 (MQ=255)                                                gagTTAAACATCTCCACCCAACAGTATTCCTATATCATCGCAGCCTATTCTGCTGCTTATACGGTCATGCa                        <  1:2241555/71‑1 (MQ=255)                                                gagTTAAACATCTCCACCCAACAATATTCCTATATCATCGCAGCCTATTCTGCTGCTTATACGGTCAt                           <  1:1248732/68‑1 (MQ=255)                                                  gTTAAACATCTCCACCCAACAGTATTCCTATATCATCGCAGCCTATTCTGCTGCTTATACGGTCATGCaa                       <  1:2880581/70‑1 (MQ=255)                                                     aaaCATCTCCACCCAACAGTATTCCTATATCATCGCAGcc                                                  >  1:2710079/1‑40 (MQ=255)                                                     aaaCATCTCCACCCAACAGTATTCCTATATCATCGCAGcc                                                  >  1:3372069/1‑40 (MQ=255)                                                     aaaCATCTCCACCCAACAGTATTCCTATATCATCGCAGCCTATTCTGCTGCTta                                    >  1:12642/1‑54 (MQ=255)                                                     aaaCATCTCCACCCAACAGTATTCCTATATCATCGCAGCCTATTCTGCTGCTTATACGGTCATGCAACCgg                   >  1:1073533/1‑71 (MQ=255)                                                     aaaCATCTCCACCCAACAGTATTCCTATATCATCGCAGCCTATTCTGCTGCTTATACGGTCATGCAACCgg                   >  1:2082348/1‑71 (MQ=255)                                                      aaCATCTCCACCCAACAGTATTCCTATATCATCGCAGCCTATTCTGCTGCTTATACGGTCATGCAACCgg                   <  1:3483823/70‑1 (MQ=255)                                                      aaCATCTCCACCCAACAGTATTCCTATATCATCGCAGCCTATTCTGCTGCTTATACGGTCATGCAACCgg                   <  1:2494542/70‑1 (MQ=255)                                                      aaCATCTCCACCCAACAGTATTCCTATATCATCGCAGCCTATTCTGCTGCTTATACGGTCATGCAACCgg                   >  1:2856166/1‑70 (MQ=255)                                                       aCATCTGCACCCAACAGTATTCCTATATCATCGCAGCCTATTCTGCTGCTTATACGGTCATGCAACCgg                   >  1:505006/1‑69 (MQ=255)                                                         aTCTCCACCCAACAGTATTCCTATATCATCGCAGCCTATTCTGCTGCTTATACGGTCATGCAACCGGTAg                >  1:485491/1‑70 (MQ=255)                                                           ctcCACCCAACAGTATTCCTATATCATCGCAGCCTATTctgct                                         <  1:1493532/43‑1 (MQ=255)                                                            tcCACCCAACAGTATTCCTATATCATCGCAGCCTAtt                                              >  1:3541953/1‑37 (MQ=255)                                                             ccACCCAACAGTATTCCTATATCATCGCAGCCTATTCtgctgc                                       >  1:1971261/1‑43 (MQ=255)                                                              cACCCAACAGTATTCCTATATCATCGCAGCCTATTCTGCTGCTTATACGGTCATGCAACCGGTAGCAGGt           <  1:1955568/70‑1 (MQ=255)                                                                cccAACAGTATTCCTATATCATCGCAGCCTATTCTGCTGCTTATACGGTCATGCAACCGGTAGCAGGTTa         >  1:99442/1‑70 (MQ=255)                                                                cccAACAGTATTCCTATATCATCGCAGCCTATTCTGCTGCTTATACGGTCATGCAACCGGTAGCAGGTTAt        <  1:1449690/71‑1 (MQ=255)                                                                  cAACAGTATTCCTATATCATCGCAGCCTATTCTGCTGCTTATACGGTCATGCAACCGGTAGCAGGTTATg       >  1:228150/1‑70 (MQ=255)                                                                   aaCAGTATTCCTATATCATCGCAGCCTATTCTGCTGCTTATACGGTCATGCAACCGGTAGc               <  1:1679900/61‑1 (MQ=255)                                                                   aaCAGTATTCCTATATCATCGCAGCCTATTCTGCTGCTTATACGGTCATGCAACCGGTAGCAGGTTATGTa     >  1:1617313/1‑71 (MQ=255)                                                                   aaCAGTATTCCTATATCATCGCAGCCTATTCTGCTGCTTATACGGTCATGCAACCGGTAGCAGGTTATGTa     <  1:544711/71‑1 (MQ=255)                                                                   aaCAGTATTCCTATATCATCGCAGCCTATTCTGCTGCTTATACGGTCATGCAACCGGTAGCAGGTTATGTa     >  1:3932903/1‑71 (MQ=255)                                                                   aaCAATATTCCTATATCATCGCAGCCTATTCTGCTGCTTATACGGTCATGCAAc                      <  1:41171/54‑1 (MQ=255)                                                                     cAGTATTCCTATATCATCGCAGCCTATTCTGCTGCTTATACGGTc                             >  1:3301219/1‑45 (MQ=255)                                                                     cAGTATTCCTATATCATCGCAGCCTATTCTGCTGCTTATACGGTCATGCAACCGGTAGCAGGTTATGTACt   >  1:3202040/1‑71 (MQ=255)                                                                     cAGTATTCCTATATCATCGCAGCCTATTCTGCTGCTTATACGGTCATGCAACCGGTAGCAGGTTATGTACt   >  1:1264853/1‑71 (MQ=255)                                                                     cAGTATTCCTATATCATCGCAGCCTATTCTGCTGCTTATACGGTCATGCAACCGGTAGCAGGTTATGTACt   >  1:2216457/1‑71 (MQ=255)                                                                     cAGTATTCCTATATCATCGCAGCCTATTCTGCTGCTTATACGGTCATGCAACCGGTAGCAGGTTATGTACt   >  1:381967/1‑71 (MQ=255)                                                                      aGTATTCCTATATCATCGCAGCCTATTCTGCTGCTTATACGGTCATGCAACCGGTAGCAgg            <  1:296771/61‑1 (MQ=255)                                                                      aGTATTCCTATATCATCGCAGCCTATTCTGCTGCTTATACGGTCATGCAACCGGTAGCAgg            <  1:2913182/61‑1 (MQ=255)                                                                      aGTATTCCTATATCATCGCAGCCTATTCTGCTGCTTATACGGTCATGCAACCGGTAGCAGGTTATGTACTg  <  1:196223/71‑1 (MQ=255)                                                                       gTATTCCTATATCATCGCAGCCTATTCTGCTGCTTATACGGTCATGCAACCGGTAGc               <  1:2320495/57‑1 (MQ=255)                                                                       |                                                                      ACCTGACGCGTAACACTGTGGCGGCAGCTGCGCCAACTCTGATGGAAGAGTTAAACATCTCCACCCAACAGTATTCCTATATCATCGCAGCCTATTCTGCTGCTTATACGGTCATGCAACCGGTAGCAGGTTATGTACTG  >  minE/2080493‑2080632

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