Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	303,058	G→A	32.0%	D191D (GAC→GAT)	tesB ←	acyl‑CoA thioesterase II

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	303,058	0	G	A	32.0%	74.3 / 40.1	84	D191D (GAC→GAT)	tesB	acyl‑CoA thioesterase II
Reads supporting (aligned to +/- strand):  ref base G (24/33);  new base A (12/15);  total (36/48)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CAGGAAGTTAAGATCAGAAGCGTAACCGAGCAGATACTGATGAACGCGCAGGTCATCCGGCACGCTACCATTTGCGCGGATCCACACCTGACGATGTGGTTCTGCGACGTGAC  >  minE/303007‑303119                                                    |                                                              cAGGAAGTTAAGATCAGAAGCGTAACCGAGCAGATACTGATGAACGCGCAGGTCATCCGGCACGCTAcc                                              <  1:2867206/69‑1 (MQ=255)         tAAGATCAGAAGCGTAACCGAGCAGATACTGATGAACGCGCAGGTCATCCGGCACGCTACCATTTgcg                                       >  1:968107/1‑68 (MQ=255)         tAAGATCAGAAGCGTAACCGAGCAGATACTGATGAACGCGCAGGTCATCCGGCACGCTACCATTTgcg                                       >  1:861970/1‑68 (MQ=255)         tAAGATCAGAAGCGTAACCGAGCAGATACTGATGAACGCGCAGGTCATCCGGCACGCTACCATTTgcg                                       >  1:491556/1‑68 (MQ=255)         tAAGATCAGAAGCGTAACCGAGCAGATACTGATGAACGCGCAGGTCATCCGGCACGCTACCATTTgcg                                       >  1:419406/1‑68 (MQ=255)         tAAGATCAGAAGCGTAACCGAGCAGATACTGATGAACGCGCAGGTCATCCGGCACGCTACCATTTgcg                                       >  1:369033/1‑68 (MQ=255)         tAAGATCAGAAGCGTAACCGAGCAGATACTGATGAACGCGCAGGTCATCCGGCACGCTACCATTTgcg                                       >  1:3228233/1‑68 (MQ=255)         tAAGATCAGAAGCGTAACCGAGCAGATACTGATGAACGCGCAGGTCATCCGGCACGCTACCATTTgcg                                       >  1:1792018/1‑68 (MQ=255)         tAAGATCAGAAGCGTAACCGAGCAGATACTGATGAACGCGCAGGTCATCCGGCACGCTACCATTTgcg                                       >  1:2088219/1‑68 (MQ=255)         tAAGATCAGAAGCGTAACCGAGCAGATACTGATGAACGCGCAGGTCATCCGGCACGCTACCATTTgcg                                       >  1:2696236/1‑68 (MQ=255)         tAAGATCAGAAGCGTAACCGAGCAGATACTGATGAACGCGCAGGTCATCCGGCACGCTACCATTTgcg                                       >  1:1506718/1‑68 (MQ=255)         tAAGATCAGAAGCGTAACCGAGCAGATACTGATGAACGCGCAGGTCATCCGGCACGCTACCATTTgcg                                       >  1:2280196/1‑68 (MQ=255)                              gCAGATACTGATGAACGCGCAGATCATCCGGCACGCTACCAt                                            >  1:746174/1‑42 (MQ=255)                              gCAGATACTGATGAACGCGCAGATCATCCGGCACGCTACCATTTGCGCGGATCCACACCTGACGAtgcg                 >  1:2035747/1‑67 (MQ=255)                              gCAGATACTGATGAACGCGCAGATCATCCGGCACGCTACCATTTGCGCGGATCCACACCTGACGAtgcg                 >  1:1093473/1‑67 (MQ=255)                              gCAGATACTGATGAACGCGCAGATCATCCGGCACGCTACCATTTGCGCGGATCCACACCTGACGAtgcg                 >  1:2185183/1‑67 (MQ=255)                              gCAGATACTGATGAACGCGCAGATCATCCGGCACGCTACCATTTGCGCGGATCCACACCTGACGAtgcg                 >  1:2402553/1‑67 (MQ=255)                              gCAGATACTGATGAACGCGCAGATCATCCGGCACGCTACCATTTGCGCGGATCCACACCTGACGAtgcg                 >  1:2540670/1‑67 (MQ=255)                              gCAGATACTGATGAACGCGCAGATCATCCGGCACGCTACCATTTGCGCGGATCCACACCTGACGAtgcg                 >  1:2993874/1‑67 (MQ=255)                              gCAGATACTGATGAACGCGCAGATCATCCGGCACGCTACCATTTGCGCGGATCCACACCTGACGAtgcg                 >  1:3118205/1‑67 (MQ=255)                              gCAGATACTGATGAACGCGCAGATCATCCGGCACGCTACCATTTGCGCGGATCCACACCTGACGAtgcg                 >  1:3237213/1‑67 (MQ=255)                              gCAGATACTGATGAACGCGCAGATCATCCGGCACGCTACCATTTGCGCGGATCCACACCTGACGAtgcg                 >  1:1101636/1‑67 (MQ=255)                              gCAGATACTGATGAACGCGCAGATCATCCGGCACGCTACCATTTGCGCGGATCCACACCTGACGAtgcg                 >  1:1420078/1‑67 (MQ=255)                              gCAGATACTGATGAACGCGCAGATCATCCGGCACGCTACCATTTGCGCGGATCCACACCTGACGAtgcg                 >  1:1622812/1‑67 (MQ=255)                                aGATACTGATGAACGCGCAGGTCATCCGGCACGCTACCATTTGCGCGGATCCACACCTGACGATGTGGt               <  1:2942130/69‑1 (MQ=255)                                aGATACTGATGAACGCGCAGGTCATCCGGCACGCTACCATTTGCGCGGATCCACACCTGACGATGTGGt               <  1:1068201/69‑1 (MQ=255)                                aGATACTGATGAACGCGCAGGTCATCCGGCACGCTACCATTTGCGCGGATCCACACCTGACGATGTGGt               <  1:1112356/69‑1 (MQ=255)                                aGATACTGATGAACGCGCAGGTCATCCGGCACGCTACCATTTGCGCGGATCCACACCTGACGATGTGGt               <  1:3147161/69‑1 (MQ=255)                                aGATACTGATGAACGCGCAGGTCATCCGGCACGCTACCATTTGCGCGGATCCACACCTGACGATGTGGt               <  1:2945936/69‑1 (MQ=255)                                aGATACTGATGAACGCGCAGGTCATCCGGCACGCTACCATTTGCGCGGATCCACACCTGACGATGTGGt               <  1:2906028/69‑1 (MQ=255)                                aGATACTGATGAACGCGCAGGTCATCCGGCACGCTACCATTTGCGCGGATCCACACCTGACGATGTGGt               <  1:2713862/69‑1 (MQ=255)                                aGATACTGATGAACGCGCAGGTCATCCGGCACGCTACCATTTGCGCGGATCCACACCTGACGATGTGGt               <  1:2177785/69‑1 (MQ=255)                                    aCTGATGAACGCGCAGGTCATCCGGCACGCTACCATTTGCGCGGATCCACACCTGACGATGTGGTTCt            >  1:2223021/1‑68 (MQ=255)                                    aCTGATGAACGCGCAGGTCATCCGGCACGCTACCATTTGCGCGGATCCACACCTGACGATGTGGTTCt            >  1:1119662/1‑68 (MQ=255)                                    aCTGATGAACGCGCAGGTCATCCGGCACGCTACCATTTGCGCGGATCCACACCTGACGATGTGGTTCt            >  1:1240142/1‑68 (MQ=255)                                    aCTGATGAACGCGCAGGTCATCCGGCACGCTACCATTTGCGCGGATCCACACCTGACGATGTGGTTCt            >  1:3105058/1‑68 (MQ=255)                                    aCTGATGAACGCGCAGGTCATCCGGCACGCTACCATTTGCGCGGATCCACACCTGACGATGTGGTTCt            >  1:2878620/1‑68 (MQ=255)                                    aCTGATGAACGCGCAGGTCATCCGGCACGCTACCATTTGCGCGGATCCACACCTGACGATGTGGTTCt            >  1:2061708/1‑68 (MQ=255)                                    aCTGATGAACGCGCAGGTCATCCGGCACGCTACCATTTGCGCGGATCCACACCTGACGATGTGGTTCt            >  1:150081/1‑68 (MQ=255)                                    aCTGATGAACGCGCAGGTCATCCGGCACGCTACCATTTGCGCGGATCCACACCTGACGATGTGGTTCt            >  1:2606718/1‑68 (MQ=255)                                    aCTGATGAACGCGCAGGTCATCCGGCACGCTACCATTTGCGCGGATCCACACCTGACGATGTGGTTCt            >  1:2631960/1‑68 (MQ=255)                                    aCTGATGAACGCGCAGGTCATCCGGCACGCTACCATTTGCGCGGATCCACACCTGACGATGTGGTTCt            >  1:949379/1‑68 (MQ=255)                                          gaACGCGCAGGTCATCCGGCACGCTACCATTTGCGCGGATCCACACCTGACGATGTGGTTCTGCTACGt     <  1:2720363/69‑1 (MQ=255)                                          gaACGCGCAGGTCATCCGGCACGCTACCATTTGCGCGGATCCACACCTGACGATGTGGTTCTGCGACGt     <  1:894715/69‑1 (MQ=255)                                          gaACGCGCAGGTCATCCGGCACGCTACCATTTGCGCGGATCCACACCTGACGATGTGGTTCTGCGACGt     <  1:367495/69‑1 (MQ=255)                                          gaACGCGCAGGTCATCCGGCACGCTACCATTTGCGCGGATCCACACCTGACGATGTGGTTCTGCGACGt     <  1:433660/69‑1 (MQ=255)                                          gaACGCGCAGGTCATCCGGCACGCTACCATTTGCGCGGATCCACACCTGACGATGTGGTTCTGCGACGt     <  1:2137572/69‑1 (MQ=255)                                          gaACGCGCAGATCATCCGGCACGCTACCATTTGCGCGGATCCACACCTGACGATGCGGTTCTGCGACGt     <  1:1130791/69‑1 (MQ=255)                                          gaACGCGCAGATCATCCGGCACGCTACCATTTGCGCGGATCCACACCTGACGATGCGGTTCTGCGACGt     <  1:1662781/69‑1 (MQ=255)                                          gaACGCGCAGATCATCCGGCACGCTACCATTTGCGCGGATCCACACCTGACGATGCGGTTCTGCGACGt     <  1:43773/69‑1 (MQ=255)                                          gaACGCGCAGATCATCCGGCACGCTACCATTTGCGCGGATCCACACCTGACGATGCGGTTCTGCGACGt     <  1:94334/69‑1 (MQ=255)                                          gaACGCGCAGATCATCCGGCACGCTACCATTTGCGCGGATCCACACCTGACGATGCGGTTCTGCGACGt     <  1:61664/69‑1 (MQ=255)                                          gaACGCGCAGATCATCCGGCACGCTACCATTTGCGCGGATCCACACCTGACGATGCGGTTCTGCGACGt     <  1:2920777/69‑1 (MQ=255)                                          gaACGCGCAGATCATCCGGCACGCTACCATTTGCGCGGATCCACACCTGACGATGCGGTTCTGCGACGt     <  1:220760/69‑1 (MQ=255)                                          gaACGCGCAGATCATCCGGCACGCTACCATTTGCGCGGATCCACACCTGACGATGCGGTTCTGCGACGt     <  1:2408064/69‑1 (MQ=255)                                          gaACGCGCAGATCATCCGGCACGCTACCATTTGCGCGGATCCACACCTGACGATGCGGTTCTGCGACGt     <  1:2760041/69‑1 (MQ=255)                                          gaACGCGCAGATCATCCGGCACGCTACCATTTGCGCGGATCCACACCTGACGATGCGGTTCTGCGACGt     <  1:2754100/69‑1 (MQ=255)                                           aaCGCGCAGGTCATCCGGCACGCTACCATTTGCGCGGATCCACACCTGACGATGTGGTTCTGCGACGTg    <  1:496739/69‑1 (MQ=255)                                           aaCGCGCAGGTCATCCGGCACGCTACCATTTGCGCGGATCCACACCTGACGATGTGGTTCTGCGACGTg    <  1:1115361/69‑1 (MQ=255)                                           aaCGCGCAGGTCATCCGGCACGCTACCATTTGCGCGGATCCACACCTGACGATGTGGTTCTGCGACGTg    <  1:771355/69‑1 (MQ=255)                                           aaCGCGCAGGTCATCCGGCACGCTACCATTTGCGCGGATCCACACCTGACGATGTGGTTCTGCGACGTg    <  1:483842/69‑1 (MQ=255)                                           aaCGCGCAGGTCATCCGGCACGCTACCATTTGCGCGGATCCACACCTGACGATGTGGTTCTGCGACGTg    <  1:836672/69‑1 (MQ=255)                                           aaCGCGCAGGTCATCCGGCACGCTACCATTTGCGCGGATCCACACCTGACGATGTGGTTCTGCGACGTg    >  1:858129/1‑69 (MQ=255)                                           aaCGCGCAGGTCATCCGGCACGCTACCATTTGCGCGGATCCACACCTGACGATGTGGTTCTGCGACGTg    <  1:2232916/69‑1 (MQ=255)                                           aaCGCGCAGGTCATCCGGCACGCTACCATTTGCGCGGATCCACACCTGACGATGTGGTTCTGCGACGTg    <  1:884585/69‑1 (MQ=255)                                           aaCGCGCAGGTCATCCGGCACGCTACCATTTGCGCGGATCCACACCTGACGATGTGGTTCTGCGACGTg    <  1:1415039/69‑1 (MQ=255)                                           aaCGCGCAGGTCATCCGGCACGCTACCATTTGCGCGGATCCACACCTGACGATGTGGTTCTGCGACGTg    <  1:2293220/69‑1 (MQ=255)                                           aaCGCGCAGGTCATCCGGCACGCTACCATTTGCGCGGATCCACACCTGACGATGTGGTTCTGCGACGTg    >  1:2948861/1‑69 (MQ=255)                                           aaCGCGCAGGTCATCCGGCACGCTACCATTTGCGCGGATCCACACCTGACGATGTGGTTCTGCGACGTg    >  1:1708104/1‑69 (MQ=255)                                           aaCGCGCAGGTCATCCGGCACGCTACCATTTGCGCGGATCCACACCTGACGATGTGGTTCTGCGACGTg    <  1:189847/69‑1 (MQ=255)                                           aaCGCGCAGGTCATCCGGCACGCTACCATTTGCGCGGATCCACACCTGACGATGTGGTTCTGCGACGTg    <  1:2911532/69‑1 (MQ=255)                                           aaCGCGCAGGTCATCCGGCACGCTACCATTTGCGCGGATCCACACCTGACGATGTGGTTCTGCGACGTg    <  1:1470764/69‑1 (MQ=255)                                           aaCGCGCAGGTCATCCGGCACGCTACCATTTGCGCGGATCCACACCTGACGATGTGGTTCTGCGACGTg    <  1:2737515/69‑1 (MQ=255)                                           aaCGCGCAGGTCATCCGGCACGCTACCATTTGCGCGGATCCACACCTGACGATGTGGTTCTGCGACGTg    <  1:2213932/69‑1 (MQ=255)                                           aaCGCGCAGATCATCCGGCACGCTACCATTTGCGCGGATCCACACCTGACGATGCGGTTCTGCGACGt     <  1:3316344/68‑1 (MQ=255)                                           aaCGCGCAGATCATCCGGCACGCTACCATTTGCGCGGATCCACACCTGACGATGCGGTTCTGCGACGt     <  1:806696/68‑1 (MQ=255)                                             cgcgCAGGTCATCCGGCACGCTACCATTTGCGCGGATCCACACCTGACGATGTGGTTCTGCGACGTGAc  <  1:1882526/69‑1 (MQ=255)                                             cgcgCAGGTCATCCGGCACGCTACCATTTGCGCGGATCCACACCTGACGATGTGGTTCTGCGACGTGAc  <  1:1624051/69‑1 (MQ=255)                                             cgcgCAGGTCATCCGGCACGCTACCATTTGCGCGGATCCACACCTGACGATGTGGTTCTGCGACGTGAc  <  1:903075/69‑1 (MQ=255)                                             cgcgCAGGTCATCCGGCACGCTACCATTTGCGCGGATCCACACCTGACGATGTGGTTCTGCGACGTGAc  <  1:1735277/69‑1 (MQ=255)                                             cgcgCAGGTCATCCGGCACGCTACCATTTGCGCGGATCCACACCTGACGATGTGGTTCTGCGACGTGAc  <  1:2998196/69‑1 (MQ=255)                                             cgcgCAGATCATCCGGCACGCTACCATTTGCGCGGATCCACACCTGACGATGCGGTTCTGCGACGTgcc  <  1:2212479/69‑3 (MQ=255)                                             cgcgCAGATCATCCGGCACGCTACCATTTGCGCGGATCCACACCTGACGATGCGGTTCTGCGACGTgcc  <  1:2801069/69‑3 (MQ=255)                                             cgcgCAGATCATCCGGCACGCTACCATTTGCGCGGATCCACACCTGACGATGCGGTTCTGCGACGTgcc  <  1:1361838/69‑3 (MQ=255)                                                    |                                                              CAGGAAGTTAAGATCAGAAGCGTAACCGAGCAGATACTGATGAACGCGCAGGTCATCCGGCACGCTACCATTTGCGCGGATCCACACCTGACGATGTGGTTCTGCGACGTGAC  >  minE/303007‑303119

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