Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	2,094,006	T→C	23.4%	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	23.4%	132.9 / 40.1	77	V282V (GTA→GTG)	tdcG	L‑serine dehydratase 3
Reads supporting (aligned to +/- strand):  ref base T (36/23);  new base C (14/4);  total (50/27)
Fisher's exact test for biased strand distribution p-value = 2.63e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.38e-01

GACGACGGAACTTATCGTAATAAGCCAGTACTGCCGGAATAATGCCGCACGCACCGTTAGTCGGTGCCGTTACCACGCGCCCGCCAGCTGCGTTTTCTTCACTAACCGCCAGCGCGTACATGTTGATCCAGTCGATGACA  >  minE/2093936‑2094075                                                                       |                                                                      gacgacGGAACTTATCGTAATAAGCCAGTACTGCCGGAATAATGCCGCACGCACCGTTAGTCGGTGCCGtt                                                                       <  1:3766914/71‑1 (MQ=255)  acgacgGAACTTATCGTAATAAGCCAGTACTGCCGGAATAATGCCGCACGCACCGTTAGTCGGTGCCGTTa                                                                      <  1:4912529/71‑1 (MQ=255)    gacgGAACTTATCGTAATAAGCCAGTACTGCCGGAATAATGCCGCACGCACCGTTAGTCGGTGCCGTTAcc                                                                    <  1:330598/71‑1 (MQ=255)    gacgGAACTTATCGTAATAAGCCAGTACTGCCGGAATAATGCCGCACGCACCGTTAGTCGGTGCCGTTAcc                                                                    <  1:5573537/71‑1 (MQ=255)       taaaCTTGTCGTAGTACGCCAGAACTGCCGGGATAATCCCGCACGCACCGTTAGTCGGCGCAGTCAc                                                                     >  1:5127424/3‑67 (MQ=255)       ggAACTTATCGTAATAAGCCAGTACTGCCGGAATAATGCCGCACGCACCGTTAGTCGGTGCCGTTACCAcg                                                                 >  1:164002/1‑71 (MQ=255)          aCTTATCGTAATAAGCCAGTACTGCCGGAATAATGCCGCACGCACCGTTAGTCGGTGCCGTTACCACgcgc                                                              >  1:3708721/1‑71 (MQ=255)          aCTTATCGTAATAAGCCAGTACTGCCGGAATAATGCCGCACGCACCGTTAGTCGGTGCCGTTACCACgcgc                                                              >  1:4843/1‑71 (MQ=255)            ttgtCGTAGTACGCCAGAACTGCCGGGATAATCCCGCACGCACCGTTAGTCGGCGCAGTCAc                                                                     <  1:3702030/59‑1 (MQ=255)            ttgtCGTAGTACGCCAGAACTGCCGGGATAATCCCGCACGCACCGTTAGTCGGCGCAGTCACCACGcgacc                                                            >  1:938376/4‑71 (MQ=255)            ttgtCGTAGTACGCCAGAACTGCCGGGATAATCCCGCACGCACCGTTAGTCGGCGCAGTCACCACGcgacc                                                            >  1:5468454/4‑71 (MQ=255)            ttgtCGTAGTACGCCAGAACTGCCGGGATAATCCCGCACGCACCGTTAGTCGGCGCAGTCACCACGcgacc                                                            >  1:2066948/4‑71 (MQ=255)             tATCGTAATAAGCCAGTACTGCCGGAATAATGCCGCACGCACCGTTAGTCGGTGCCGTTACCAcgc                                                                >  1:217407/1‑66 (MQ=255)             tATCGTAATAAGCCAGTACTGCCGGAATAATGCCGCACGCACCGTTAGTCGGTGCCGTTACCAcgc                                                                >  1:3890246/1‑66 (MQ=255)             tATCGTAATAAGCCAGTACTGCCGGAATAATGCCGCACGCACCGTTAGTCGGTGCCGTTACCAcgc                                                                >  1:5502257/1‑66 (MQ=255)              gtCGTAGTACGCCAGAACTGCCGGGATAATCCCGCACGCACCGTTAGTCGGCGCAGTCACCa                                                                   >  1:3176502/2‑62 (MQ=255)              gtCGTAGTACGCCAGAACTGCCGGGATAATCCCGCACGCACCGTTAGTCGGCGCAGTCACCa                                                                   >  1:3169996/2‑62 (MQ=255)              aTCGTAATAAGCCAGTACTGCCGGAATAATGCCGCACGCACCGTTAGTCGGTGCCGTTa                                                                      >  1:4688323/1‑59 (MQ=255)              aTCGTAATAAGCCAGTACTGCCGGAATAATGCCGCACGCACCGTTAGTCGGTGCCGTTACCACGcgcccgc                                                          >  1:3695590/1‑71 (MQ=255)              aTCGTAATAAGCCAGTACTGCCGGAATAATGCCGCACGCACCGTTAGTCGGTGCCGTTACCACGcgcccgc                                                          >  1:2676583/1‑71 (MQ=255)              aTCGTAATAAGCCAGTACTGCCGGAATAATGCCGCACGCACCGTTAGTCGGTGCCGTTACCACGcgcccgc                                                          >  1:4637198/1‑71 (MQ=255)              aTCGTAATAAGACAGTACTGCCGGAATAATGCCGCACGCACCGTTAGTCGGTGCCGTTa                                                                      >  1:168245/1‑59 (MQ=255)               tCGTAGTACGCCAGAACTGCCGGGATAATCCCGCACGCACCGTTAGTCGGCGCAGTcaccac                                                                  >  1:3411282/1‑62 (MQ=255)               tCGTAGTACGCCAGAACTGCCGGGATAATCCCGCACGCACCGTTAGTCGGCGCAGTcaccac                                                                  >  1:1527158/1‑62 (MQ=255)                 gtagtaCGCCAGAACTGCCGGGATAATCCCGCACGCACCGTTAGTCGGCGCAGTCACCACGCGACCGCCAg                                                       >  1:806680/1‑71 (MQ=255)                    ataaGCCAGTACTGCCGGAATAATGCCGCACGCACCGTTAGTCGGTGCCGTTACCACGCGCCCGCCAGCTg                                                    >  1:4459363/1‑71 (MQ=255)                      aaGCCAGTACTGCCGGAATAATGCCGCACGCACCGTTAGTCGGTGCCGTTACCa                                                                   >  1:3661131/1‑54 (MQ=255)                       aGCCAGTACTGCCGGAATAATGCCGCACGCACCGTTAGTCGGTGCCGTTACCACGCGCCCGCCAGCTGCg                                                  <  1:659736/70‑1 (MQ=255)                       aGCCAGTACTGCCGGAATAATGCCGCACGCACCGTTAGTCGGTGCCGTTACCACGCGCCCGCCAGCTGCGt                                                 >  1:2956915/1‑71 (MQ=255)                        gCCAGTACTGCCGGAATAATGCCGCACGCACCGTTAGTCGGTGCCGTTACCACGCGCCCGCCAGCTGCg                                                  <  1:3037432/69‑1 (MQ=255)                        gCCAGTACTGCCGGAATAATGCCGCACGCACCGTTAGTCGGTGCCGTTACCACGCGCCCGCCAGCTGCGtt                                                <  1:5166382/71‑1 (MQ=255)                        gCCAGTACTGCCGGAATAATGCCGCACGCACCGTTAGTCGGTGCCGTTACCACGCGCCCGCCAGCTGCGtt                                                <  1:465832/71‑1 (MQ=255)                         ccAGAACTGCCGGGATAATCCCGCACGCACCGTTAGTCGGCGCAGTCACCACGCGACCGCCAGCAGCGttc                                               >  1:2181308/1‑70 (MQ=255)                         ccAGAACTGCCGGGATAATCCCGCACGCACCGTTAGTCGGCGCAGTCACCACGCGACCGCCAGCAGCGttc                                               >  1:2620699/1‑70 (MQ=255)                         ccAGAACTGCCGGGATAATCCCGCACGCACCGTTAGTCGGCGCAGTCACCACGCGACCGCCAGCAGCGttc                                               <  1:2552742/71‑2 (MQ=255)                          cAGTACTGCCGGAATAATGCCGCACGCACCGTTAGTCGGTGCCGTTACCACGCGCCCGCCAGCTGCGtttt                                              <  1:182662/71‑1 (MQ=255)                          cAGTACTGCCGGAATAATGCCGCACGCACCGTTAGTCGGTGCCGTTACCACGCGCCCGCCAGCTGCGtttt                                              <  1:2560447/71‑1 (MQ=255)                          cAGTACTGCCGGAATAATGCCGCACGCACCGTTAGTCGGTGCCGTTACCACGCGCCCGCCAGCTGCGtttt                                              >  1:1294979/1‑71 (MQ=255)                           aGTACTGCCGGAATAATGCCGCACGCACCGTTAGTCGGTGCCGTTACCACGCGCCCGCCAGCTGCGTTttc                                             >  1:3581412/1‑71 (MQ=255)                           aGTACTGCCGGAATAATGCCGCACGCACCGTTAGTCGGTGCCGTTACCACGCGCCCGCCAGCTGCGTTttc                                             >  1:58555/1‑71 (MQ=255)                             tACTGCCGGAATAATGCCGCACGCACCGTTAGTCGGTGCCGTTACCACGCGCCCGCCAGCTGCGTTttctt                                           <  1:3112079/71‑1 (MQ=255)                             tACTGCCGGAATAATGCCGCACGCACCGTTAGTCGGTGCCGTTACCACGCGCCCGCCAGCTGCGTTttctt                                           <  1:4390054/71‑1 (MQ=255)                               cTGCCGGAATAATGCCGCACGCACCGTTAGTCGGTGCCGtt                                                                       >  1:3661615/1‑41 (MQ=255)                               cTGCCGGAATAATGCCGCACGCACCGTTAGTCGGTGCCGtt                                                                       >  1:3664791/1‑41 (MQ=255)                               cTGCCGGAATAATGCCGCACGCACCGTTAGTCGGTGCCGtt                                                                       >  1:3668176/1‑41 (MQ=255)                               cTGCCGGAATAATGCCGCACGCACCGTTAGTCGGTGCCGTTACCACGCGCCCGCCAGCTGCGTTttctt                                           >  1:1208784/1‑69 (MQ=255)                                tGCCGGAATAATGCCGCACGCACCGTTAGTCGGTGCCGTTa                                                                      <  1:3043430/41‑1 (MQ=255)                                tGCCGGAATAATGCCGCACGCACCGTTAGTCGGTGCCGTTa                                                                      <  1:556242/41‑1 (MQ=255)                                 gCCGGAATAATGCCGCACGCACCGTTAGTCGGTGCCGTTACCACGCGCCCGCCAGCTGCGtt                                                <  1:1173057/62‑1 (MQ=255)                                 gCCGGAATAATGCCGCACGCACCGTTAGTCGGTGCCGTTACCACGCGCCCGCCAGCTGCGTTTTCTTCAc                                        >  1:3455443/1‑70 (MQ=255)                                  ccGGAATAATGCCGCACGCACCGTTAGTCGGTGCCGTTa                                                                      <  1:2388891/39‑1 (MQ=255)                                   cGGAATAATGCCGCACGCACCGTTAGTCGGTGCCGTTACCACGCGCCCGCCAGCTGCGTTTTCTTCACTaa                                     <  1:384785/71‑1 (MQ=255)                                    ggAATAATGCCGCACGCACCGTTAGTCGGTGCCGTTACCACGCGCCCGCCAGCTGCGTTTTCTTCACTAAc                                    >  1:3191954/1‑71 (MQ=255)                                     ggataatCCCGCACGCACCGTTAGTCGGCGCAGTCACCACGCGACCGCCAGCAGCGtt                                                <  1:5387372/56‑1 (MQ=255)                                     gAATAATGCCGCACGCACCGTTAGTCGGTGCCGTTACCACGCGCCCGCCAGCTGCGTTTTCTTCACTaa                                     >  1:4677484/1‑69 (MQ=255)                                       ataatGCCGCACGCACCGTTAGTCGGTGCCGTTACCACGCGCCCGCCAGCTGCGTTTTCTTCACTAACCGc                                 >  1:4308065/1‑71 (MQ=255)                                       ataatGCCGCACGCACCGTTAGTCGGTGCCGTTACCACGCGCCCGCCAGCTGCGTTTTCTTCACTAACCGc                                 >  1:5546509/1‑71 (MQ=255)                                         aatGCCGCACGCACCGTTAGTCGGTGCCGTTACCACGCGCCCGCCAGCTGCGTTttct                                            <  1:4361282/58‑1 (MQ=255)                                         aatGCCGCACGCACCGTTAGTCGGTGCCGTTACCACGCGCCCGCCAGCTGCGTTttct                                            <  1:2004588/58‑1 (MQ=255)                                         aatGCCGCACGCACCGTTAGTCGGTGCCGTTACCACGCGCCCGCCAGCTGCGTTttct                                            <  1:604785/58‑1 (MQ=255)                                           tcccGCACGCACCGTTAGTCGGCGCAGTCACCACGCGACCGCCAGCAGCGTTCTCTTCGTTCACTGCCAgt                             >  1:1552554/3‑70 (MQ=255)                                              cgcacgcaCCGTTAGTCGGTGCCGTTACCACGCGCCCGCCAGCTGCGTTTTCTTCACTAACCGCCcgcgcg                          >  1:1346859/1‑71 (MQ=255)                                               gcacgcacCGTTAGTCGGTGCCGTTACCACGCGCCCGCCAGCTGCGTTTTCTTCACTAACCGCCAGcgcg                          >  1:3312505/1‑70 (MQ=255)                                               gcacgcacCGTTAGTCGGCGCAGTCACCACGCGACCGCCAGCAGCGTTCTCTTCGTTCACTGCCAgtgca                          >  1:2334490/1‑69 (MQ=255)                                               gcacgcacCGTTAGTCGGCGCAGTCACCACGCGACCGCCAGCAGCGTTCTCTTCGTTCACTGCCAgtgca                          <  1:1044766/70‑2 (MQ=255)                                                cacgcacCGTTAGTCGGCGCAGTCACCACGCGACCGCCAGCAGCGTTCTCTTCGTTCACTGCCAgtgcaaa                        >  1:620203/1‑68 (MQ=255)                                                  cgcacCGTTAGTCGGTGCCGTTACCACGCGCCCGCCAGCTGCGTTTTCTTCACTAACCGCCAGCGCGTACa                      >  1:5409053/1‑71 (MQ=255)                                                  cgcacCGTTAGTCGGTGCCGTTACCACGCGCCCGCCAGCTGCGTTTTCTTCACTAACCGCCAGCGCGTACa                      >  1:2567839/1‑71 (MQ=255)                                                       cGTTAGTCGGTGCCGTTACCACGCGCCCGCCAGCTGCGTTTTCTTCACTAACCGCCAGCGCGTACATGTTg                 <  1:2775499/71‑1 (MQ=255)                                                       cGTTAGTCGGCGCAGTCACCACGCGACCGCCAGCAGCGTTCTCTTCGTTCACTGCCAGTGCAAACATGt                   >  1:5522383/1‑69 (MQ=255)                                                        gTTAGTCGGTGCCGTTACCACGCGCCCGCCAGCTGCGTTTTCTTCACTAACCGCCAgc                             >  1:2495112/1‑58 (MQ=255)                                                        gTTAGTCGGTGCCGTTACCACGCGCCCGCCAGCTGCGTTTTCTTCACTAACCGCCAGcgcg                          <  1:5421752/61‑1 (MQ=255)                                                        gTTAGTCGGTGCCGTTACCACGCGCCCGCCAGCTGCGTTTTCTTCACTAACCGCCAGCGCGTACATGTTGa                >  1:3041413/1‑71 (MQ=255)                                                        gTTAGTCGGTGCCGTTACCACGCGCCCGCCAGCTGCGTTTTCTTCACTAACCGCCAGCGCGTACATGTTGa                >  1:925203/1‑71 (MQ=255)                                                         ttAGTCGGTGCCGTTACCACGCGCCCGCCAGCTGCGTTTTCTTCACTAACCGCCAGCGCGTACATGTTGAt               >  1:5399221/1‑71 (MQ=255)                                                         ttAGTCGGTGCCGTTACCACGCGCCCGCCAGCTGCGTTTTCTTCACTAACCGCCAGCGCGTACATGTTGAt               >  1:466892/1‑71 (MQ=255)                                                          tAGTCGGTGCCGTTACCACGCGCCCGCCAGCTGCGTTTTCTTCACTAACCGCCAGCGCGTACATGTTGa                <  1:3601001/69‑1 (MQ=255)                                                           aGTCGGCGCAGTCACCACGCGACCGCCAGCAGCGTTCTCTTCGTTCACTGCCAGTGCAAACATGTTGATc              >  1:3908176/1‑70 (MQ=255)                                                             tCGGTGCCGTTACCACGCGCCCGCCAGCTGCGTTttctt                                           <  1:5136975/39‑1 (MQ=255)                                                               ggTGCCGTTACCACGCGCCCGCCAGCTGCGTTTTCTTCACTAACCGCCAGCGCGTACATGTTGATCCAGTc         >  1:2580445/1‑71 (MQ=255)                                                               ggTGCCGTTACCACGCGCCCGCCAGCTGCGTTTTCTTCACTAACCGCCAGCGCGTACATGTTGATCCAGTc         >  1:749459/1‑71 (MQ=255)                                                                   ccGTTACCACGCGCCCGCCAGCTGCGTTTTCTTCACTAACCGCCAGCGCGTACATGTTGATCCAGTCGAt      >  1:2698284/1‑70 (MQ=255)                                                                     gTTACCACGCGCCCGCCAGCTGCGTTTTCTTCACTAACCGCCAGCGCGTACATGTTGATCCAGTCGa       >  1:40654/1‑67 (MQ=255)                                                                      ttACCACGCGCCCGCCAGCTGCGTTTTCTTCACTAACCGCCAGCGCGTACATGTTGATCCAGTCGATGACa  <  1:977343/71‑1 (MQ=255)                                                                       |                                                                      GACGACGGAACTTATCGTAATAAGCCAGTACTGCCGGAATAATGCCGCACGCACCGTTAGTCGGTGCCGTTACCACGCGCCCGCCAGCTGCGTTTTCTTCACTAACCGCCAGCGCGTACATGTTGATCCAGTCGATGACA  >  minE/2093936‑2094075

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