Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	454,658	T→G	12.1%	M143L (ATG→CTG)	ybgH ←	predicted transporter

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	454,658	0	T	G	12.1%	144.4 / 16.5	66	M143L (ATG→CTG)	ybgH	predicted transporter
Reads supporting (aligned to +/- strand):  ref base T (24/34);  new base G (6/2);  total (30/36)
Fisher's exact test for biased strand distribution p-value = 1.28e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CTGTACTCTTCCTGGGCGTAACCACAGGCGATAGGTGCGATAATAGACCCCACGTTACCCGCCGCATACATCAGCGAGAAACCGCCATCACGACGCGGATCGGTTGGCTCATACAGCTCGCCGAGCAGACAGCTTACGTT  >  minE/454588‑454727                                                                       |                                                                      cTGTACTCTTCCTGGGCGTAACCACAGGCGATAGGTGCGATAATAGACCCCACGTTACCCGCCGCATACAt                                                                       >  1:4451589/1‑71 (MQ=255)    tACTCTTCCTGGGCGTAACCACAGGCGATAGGTGCGATAATAGACCCCACGTTACCCGCCGCATACATCAg                                                                    >  1:4705219/1‑71 (MQ=255)          tCCTGGGCGTAACCACAGGCGATAGGTGCGATAATAGACCCCACGTTACCCGCCGCATACATCAGCgaga                                                               >  1:1479483/1‑70 (MQ=255)          tCCTGGGCGTAACCACAGGCGATAGGTGCGATAATAGACCCCACGTTACCCGCCGCATACATCAGCgaga                                                               >  1:3621723/1‑70 (MQ=255)          tCCTGGGCGTAACCACAGGCGATAGGTGCGATAATAGACCCCACGTTACCCGCCGCATACATCAGCgaga                                                               >  1:2600754/1‑70 (MQ=255)           ccTGGGCGTAACCACAGGCGATAGGTGCGATAATAGACCCCACGTTACCCGCCGCATACATCAGCgaga                                                               <  1:5217862/69‑1 (MQ=255)           ccTGGGCGTAACCACAGGCGATAGGTGCGATAATAGACCCCACGTTACCCGCCGCATACATCAGCgaga                                                               <  1:4953484/69‑1 (MQ=255)           ccTGGGCGTAACCACAGGCGATAGGTGCGATAATAGACCCCACGTTACCCGCCGCATACATCAGCgaga                                                               <  1:1668246/69‑1 (MQ=255)           ccTGGGCGTAACCACAGGCGATAGGTGCGATAATAGACCCCACGTTACCCGCCGCATACATCAGCgag                                                                <  1:1013263/68‑1 (MQ=255)            cTGGGCGTAACCACAGGCGATAGGTGCGATAATAGACCCCACGTTACCCGCCGCATACATCAGCGAGAAAc                                                            >  1:3839019/1‑71 (MQ=255)            cTGGGCGTAACCACAGGCGATAGGTGCGATAATAGACCCCACGTTACCCGCCGCATACATCAGCGAGAAAc                                                            <  1:4596040/71‑1 (MQ=255)            cTGGGCGTAACCACAGGCGATAGGTGCGATAATAGACCCCACGTTACCCGCCGCATACATCAGCGAGAAAc                                                            >  1:340751/1‑71 (MQ=255)                  gTAACCACAGGCGATAGGTGCGATAATAGACCCCACGTTACCCGCCGCATACATCAGCGAGAAAc                                                            <  1:779128/65‑1 (MQ=255)                  gTAACCACAGGCGATAGGTGCGATAATAGACCCCACGTTACCCGCCGCATACATCAGCGAGAAAc                                                            <  1:3825161/65‑1 (MQ=255)                  gTAACCACAGGCGATAGGTGCGATAATAGACCCCACGTTACCCGCCGCATACATCAGCGAGAAACCGCCAt                                                      <  1:5449280/71‑1 (MQ=255)                   tAACCACAGGCGATAGGTGCGATAATAGACCCCACGTTACCCGCCGCATACATCAGCGAGAAACCGCCAt                                                      >  1:5243062/1‑70 (MQ=255)                   tAACCACAGGCGATAGGTGCGATAATAGACCCCACGTTACCCGCCGCATACATCAGCGAGAAACCGCCAt                                                      >  1:934841/1‑70 (MQ=255)                     aCCACAGGCGATAGGTGCGATAATAGACCCCACGTTACCCGCCGCATACATCAGCGAGAAACCGCCAt                                                      <  1:654870/68‑1 (MQ=255)                         caGGCGATAGGTGCGATAATAGACCCCACTTTACCCGCCGCATACATCAGCGAGAAACCGCCATCacgacg                                               <  1:1461168/71‑1 (MQ=255)                               aTAGGTGCGATAATAGACCCCACGTTACCCGCCGCATACATCAGCGAGAAACCGCCATCACGACGCGGATc                                         <  1:504591/71‑1 (MQ=255)                                   gTGCGATAATAGACCCCACGTTACCCGCCGCATACATCAGCGAGAAACCGCCATCACGACGCGGATCGGtt                                     >  1:2602373/1‑71 (MQ=255)                                        ataataGACCCCACGTTACCCGCCGCATACATCAGCGAGAAACCGcc                                                        <  1:1588777/47‑1 (MQ=255)                                        ataataGACCCCACGTTACCCGCCGCATACATCAGCGAGAAACCGcc                                                        <  1:5299737/47‑1 (MQ=255)                                        ataataGACCCCACGTTACCCGCCGCATACATCAGCGAGAAACCGcc                                                        <  1:201918/47‑1 (MQ=255)                                         taataGACCCCACGTTACCCGCCGCATACATCAGCGAGAAACCGCCATCACGACGCGGATCGGTTGGCTCa                               >  1:4584737/1‑71 (MQ=255)                                         taataGACCCCACGTTACCCGCCGCATACATCAGCGAGAAACCGCCATCACGACGCGGATCGGTTGGCTCa                               <  1:1288495/71‑1 (MQ=255)                                         taataGACCCCACGTTACCCGCCGCATACATCAGCGAGAAACCGCCATCACGACGCGGATCGGTTGGCTCa                               >  1:4852968/1‑71 (MQ=255)                                          aataGACCCCACGTTACCCGCCGCATACATCAGCGAGaaa                                                             >  1:2647300/1‑40 (MQ=255)                                          aataGAACCCACGTTACCCGCCGCATACATCAGCGAGAAACCGCCATCACGCCGCGGATCGGt                                      >  1:733942/1‑63 (MQ=255)                                           ataGACCCCACGTTACCCGCCGCATACATCAGCGAGAAACCGCCATCACGACGCGGATCGGTTGGCTCAt                              <  1:5206795/70‑1 (MQ=255)                                              gACCCCACGTTACCCGCCGCATACATCAGCGAGAAACCGCCATCACGACGCGGATCGGTTGGCTCATACAg                          >  1:1597261/1‑71 (MQ=255)                                                 cccACGTTACCCGCCGCATACATCAGCGAGAAACCGCCATCac                                                   <  1:677528/43‑1 (MQ=255)                                                 cccACGTTACCCGCCGCATACATCAGCGAGAAACCGCCATCACGACGCGGATCGGTTGGc                                  <  1:2548731/60‑1 (MQ=255)                                                 cccACGTTACCCGCCGCATACATCAGCGAGAAACCGCCATCACGACGCGGATCGGTTGGCTCATACAGc                         <  1:3234855/69‑1 (MQ=255)                                                  ccACGTTACCCGCCGCATACATCAGCGAGAAACCGCCATCACGACGCgg                                            >  1:4515173/1‑49 (MQ=255)                                                  ccACGTTACCCGCCGCATACATCAGCGAGAAACCGCCATCACGACGCGGATCGGTTGGCTCATACAGCTc                       <  1:5219523/70‑1 (MQ=255)                                                   cACGTTACCCGCCGCATACATCAGCGAGAAACCGCCATCACGACGCGGAt                                          >  1:5310487/1‑50 (MQ=255)                                                     cGTTACCCGCCGCATACATCAGCGAGAAACCGCCATCACGACGCGGATCGGTTGGCTCATACAGCTCGcc                    >  1:347294/1‑70 (MQ=255)                                                     cGTTACCCGCCGCATACATCAGCGAGAAACCGCCATCACGACGCGGATCGGTTGGCTCATACAGCTCGCCg                   <  1:5280789/71‑1 (MQ=255)                                                      gTTACCCGCCGCATACATCAGCGAGAAACCGCCATCACGACGCGGATCGGtt                                     <  1:1217606/52‑1 (MQ=255)                                                       ttACCCGCCGCATACATCAGCGAGAAACCGCCATCACGACGCGGATCGGTTGGCTCa                               <  1:4590272/57‑1 (MQ=255)                                                       ttACCCGCCGCATACATCAGCGAGAAACCGCCATCACGACGCGGATCGGTTGGCTCATACAGCTCGCCGAg                 <  1:2071261/71‑1 (MQ=255)                                                        tACCCGCCGCATACATCAGCGAGAAACCGCCATCACGACGCGGATCGGTTg                                    >  1:98821/1‑51 (MQ=255)                                                        tACCCGCCGCATACATCAGCGAGAAACCGCCATCACGACGCGGATCGGTTGGCTCATa                             >  1:309591/1‑58 (MQ=255)                                                        tACCCGCCGCATACATCAGCGAGAAACCGCCATCACGACGCGGATCGGTTGGCTCATACAGCTCGCCGAg                 >  1:261339/1‑70 (MQ=255)                                                        tACCCGCCGCATACATCAGCGAGAAACCGCCATCACGACGCGGATCGGTTGGCTCATACAGCTCGCCGAg                 >  1:4880521/1‑70 (MQ=255)                                                         gcccGCAGCATACAGCAGCGAAAAACCGCCATCACGTCGATGATCGTTCTCGTCGTAGAGCTCGCCAAGc                >  1:1873060/2‑70 (MQ=255)                                                          cccGCCGCATACATCAGCGAGAAACCGCCATCACGACGCGGATCGGTTGGCTCATACAGCTCGCCGAGCa               <  1:2095070/70‑1 (MQ=255)                                                          cccGCCGCATACATCAGCGAGAAACCGCCATCACGACGCGGATCGGTTGGCTCATACAGCTCGCCGAGCa               <  1:3502401/70‑1 (MQ=255)                                                          cccGCCGCATACATCAGCGAGAAACCGCCATCACGACGCGGATCGGTTGGCTCATACAGCTCGCCGAGCa               <  1:3495535/70‑1 (MQ=255)                                                          cccGCAGCATACAGCAGCGAAAAACCGCCATCACGTCGATGATCGTTCTCGTCGTAGAGCTCGCCAAGCa               <  1:4836805/70‑1 (MQ=255)                                                           ccgccgCATACATCAGCGAGAAACCGCCATCACGACGCGGATCGGTTGGCTCa                               <  1:3817972/53‑1 (MQ=255)                                                               agcATACAGCAGCGAAAAACCGCCATCACGTCGATGATCGTTCTCGTCGTAGAGCTCGCCAAGCAAAc            <  1:5091729/67‑1 (MQ=255)                                                               agcATACAGCAGCGAAAAACCGCCATCACGTCGATGATCGTTCTCGTCGTAGAGCTCGCCAAGCAAACAg          >  1:4744271/2‑70 (MQ=255)                                                                 cATACAGCAGCGAAAAACCGCCATCACGTCGATGATCGTTCTCGTCGTAGAGCTCGCCAAGCAAACAGCt        >  1:4100189/1‑70 (MQ=255)                                                                  aTACAGCAGCGAAAAACCGCCATCACGTCGATGATCGTTCTCGTCGTAGAGCTCGCCAAGCAAACAGCtga      >  1:49635/1‑69 (MQ=255)                                                                   tACATCAGCGAGAAACCGCCATCACGACGCGGATCGGTTGGCTCa                               <  1:2830171/45‑1 (MQ=255)                                                                   tACATCAGCGAGAAACCGCCATCACGACGCGGATCGGTTGGCTCa                               <  1:4205409/45‑1 (MQ=255)                                                                   tACATCAGCGAGAAACCGCCATCACGACGCGGATCGGTTGGCTCATACAGCTCGCCGAg                 <  1:3852371/59‑1 (MQ=255)                                                                   tACATCAGCGAGAAACCGCCATCACGACGCGGATCGGTTGGCTCATACAGCTCGCCGAg                 <  1:5309415/59‑1 (MQ=255)                                                                   tACAGCAGCGAAAAACCGCCATCACGTCGATGATCGTTCTCGTCGTAGAGCTCGCCAAGCAAACAGCtgat     >  1:296372/1‑68 (MQ=255)                                                                   tACAGCAGCGAAAAACCGCCATCACGTCGATGATCGTTCTCGTCGTAGAGCTCGCCAAGCAAACAGCtgat     >  1:3753874/1‑68 (MQ=255)                                                                     cATCAGCGAGAAACCGCCATCACGACGCGGATCGGTTGGCTCATACAGCTCGCCGAGCa               >  1:1126806/1‑59 (MQ=255)                                                                     cATCAGCGAGAAACCGCCATCACGACGCGGATCGGTTGGCTCATACAGCTCGCCGAGCAGACAGCTTACGt   >  1:4004165/1‑71 (MQ=255)                                                                     cATCAGCGAGAAACCGCCATCACGACGCGGATCGGTTGGCTCATACAGCTCGCCGAGCAGACAGCTTACGt   >  1:1176694/1‑71 (MQ=255)                                                                      aTCAGCGAGAAACCGCCATCACGACGCGGATCGGTTGGCTCATACAGCTCGCCGAGCAGACAGCTTACGtt  <  1:3866004/71‑1 (MQ=255)                                                                      aTCAGCGAGAAAACGCCATCACGACGCGGATCGGTTGGCTCa                               <  1:5260130/42‑1 (MQ=255)                                                                       tCAGCGAGAAACCGCCATCACGACGCGGATCGGTTGGCTCATACAGCTCGCCGAGCa               >  1:5491708/1‑57 (MQ=255)                                                                       |                                                                      CTGTACTCTTCCTGGGCGTAACCACAGGCGATAGGTGCGATAATAGACCCCACGTTACCCGCCGCATACATCAGCGAGAAACCGCCATCACGACGCGGATCGGTTGGCTCATACAGCTCGCCGAGCAGACAGCTTACGTT  >  minE/454588‑454727

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