Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	2,094,015	C→A	18.2%	G279G (GGG→GGT)	tdcG ←	L‑serine dehydratase 3

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	2,094,015	0	C	A	18.2%	147.5 / 13.1	54	G279G (GGG→GGT)	tdcG	L‑serine dehydratase 3
Reads supporting (aligned to +/- strand):  ref base C (24/20);  new base A (4/6);  total (28/26)
Fisher's exact test for biased strand distribution p-value = 4.94e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CTTATCGTAATAAGCCAGTACTGCCGGAATAATGCCGCACGCACCGTTAGTCGGTGCCGTTACCACGCGCCCGCCAGCTGCGTTTTCTTCACTAACCGCCAGCGCGTACATGTTGATCCAGTCGATGACATTCATCGG  >  minE/2093946‑2094083                                                                      |                                                                     cTTATCGTAATAAGCCAGTACTGCCGGAATAATGCCGCACGCACCGTTAGTCGGTGCCGTTACCACgcgc                                                                      >  1:389199/1‑70 (MQ=255) cTTATCGTAATAAGCCAGTACTGCCGGAATAATGCCGCACGCACCGTTAGTCGGTGCCGTTACCACGcgcc                                                                     >  1:2515975/1‑71 (MQ=255) cTTATCGTAATAAGCCAGTACTGCCGGAATAATGCCGCACGCACCGTTAGTCGGTGCCGTTACCACGcgcc                                                                     >  1:2520396/1‑71 (MQ=255)   tgtCGTAGTACGCCAGAACTGCCGGGATAATCCCGCACGCACCGTTAGTCGGCGCAGTCACCACGcgaccg                                                                   >  1:3615862/3‑71 (MQ=255)   tgtCGTAGTACGCCAGAACTGCCGGGATAATCCCGCACGCACCGTTAGTCGGCGCAGTCACCACGcgacc                                                                    >  1:3278248/3‑70 (MQ=255)    aTCGTAATAAGCCAGTACTGCCGGAATAATGCCGCACGCACCGTTAGTCGGTGCCGTTACCACGcgcccgc                                                                  >  1:73143/1‑71 (MQ=255)    aTCGTAATAAGCCAGTACTGCCGGAATAATGCCGCACGCACCGTTAGTCGGTGCCGTTACCACGcgcccgc                                                                  >  1:1539540/1‑71 (MQ=255)    aTCGTAATAAGCCAGTACTGCCGGAATAATGCCGCACGCACCGTTAGTCGGTGCCGTTACCACGcgcccgc                                                                  >  1:1591342/1‑71 (MQ=255)    aTCGTAATAAGCCAGTACTGCCGGAATAATGCCGCACGCACCGTTAGTCGGTGCCGTTACCACGcgcccgc                                                                  >  1:1624712/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)       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)                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:2554477/71‑1 (MQ=255)                         cGGAATAATGCCGCACGCACCGTTAGTCGGTGCCGTTACCACGCGCCCGCCAGCTGCGTTTTCTTCACTaa                                             <  1:421292/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)                                     gcacgcacCGTTAGTCGGTGCGGTTACCACGCGCCCGCCAGCTGCGTTTTCTTCACTAACCGCCAGcgcg                                  >  1:777986/1‑70 (MQ=255)                                     gcacgcacCGTTAGTCGGTGCCGTTACCACGCGCCCGCCAGCTGCGTTTTCTTCACTAACCGCCAgcgc                                   >  1:1339390/1‑69 (MQ=255)                                      cacgcacCGTTAGTCGGTGCCGTTACCACGCGCCCGCCAGCTGCGTTTTCTTCACTAACCGCCAGCGCGt                                 >  1:2436712/1‑70 (MQ=255)                                         gcacCGTTAGTCGGTGCCGTTACCACGCGCCCGCCAGCTGCGTTTTCTTCACTAACCGc                                         <  1:1597629/59‑1 (MQ=255)                                         gcacCGTTAGTCGGTGCCGTTACCACGCGCCCGCCAGCTGCGTTTTCTTCACTAACCGCCAgcg                                    <  1:1285285/64‑1 (MQ=255)                                         gcacCGTTAGTCGGTGCCGTTACCACGCGCCCGCCAGCTGCGTTTTCTTCACTAACCGCCAGCGCGTACAt                             >  1:1447402/1‑71 (MQ=255)                                          cacCGTTAGTCGGCGCAGTCACCACGCGACCGCCAGCAGCGTTCTCTTCGTTCACTGCCAGTGCAAACATg                            <  1:1704461/71‑1 (MQ=255)                                             cGTTAGTCGGTGCCGTGACCACGCGCCCGCCAGCTGCGTTttct                                                    >  1:56041/1‑44 (MQ=255)                                              gTTAGTCGGTGCCGTTACCACGCGCCCGCCAGCTGCGTTTTCTTCACt                                               <  1:1731682/48‑1 (MQ=255)                                              gTTAGTCGGTGCCGTTACCACGCGCCCGCCAGCTGCGTTTTCTTCACTAACCTCCAGCGCGTACATGTTGa                        <  1:3195049/71‑1 (MQ=255)                                                tAGTCGGTGCCGTTACCACGCGCCCGCCAGCTGCGTTTTCTTCACTAACCGCCAGCGCGTACATGTTGATc                      >  1:3610100/1‑71 (MQ=255)                                                 aGTCGGTGCCGTTACCACGCGCCCGCCAGCTGCGTTTTCTTCACTAACCGCCAGCGCGTACATGTTGATc                      <  1:745603/70‑1 (MQ=255)                                                 aGTCGGTGCCGTTACCACGCGCCCGCCAGCTGCGTTTTCTTCACTAACCGCCAGCGCGTACATGTTGATc                      <  1:2531017/70‑1 (MQ=255)                                                  gTCGGTGCCGTTACCACGCGCCCGCCAGCTGCGTTTTCTTCACTAACCGCCAGCGCGTACATGTTGATCCa                    <  1:1597341/71‑1 (MQ=255)                                                   tCGGTGCCGTTACCACGCGCCCGCCAGCTGCGTTTTCTTCACTAACCGCCAGCGCGTACATGtt                          <  1:844842/64‑1 (MQ=255)                                                         ccGTTACCACGCGCCCGCCAGCTGCGTTTTCTTCACTAACCGCCAGCGCGTACATGTTGATCCAGTCGAt              <  1:1377253/70‑1 (MQ=255)                                                         ccGTTACCACGCGCCCGCCAGCTGCGTTTTCTTCACTAACCGCCAGCGCGTACATGTTGATCCAGTCGAt              <  1:2473296/70‑1 (MQ=255)                                                         ccGTTACCACGCGCCCGCCAGCTGCGTTTTCTTCACTAACCGCCAGCGCGTACATGTTGATCCAGTCGAt              >  1:291964/1‑70 (MQ=255)                                                          cGTTACGACGCGCCCGCCAGCTGCGTTTTCTTCACTAACCGCCAGCGCGTACATGTTGAt                       >  1:3359/1‑60 (MQ=255)                                                          cGTTACCACGCGCCCGCCAGCTGCGTTTTCTTCACTAACCGCCAGCGCGTACATGTTGAt                       >  1:1839150/1‑60 (MQ=255)                                                           gTTACCACGCGCCCGCCAGCTGCGTTTTCTTCACTAACCGCCAGCGCGTACATGTTGa                        <  1:189948/58‑1 (MQ=255)                                                            ttACCACGCGCCCGCCAGCTGCGTTTTCTTCACTAACCGCCAGCGCGTACATGTTGATCCAGTCGATGACa          >  1:3949123/1‑71 (MQ=255)                                                             tACCACGCGCCCGCCAGCTGCGTTTTCTTCACTAACCGCCAGCGCGTACATGTTGATCCAGTCGATGACAt         >  1:1978962/1‑71 (MQ=255)                                                                cACGCGCCCGCCAGCTGCGTTTTCTTCACTAACCGCCAGCGCGTACAt                             <  1:773673/48‑1 (MQ=255)                                                                cACGCGCCCGCCAGCTGCGTTTTCTTCACTAACCGCCAGCGCGTACATGTTGATCCAGTCGATGACATTCa      >  1:3487004/1‑71 (MQ=255)                                                                cACGCGCCCGCCAGCTGCGTTTTCTTCACTAACCGCCAGCGCGTACATGTTGATCCAGTCGATGACATTCa      <  1:963412/71‑1 (MQ=255)                                                                cACGCGCCCGCCAGCTGCGTTTTCTTCACTAACCGCCAGCGCGTACATGTTGATCCAGTCGATGACATTCa      <  1:3628263/71‑1 (MQ=255)                                                                 aCGCGCCCGCCAGCTGCGTTTTCTTCACTAACCGCCAGCGCGTACATGTTGATCCAGTCGATGACAtt        >  1:1038050/1‑68 (MQ=255)                                                                  cgcgCCCGCCAGCTGCGTTTTCTTCACTAACCGCCAGCGCGTACATGTTGATCCAGTCGATGACATTCATc    >  1:2447557/1‑71 (MQ=255)                                                                  cgcgCCCGCCAGCTGCGTTTTCTTCACTAACCGCCAGCGCGTACATGTTGATCCAGTCGATGACATTCATc    >  1:1682194/1‑71 (MQ=255)                                                                  cgcgCCCGCCAGCTGCGTTTTCTTCACTAACCGCCAGCGCGTACATGTTGATCCAGTCGATGACATTCATc    >  1:2627380/1‑71 (MQ=255)                                                                  cgcgCCCGCCAGCTGCGTTTTCTTCACTAACCGCCAGCGCGTACATGTTGATCCAGTCGAAGACATTCATc    >  1:535073/1‑71 (MQ=255)                                                                   gcgcCCGCCAGCTGCGTTTTCTTCACTAACCGCCAGCGCGTACATGTTGATCCAGt                  >  1:3897946/1‑56 (MQ=255)                                                                     gcccgccAGCTGCGTTTTCTTCACTAACCGCCAGCGCGTACATGTTGATCCAGTCGATGACATTCATCgg  <  1:2891269/70‑1 (MQ=255)                                                                      |                                                                     CTTATCGTAATAAGCCAGTACTGCCGGAATAATGCCGCACGCACCGTTAGTCGGTGCCGTTACCACGCGCCCGCCAGCTGCGTTTTCTTCACTAACCGCCAGCGCGTACATGTTGATCCAGTCGATGACATTCATCGG  >  minE/2093946‑2094083

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