Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	2,093,752	T→A	16.3%	N367I (AAT→ATT) ‡	tdcG ←	L‑serine dehydratase 3

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	2,093,752	0	T	A	16.3%	164.6 / 22.9	67	N367I (AAT→ATT) ‡	tdcG	L‑serine dehydratase 3
Reads supporting (aligned to +/- strand):  ref base T (31/25);  new base A (2/9);  total (33/34)
Fisher's exact test for biased strand distribution p-value = 4.46e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.49e-01

ATTTGTACCTGTCCGGCAACCGGATCGCAGGTCAGCCCAAGGTTATGCTCCATCGCGATTTCCGCCGCATTGCATACCTGCGCCGGACTGCCGCCCAGTAGTTCAGTTAACCCTGCCGCCGCCATTGAACAGGCCACGC  >  minE/2093683‑2093821                                                                      |                                                                      aTTTGTACCTGTCCGGCAACCGGATCGCAGGTCAGCCCAAGGTTATGCTCCATCGCGATTTCCGCCGCAtt                                                                      <  1:755610/71‑1 (MQ=255) aTTTGTACCTGTCCGGCAACCGGATCGCAGGTCAGCCCAAGGTTATGCTCCATCGCGATTTCCGCCGCAtt                                                                      >  1:1312345/1‑71 (MQ=255) aTTTGTACCTGTCCGGCAACCGGATCGCAGGTCAGCCCAAGGTTATGCTCCATCGCGATTTACGCCGCAtt                                                                      <  1:3689283/71‑1 (MQ=255)   ttGTACCTGTCCGGCAACCGGATCGCAGGTCAGCCCAAGGTTATGCTCCATCGCGATTTCCGCCGCATTg                                                                     <  1:2937035/70‑1 (MQ=255)    tGTACCTGTCCGGCAACCGGATCGCAGGTCAGCCCAAGGTTATGCTCCATCGCGATTTCCGCCGCATTg                                                                     <  1:2619489/69‑1 (MQ=255)      tACCTGTCCGGCGACCGGGTCACACGTCAGACCGAGGTTGTGCTCCATGGCGATTTCCGCCGCGATGCac                                                                  <  1:2748418/70‑2 (MQ=255)      tACCTGTCCGGCAACCGGATCGCAGGTCAGCCCAAGGTTATGCTCCATCGCGATTTCCGCCGCATTGCAt                                                                  <  1:660514/70‑1 (MQ=255)      tACCTGTCCGGCAACCGGATCGCAGGTCAGCCCAAGGTTATGCTCCATCGCGATTTCCGCCGCATTGCAt                                                                  <  1:665244/70‑1 (MQ=255)      tACCTGTCCGGCAACCGGATCGCAGGTCAGCCCAAGGTTATGCTCCATCGCGATTTCCGCCGCATTGCAt                                                                  <  1:1846827/70‑1 (MQ=255)        ccTGTCCGGCGACCGGGTCACACGTCAGACCGAGGTTGTGCTCCATGGCGATTTCCGCCGCGATGCACAcc                                                               <  1:3773520/71‑1 (MQ=255)        ccTGTCCGGCGACCGGGTCACACGTCAGACCGAGGTTGTGCTCCATGGCGATTTCCGCCGCGATGCACAcc                                                               <  1:621551/71‑1 (MQ=255)         cTGTCCGGCAACCGGATCGCAGGTCAGCCCAAGGTTATGCTCCATCGCGATTTCCGCCGCATTGCATACCt                                                              <  1:3419973/71‑1 (MQ=255)         cTGTCCGGCAACCGGATCGCAGGTCAGCCCAAGGTTATGCTCCATCGCGATTTCCGCCGCATTGCATACCt                                                              <  1:3625457/71‑1 (MQ=255)            tCCGGCGACCGGGTCACACGTCAGACCGAGGTTGTGCTCCATGGCGATTTCCGCCGCGATGCACACCTg                                                             <  1:3655475/69‑1 (MQ=9)              cGGCAACCGGATCGCAGGTCAGCCCAAGGTTATGCTCCATCGCGATTTCCGCCGCATTGCATACCTgcgc                                                          <  1:2887763/70‑1 (MQ=255)                gCAACCGGATCGCAGGTCAGCCCAAGGTTATGCTCAATCGCGATTTCCGCCGCATTGCATACCTGCGCCg                                                        >  1:3141530/1‑70 (MQ=255)                 cgaCCGGGTCACACGTCAGACCGAGGTTGTGCTCCATGGCGATTTCCGCCGCGATGCACACCTGCGCCggg                                                      <  1:1670782/69‑2 (MQ=255)                 cgaCCGGGTCACACGTCAGACCGAGGTTGTGCTCCATGGCGATTTCCGCCGCGATGCACACCTGCGCCggg                                                      <  1:2687725/69‑2 (MQ=255)                  aaCCGGATCGCAGGTCAGCCCAAGGTTATGCTCCATCGCGATTTCCGCCGCATTGCATACCTGCGCCGGAc                                                     <  1:1732695/71‑1 (MQ=255)                  aaCCGGATCGCAGGTCAGCCCAAGGTTATGCTCCATCGCGATTTCCGCCGCATTGCATACCTGCGCCGGAc                                                     >  1:513882/1‑71 (MQ=255)                  aaCCGGATCGCAGGTCAGCCCAAGGTTATGCTCCATCGCGATTTCCGCCGCATTGCATACCTGCGCCGGAc                                                     <  1:3014455/71‑1 (MQ=255)                   aCCGGATCGCAGGTCAGCCCAAGGTTATGCTCCATCGCGATTTCCGCCGCATTGCATACCTGCGCCGGACt                                                    >  1:1686890/1‑71 (MQ=255)                   aCCGGATCGCAGGTCAGCCCAAGGTTATGCTCCATCGCGATTTCCGCCGCATTGCATACCTGCGCCGGACt                                                    >  1:2419385/1‑71 (MQ=255)                    ccGGGTCACACGTCAGACCGAGGTTGTGCTCCATGGCGATTTCCGCCGCGATGCACACCTGCGCCGGGCtt                                                   >  1:2433961/1‑70 (MQ=255)                     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)                                    cccAAGGTTATGGTCCATCGCGATTTCCGCCGCATTGCATACCTgcgc                                                          >  1:2891604/1‑48 (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)                                        aGGTTATGCTCCATCGCGATTTCCGCCGCATTGCATACCTgcgc                                                          >  1:833009/1‑44 (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)                                               gCTCCATCGCGATTTCCGCCGCATTGCATACCTGCGcc                                                         <  1:2275269/38‑1 (MQ=255)                                               gCTCCATCGCGATTTCCGCCGCATTGCATACCTGCGCCGGACTGCCGCCCAGTAGTTCAg                                   >  1:2008961/1‑60 (MQ=255)                                                 tCCATCGCGATTTCCGCCGCATTGCATACCTGCGCCGGACTGCCGCCCAGTAGTTCAGTTAACCCTgccgc                      >  1:2338377/1‑71 (MQ=255)                                                  ccATCGCGATTTCCGCCGCATTGCATACCTGCGCCGGACTGCCGCCCAGTAGTTCAGTTAACCCTgccgcc                     >  1:902015/1‑71 (MQ=255)                                                  ccATCGCGATTTCCGCCGCATTGCATACCTGCGCCGGACTGCCGCCCAGTAGTTCAGTTAACCCTgccgcc                     >  1:3727431/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:1352718/1‑71 (MQ=255)                                                       gcgATTTCCGCCGCATTGCATACCTGCGCCGGACTGCCGCCCAGTAGTTCAGTTAACCCTGCCGCCGCCAt                >  1:790937/1‑71 (MQ=255)                                                       gcgATTTCCGCCGCATTGCATACCTGCGCCGGACTGCCGCCCAGTAGTTCAGTTAACCCTGCCGCCGCCAt                >  1:2618907/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:2260235/1‑69 (MQ=255)                                                             tCCGCCGCATTGCATACCTGCGCCGGACTGCCGCCCAGTAGTTCAGTTAACCCTGCCGCCGCCATTGaa            >  1:3228860/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)                                                                      |                                                                      ATTTGTACCTGTCCGGCAACCGGATCGCAGGTCAGCCCAAGGTTATGCTCCATCGCGATTTCCGCCGCATTGCATACCTGCGCCGGACTGCCGCCCAGTAGTTCAGTTAACCCTGCCGCCGCCATTGAACAGGCCACGC  >  minE/2093683‑2093821

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