Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	2,093,979	G→C	18.7%	G291G (GGC→GGG)	tdcG ←	L‑serine dehydratase 3

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	2,093,979	0	G	C	18.7%	175.1 / 31.1	80	G291G (GGC→GGG)	tdcG	L‑serine dehydratase 3
Reads supporting (aligned to +/- strand):  ref base G (37/28);  new base C (8/7);  total (45/35)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GCAATTGACCGCTCGTTTACCGGACGACGGAACTTATCGTAATAAGCCAGTACTGCCGGAATAATGCCGCACGCACCGTTAGTCGGTGCCGTTACCACGCGCCCGCCAGCTGCGTTTTCTTCACTAACCGC  >  minE/2093914‑2094044                                                                  |                                                                  gCAATTGACCGCTCGTTTACCGGACGACGGAACTTATCGTAATAAGCCAGTACTGCCGGAATAATGCcgca                                                              >  1:1548272/1‑71 (MQ=255)   aaTTGACCGCTCGTTTACCGGACGACGGAACTTATCGTAATAAGCCAGTACTGCCGGAATAATGCcgcacg                                                            >  1:67203/1‑71 (MQ=255)   aaTTGACCGCTCGTTTACCGGACGACGGAACTTATCGTAATAAGCCAGTACTGCCGGAATAATGCcgcacg                                                            >  1:3103246/1‑71 (MQ=255)      tGACCGCTCGTTTACCGGACGACGGAACTTATCGTAATAAGCCAGTACTGCCGGAATAATGCCgcacgcac                                                         >  1:2653440/1‑71 (MQ=255)       gACCGCTCGTTTACCGGACGACGGAACTTATCGTAATAAGCCAGTACTGCCGGAATAATGCCGCACGCAcc                                                        >  1:1605929/1‑71 (MQ=255)        aCCGCTCGTTTACCGGACGACGGAACTTATCGTAATAAGCCAGTACTGCCGGAATAATGCCGCACGCACCg                                                       <  1:2274538/71‑1 (MQ=255)         ccGCTCGTTTACCGGACGACGGAACTTATCGTAATAAGCCAGTACTGCCGGAATAATGCCGCACGCACCGt                                                      <  1:1093574/71‑1 (MQ=255)          cGCTCGTTTACCGGACGACGGAACTTATCGTAATAAGCCAGTACTGCCGGAATAATGCCGCACGCACCGt                                                      >  1:2007860/1‑70 (MQ=255)          cGCTCGTTTACCGGACGACGGAACTTATCGTAATAAGCCAGTACTGCCGGAATAATGCCGCACGCACCGt                                                      >  1:2012393/1‑70 (MQ=255)             tCGTTTACCGGACGACGGAACTTATCGTAATAAGCCAGTACTGCCGGAATAATGCCGCACGCACCGTTAGt                                                  >  1:3372359/1‑71 (MQ=255)               gTTTACCGGACGACGGAACTTATCGTAATAAGCCAGTACTGCCGGAATAATGCCGCACGCACCGTTAGTc                                                 <  1:310760/70‑1 (MQ=255)               gTTTACCGGACGACGGAACTTATCGTAATAAGCCAGTACTGCCGGAATAATGC‑‑CACGCACCGTTAGTCg                                                <  1:1137838/69‑1 (MQ=255)               gTTTACCGGACGACGGAACTTATCGTAATAAGCCAGTAATGCCGGAATAATGCCGCACGCACCGTTAGTCg                                                >  1:3272342/1‑71 (MQ=255)                tttACCGGACGACGGAACTTATCGTAATAAGCCAGTACTGCCGGAATAATGCCGCACGCACCg                                                       <  1:2611566/63‑1 (MQ=255)                  tACCGGACGACGGAACTTATCGTAATAAGCCAGTACTGCCGGAATAATGCCGCACGCACCGTTAGTCGGt                                              <  1:746612/70‑1 (MQ=255)                    ccGGACGACGGAACTTATCGTAATAAGCCAGTACTGCCGGAATAATGCCGCACGCACCGTTAGTCGGTGcc                                           <  1:1854728/71‑1 (MQ=255)                    ccGGACGACGGAACTTATCGTAATAAGCCAGTACTGCCGGAATAATGCCGCACGCACCGTTAGTCGGTGcc                                           <  1:3129887/71‑1 (MQ=255)                     cGGACGACGGAACTTATCGTAATAAGCCAGTACTGCCGGAATAATGCCGCACGCACCGTTAGTCGGTGCCg                                          >  1:893794/1‑71 (MQ=255)                     cGGACGACGGAACTTATCGTAATAAGCCAGTACTGCCGGAATAATGCCGCACGCACCGTTAGTCGGTGCCg                                          <  1:1144772/71‑1 (MQ=255)                     cGGACGACGGAACTTATCGTAATAAGCCAGTACTGCCGGAATAATGCCGCACGCACCGTTAGTCGGTGCCg                                          >  1:3357418/1‑71 (MQ=255)                     cGGACGACGGAACTTATCGTAATAAGCCAGTACTGCCGGAATAATGCCGCACGCACCGTTAGTCGGTGCCg                                          >  1:2878945/1‑71 (MQ=255)                      ggACGACGGAACTTATCGTAATAAGCCAGTACTGCCGGAATAATGCCGCACGCACCGTTAGTCGGTGCCGt                                         <  1:989600/71‑1 (MQ=255)                       gacgacGGAACTTATCGTAATAAGCCAGTACTGCCGGAATAATGCCGCACGCACCGTTAGTCGGTGCCGt                                         >  1:3551870/1‑70 (MQ=255)                       gacgacGGAACTTATCGTAATAAGCCAGTACTGCCGGAATAATGCCGCACGCACCGTTAGTCGGTGCCGt                                         >  1:409085/1‑70 (MQ=255)                         cgacgGAACTTATCGTAATAAGCCAGTACTGCCGGAATAATGCCGCACGCACCGTTAGTCGGTGCCGTTAc                                      <  1:1173665/71‑1 (MQ=255)                         cgacgGAACTTATCGTAATAAGCCAGTACTGCCGGAATAATGCCGCACGCACCGTTAGTCGGTGCCGTTAc                                      <  1:3261656/71‑1 (MQ=255)                          gacgGAACTTATCGTAATAAGCCAGTACTGCCGGAATAATGCCGCACGCACCGTTAGTCGGTGCCGTTa                                       <  1:3471859/69‑1 (MQ=255)                          gacgGAACTTATCGTAATAAGCCAGTACTGCCGGAATAATGCCGCACGCACCGTTAGTCGGTGCCGTTAcc                                     <  1:1039982/71‑1 (MQ=255)                          gacgGAACTTATCGTAATAAGCCAGTACTGCCGGAATAATGCCGCACGCACCGTTAGTCGGTGCCGTTAcc                                     <  1:3532993/71‑1 (MQ=255)                          gacgGAACTTATCGTAATAAGCCAGTACTGCCGGAATAATGCCGCACGCACCGTTAGTCGGTGCCGTTAcc                                     <  1:84647/71‑1 (MQ=255)                             taaaCTTGTCGTAGTACGCCAGAACTGCCGGGATAATCCCGCACGCACCGTTAGTCGGCGCAGTcaccac                                   <  1:711265/68‑1 (MQ=255)                              gAACTTATCGTAATAAGCCAGTACTGCCGGAATAATGCCGCACGCACCGTTAGTCGGTGCCGTTACCAcgc                                 <  1:186443/71‑1 (MQ=255)                               aaCTTGTCGTAGTACGCCAGAACTGCCGGGATAATCCCGCACGCACCGTTAGTCGGCGCAGTCACCAcgc                                 <  1:163734/70‑1 (MQ=255)                               aaCTTATCGTAATAAGCCAGTACTGCCGGAATAATGCCGCACGCACCGTTAGTCGGTGCCGTTACCAcgcg                                >  1:99484/1‑71 (MQ=255)                               aaCTTATCGTAATAAGCCAGTACTGCCGGAATAATGCCGCACGCACCGTTAGTCGGTGCCGTTACCAcgcg                                >  1:279286/1‑71 (MQ=255)                                aCTTATCGTAATAAGCCAGTACTGCCGGAATAATGCcgcacgc                                                           <  1:3901322/43‑1 (MQ=255)                                 cTTATCGTAATAAGCCAGTACTGCCGGAATAATGCCGCACGCACCGt                                                      >  1:455183/1‑47 (MQ=255)                                 cTTATCGTAATAAGCCAGTACTGCCGGAATAATGCCGCACGCACCGTTa                                                    <  1:653536/49‑1 (MQ=255)                                 cTTATCGTAATAAGCCAGTACTGCCGGAATAATGCCGCACGCACCGTTAGTCGGTGCCGTTACCACgcgc                               >  1:389199/1‑70 (MQ=255)                                 cTTATCGTAATAAGCCAGTACTGCCGGAATAATGCCGCACGCACCGTTAGTCGGTGCCGTTACCACGcgcc                              >  1:2520396/1‑71 (MQ=255)                                 cTTATCGTAATAAGCCAGTACTGCCGGAATAATGCCGCACGCACCGTTAGTCGGTGCCGTTACCACGcgcc                              >  1:2515975/1‑71 (MQ=255)                                  ttATCGTAATAAGCCAGTACTGCCGGAATAATGCCGCACGCACCGTTAGTCGGTGCCGTTACCAcg                                  <  1:3082986/66‑1 (MQ=255)                                   tgtCGTAGTACGCCAGAACTGCCGGGATAATCCCGCACGCACCGTTAGTCGGCGCAGTCACCACGcgaccg                            >  1:3615862/3‑71 (MQ=255)                                   tgtCGTAGTACGCCAGAACTGCCGGGATAATCCCGCACGCACCGTTAGTCGGCGCAGTCACCACGcgacc                             >  1:3278248/3‑70 (MQ=255)                                    aTCGTAATAAGCCAGTACTGCCGGAATAATGCCGCACGCACCGTTAGTc                                                 >  1:3603563/1‑49 (MQ=255)                                    aTCGTAATAAGCCAGTACTGCCGGAATAATGCCGCACGCACCGTTAGTCGGTGCCGTTACCACGcgcccgc                           >  1:1591342/1‑71 (MQ=255)                                    aTCGTAATAAGCCAGTACTGCCGGAATAATGCCGCACGCACCGTTAGTCGGTGCCGTTACCACGcgcccgc                           >  1:1624712/1‑71 (MQ=255)                                    aTCGTAATAAGCCAGTACTGCCGGAATAATGCCGCACGCACCGTTAGTCGGTGCCGTTACCACGcgcccgc                           >  1:73143/1‑71 (MQ=255)                                    aTCGTAATAAGCCAGTACTGCCGGAATAATGCCGCACGCACCGTTAGTCGGTGCCGTTACCACGcgcccgc                           >  1:1539540/1‑71 (MQ=255)                                     tCGTAATAAGCCAGTACTGCCGGAATAATGCCGCACGCACCGTTAGTCGGTGCCGTTACCACGcgcccgcc                          <  1:2693478/71‑1 (MQ=255)                                     tCGTAATAAGCCAGTACTGCCGGAATAATGCCGCACGCACCGTTAGTCGGTGCCGTTACCACGcgcccgc                           <  1:2267458/70‑1 (MQ=255)                                      cGTAATAAGCCAGTACTGCCGGAATAATGCCGCACGCACCGTTAGTCGGTGCCGTTACCACGcgcccgc                           >  1:1069282/1‑69 (MQ=255)                                      cGTAATAAGCCAGTACTGCCGGAATAATGCCGCACGCACCGTTAGTCGGTGCCGTTACCACGcgcccgc                           >  1:3070272/1‑69 (MQ=255)                                       gTAATAAGCCAGTACTGCCGGAATAATGCCGCACGCACCgac                                                     <  1:151593/42‑3 (MQ=255)                                       gTAATAAGCCAGTACTGCCGGAATAATGCCGCACGCACCGTTAGTCGGTGCCGTTa                                       >  1:1458611/1‑56 (MQ=255)                                       gTAATAAGCCAGTACTGCCGGAATAATGCCGCACGCACCGTTAGTCGGTGCCGTTACCACGCGCCCGCCAg                        <  1:885111/71‑1 (MQ=255)                                       gTAATAAGCCAGTACTGCCGGAATAATGCCGCACGCACCGTTAGTCGGTGCCGTTACCACGCGCCCGCCAg                        >  1:2664009/1‑71 (MQ=255)                                        taataaGCCAGTACTGCCGGAATAATGCCGCACGCACCGTTAGTCGGTGCCGTTACCACGCGCCCGCCAGc                       >  1:2931244/1‑71 (MQ=255)                                         aataaGCCAGTACTGCCGGAATAATGCCGCACGCACCGTTAGTCGGTGCCGTTACCACGcgccc                             >  1:103838/1‑64 (MQ=255)                                             cgCCAGAACTGCCGGGATAATCCCGCACGCACCGTTAGTCGGCGCAGTCACCACGCGACCGCCAGCAGCGt                  <  1:878775/70‑1 (MQ=255)                                              gCCAGTACTGCCGGAATAATGCCGCACGCACCGTTAGTCGGTGCCGTTACCACGCGCCCGCCAGCTGCGtt                 <  1:3639644/71‑1 (MQ=255)                                              gCCAGAACTGCCGGGATAATCCCGCACGCACCGTTAGTCGGCTCAGTCACCACGCGACCGCCAGCAGCGtt                 <  1:1300937/71‑1 (MQ=255)                                               ccAGAACTGCCGGGATAATCCCGCACGCACCGTTAGTCGGCGCAGTCACCACGCGACCGCCAGCAGCGttc                <  1:344529/71‑2 (MQ=255)                                                cAGAACTGCCGGGATAATCCCGCACGCACCGTTAGTCGGCGCAGTCACCAcgcgac                              >  1:1228213/1‑54 (MQ=255)                                                cAGAACTGCCGGGATAATCCCGCACGCACCGTTAGTCGGCGCAGTCACCACGCGACCGCCAGCAGCGttct               >  1:3765492/1‑69 (MQ=255)                                                  gaaCTGCCGGGATAATCCCGCACGCACCGTTAGTCGGCGCAGTCACCACGCGACCGCCAg                        >  1:636282/3‑60 (MQ=255)                                                  gTACTGCCGGAATAATGCCGCACGCACCGTTAGTCGGTGCCGTTACCACGCGCCCGCCAGCTGCGTTttct             >  1:3707312/1‑71 (MQ=255)                                                  gTACTGCCGGAATAATGCCGCACGCACCGTTAGTCGGTGCCGTTACCACGCGCCCGCCAGCTGCGTTttct             >  1:2667918/1‑71 (MQ=255)                                                   tACTGCCGGAATAATGCCGCACGCACCGTTAGTCGGTGCCGTTACCACGCGCCCGCCAGCTGCGTTttctt            <  1:3127909/71‑1 (MQ=255)                                                   aaCTGCCGGGATAATCCCGCACGCACCGTTAGTCGGCGCAGTCACCACGCGACCGCCAGCAGCGTTctctt            >  1:3484239/2‑71 (MQ=255)                                                    aCTGCCGGGATAATCCCGCACGCACCGTTAGTCGGCGCAGTCACCACGcgacc                             >  1:3847039/1‑53 (MQ=255)                                                     cTGCCGGGATAATCCCGCACGCACCGTTAGTCGGCGCAGTCACCACGCGACCGCCAGCAGCGTTctcttcg          <  1:1568880/71‑2 (MQ=255)                                                     cTGCCGGAATAATGCCGCACGCACCGTTAGTCGGTGCCGTTACCACGCGCCCGCCAGCTGCGTTttcttc           >  1:1453895/1‑70 (MQ=255)                                                       gCCGGGATAATCCCGCACGCACCGTTAGTCGGCGCAGTCACCACGCGACCGCCAGCAGCGTtctc              <  1:346085/65‑1 (MQ=255)                                                       gCCGGAATAATGCCGCACGCACCGTTAGTCGGTGCCGTTACCACGCGCCCGCCAGCTGCGTTttctt            >  1:1656806/1‑67 (MQ=255)                                                        ccGGGATAATCCCGCACGCACCGTTAGTCGGCGCAGTCACCACGCGACCGCCAGCAGCGtt                 >  1:986529/1‑61 (MQ=255)                                                         cGGAATAATGCCGCACGCACCGTTAGTCGGTGCCGTTACCACGCGCCCGCCAGCTGCGTTTTCTTCACTaa      <  1:421292/71‑1 (MQ=255)                                                         cGGAATAATGCCGCACGCACCGTTAGTCGGTGCCGTTACCACGCGCCCGCCAGCTGCGTTTTCTTCACTaa      <  1:2554477/71‑1 (MQ=255)                                                          ggAATAATGCCGCACGCACCGTTAGTCGGTGCCGTTACCACGCGCCCGCCAGCTGCGTTTTCTTCACTAAc     >  1:3107302/1‑71 (MQ=255)                                                             ataatGCCGCACGCACCGTTAGTCGGTGCCGTTACCACGcgcccgcc                          >  1:3078741/1‑47 (MQ=255)                                                                  gCCGCACGCACCGTTAGTCGGTGCCGTTACCACGCGCCCGCCAGCTGCGTTTTCTTCACTAACCGc  <  1:1832254/66‑1 (MQ=255)                                                                  |                                                                  GCAATTGACCGCTCGTTTACCGGACGACGGAACTTATCGTAATAAGCCAGTACTGCCGGAATAATGCCGCACGCACCGTTAGTCGGTGCCGTTACCACGCGCCCGCCAGCTGCGTTTTCTTCACTAACCGC  >  minE/2093914‑2094044

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