Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	2,093,787	A→T	6.3%	T355T (ACT→ACA) ‡	tdcG ←	L‑serine dehydratase 3

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	2,093,787	0	A	T	6.3%	177.8 / 3.5	63	T355T (ACT→ACA) ‡	tdcG	L‑serine dehydratase 3
Reads supporting (aligned to +/- strand):  ref base A (37/22);  new base T (2/2);  total (39/24)
Fisher's exact test for biased strand distribution p-value = 6.32e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.71e-01

CAAGGTTATGCTCCATCGCGATTTCCGCCGCATTGCATACCTGCGCCGGACTGCCGCCCAGTAGTTCAGTTAACCCTGCCGCCGCCATTGAACAGGCCACGCCAATCTCCCCCTGACAGCCGACTTCCGCGCCAGA  >  minE/2093720‑2093855                                                                    |                                                                     cAAGGTTATGCTCCATCGCGATTTCCGCCGCATTGCATACCTGCGCCGGACTGCCGCCCAGTAGTTCAGt                                                                    >  1:556013/1‑70 (MQ=255)  aaGGTTATGCTCCATCGCGATTTCCGCCGCATTGCATACCTGCGCCGGACTGCCGCCCAGTAGTTCAGTTa                                                                  <  1:1068627/71‑1 (MQ=255)  aaGGTTATGCTCCATCGCGATTTCCGCCGCATTGCATACCTGCGCCGGACTGCCGCCCAGTAGTTCAGTTa                                                                  <  1:1061863/71‑1 (MQ=255)   aGGTTATGCTCCATCGCGATTTCCGCCGCATTGCATACCTGCGCCGGACTGCCGCCCAGTAGTTCAGTTaa                                                                 <  1:3242207/71‑1 (MQ=255)    ggTTATGCTCCATCGCGATTTCCGCCGCATTGCATACCTGCGCCGGACTGCCGCCCAGTAGTTCAg                                                                     >  1:2795890/1‑66 (MQ=255)    ggTTATGCTCCATCGCGATTTCCGCCGCATTGCATACCTGCGCCGGACTGCCGCCCAGTAGTTCAg                                                                     >  1:80479/1‑66 (MQ=255)    ggTTATGCTCCATCGCGATTTCCGCCGCATTGCATACCTGCGCCGGACTGCCGCCCAGTAGTTCAg                                                                     >  1:2639012/1‑66 (MQ=255)    ggTTATGCTCCATCGCGATTTCCGCCGCATTGCATACCTGCGCCGGACTGCCGCCCAGTAGTTCAg                                                                     >  1:1757215/1‑66 (MQ=255)     gTTATGCTCCATCGCGATTTCCGCCGCATTGCATACCTGCGCCGGACTGCCGCCCAGTAGTTCAGTTAAcc                                                               <  1:829348/71‑1 (MQ=255)     gTTATGCTCCATCGCGATTTCCGCCGCATTGCATACCTGCGCCGGACTGCCGCCCAGTAGTTCAGTTAAcc                                                               <  1:4757352/71‑1 (MQ=255)              cATCGCGATTTCCGCCGCATTGCATACCTGCGCCGGACTGCCGCCCAGTAGTTCAGTTAACCCTgccgccg                                                      >  1:5469700/1‑71 (MQ=255)               aTCGCGATTTCCGCCGCATTGCATACCTGCGCCGGACTGCCGCCCAGTAGTTCAGTTAACCCTgccgccgc                                                     >  1:3856781/1‑71 (MQ=255)                tCGCGATTTCCGCCGCATTGCATACCTGCGCCGGACTGCCGCCCAGTAGTTCAGTTAACCCTgccgccgc                                                     >  1:2843138/1‑70 (MQ=255)                  gcgATTTCCGCCGCATTGCATACCTGCGCCGGACTGCCGCCCAGTAGTTCAGtt                                                                   >  1:682396/1‑54 (MQ=255)                  gcgATTTCCGCCGCATTGCATACCTGCGCCGGACTGCCGCCCAGTAGTTCAGTTAAcc                                                               >  1:2489800/1‑58 (MQ=255)                  gcgATTTCCGCCGCATTGCATACCTGCGCCGGACTGCCGCCCAGTAGTTCAGTTAAcc                                                               >  1:5362224/1‑58 (MQ=255)                  gcgATTTCCGCCGCATTGCATACCTGCGCCGGACTGCCGCCCAGTAGTTCAGTTAAcc                                                               >  1:1375060/1‑58 (MQ=255)                  gcgATTTCCGCCGCATTGCATACCTGCGCCGGACTGCCGCCCAGTAGTTCAGTTAACCCTGCCGCCGCCAt                                                  >  1:849638/1‑71 (MQ=255)                    gATTTCCGCCGCATTGCATACCTGCGCCGGACTGCCGCCCAGTAGTTCAGTTAACCCTGCCGCCGCCAtt                                                 >  1:4173286/1‑70 (MQ=255)                    gATTTCCGCCGCATTGCATACCTGCGCCGGACTGCCGCCCAGTAGTTCAGTTAACCCTGCCGCCGCCATTg                                                <  1:4048256/71‑1 (MQ=255)                      tttCCGCCGCATTGCATACCTGCGCCGGACTGCCGCCCAGTAGTTCAGTTAACCCTgccgccgc                                                     <  1:2904729/64‑1 (MQ=255)                         ccgccgCGATGCACACCTGCGCCGGGCTTGCGCCTAATAGTTCTGCCAGACCAGCCGCCGCCATTGagc                                             <  1:4323993/69‑3 (MQ=255)                         ccgccgCATTGCATACCTGCGCCGGACTGCCGCCCAGTAGTTCAGTTAACCCTGCCGCCGCCATTGAACa                                            <  1:161852/70‑1 (MQ=255)                          cgccgcGATGCACACCTGCGCCGGGCTTGCGCCTAATAGTTCTGCCAGACCAGCCGCCGCCATTGAGCacg                                          >  1:4570746/1‑69 (MQ=255)                               cATTGCATACCTGCGCCGGACTGCCGCCCAGTAGTTCAGTTAACCCTGCCGCCGCCa                                                   >  1:2495466/1‑57 (MQ=255)                                 ttGCATACCTGCGCCGGACTGCCGCCCAGTAGTTCAGTTAACCCTGCCGCCGCCATTGAACAGGCCACGc                                    >  1:4963297/1‑70 (MQ=255)                                  tGCACACCTGCGCCGGGCTTGCGCCTAATAGTTCTGCCAGACCAGCCGCCGCCATTGAGCAcgccacgcca                                  <  1:543722/71‑1 (MQ=255)                                   gCATACCTGCGCCGGACTGCCGCCCAGTAGTTCAg                                                                     >  1:2572045/1‑35 (MQ=255)                                   gCATACCTGCGCCGGACTGCCGCCCAGTAGTTCAGTTAACCCTGCCGCCGCCATTGAAc                                             >  1:1229003/1‑59 (MQ=255)                                   gCATACCTGCGCCGGACTGCCGCCCAGTAGTTCAGTTAACCCTGCCGCCGCCATTGAACAGGCCACGCCa                                  <  1:4730031/70‑1 (MQ=255)                                   gCATACCTGCGCCGGACTGCCGCCCAGTAGTTCAGTTAACCCTGCCGCCGCCATTGAACAGGCCACGCCa                                  <  1:4668797/70‑1 (MQ=255)                                       aCCTGCGCCGGGCTTGCGCCTAATAGTTCTGCCAGACCAGCCGCCGCCATTGAGCACGCCACGCCaactt                              >  1:3027797/1‑67 (MQ=255)                                         cTGCGCCGGACTGCCGCCCAGTAGTTCAGTTAACCCTGCCGCCGCCATTGAACAGGCCACGCCAATCTccc                           >  1:4067880/1‑71 (MQ=255)                                            cgcCGGACTGCCGCCCAGTAGTTCAGTTAACCCTGCCGCCGCCATTGAACAGGCCACGCCAAtc                               >  1:5441428/1‑64 (MQ=255)                                            cgcCGGACTGCCGCCCAGTAGTTCAGTTAACCCTGCCGCCGCCATTGAACAGGCCACGCCAATCTCCCCCt                        <  1:4555033/71‑1 (MQ=255)                                            cgcCGGACTGCCGCCCAGTAGTTCAGTTAACCCTGCCGCCGCCATTGAACAGGCCACGCCAATCTCCCCCt                        <  1:756842/71‑1 (MQ=255)                                              ccGGACTGCCGCCCAGTAGTTCAGTTAACCCTGCCGCCGCCATTGAACAGGCCACGCCAATCTCCCCCTGa                      <  1:344475/71‑1 (MQ=255)                                              ccGGACTGCCGCCCAGTAGTTCAGTTAACCCTGCCGCCGCCATTGAACAGGCCACGCCAATCTCCCCCTGa                      <  1:663082/71‑1 (MQ=255)                                                ggACTGCCGCCCAGTAGTTCAGTTAACCCTGCCGCCGCCATTGaa                                              >  1:3686379/1‑45 (MQ=255)                                                ggACTGCCGCCCAGTAGTTCAGTTAACCCTGCCGCCGCCATTGAACAGGCCACGCCAATCTCCCCCTGACa                    <  1:1509830/71‑1 (MQ=255)                                                 gACTGCCGCCCAGTAGTTCAGTTAACCCTGCCGCCGCCATTGAACAGGCCACGCCAATCTCCCCCTGACAg                   >  1:258567/1‑71 (MQ=255)                                                     gccgccCAGTAGTTCAGTTAACCCTGCCGCCGCCATTGAACAGGCCACGCCAATCTCCCCCTGACAGCCg                >  1:2350007/1‑70 (MQ=255)                                                      ccgccCAGTAGTTCAGTTAACCCTGCCGCCGCCATTGAACAGGCCACGCCAATCTCCCCCTGACAGCCGa               <  1:384699/70‑1 (MQ=255)                                                      ccgccCAGTAGTTCAGTTAACCCTGCCGCCGCCATTGAACAGGCCACGCCAATCTCCCCCTGACAGCCGa               >  1:2655589/1‑70 (MQ=255)                                                          ccAGTAGTTCAGTTAACCCTGCCGCCGCCATTGAACAGGCCACGCCAATCTCCCCCTGACAGCCGACTTcc          <  1:995537/71‑1 (MQ=255)                                                           cAGTAGTTCAGTTAACCCTGCCGCCGCCATTGAACAGGCCACGc                                    >  1:3114831/1‑44 (MQ=255)                                                           cAGTAGTTCAGTTAACCCTGCCGCCGCCATTGAACAGGCCACGc                                    >  1:552915/1‑44 (MQ=255)                                                               agtTCAGTTAACCCTGCCGCCGCCATTGAACAGGCCACGCCa                                  <  1:4597969/42‑1 (MQ=255)                                                                gtTCAGTTAACCCTGCCGCCGCCATTGAACAGGCCAc                                      >  1:1670102/1‑37 (MQ=255)                                                                gtTCAGTTAACCCTGCCGCCGCCATTGAACAGGCCAc                                      >  1:998413/1‑37 (MQ=255)                                                                gtTCAGTTAACCCTGCCGCCGCCATTGAACAGGCCAc                                      >  1:3106350/1‑37 (MQ=255)                                                                gtTCAGTTAACCCTGCCGCCGCCATTGAACAGGCCACGCCAAt                                >  1:5150474/1‑43 (MQ=255)                                                                gtTCAGTTAACCCTGCCGCCGCCATTGAACAGGCCACGCCAAt                                >  1:756109/1‑43 (MQ=255)                                                                gtTCAGTTAACCCTGCCGCCGCCATTGAACAGGCCACGCCAATCTCCCCCTGACAGCCGACTTCCgcgc      >  1:3478369/1‑69 (MQ=255)                                                                gtTCAGTTAACACTGCCGCCGCCATTGAACAGGCCAc                                      >  1:5085282/1‑37 (MQ=255)                                                                 ttCAGTTAACCCTGCCGCCGCCATTGAACAGGCCACGCCAATCTCCCCCTGACAGCCGa               <  1:2961925/59‑1 (MQ=255)                                                                 ttCAGTTAACCCTGCCGCCGCCATTGAACAGGCCACGCCAATCTCCCCCTGACAGCCGACTTc           <  1:1246004/63‑1 (MQ=255)                                                                  tCAGTTAACCCTGCCGCCGCCATTGAACAGGCCACGCCAATCTCCCCCTGACa                    <  1:4327679/53‑1 (MQ=255)                                                                   cAGTTAACCCTGCCGCCGCCATTGAACAGGCCACGCCAATCTCCCCCTGACAGCCGACTTCCGCGCCaga  <  1:155535/70‑1 (MQ=255)                                                                   cAGTTAACCCTGCCGCCGCCATTGAACAGGCCACGCCAATCTCCCCCTGACAGCCGACTTCCGCGCCa    >  1:5442082/1‑68 (MQ=255)                                                                   cAGTTAACCCTGCCGCCGCCATTGAACAGGCCACGCCAATCTCCCCCTGACAGCCGACTTCCGCGCCa    >  1:1812151/1‑68 (MQ=255)                                                                   cAGTTAACCCTGCCGCCGCCATTGAACAGGCCACGCCAATCTCCCCCTGACAGCCGACTTCCGCGCCa    >  1:4780405/1‑68 (MQ=255)                                                                   cAGTTAACCCTGCCGCCGCCATTGAACAGGCCACGCCAATCTCCCCCTGACAGCCGACTTCCGCGCCa    >  1:4659274/1‑68 (MQ=255)                                                                    |                                                                     CAAGGTTATGCTCCATCGCGATTTCCGCCGCATTGCATACCTGCGCCGGACTGCCGCCCAGTAGTTCAGTTAACCCTGCCGCCGCCATTGAACAGGCCACGCCAATCTCCCCCTGACAGCCGACTTCCGCGCCAGA  >  minE/2093720‑2093855

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