Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	2,093,772	G→T	12.7%	G360G (GGC→GGA)	tdcG ←	L‑serine dehydratase 3

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	2,093,772	0	G	T	12.7%	142.8 / 5.8	47	G360G (GGC→GGA)	tdcG	L‑serine dehydratase 3
Reads supporting (aligned to +/- strand):  ref base G (25/16);  new base T (2/4);  total (27/20)
Fisher's exact test for biased strand distribution p-value = 3.79e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.94e-01

CGGATCGCAGGTCAGCCCAAGGTTATGCTCCATCGCGATTTCCGCCGCATTGCATACCTGCGCCGGACTGCCGCCCAGTAGTTCAGTTAACCCTGCCGCCGCCATTGAACAGGCCACGCCAATCTCCCCCTGACAG  >  minE/2093703‑2093838                                                                      |                                                                   cGGATCGCAGGTCAGCCCAAGGTTATGCTCCATCGCGATTTCCGCCGCATTGCATACCTGCGCCGGACTgc                                                                   >  1:613875/1‑71 (MQ=255)       gCAGGTCAGCCCAAGGTTATGCTCCATCGCGATTTCCGCCGCATTGCATACCTGCGCCGGACTGCCGccc                                                              <  1:838594/70‑1 (MQ=255)        cAGGTCAGCCCAAGGTTATGCTCCATCGCGATTTCCGCCGCATTGCATACCTGCGCCGGACTGCCGCCCag                                                            >  1:2013484/1‑71 (MQ=255)            tCAGCCCAAGGTTATGCTCCATCGCGATTTCCGCCGCATTGCATACCTGCGCCGGACTGCCGCCCagtagt                                                        <  1:103575/71‑1 (MQ=255)              aGCCCAAGGTTATGCTCCATCGCGATTTCCGCCGCATTTCATACCTGCGCCGGACTGCCGCCCAGTAGTTc                                                      <  1:1713055/71‑1 (MQ=255)              aGCCCAAGGTTATGCTCCATCGCGATTTCCGCCGCATTGCATACCTGCGCCGGACTGCCGCCCAGTAGtt                                                       >  1:328614/1‑70 (MQ=255)               gCCCAAGGTTATGCTCCATCGCGATTTCCGCCGCATTGCATACCTGCGCCGGACTGCCGCCCAGTAGTTCa                                                     <  1:1277156/71‑1 (MQ=255)                  cAAGGTTATGCTCCATCGCGATTTCCGCCGCATTGCATACCTGCGCCGGACTGCCGCCCAGTAGTTCa                                                     <  1:3106494/68‑1 (MQ=255)                  cAAGGTTATGCTCCATCGCGATTTCCGCCGCATTGCATACCTGCGCCGGACTGCCGCCCAGTAGTTCAGtt                                                  <  1:2743762/71‑1 (MQ=255)                  cAAGGTTATGCTCCATCGCGATTTCCGCCGCATTGCATACCTGCGCCGGACTGCCGCCCAGTAGTTCAGt                                                   >  1:1230071/1‑70 (MQ=255)                   aaGGTTATGCTCCATCGCGATTTCCGCCGCATTGCATACCTGCGCCGGACTGCCGCCCagt                                                           <  1:3016703/61‑1 (MQ=255)                    aGGTTATGCTCCATCGCGATTTCCGCCGCATTGCATACCTGCGCCGGACTGCCGCCCAGTAGTTCAGTTaa                                                >  1:143060/1‑71 (MQ=255)                    aGGTTATGCTCCATCGCGATTTCCGCCGCATTGCATACCTGCGCCGGACTGCCGCCCAGTAGTTCAGTTa                                                 >  1:76089/1‑70 (MQ=255)                     ggTTATGCTCCATCGCGATTTCCGCCGCATTGCATACCTGCGCCGGACTGCCGCCCAGTAGTTCAGTTAAc                                               <  1:265030/71‑1 (MQ=255)                      gTTGTGCTCCATGGCGATTTCCGCCGCGATGCACACCTGCGCCGGGCTTGCGCCTAATAGtt                                                       <  1:3850049/62‑1 (MQ=255)                      gTTATGCTCCATCGCGATTTCCGCCGCATTGCATACCTGCGCCGGACTg                                                                    >  1:1920682/1‑49 (MQ=255)                       ttATGCTCCATCGCGATTTCCGCCGCATTGCATACCTGCGCCGGACTGCCGCCCAGTAGTTCAGTTAAccc                                             <  1:306517/71‑1 (MQ=255)                         aTGCTCCATCGCGATTTCCGCCGCATTGCATACCTGCGCCGGACTGCCGCCCAGTAGTTCAGTTAACCCt                                            <  1:2570331/70‑1 (MQ=255)                          tGCTCCATCGCGATTTCCGCCGCATTGCATACCTGCGCCGGACTGCCGCCCAGTAGTTCAGTTAACCCTg                                           <  1:995383/70‑1 (MQ=255)                           gCTCCATCGCGATTTCCGCCGCATTGCATACCTGCGCCGGACTGCCGCCCAGTAGTTCAg                                                    >  1:2008961/1‑60 (MQ=255)                             tCCATCGCGATTTCCGCCGCATTGCATACCTGCGCCGGACTGCCGCCCAGTAGTTCAGTTAACCCTgccgc                                       >  1:2338377/1‑71 (MQ=255)                              ccATCGCGATTTCCGCCGCATTGCATACCTGCGCCGGACTGCCGCCCAGTAGTTCAGTTAACCCTgccgcc                                      >  1:3727431/1‑71 (MQ=255)                              ccATCGCGATTTCCGCCGCATTGCATACCTGCGCCGGACTGCCGCCCAGTAGTTCAGTTAACCCTgccgcc                                      >  1:902015/1‑71 (MQ=255)                               cATCGCGATTTCCGCCGCATTGCATACCTGCGCCGGACTGCCGCCCAGTAGTTCAGTTAACCCTgccgccg                                     >  1:115773/1‑71 (MQ=255)                                aTCGCGATTTCCGCCGCATTGCATACCTGCGCCGGACTGCCGCCCAGTAGTTCAGTTAACCCTgcc                                         >  1:1887974/1‑66 (MQ=255)                                 tCGCGATTTCCGCCGCATTGCATACCTGCGCCGGACTGCCGCCCAGTAGTTCAGTTAACCCTgccgccgcc                                   <  1:2749942/71‑1 (MQ=255)                                   gcgATTTCCGCCGCATTGCATACCTGCGCCGGACTGCCGCCCAGTAGTTCAGTTAACCCTGCCGCCGCCa                                  <  1:3865026/70‑1 (MQ=255)                                   gcgATTTCCGCCGCATTGCATACCTGCGCCGGACTGCCGCCCAGTAGTTCAGTTAACCCTGCCGCCGCCAt                                 >  1:2618907/1‑71 (MQ=255)                                   gcgATTTCCGCCGCATTGCATACCTGCGCCGGACTGCCGCCCAGTAGTTCAGTTAACCCTGCCGCCGCCAt                                 >  1:790937/1‑71 (MQ=255)                                   gcgATTTCCGCCGCATTGCATACCTGCGCCGGACTGCCGCCCAGTAGTTCAGTTAACCCTGCCGCCGCCAt                                 >  1:1352718/1‑71 (MQ=255)                                    cgATTTCCGCCGCGATGCACACCTGCGCCGGGCTTGCGCCTAATAGTTCTGCCAGACCAGCCGCCGCCAtt                                <  1:31148/71‑1 (MQ=255)                                     gATTTCCGCCGCATTGCATACCTGCGCCGGACTGCCGCCCAGTAGtt                                                       >  1:869469/1‑47 (MQ=255)                                      aTTTCCGCCGCATTGCATACCTGCGCCGGACTGCCGCCCAGTAGTTCAGTTAACCCTGCCGCCGCCATTg                               >  1:1832840/1‑70 (MQ=255)                                       tttCCGCCGCGATGCACACCTGCGCCGGGCTTGCGCCTAATAGTTCTGCCAGACCAGCCGCCGCCAt                                 <  1:546346/67‑1 (MQ=255)                                        ttCCGCCGCGATGCACACCTGCGCCGGGCTTGCGCCTAATAGTTCTGCCAGACCAGCCGCCGCCATTGag                             >  1:3245122/1‑69 (MQ=255)                                        ttCCGCCGCATTGCATACCTGCGCCGGACTGCCGCCCAGTAGTTCAGtt                                                  <  1:3523834/49‑1 (MQ=255)                                        ttCCGCCGCATTGCATACCTGCGCCGGACTGCCGCCCAGTAGTTCAGTTaa                                                >  1:3488973/1‑51 (MQ=255)                                         tCCGCCGCATTGCATACCTGCGCCGGACTGCCGCCCAGTAGTTCAGTTaa                                                >  1:319873/1‑50 (MQ=255)                                         tCCGCCGCATTGCATACCTGCGCCGGACTGCCGCCCAGTAGTTCAGTTAACCCTGCCGCCGCCATTGaa                             >  1:3228860/1‑69 (MQ=255)                                         tCCGCCGCATTGCATACCTGCGCCGGACTGCCGCCCAGTAGTTCAGTTAACCCTGCCGCCGCCATTGaa                             >  1:2260235/1‑69 (MQ=255)                                             ccgcATTGCATACCTGCGCCGGACTGCCGCCCAGTAGtt                                                       >  1:2083737/1‑39 (MQ=255)                                             ccgcATTGCATACCTGCGCCGGACTGCCGCCCAGTAGtt                                                       >  1:2461100/1‑39 (MQ=255)                                              cgcATTGCATACCTGCGCCGGACTTCCGCCCAGTAGTTCAGTTAACCCTGCCGCCGCCATTGAACAGGCCa                      >  1:1650097/1‑71 (MQ=255)                                                  ttGCATACCTGCGCCGGACTGCCGCCCAGTAGTTCAGTTAACCCTGCCGCCGCCATTGAACAGGCCACGc                   >  1:3167169/1‑70 (MQ=255)                                                  atGCACACCTGCGCCGGGCTTGCGCCTAATAGTTCTGCCAGACCAGCCGCCGCCATTGAGCacgccacgcc                  <  1:788451/70‑1 (MQ=255)                                                    gCATACCTGCGCCGGACTGCCGCCCAGTAGTTCAGTTAACCCTg                                           >  1:1979594/1‑44 (MQ=255)                                                    gCATACCTGCGCCGGACTGCCGCCCAGTAGTTCAGTTAACCCTGCCGCCGCCATTGAACAGGCCACGCCaa                >  1:3864670/1‑71 (MQ=255)                                                             cgcCGGACTGCCGCCCAGTATTTCAGTTAACCCTGCCGCCGCCATTGAACAGGCCACGCCAATCTCCCCCt       <  1:1525452/71‑1 (MQ=255)                                                                  gACTGCCGCCCAGTAGTTCAGTTAACCCTGCCGCCGCCATTGAACAGGCCACGCCAATCTCCCCCTGACAg  <  1:1332275/71‑1 (MQ=255)                                                                      |                                                                   CGGATCGCAGGTCAGCCCAAGGTTATGCTCCATCGCGATTTCCGCCGCATTGCATACCTGCGCCGGACTGCCGCCCAGTAGTTCAGTTAACCCTGCCGCCGCCATTGAACAGGCCACGCCAATCTCCCCCTGACAG  >  minE/2093703‑2093838

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑