Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	2,094,792	C→A	10.3%	V20V (GTG→GTT)	tdcG ←	L‑serine dehydratase 3

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	2,094,792	0	C	A	10.3%	240.7 / 8.6	77	V20V (GTG→GTT)	tdcG	L‑serine dehydratase 3
Reads supporting (aligned to +/- strand):  ref base C (26/43);  new base A (6/2);  total (32/45)
Fisher's exact test for biased strand distribution p-value = 6.06e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.53e-01

CGTCGCGGTTAATAAGCCGCTACTTTCCAGCCGATCAATAAAACTTTTTCCGGCATTCATTGGCCCCACGGTATGCGAACTGGAGGGACCAATCCCAATTTTGAAAATATCGAATGCACTAATCATATCCACACCCT  >  minE/2094726‑2094862                                                                   |                                                                       cgtcgCGGTTAATAAGCCGCTACTTTCCAGCCGATCAATAAAACTTTTTCCGGCATTCATTGGCCCCACgt                                                                    >  1:4200169/1‑70 (MQ=255)  gtcgCGGTTAATAAGCCGCTACTTTCCAGCCGATCAATAAAACTTTTTCCGGCATTCATTGGCCCCACgg                                                                    >  1:3269073/1‑70 (MQ=255)   tcgCGGTTAATAAGCCGCTACTTTCCAGCCGATCAATAAAACTTTTTCCGGCATTCATTGGCCCCACGGt                                                                   <  1:853704/70‑1 (MQ=255)    cgcgGTTAATAAGCCGCTACTTTCCAGCCGATCAATAAAACTTTTTCCGGCATTCATTGGCCCCACGGTAt                                                                 >  1:1860267/1‑71 (MQ=255)       ggTTAATAAGCCGCTACTTTCCAGCCGATCAATAAAACTTTTTCCGGCATTCATTGGCCCCACGGt                                                                   >  1:696469/1‑66 (MQ=255)       ggTTAATAAGCCGCTACTTTCCAGCCGATCAATAAAACTTTTTCCGGCATTCATTGGCCCCACGGTATGCg                                                              <  1:1441762/71‑1 (MQ=255)        gTTAATAAGCCGCTACTTTCCAGCCGATCAATAAAACTTTTTCCGGCATTCATTGGCCCCACGGTATGCGa                                                             <  1:699096/71‑1 (MQ=255)            ataaGCCGCTACTTTCCAGCCGATCAATAAAACTTTTTCCGGCATTCATTGGCCCCACGGTATGCGAACTg                                                         <  1:1327468/71‑1 (MQ=255)              aaGCCGCTACTTTCCAGCCGATCAATAAAACTTTTTCCGGCATTCATTGGCCCCACGGTATGCGAACTg                                                         >  1:4823722/1‑69 (MQ=255)               aGCCGCTACTTTCCAGCCGATCAATAAAACTTTTTCCGGCATTCATTGGCCCCACGGTATGCGAACTGGa                                                       <  1:3861389/70‑1 (MQ=255)               aGCCGCTACTTTCCAGCCGATCAATAAAACTTTTTCCGGCATTCATTGGCCCCACGGTATGCGAACTGGa                                                       <  1:1862327/70‑1 (MQ=255)               aGCCGCTACTTT‑CAGCCGATCAATAAAACTTTTTCCGGCATTCATTGGCCCCa                                                                       >  1:70190/1‑53 (MQ=255)                       cTTTCCAGCCGATCAATAAAACTTTTTCCGGCATTCATTGGCCCCACGGTATGCGAACTGGAGGGACCaa                                               <  1:4652208/70‑1 (MQ=255)                       cTTTCCAGCCGATCAATAAAACTTTTTCCGGCATTCATTGGCCCCACGGTATGCGAACTGGAGGGACCAAt                                              <  1:4763655/71‑1 (MQ=255)                       cTTTCCAGCCGATCAATAAAACTTTTTCCGGCATTCATTGGCCCCACGGTATGCGAACTGGAGGGACCAAt                                              >  1:5313544/1‑71 (MQ=255)                         ttCCAGCCGATCAATAAAACTTTTTCCGGCATTCATTGGCCCCACGGTATGCGAACTGGAGGGACCAATcc                                            >  1:5114875/1‑71 (MQ=255)                          tCCAGCCGATCAATAAAACTTTTTCCGGCATTCATTGGCCCCACGGTATGCGAACTGGAGGGACCAATcc                                            <  1:1358814/70‑1 (MQ=255)                          tCCAGCCGATCAATAAAACTTTTTCCGGCATTCATTGGCCCCACGGTATGCGAACTGGAGGGACCAATcc                                            >  1:5202564/1‑70 (MQ=255)                          tCCAGCCGATCAATAAAACTTTTTCCGGCATTCATTGGCCCCACGGTATGCGAACTGGAGGAACCAATcc                                            >  1:4056668/1‑70 (MQ=255)                          tCCAGCCGATCAATAAAACTTTTTCCGGCATTCATTGGCCACACGGTATGCGAACTGGAGGGACCAATcc                                            >  1:491368/1‑70 (MQ=255)                           ccAGCCGATCAATAAAACTTTTTCCGGCATTCATTGGCCCCACGGTATGCGAACt                                                          <  1:2402142/55‑1 (MQ=255)                            cAGCCGATCAATAAAACTTTTTCCGGCATTCATTGGCCCCACGGTATGc                                                               >  1:5368591/1‑49 (MQ=255)                                  aTCAATAAAACTTTTTCCGGCATTCATTGGCCCCACGGTATGCg                                                              <  1:2889229/44‑1 (MQ=255)                                   tCAATAAAACTTTTTCCGGCATTCATTGGCCCCACGGTATGCGAACTGGAGGGACCAATCCCAATTTTGaa                                  <  1:4121931/71‑1 (MQ=255)                                   tCAATAAAACTTTTTCCGGCATTCATTGGCCCCACGGTATGCGAACTGGAGGGACCAATCCCAATTTTGa                                   >  1:3134508/1‑70 (MQ=255)                                     aaTAAAACTTTTTCCGGCATTCATTGGCCCCACGGTATGCGAACTGGAGGGACCAATCCCAATTTTGaaaa                                <  1:4592910/71‑1 (MQ=255)                                         aaaCTTTTTCCGGCATTCATTGGCCCCACGGTATGCg                                                              <  1:3587636/37‑1 (MQ=255)                                         aaaCTTTTTCCGGCATTCATTGGCCCCACGGTATGCg                                                              <  1:5395617/37‑1 (MQ=255)                                         aaaCTTTTTCCGGCATTCATTGGCCCCACGGTATGCGAACTGGa                                                       >  1:1471810/1‑44 (MQ=255)                                         aaaCTTTTTCCGGCATTCATTGGCCCCACGGTATGCGAACTGGAGGGACCAATCCCAATTTTGAAAATATc                            <  1:5076281/71‑1 (MQ=255)                                          aaCTTTTTCCGGCATTCATTGGCCCCACGGTATGCGAACTGGAGGGACCAATCCCAATTTTGAAAATATCg                           <  1:1802592/71‑1 (MQ=255)                                          aaCTTTTTCCGGCATTCATTGGCCCCACGGTATGCGAACTGGAGGGACCAATCCCAATTTTGAAAATATCg                           <  1:1765634/71‑1 (MQ=255)                                          aaCTTTTTCCGGCATTCATTGGCCCCACGGTATGCGAACTGGAGGGACCAATCCCAATTTTGAAAATATCg                           <  1:3942990/71‑1 (MQ=255)                                           aCTTTTTCCGGCATTCATTGGCCCCACGGTATGCGAACTGGAGGGACCAATCCCAATTTTGAAAATATCg                           <  1:4506995/70‑1 (MQ=255)                                           aCTTTTTCCGGCATTCATTGGCCCCACGGTATGCGAACTGGAGGGACCAATCCCAATTTTGAAAATATCGa                          <  1:218863/71‑1 (MQ=255)                                             tttttCCGGCATTCATTGGCCCCACGGTATGCGAACTGGAGGGACCAAt                                              >  1:4217547/1‑49 (MQ=255)                                              ttttCCGGCATTCATTGGCCCCACGGTATGCGAACTGGAGGGACCAATCCCAATTTTGAAAATATCGAATg                       >  1:3337653/1‑71 (MQ=255)                                              ttttCCGGCATTCATTGGCCCCACGGTATGCGAACTGGAGGGACCAATCCCAATTTTGAAAATATCGAATg                       >  1:2724242/1‑71 (MQ=255)                                              ttttCCGGCATTCATTGGCCCCACGGTATGCGAACTGGAGGGACCAATCCCAATTTTGAAAATATCGAATg                       >  1:5202255/1‑71 (MQ=255)                                               tttCCGGCATTCATTGGCCCCACGGTATGCGAACTGGAGGGACCAATCCCAAt                                        >  1:3459377/1‑53 (MQ=255)                                               tttCCGGCATTCATTGGCCCCACGGTATGCGAACTGGAGGGACCAATCCCAAt                                        >  1:5331385/1‑53 (MQ=255)                                                   cGGCATTCATTGGCCCCACGGTATGCGAACTGGAGGGACCAATCCCAATTTTGAAAATATCGAATGCACt                   <  1:4448373/70‑1 (MQ=255)                                                   cGGCATTCATTGGCCCCACGGTATGCGAACTGGAGGGACCAATCCCAATTTTGAAAATATCGAATGCACt                   <  1:3990211/70‑1 (MQ=255)                                                    ggCATTCATTGGCCCCACGGTATGCGAACTGGAGGGAc                                                  >  1:249946/1‑38 (MQ=255)                                                    ggCATTCATTGGCCCCACGGTATGCGAACTGGAGGGACCAATCCCAATTTTGAAAATATCGAATGCACTaa                 <  1:3967628/71‑1 (MQ=255)                                                     gctttcattGGTCCAACGGTATGAGAACTGGAAGGGCCAATGCCGATTTTGAAAATATCGAATACGCTAAt                >  1:3441019/4‑71 (MQ=255)                                                      ctttcattGGTCCAACGGTATGAGAACTGGAAGGGCCAATGCCGATTTTGAAAATATCTAAta                       >  1:569464/3‑62 (MQ=255)                                                      ctttcattGGTCCAACGGTATGAGAACTGGAAGGGCCAATGCCGATTTTGAAAATATCGAAta                       >  1:4994076/3‑62 (MQ=255)                                                      ctttcattGGTCCAACGGTATGAGAACTGGAAGGGCAAATGCCGATTTTGAAAata                              >  1:2413606/3‑56 (MQ=255)                                                       tttcattGGTCCAACGGTATGAGAACTGGAAGGGCCAATGCCGATTTTGAAAATATCGAAta                       >  1:2722853/2‑61 (MQ=255)                                                       attcattGGCCCCACGGTATGCGAACTGGAGGGACCAATCCCAATTTTGAAAATATCGAATGCACTAATCa              >  1:1953523/1‑71 (MQ=255)                                                        ttcattGGCCCCACGGTATGCGAACTGGAGGGACCAATCCCAATTTTGaa                                  <  1:975789/50‑1 (MQ=255)                                                        ttcattGGCCCCACGGTATGCGAACTGGAGGGACCAATCCCAATTTTGaa                                  <  1:2567559/50‑1 (MQ=255)                                                        ttcattGGCCCCACGGTATGCGAACTGGAGGGACCAATCCCAATTTTGaa                                  <  1:3304151/50‑1 (MQ=255)                                                        ttcattGGCCCCACGGTATGCGAACTGGAGGGACCAATCCCAATTTTGAAAATATCg                           <  1:1480598/57‑1 (MQ=255)                                                        ttcattGGCCCCACGGTATGCGAACTGGAGGGACCAATCCCAATTTTGAAAATATCg                           <  1:4444206/57‑1 (MQ=255)                                                        ttcattGGCCCCACGGTATGCGAACTGGAGGGACCAATCCCAATTTTGAAAATATCGAATg                       <  1:4680636/61‑1 (MQ=255)                                                        ttcattGGCCCCACGGTATGCGAACTGGAGGGACCAATCCCAATTTTGAAAATATCGAATg                       <  1:4937353/61‑1 (MQ=255)                                                        ttcattGGCCCCACGGTATGCGAACTGGAGGGACCAATCCCAATTTTGAAAATATCGAATg                       <  1:4967718/61‑1 (MQ=255)                                                        ttcattGGCCCCACGGTATGCGAACTGGAGGGACCAATCCCAATTTTGAAAATATCGAATg                       <  1:4988320/61‑1 (MQ=255)                                                        ttcattGGCCCCACGGTATGCGAACTGGAGGGACCAATCCCAATTTTGAAAATATCGAATg                       <  1:2780132/61‑1 (MQ=255)                                                        ttcattGGCCCCACGGTATGCGAACTGGAGGGACCAATCCCAATTTTGAAAATATCGAATg                       <  1:5178284/61‑1 (MQ=255)                                                        ttcattGGCCCCACGGTATGCGAACTGGAGGGACCAATCCCAATTTTGAAAATATCGAATGCACTAATCat             <  1:1549673/71‑1 (MQ=255)                                                         tcattGGCCCCACGGTATGCGAACTGGAGGGACCAATCCCAATTTTGAAAATATc                            <  1:201504/55‑1 (MQ=255)                                                         tcattGGCCCCACGGTATGCGAACTGGAGGGACCAATCCCAATTTTGAAAATATCGAATGCACTAATCata            >  1:4600281/1‑71 (MQ=255)                                                         tcattGGCCCCACGGTATGCGAACTGGAGGGACCAATCCCAATTTTGAAAATATCGAATGCACTAATCata            >  1:185753/1‑71 (MQ=255)                                                          cattGGTCCAACGGTATGAGAACTGGAAGGGCCAATGCCGATTTTGaaaa                                >  1:2021963/1‑50 (MQ=255)                                                           attGGCCCCACGGTATGCGAACTGGAGGGACCAATCCCAATTTTGAAAATATCGAATGCACTAATCatat           <  1:3790774/70‑1 (MQ=255)                                                           attGGCCCCACGGTATGCGAACTGGAGGGACCAATCCCAATTTTGAAAATATCGAATGCACTAATCatat           <  1:3673656/70‑1 (MQ=255)                                                           attGGCCCCACGGTATGCGAACTGGAGGGACCAATCCCAATTTTGAAAATATCGAATGCACTAATCATATc          <  1:2866449/71‑1 (MQ=255)                                                           attGGCCCCACGGTATGCGAACTGGAGGGACCAATCCCAATTTTGAAAATATCGAATGCACTAATCATATc          <  1:3152344/71‑1 (MQ=255)                                                            ttGGTCCAACGGTATGAGAACTGGAAGGGCCAATGCCGATTTTGAAAATATCGAATACGCTAATCatagga         <  1:1323247/71‑4 (MQ=255)                                                            ttGGCCCCACGGTATGCGAACTGGAGGGACCAATCCCAATTTTGAAAATATCGAATGCACTAATCATATc          <  1:4339447/70‑1 (MQ=255)                                                             tGGTCCAACGGTATGAGAACTGGAAGGGCCAATGCCGATTTTGAAAATATCGAATACGCTAATCatagga         <  1:3516884/70‑4 (MQ=255)                                                              ggCCCCACGGTATGCGAACTGGAGGGACCAATCCCAATTTTGAAAatat                             >  1:335444/1‑49 (MQ=255)                                                              ggCCCCACGGTATGCGAACTGGAGGGACCAATCCCAATTTTGAAAat                               >  1:3137799/1‑47 (MQ=255)                                                              ggCCCCACGGTATGCGAACTGGAGGGACCAATCCCAATTTTGAAAATATCGAATGCACTAATCATATCca        <  1:1537781/70‑1 (MQ=255)                                                                  ccACGGTATGCGAACTGGAGGGACCAATCCCAATTTTGa                                   <  1:697124/39‑1 (MQ=255)                                                                  ccACGGTATGCGAACTGGAGGGACCAATCCCAATTTTGAAAATATCg                           <  1:4611387/47‑1 (MQ=255)                                                                  ccACGGTATGCGAACTGGAGGGACCAATCCCAATTTTGAAAATATCGAATGCACTAATCATATCCACAccc   >  1:5661671/1‑71 (MQ=255)                                                                  ccACGGTATGCGAACTGGAGGGACCAATCCCAATTTTGAAAATATCGAATGCACTAATCATATCCACAcc    >  1:4446386/1‑70 (MQ=255)                                                                   cACGGTATGCGAACTGGAGGGACCAATCCCAATTTTg                                    <  1:5688508/37‑1 (MQ=255)                                                                   cACGGTATGCGAACTGGAGGGACCAATCCCAATTTTg                                    <  1:4752548/37‑1 (MQ=255)                                                                   cACGGTATGCGAACTGGAGGGACCAATCCCAATTTTGAAAat                               >  1:5219107/1‑42 (MQ=255)                                                                   cACGGTATGCGAACTGGAGGGACCAATCCCAATTTTGAAAATATCGAATGCACTAATCATATCCACACCCt  >  1:3152767/1‑71 (MQ=255)                                                                   |                                                                       CGTCGCGGTTAATAAGCCGCTACTTTCCAGCCGATCAATAAAACTTTTTCCGGCATTCATTGGCCCCACGGTATGCGAACTGGAGGGACCAATCCCAATTTTGAAAATATCGAATGCACTAATCATATCCACACCCT  >  minE/2094726‑2094862

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑