Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	887,945	C→T	46.2%	E373K (GAG→AAG)	yciW ←	predicted oxidoreductase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	887,945	0	C	T	46.2%	13.2 / 79.8	63	E373K (GAG→AAG)	yciW	predicted oxidoreductase
Reads supporting (aligned to +/- strand):  ref base C (20/13);  new base T (14/15);  total (35/28)
Fisher's exact test for biased strand distribution p-value = 4.44e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.72e-01

TTTTGGTGGCCATTAGAGCGGCTGACAGTTTTACGCGAATCTGTCTGACGCGGCAAGGTTAATATGTCTCACCCAACGCGATTTTTAAGCGATTTATCCAGCCACACAGCCCGCTCCATACCAGCAGATTAATAATC  >  minE/887876‑888012                                                                      |                                                                    ttttGGTGGCCATTAGAGCGGCTGACAGTTTTACGCGAATCTGTCTGACGCGGCAAGGTTAATATGTCTCa                                                                    >  1:3508053/1‑71 (MQ=255)      gtggCCATTAGAGCGGCTGACAGTTTTACGCGAATCTGTCTGACGCGGCAAGGTTAATATGTCTTACCCaa                                                               <  1:4630952/71‑1 (MQ=255)      gtggCCATTAGAGCGGCTGACAGTTTTACGCGAATCTGTCTGACGCGGCAAGGTTAATATGTCTTACCCa                                                                <  1:438418/70‑1 (MQ=255)         gCCATTAGAGCGGCTGACAGTTTTACGCGAATCTGTCTGACGCGGCAAGGTTAATATGTCTCACCCAAcgc                                                            <  1:4651169/71‑1 (MQ=255)         gCCATTAGAGCGGCTGACAGTTTTACGCGAATCTGTCTGACGCGGCAAGGTTAATATGTCTCACCCAAcgc                                                            <  1:1335718/71‑1 (MQ=255)              tAGAGCGGCTGACAGTTTTACGCGAATCTGTCTGACGCGGCAAGGTTAATATGTCTTACCCAACGCGAttt                                                       >  1:4660980/1‑71 (MQ=255)              tAGAGCGGCTGACAGTTTTACGCGAATCTGTCTGACGCGGCAAGGTTAATATGTCTCACCCAACGCGAttt                                                       >  1:2050289/1‑71 (MQ=255)              tAGAGCGGCTGACAGTTTTACGCGAATCTGTCTGACGCGGCAAGGTTAATATGTCTCACCCAACGCGAttt                                                       >  1:2200646/1‑71 (MQ=255)               agagCGGCTGACAGTTTTACGCGAATCTGTCTGACGCGGCAAGGTTAATATGTCTCACCCAACGCGAtttt                                                      >  1:2963786/1‑71 (MQ=255)               agagCGGCTGACAGTTTTACGCGAATCTGTCTGACGCGGCAAGGTTAATATGTCTCACCCAACGCGAtttt                                                      >  1:3788940/1‑71 (MQ=255)                gagCGGCTGACAGTTTTACGCGAATCTGTCTGACGCGGCAAGGTTAATATGTCTCACCCAACGCGAttttt                                                     <  1:1187614/71‑1 (MQ=255)                   cGGCTGACAGTTTTACGCGAATCTGTCTGACGCGGCAAGGTTAATATGTCTCACCCAACGCGATTTTTAAg                                                  >  1:2130468/1‑71 (MQ=255)                     gCTGACAGTTTTACGCGAATCTGTCTGACGCGGCAAGGTTAATATGTCTTACCCAACGCGATTTTTAAGCg                                                <  1:124988/71‑1 (MQ=255)                     gCTGACAGTTTTACGCGAATCTGTCTGACGCGGCAAGGTTAATATGTCTCACCCAACGCGATTTTTAAGc                                                 >  1:261562/1‑70 (MQ=255)                         aCAGTTTTACGCGAATCTGTCTGACGCGGCAAGGTTAATATGTCTTACCCAACGCGATTTTTAAGCGAtt                                             >  1:194227/1‑70 (MQ=255)                         aCAGTTTTACGCGAATCTGTCTGACGCGGCAAGGTTAATATGTCTCACCCAACGCGATTTTTAAGCGAtt                                             <  1:365835/70‑1 (MQ=255)                          cAGTTTTACGCGAATCTGTCTGACGCGGCAAGGTTAATATGTCTTACCCAACGCGATTTTTAAGCGATTTa                                           <  1:2383645/71‑1 (MQ=255)                           aGTTTTACGCGAATCTGTCTGACGCGGCAAGGTTAATATGTCTTACCCAACGCGATTTTTAAGCGATTTAt                                          >  1:1280302/1‑71 (MQ=255)                             ttttACGCGAATCTGTCTGACGCGGCAAGGTTAATATGTCTGACCCAACGCGATTTTTAAGCGAt                                              >  1:1667967/1‑65 (MQ=255)                             ttttACGCGAATCTGTCTGACGCGGCAAGGTTAATATGTCTCACCCAACGCGATTTTTAAGCGATTTATc                                         >  1:3716983/1‑70 (MQ=255)                              tttACGCGAATCTGTCTGACGCGGCAAGGTTAATATGTCTTACCCAACGCGATTTTTAAGCGa                                               >  1:4313179/1‑63 (MQ=255)                               ttACGCGAATCTGTCTGACGCGGCAAGGTTAATATGTCTTACCCAACGCGATTTTTAAGCGATTTATCCa                                       <  1:3710660/70‑1 (MQ=255)                               ttACGCGAATCTGTCTGACGCGGCAAGGTTAATATGTCTTACCCAACGCGATTTTTAAGCGATTTATCCAg                                      >  1:3105316/1‑71 (MQ=255)                               ttACGCGAATCTGTCTGACGCGGCAAGGTTAATATGTCTCACCCAACGCGATTTTTAAGCGAt                                              <  1:3757538/63‑1 (MQ=255)                                tACGCGAATCTGTCTGACGCGGCAAGGTTAATATGTCTTACCCAACGCGATTTTTAAg                                                  >  1:4418720/1‑58 (MQ=255)                                  cgcgAATCTGTCTGACGCGGCAAGTTTAATATGTCTTACCCAACGAGATTTTTAAGCGATTTATcc                                        >  1:2108165/1‑66 (MQ=255)                                  cgcgAATCTGTCTGACGCGGCAAGGTTAATATGTCTCACCCAACGCGAttttt                                                     >  1:3169484/1‑53 (MQ=255)                                  cgcgAATCTGTCTGACGCGGCAAGGTTAATATGTCTCACCCAACGCGATTTTTAAGCGATTTATCCAGCca                                   <  1:4239221/71‑1 (MQ=255)                                    cgAATCTGTCTGACGCGGCAAGGTTAATATGTCTCACCCAACGCGATTTTTAAGCg                                                <  1:3964753/56‑1 (MQ=255)                                    cgAATCTGTCTGACGCGGCAAGGTTAATATGTCTCACCCAACGCGATTTTTAAGCGATTTATCCAGCcaca                                 <  1:288682/71‑1 (MQ=255)                                     gAATCTGTCTGACGCGGCAAGGTTAATATGTCTCACCCAACGCGATTTTTAAGCGATTTATc                                         >  1:1576102/1‑62 (MQ=255)                                     gAATCTGTCTGACGCGGCAAGGTTAATATGTCTCACCCAACGCGATTTTTAAGCGATTTATCCAGCcaca                                 >  1:1764579/1‑70 (MQ=255)                                      aaTCTGTCTGACGCGGCAAGGTTAATATGTCTCACCCAACGCGATTTTTAAGCGATTTATc                                         >  1:382202/1‑61 (MQ=255)                                         ctgtctgACGCGGCAAGGTTAATATGTCTTACCCaa                                                               >  1:105282/1‑36 (MQ=255)                                         ctgtctgACGCGGCAAGGTTAATATGTCTTACCCaa                                                               >  1:3224012/1‑36 (MQ=255)                                           gtctgACGCGGCAAGGTTAATATGTCTCACCCAACGCGATTTTTa                                                    >  1:1968177/1‑45 (MQ=255)                                           gtctgACGCGGCAAGGTTAATATGTCTCACCCAACGCGATTTTTa                                                    >  1:4137876/1‑45 (MQ=255)                                            tctgACGCGGCAAGGTTAATATGTCTTACCCAACGCGATTTTTAAGCGATTTATCCAGCCACACAGcc                            <  1:2748106/68‑1 (MQ=255)                                            tctgACGCGGCAAGGTTAATATGTCTTACCCAACGCGATTTTTAAGCGATTTATCCAGCCACACAGCCCGc                         >  1:3666520/1‑71 (MQ=255)                                            tctgACGCGGCAAGGTTAATATGTCTCACCCAACGCGATTTTTAAGCGATTTATCCAGCCACACAGcc                            <  1:1899710/68‑1 (MQ=255)                                            tctgACGCGGCAAGGTTAATATGTCTCACCCAACGCGATTTTTAAGCGATTTATCCAGCCACACAGCCCg                          >  1:2133621/1‑70 (MQ=255)                                              tgACGCGGCAAGGTTAATATGTCTCACCCAACGCGATTTTTAAGCGATTTATCCAGCCACACAGCCCGCTc                       >  1:1893574/1‑71 (MQ=255)                                                aCGCGGCAAGGTTAATATGTCTTACCCAACGCGATTTTTAAGCGATTTATCCAGCCACACAGCCCGCTc                       <  1:3227601/69‑1 (MQ=255)                                                aCGCGGCAAGGTTAATATGTCTCACCCAACGCGATTTTTAAGCGATTTATCCAGCCACACAGCCCGCTc                       <  1:4100634/69‑1 (MQ=255)                                                       aaGGTTAATATGTCTCACCCAACGCGATTTTTAAGCGa                                               <  1:2342254/38‑1 (MQ=255)                                                         ggTTAATATGTCTTACCCAACGCGATTTTTAAGCGATTTATCCAGCCACACAGCCCGCTCCATACcagcag            <  1:653087/71‑1 (MQ=255)                                                          gTTAATATGTCTTACCCAACGCGATTTTTAAGCGATTTATCCAGCCACACAg                              >  1:3395621/1‑52 (MQ=255)                                                           ttAATATGTCTTACCCAACGCGATTTTTAAGCGATTTATCCAGCCACACAGCCCGCTCCATACCAGCAGAt          <  1:2992031/71‑1 (MQ=255)                                                           ttAATATGTCTTACCCAACGCGATTTTTAAGCGATTTATCCAGCCACACAGCCCGCTCCATACCAGCAGAt          <  1:597668/71‑1 (MQ=255)                                                           ttAATATGTCTTACCCAACGCGATTTTTAAGCGATTTATCCAGCCACACAGCCCGCTCCATACCAGCAGAt          <  1:1907405/71‑1 (MQ=255)                                                           ttAATATGTCTTACCCAACGCGATTTTTAAGCGATTTATCCAGCCACACAGCCCGCTCCATACCAGCAGAt          >  1:1061517/1‑71 (MQ=255)                                                            tAATATGTCTTACCCAACGCGATTTTTAAGCGATTTATCCAGCCACACAGCCCGCTcc                      >  1:3760677/1‑58 (MQ=255)                                                            tAATATGTCTTACCCAACGCGATTTTTAAGCGATTTATCCAGCCACACAGCCCGCTCCa                     >  1:2759757/1‑59 (MQ=255)                                                            tAATATGTCTCACCCAACGCGATTTTTAAGCGATTTATCCAGCCACACAGCCCGCTCCATACCAGCAGAt          >  1:2870514/1‑70 (MQ=255)                                                            tAATATGTCTCACCCAACGCGATTTTTAAGCGATTTATCCAGCCACACAGCCCGCTCCATACCAGCAGAt          >  1:490589/1‑70 (MQ=255)                                                              atatGTCTTACCCAACGCGATTTTTAAGCGATTTATCCAGCCACACAGCCCGCTCCATACCAGCAGATtaa       <  1:1118706/71‑1 (MQ=255)                                                              atatGTCTTACCCAACGCGATTTTTAAGCGATTTATCCAGCCACACAGCCCGCTCCATACCAGCAGATta        <  1:4986/70‑1 (MQ=255)                                                                 tGTCTTACCCAACGCGATTTTTAAGCGATTTATCCAGCCACACAGCCCGCTCCATACCAGCAGATtaata     <  1:3624682/70‑1 (MQ=255)                                                                 tGTCTTACCCAACGCGATTTTTAAGCGATTTATCCAGCCACACAGCCCGCTCCATACCAGCAGATtaata     <  1:3379701/70‑1 (MQ=255)                                                                 tGTCTCACCCAACGCGATTTTTAAGCGATTTATCCAGCcacaca                               <  1:1986861/44‑1 (MQ=255)                                                                 tGTCTCACCCAACGCGATTTTTAAGCGATTTATCCAGCCACACAGCCCGCTCCATACCAGCAGATtaata     >  1:3053859/1‑70 (MQ=255)                                                                 tGTCTCACCCAACGCGATTTTTAAGCGATTTATCCAGCCACACAGCCCGCTCCATACCAGCAGATtaata     >  1:1555172/1‑70 (MQ=255)                                                                 tGTCTCACCCAACGCGATTTTTAAGCGATTTATCCAGCCACACAGCCCGCTCCATACCAGCAGATtaa       <  1:2952184/68‑1 (MQ=255)                                                                   tctcACCCAACGCGATTTTTAAGCGATTTATCCAGCCACACAGCCCGCTCCATACCAGCAGATTAATAATc  <  1:1317276/71‑1 (MQ=255)                                                                      |                                                                    TTTTGGTGGCCATTAGAGCGGCTGACAGTTTTACGCGAATCTGTCTGACGCGGCAAGGTTAATATGTCTCACCCAACGCGATTTTTAAGCGATTTATCCAGCCACACAGCCCGCTCCATACCAGCAGATTAATAATC  >  minE/887876‑888012

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