Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	2,094,003	C→A	23.1%	T283T (ACG→ACT)	tdcG ←	L‑serine dehydratase 3

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	2,094,003	0	C	A	23.1%	194.6 / 34.5	86	T283T (ACG→ACT)	tdcG	L‑serine dehydratase 3
Reads supporting (aligned to +/- strand):  ref base C (41/25);  new base A (16/4);  total (57/29)
Fisher's exact test for biased strand distribution p-value = 1.81e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CCGGACGACGGAACTTATCGTAATAAGCCAGTACTGCCGGAATAATGCCGCACGCACCGTTAGTCGGTGCCGTTACCACGCGCCCGCCAGCTGCGTTTTCTTCACTAACCGCCAGCGCGTACATGTTGATCCAGTCGAT  >  minE/2093933‑2094071                                                                       |                                                                     ccGGACGACGGAACTTATCGTAATAAGCCAGTACTGCCGGAATAATGCCGCACGCACCGTTAGTCGGTGcc                                                                      <  1:5308122/71‑1 (MQ=255)   ggACGACGGAACTTATCGTAATAAGCCAGTACTGCCGGAATAATGCCGCACGCACCGTTAGTCGGTGCCGt                                                                    >  1:936761/1‑71 (MQ=255)   ggACGACGGAACTTATCGTAATAAGCCAGTACTGCCGGAATAATGCCGCACGCACCGTTAGTCGGTGCCGt                                                                    >  1:815046/1‑71 (MQ=255)   ggACGACGGAACTTATCGTAATAAGCCAGTACTGCCGGAATAATGCCGCACGCACCGTTAGTCGGTGCCGt                                                                    >  1:1516816/1‑71 (MQ=255)   ggACGACGGAACTTATCGTAATAAGCCAGTACTGCCGGAATAATGCCGCACGCACCGTTAGTCGGTGCCGt                                                                    <  1:351582/71‑1 (MQ=255)    gacgacGGAACTTATCGTAATAAGCCAGTACTGCCGGAATAATGCCGCACGCACCGTTAGTCGGTGCCGtt                                                                   <  1:3766914/71‑1 (MQ=255)    gacgacGGAACTTATCGTAATAAGCCAGTACTGCCGGAATAATGCCGCACGCACCGTTAGTCGGTGCCGt                                                                    >  1:1967858/1‑70 (MQ=255)     acgacgGAACTTATCGTAATAAGCCAGTACTGCCGGAATAATGCCGCACGCACCGTTAGTCGGTGCCGTTa                                                                  <  1:4912529/71‑1 (MQ=255)       gacgGAACTTATCGTAATAAGCCAGTACTGCCGGAATAATGCCGCACGCACCGTTAGTCGGTGCCGTTAcc                                                                <  1:5573537/71‑1 (MQ=255)       gacgGAACTTATCGTAATAAGCCAGTACTGCCGGAATAATGCCGCACGCACCGTTAGTCGGTGCCGTTAcc                                                                <  1:330598/71‑1 (MQ=255)          taaaCTTGTCGTAGTACGCCAGAACTGCCGGGATAATCCCGCACGCACCGTTAGTCGGCGCAGTCAc                                                                 >  1:5127424/3‑67 (MQ=255)          ggAACTTATCGTAATAAGCCAGTACTGCCGGAATAATGCCGCACGCACCGTTAGTCGGTGCCGTTACCAcg                                                             >  1:164002/1‑71 (MQ=255)             aCTTATCGTAATAAGCCAGTACTGCCGGAATAATGCCGCACGCACCGTTAGTCGGTGCCGTTACCACgcgc                                                          >  1:4843/1‑71 (MQ=255)             aCTTATCGTAATAAGCCAGTACTGCCGGAATAATGCCGCACGCACCGTTAGTCGGTGCCGTTACCACgcgc                                                          >  1:3708721/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:5502257/1‑66 (MQ=255)                tATCGTAATAAGCCAGTACTGCCGGAATAATGCCGCACGCACCGTTAGTCGGTGCCGTTACCAcgc                                                            >  1:3890246/1‑66 (MQ=255)                tATCGTAATAAGCCAGTACTGCCGGAATAATGCCGCACGCACCGTTAGTCGGTGCCGTTACCAcgc                                                            >  1:217407/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:2676583/1‑71 (MQ=255)                 aTCGTAATAAGCCAGTACTGCCGGAATAATGCCGCACGCACCGTTAGTCGGTGCCGTTACCACGcgcccgc                                                      >  1:3695590/1‑71 (MQ=255)                 aTCGTAATAAGCCAGTACTGCCGGAATAATGCCGCACGCACCGTTAGTCGGTGCCGTTACCACGcgcccgc                                                      >  1:4637198/1‑71 (MQ=255)                 aTCGTAATAAGACAGTACTGCCGGAATAATGCCGCACGCACCGTTAGTCGGTGCCGTTa                                                                  >  1:168245/1‑59 (MQ=255)                  tCGTAGTACGCCAGAACTGCCGGGATAATCCCGCACGCACCGTTAGTCGGCGCAGTcaccac                                                              >  1:1527158/1‑62 (MQ=255)                  tCGTAGTACGCCAGAACTGCCGGGATAATCCCGCACGCACCGTTAGTCGGCGCAGTcaccac                                                              >  1:3411282/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:465832/71‑1 (MQ=255)                           gCCAGTACTGCCGGAATAATGCCGCACGCACCGTTAGTCGGTGCCGTTACCACGCGCCCGCCAGCTGCGtt                                            <  1:5166382/71‑1 (MQ=255)                            ccAGAACTGCCGGGATAATCCCGCACGCACCGTTAGTCGGCGCAGTCACCACGCGACCGCCAGCAGCGttc                                           >  1:2620699/1‑70 (MQ=255)                            ccAGAACTGCCGGGATAATCCCGCACGCACCGTTAGTCGGCGCAGTCACCACGCGACCGCCAGCAGCGttc                                           <  1:2552742/71‑2 (MQ=255)                            ccAGAACTGCCGGGATAATCCCGCACGCACCGTTAGTCGGCGCAGTCACCACGCGACCGCCAGCAGCGttc                                           >  1:2181308/1‑70 (MQ=255)                             cAGTACTGCCGGAATAATGCCGCACGCACCGTTAGTCGGTGCCGTTACCACGCGCCCGCCAGCTGCGtttt                                          <  1:2560447/71‑1 (MQ=255)                             cAGTACTGCCGGAATAATGCCGCACGCACCGTTAGTCGGTGCCGTTACCACGCGCCCGCCAGCTGCGtttt                                          <  1:182662/71‑1 (MQ=255)                             cAGTACTGCCGGAATAATGCCGCACGCACCGTTAGTCGGTGCCGTTACCACGCGCCCGCCAGCTGCGtttt                                          >  1:1294979/1‑71 (MQ=255)                              aGTACTGCCGGAATAATGCCGCACGCACCGTTAGTCGGTGCCGTTACCACGCGCCCGCCAGCTGCGTTttc                                         >  1:58555/1‑71 (MQ=255)                              aGTACTGCCGGAATAATGCCGCACGCACCGTTAGTCGGTGCCGTTACCACGCGCCCGCCAGCTGCGTTttc                                         >  1:3581412/1‑71 (MQ=255)                                tACTGCCGGAATAATGCCGCACGCACCGTTAGTCGGTGCCGTTACCACGCGCCCGCCAGCTGCGTTttctt                                       <  1:3112079/71‑1 (MQ=255)                                tACTGCCGGAATAATGCCGCACGCACCGTTAGTCGGTGCCGTTACCACGCGCCCGCCAGCTGCGTTttctt                                       <  1:4390054/71‑1 (MQ=255)                                  cTGCCGGAATAATGCCGCACGCACCGTTAGTCGGTGCCGtt                                                                   >  1:3664791/1‑41 (MQ=255)                                  cTGCCGGAATAATGCCGCACGCACCGTTAGTCGGTGCCGtt                                                                   >  1:3668176/1‑41 (MQ=255)                                  cTGCCGGAATAATGCCGCACGCACCGTTAGTCGGTGCCGtt                                                                   >  1:3661615/1‑41 (MQ=255)                                  cTGCCGGAATAATGCCGCACGCACCGTTAGTCGGTGCCGTTACCACGCGCCCGCCAGCTGCGTTttctt                                       >  1:1208784/1‑69 (MQ=255)                                   tGCCGGAATAATGCCGCACGCACCGTTAGTCGGTGCCGTTa                                                                  <  1:556242/41‑1 (MQ=255)                                   tGCCGGAATAATGCCGCACGCACCGTTAGTCGGTGCCGTTa                                                                  <  1:3043430/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:5546509/1‑71 (MQ=255)                                          ataatGCCGCACGCACCGTTAGTCGGTGCCGTTACCACGCGCCCGCCAGCTGCGTTTTCTTCACTAACCGc                             >  1:4308065/1‑71 (MQ=255)                                            aatGCCGCACGCACCGTTAGTCGGTGCCGTTACCACGCGCCCGCCAGCTGCGTTttct                                        <  1:604785/58‑1 (MQ=255)                                            aatGCCGCACGCACCGTTAGTCGGTGCCGTTACCACGCGCCCGCCAGCTGCGTTttct                                        <  1:4361282/58‑1 (MQ=255)                                            aatGCCGCACGCACCGTTAGTCGGTGCCGTTACCACGCGCCCGCCAGCTGCGTTttct                                        <  1:2004588/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)                                                                       |                                                                     CCGGACGACGGAACTTATCGTAATAAGCCAGTACTGCCGGAATAATGCCGCACGCACCGTTAGTCGGTGCCGTTACCACGCGCCCGCCAGCTGCGTTTTCTTCACTAACCGCCAGCGCGTACATGTTGATCCAGTCGAT  >  minE/2093933‑2094071

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