Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	2,093,780	G→A	8.1%	L358L (CTG→TTG) ‡	tdcG ←	L‑serine dehydratase 3

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	2,093,780	0	G	A	8.1%	137.6 / 2.4	49	L358L (CTG→TTG) ‡	tdcG	L‑serine dehydratase 3
Reads supporting (aligned to +/- strand):  ref base G (26/19);  new base A (2/2);  total (28/21)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.78e-01

CAGGTCAGCCCAAGGTTATGCTCCATCGCGATTTCCGCCGCATTGCATACCTGCGCCGGACTGCCGCCCAGTAGTTCAGTTAACCCTGCCGCCGCCATTGAACAGGCCACGCCAATCTCCCCCTGACAGCCGACTTCC  >  minE/2093710‑2093847                                                                       |                                                                    cAGGTCAGCCCAAGGTTATGCTCCATCGCGATTTCCGCCGCATTGCATACCTGCGCCGGACTGCCGCCCag                                                                     <  1:3440753/71‑1 (MQ=255)     tCAGCCCAAGGTTATGCTCCATCGCGATTTCCGCCGCATTGCATACCTGCGCCGGACTGCCGCCCagtagt                                                                 <  1:3609381/71‑1 (MQ=255)       aGCCCAAGGTTATGCTCCATCGCGATTTCCGCCGCATTGCATACCTGCGCCGGACTGCCGCCCAGTAGtt                                                                <  1:577449/70‑1 (MQ=255)           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:80479/1‑66 (MQ=255)              ggTTATGCTCCATCGCGATTTCCGCCGCATTGCATACCTGCGCCGGACTGCCGCCCAGTAGTTCAg                                                             >  1:1757215/1‑66 (MQ=255)              ggTTATGCTCCATCGCGATTTCCGCCGCATTGCATACCTGCGCCGGACTGCCGCCCAGTAGTTCAg                                                             >  1:2639012/1‑66 (MQ=255)              ggTTATGCTCCATCGCGATTTCCGCCGCATTGCATACCTGCGCCGGACTGCCGCCCAGTAGTTCAg                                                             >  1:2795890/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)                           cgcgATTTCCGCCGCATTGCATACCTGCGCCGGACTGCCGCCCagta                                                                   >  1:4412515/1‑47 (MQ=255)                            gcgATTTCCGCCGCATTGCATACCTGCGCCGGACTGCCGCCCAGTAGTTCAGtt                                                           >  1:682396/1‑54 (MQ=255)                            gcgATTTCCGCCGCATTGCATACCTGCGCCGGACTGCCGCCCAGTAGTTCAGTTAAcc                                                       >  1:1375060/1‑58 (MQ=255)                            gcgATTTCCGCCGCATTGCATACCTGCGCCGGACTGCCGCCCAGTAGTTCAGTTAAcc                                                       >  1:5362224/1‑58 (MQ=255)                            gcgATTTCCGCCGCATTGCATACCTGCGCCGGACTGCCGCCCAGTAGTTCAGTTAAcc                                                       >  1:2489800/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)                               aTTTCCGCCGCATTGCATACCTGCGCCGGACTGCCGCCCAGTAGTTc                                                               <  1:1469940/47‑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:4668797/70‑1 (MQ=255)                                             gCATACCTGCGCCGGACTGCCGCCCAGTAGTTCAGTTAACCCTGCCGCCGCCATTGAACAGGCCACGCCa                          <  1:4730031/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)                                                                       |                                                                    CAGGTCAGCCCAAGGTTATGCTCCATCGCGATTTCCGCCGCATTGCATACCTGCGCCGGACTGCCGCCCAGTAGTTCAGTTAACCCTGCCGCCGCCATTGAACAGGCCACGCCAATCTCCCCCTGACAGCCGACTTCC  >  minE/2093710‑2093847

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