Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	2,094,006	T→C	14.3%	V282V (GTA→GTG)	tdcG ←	L‑serine dehydratase 3

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	2,094,006	0	T	C	14.3%	144.0 / 14.8	56	V282V (GTA→GTG)	tdcG	L‑serine dehydratase 3
Reads supporting (aligned to +/- strand):  ref base T (26/22);  new base C (3/5);  total (29/27)
Fisher's exact test for biased strand distribution p-value = 4.62e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CGACGGAACTTATCGTAATAAGCCAGTACTGCCGGAATAATGCCGCACGCACCGTTAGTCGGTGCCGTTACCACGCGCCCGCCAGCTGCGTTTTCTTCACTAACCGCCAGCGCGTACATGTTGATCCAGTCGATG  >  minE/2093938‑2094072                                                                     |                                                                   cgacgGAACTTATCGTAATAAGCCAGTACTGCCGGAATAATGCCGCACGCACCGTTAGTCGGTGCCGTTAc                                                                  >  1:1670121/1‑71 (MQ=255)    cgGAACTTATCGTAATAAGCCAGTACTGCCGGAATAATGCCGCACGCACCGTTAGTCGGTGCCGTTACCAc                                                               >  1:4040782/1‑71 (MQ=255)     taaaCTTGTCGTAGTACGCCAGAACTGCCGGGATAATCCCGCACGCACCGTTAGTCGGCGCAGTcaccac                                                               <  1:2549982/68‑1 (MQ=255)         cTTATCGTAATAAGCCAGTACTGCCGGAATAATGCCGCACGCACCGTTAGTCGGTGCCGTTACCACGcgcc                                                          >  1:3169939/1‑71 (MQ=255)            aTCGTAATAAGCCAGTACTGCCGGAATAATGCCGCACGCACCGTTAGTCGGTGCCGTTACCACGcgcccgc                                                       >  1:3507484/1‑71 (MQ=255)            aTCGTAATAAGCCAGTACTGCCGGAATAATGCCGCACGCACCGTTAGTCGGTGCCGTTACCACGcgcccgc                                                       >  1:1998707/1‑71 (MQ=255)                 aataaGCCAGTACTGCCGGAATAATGCCGCACGCACCGTTAGTCGGTGCCGTTACCACGCGCCCGCCa                                                     >  1:30859/1‑68 (MQ=255)                 aataaGCCAGTACTGCCGGAATAATGCCGCACGCACCGTTAGTCGGTGCCGTTACCACGCGCCCGCCAGCt                                                  >  1:2630179/1‑71 (MQ=255)                  ataaGCCAGTACTGCCGGAATAATGCCGCACGCACCGTTAGTCGGTGCCGTTACCACGCGCCCGCCAGCt                                                  >  1:868048/1‑70 (MQ=255)                  ataaGCCAGTACTGCCGGAATAATGCCGCACGCACCGTTAGTCGGTGCCGTTACCACGCGCCCGCCAGCTg                                                 >  1:1189465/1‑71 (MQ=255)                         aGTACTGCCGGAATAATGCCGCACGCACCGTTAGTCGGTGCCGTTACCACGCGCCCGCCAGCTGCGTTttc                                          <  1:262743/71‑1 (MQ=255)                           tACTGCCGGAATAATGCCGCACGCACCGTTAGTCGGTGCCGTTACCACGCGCCCGCCAGCTTCGTTttctt                                        <  1:682388/71‑1 (MQ=255)                           tACTGCCGGAATAATGCCGCACGCACCGTTAGTCGGTGCCGTTACCACGCGCCCGCCAGCTGCg                                               <  1:4121814/64‑1 (MQ=255)                           tACTGCCGGAATAATGCCGCACGCACCGTTAGTCGGTGCCGTTACCACGCGCCCGCCAGCTGCGTTttctt                                        <  1:4096712/71‑1 (MQ=255)                           tACTGCCGGAATAATGCCGCACGCACCGTTAGTCGGTGCCGTTACCACGCGCCCGCCAGCTGCGTTttctt                                        <  1:291785/71‑1 (MQ=255)                            aCTGCCGGAATAATGCCGCACGCACCGTTAGTCGGTGCCGTTACCACGCGCCCGCCAGCTGCGTTttctt                                        >  1:224093/1‑70 (MQ=255)                             cTGCCGGGATAATCCCGCACGCACCGTTAGTCGGCGCAGTCACCACGCGACCGCCAGCAGCGTTctcttcg                                      <  1:2716984/71‑2 (MQ=255)                             cTGCCGGGATAATCCCGCACGCACCGTTAGTCGGCGCAGTCACCACGCGACCGCCAGCAGCGTTctcttcg                                      <  1:327164/71‑2 (MQ=255)                             cTGCCGGAATAATGCCGCACGCACCGTTAGTCGGTGCCGTTACCACGCGCCCGCCAGCTGCGTTTTCTTCa                                      >  1:2928687/1‑71 (MQ=255)                             cTGCCGGAATAATGCCGCACGCACCGTTAGTCGGTGCCGTTACAACGCGCCCGCCAGCTGCGTTTTCTTCa                                      >  1:1124931/1‑71 (MQ=255)                               gCCGGAATAATGCCGCACGCACCGTTAGTCGGTGCCGTTACCACGCGCCCGCCAGCTGCGTTTTCTTCAc                                     >  1:361529/1‑70 (MQ=255)                               gCCGGAATAATGCCGCACGCACCGTTAGTCGGTGCCGTTACCACGCGCCCGCCAGCTGCGTTTTCTTCAc                                     >  1:3214810/1‑70 (MQ=255)                                ccGGAATAATGCCTCACGCACCGTTAGTCGGTGCCGTTACCa                                                                >  1:3215485/1‑42 (MQ=255)                                ccGGAATAATGCCGCACGCACCGTTAGTCGGTGCCGTTACCACGCGCCCGCCAGCTGCGTTTTCTTCACt                                    >  1:1250371/1‑70 (MQ=255)                                ccGGAATAATGCCGCACGCACCGTTAGTCGGTGCCGTTACCACGCGCCCGCCAGCTGCGTTTTCTTCACTa                                   <  1:139694/71‑1 (MQ=255)                                 cGGAATAATGCCGCACGCACCGTTAGTCGGTGCCGTTACCACGCGCCCGCCAGCTGCGTTTTCTTCACTaa                                  <  1:2098593/71‑1 (MQ=255)                                     ataatGCCGCACGCACCGTTAGTCGGTGCCGTTACCACGCGCCCGCCAGCTGCGTTTTCTTCACTAACCGc                              >  1:3190474/1‑71 (MQ=255)                                     ataatGCCGCACGCACCGTTAGTCGGTGCCGTTACCACGCGCCCGCCAGCTGCGTTTTCTTCACTAACCGc                              >  1:3465786/1‑71 (MQ=255)                                         tGCCGCACGCACCGTTAGTCGGTGCCGTTACCACGCGCCCGCCAGCTGCGTTTTCTTCACTAAcc                                >  1:3695737/1‑65 (MQ=255)                                         tGCCGCACGCACCGTTAGTCGGTGCCGTTACCACGCGCCCGCCAGCTGCGTTTTCTTCACTAACCGTCAgc                          <  1:424534/71‑1 (MQ=255)                                           ccGCACGCACCGTTAGTCGGCGCAGTCACCACGCGACCGCCAGCAGCGTTCTCTTCGTTCACTGCCAgtg                         >  1:915853/1‑68 (MQ=255)                                           ccGCACGCACCGTTAGTCGGCGCAGTCACCACGCGACCGCCAGCAGCGTTCTCTTCGTTCACTGCCAgtg                         >  1:2807830/1‑68 (MQ=255)                                            cgcacgcaCCGTTAGTCGGCGCAGTCACCACGCGACCGCCAGCAGCGTTCTCTTCGTTCACTGCCAgtgca                       <  1:3406803/71‑2 (MQ=255)                                            cgcacgcaCCGTTAGTCGGCGCAGTCACCACGCGACCGCCAGCAGCGTTCTCTTCGTTCACTGCCAgtgca                       <  1:1148017/71‑2 (MQ=255)                                            cgcacgcaCCGTTAGTCGGCGCAGTCACCACGCGACCGCCAGCAGCGTTCTCTTCGTTCACTGCCAgtgc                        >  1:3438212/1‑70 (MQ=255)                                             gcacgcacCGTTAGTCGGTGCCGTTACCACGCGCCCGCCAGCTGCGTTTTCTTCACTAACCGCCAGcgcg                       >  1:2686903/1‑70 (MQ=255)                                                cgcacCGTTAGTCGGTGCCGTTACCACGCGCCCGCCAGCTGCGTTTTCTTCACTAACCGCCAGCGCGTACa                   >  1:3285942/1‑71 (MQ=255)                                                 gcacCGTTAGTCGGCGCAGGCACCACGCGACCGCCAGCAGCGTTctct                                         <  1:3950309/48‑1 (MQ=16)                                                  cacCGTTAGTCGGTGCCGTTACCACGCGCCCGCCAGCTGCGTTTTCTTCACTAACCGCCAGCGCGTACATg                 <  1:520067/71‑1 (MQ=255)                                                  cacCGTTAGTCGGTGCCGTTACCACGCGCCCGCCAGCTGCGTTTTCTTCACTAACCGCCAGCGCGTACATg                 <  1:2201729/71‑1 (MQ=255)                                                   acCGTTAGTCGGTGCCGTTACCACGCGCCCGCCAGCTGCGTTTTCTTCACTAACCGCCAGCGCGTACATGt                <  1:3827713/71‑1 (MQ=255)                                                    ccGTTAGTCGGTGCCGTTACCACGCGCCCGCCAGCTGCGTTttct                                         <  1:4082012/45‑1 (MQ=255)                                                       ttAGTCGGTGCCGTTACCACGCGCCCGCCAGCTGCGTTTTCTTCACTAACCGCCAGCGCGTACATGTTGAt            >  1:1401770/1‑71 (MQ=255)                                                       ttAGTCGGTGCCGTTACCACGCGCCCGCCAGCTGCGTTTTCTTCACTAACCGCCAGCGCGTACATGTTGAt            <  1:2463316/71‑1 (MQ=255)                                                        tAGTCGGTGCCGTTACCACGCGCCCGCCAGCTGCGTTTTCTTCACTAACCGCCAGCGCGTACATGTTGAt            >  1:3042194/1‑70 (MQ=255)                                                         aGTCGGTGCCGTTACCACGCGCCCGCCAGCTGCGTTTTCTTCACTa                                   <  1:2193282/46‑1 (MQ=255)                                                          gTCGGTGCCGTTACCACGCGCCCGCCAGCTGCGTTTTCTTCACTAACCGCCAGCGCGTACATGTTGATcc          <  1:2943441/70‑1 (MQ=255)                                                          gTCGGTGCCGTTACCACGCGCCCGCCAGCTGCGTTTTCTTCACTAACCGCCAGCGCGTACATGTTGATCCa         <  1:3478794/71‑1 (MQ=255)                                                           tCGGTGCCGTTACCACGCGCCCGCCAGCTGCGTTTTCTTCACTAACCGcc                             >  1:4110362/1‑50 (MQ=255)                                                           tCGGTGCCGTTACCACGCGCCCGCCAGCTGCGTTTTCTTCACTAACCGCCAGCGCGTACATGTTGATCCAg        <  1:2559309/71‑1 (MQ=255)                                                              gTGCCGTTACCACGCGCCCGCCAGCTGCGTTTTCTTCACTAACCGCCAGCGCGTACATGTTGAt            <  1:2510980/64‑1 (MQ=255)                                                              gTGCCGTTACCACGCGCCCGCCAGCTGCGTTTTCTTCACTAACCGCCAGCGCGTACATGTTGATCCAGTCg     <  1:3659028/71‑1 (MQ=255)                                                               tGCCGTTACCACGCGCCCGCCAGCTGCGTTTTCTTCACTAACCGCCAGCGCGTACATGTTGATCCAGTCGa    >  1:1423232/1‑71 (MQ=255)                                                                 ccGTTACCACGCGCCCGCCAGCTGCGTTTTCTTCACTAACCGCCAGCGCGTACATGTTGATCCAGTCGAt   <  1:1911099/70‑1 (MQ=255)                                                                 ccGTTACCACGCGCCCGCCAGCTGCGTTTTCTTCACTAACCGCCAGCGCGTACATGTTGATCCAGTCGATg  <  1:565092/71‑1 (MQ=255)                                                                  cGTTACCACGCGCCCGCCAGCTGCGTTTTCTTCACTAACCGCCAGCGCGTACATGt                >  1:3777975/1‑56 (MQ=255)                                                                   gTTACCACGCGCCCGCCAGCTGCGTTTTCTTCACTAACCGCCAGCGCGTACATGtt               <  1:170599/56‑1 (MQ=255)                                                                    ttACCACGCGCCCGCCAGCTGCGTTTTCTTCACTAACCGcc                             >  1:3771245/1‑41 (MQ=255)                                                                     |                                                                   CGACGGAACTTATCGTAATAAGCCAGTACTGCCGGAATAATGCCGCACGCACCGTTAGTCGGTGCCGTTACCACGCGCCCGCCAGCTGCGTTTTCTTCACTAACCGCCAGCGCGTACATGTTGATCCAGTCGATG  >  minE/2093938‑2094072

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