Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	2,094,006	T→C	21.1%	V282V (GTA→GTG)	tdcG ←	L‑serine dehydratase 3

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	2,094,006	0	T	C	21.1%	139.7 / 33.6	72	V282V (GTA→GTG)	tdcG	L‑serine dehydratase 3
Reads supporting (aligned to +/- strand):  ref base T (29/27);  new base C (8/7);  total (38/34)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.88e-01

CGACGGAACTTATCGTAATAAGCCAGTACTGCCGGAATAATGCCGCACGCACCGTTAGTCGGTGCCGTTACCACGCGCCCGCCAGCTGCGTTTTCTTCACTAACCGCCAGCGCGTACATGTTGATCCAGTCGATGACAT  >  minE/2093938‑2094076                                                                     |                                                                       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:84647/71‑1 (MQ=255)  gacgGAACTTATCGTAATAAGCCAGTACTGCCGGAATAATGCCGCACGCACCGTTAGTCGGTGCCGTTAcc                                                                     <  1:3532993/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)         cTTATCGTAATAAGCCAGTACTGCCGGAATAATGCCGCACGCACCGTTAGTCGGTGCCGTTACCACgcgc                                                               >  1:389199/1‑70 (MQ=255)         cTTATCGTAATAAGCCAGTACTGCCGGAATAATGCCGCACGCACCGTTAGTCGGTGCCGTTACCACGcgcc                                                              >  1:2515975/1‑71 (MQ=255)         cTTATCGTAATAAGCCAGTACTGCCGGAATAATGCCGCACGCACCGTTAGTCGGTGCCGTTACCACGcgcc                                                              >  1:2520396/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)            aTCGTAATAAGCCAGTACTGCCGGAATAATGCCGCACGCACCGTTAGTCGGTGCCGTTACCACGcgcccgc                                                           >  1:1624712/1‑71 (MQ=255)            aTCGTAATAAGCCAGTACTGCCGGAATAATGCCGCACGCACCGTTAGTCGGTGCCGTTACCACGcgcccgc                                                           >  1:73143/1‑71 (MQ=255)            aTCGTAATAAGCCAGTACTGCCGGAATAATGCCGCACGCACCGTTAGTCGGTGCCGTTACCACGcgcccgc                                                           >  1:1591342/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)               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:2667918/1‑71 (MQ=255)                          gTACTGCCGGAATAATGCCGCACGCACCGTTAGTCGGTGCCGTTACCACGCGCCCGCCAGCTGCGTTttct                                             >  1:3707312/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:2531017/70‑1 (MQ=255)                                                         aGTCGGTGCCGTTACCACGCGCCCGCCAGCTGCGTTTTCTTCACTAACCGCCAGCGCGTACATGTTGATc               <  1:745603/70‑1 (MQ=255)                                                          gTCGGTGCCGTTACCACGCGCCCGCCAGCTGCGTTTTCTTCACTAACCGCCAGCGCGTACATGTTGATCCa             <  1:1597341/71‑1 (MQ=255)                                                           tCGGTGCCGTTACCACGCGCCCGCCAGCTGCGTTTTCTTCACTAACCGCCAGCGCGTACATGtt                   <  1:844842/64‑1 (MQ=255)                                                                 ccGTTACCACGCGCCCGCCAGCTGCGTTTTCTTCACTAACCGCCAGCGCGTACATGTTGATCCAGTCGAt       <  1:2473296/70‑1 (MQ=255)                                                                 ccGTTACCACGCGCCCGCCAGCTGCGTTTTCTTCACTAACCGCCAGCGCGTACATGTTGATCCAGTCGAt       <  1:1377253/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)                                                                     |                                                                       CGACGGAACTTATCGTAATAAGCCAGTACTGCCGGAATAATGCCGCACGCACCGTTAGTCGGTGCCGTTACCACGCGCCCGCCAGCTGCGTTTTCTTCACTAACCGCCAGCGCGTACATGTTGATCCAGTCGATGACAT  >  minE/2093938‑2094076

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