Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	2,036,045	C→T	12.2%	C59C (TGC→TGT)	ygiC →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	2,036,045	0	C	T	12.2%	237.0 / 19.4	82	C59C (TGC→TGT)	ygiC	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (31/41);  new base T (3/7);  total (34/48)
Fisher's exact test for biased strand distribution p-value = 5.11e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.74e-01

CTTACTACAAGTTGACCCTCGCCCAGGTTGAAAAGCTGGAAGAAGTCACCGCCGAACTGCACCAGATGTGCCTGAAAGTGGTGGAAAAAGTGATCGCCAGCGATGAGCTGATGACCAAATTCCGCATTCCAAAA  >  minE/2035975‑2036108                                                                       |                                                                cTTACTACAAGTTGACCCTCGCCCAGGTTGAAAAGCTGGAAGAAGTCACCGCCGAACTGCACCAGATGTGc                                                                 <  1:2889722/71‑1 (MQ=255) cTTACTACAAGTTGACCCTCGCCCAGGTTGAAAAGCTGGAAGAAGTCACCGCCGAACTGCACCAGATGTGc                                                                 <  1:1596260/71‑1 (MQ=255) cTTACTACAAGTTGACCCTCGCCCAGGTTGAAAAGCTGGAAGAAGTCACCGCCGAACTGCACCAGATGTGc                                                                 >  1:1710594/1‑71 (MQ=255)   tactacAAGTTGACCCTCGCCCAGGTTGAAAAGCTGGAAGAAGTCACCGCCGAACTGCACCAGATGTGCCt                                                               >  1:2768699/1‑71 (MQ=255)     ctacAAGTTGACCCTCGCCCAGGTTGAAAAGCTGGAAGAAGTCACCGCCGAACTGCACCAGATGTGCCTg                                                              >  1:3062603/1‑70 (MQ=255)     ctacAAGTTGACCCTCGCCAAGGTTGAAAAGCTGGAAGAAGTCACCGCCGAACTGCACCAGATGTGCCTg                                                              >  1:3175109/1‑70 (MQ=255)        cAAGTTGACCCTCGCCCAGGTTGAAAAGCTGGAAGAAGTCACCGCCGAACTGCACCAGATGTGCCTg                                                              >  1:2838755/1‑67 (MQ=255)        cAAGTTGACCCTCGCCCAGGTTGAAAAGCTGGAAGAAGTCACCGCCGAACTGCACCAGATGTGCCTGaa                                                            <  1:3120733/69‑1 (MQ=255)        cAAGTTGACCCTCGCCCAGGTTGAAAAGCTGGAAGAAGTCACCGCCGAACTGCACCAGATGTGCCTGa                                                             >  1:852201/1‑68 (MQ=255)        cAAGTTGACCCTCGCCCAGGTTGAAAAGCTGGAAGAAGTCACCGCCGAACTGCACCAGATGTGCCTGa                                                             >  1:3176440/1‑68 (MQ=255)        cAAGTTGACCCTCGCCCAGGTTGAAAAGCTGGAAGAAGTCACCGCCGAACTGCACCAGATGTGCCTGAAAg                                                          >  1:873682/1‑71 (MQ=255)        cAAGTTGACCCTCGCCCAGGTTGAAAAGCTGGAAGAAGTCACCGCCGAACTGCACCAGATGTGCCTGAAAg                                                          <  1:834943/71‑1 (MQ=255)        cAAGTTGACCCTCGCCCAGGTTGAAAAGCTGGAAGAAGTCACCGCCGAACTGCACCAGATGTGCCTGAAAg                                                          >  1:2064234/1‑71 (MQ=255)        cAAGTTGACCCTCGCCCAGGTTGAAAAGCTGGAAGAAGTCACCGCCGAACTGCACCAGATGTGCCTGAAAg                                                          >  1:783836/1‑71 (MQ=255)             tGACCCTCGCCCAGGTTGAAAAGCTGGAAGAAGTCACCGCCGAACTGCACCAGATGTGCCTGAAAgtggtg                                                     >  1:2896403/1‑71 (MQ=255)             tGACCCTCGCCCAGGTTGAAAAGCTGGAAGAAGTCACCGCCGAACTGCACCAGATGTGCCTGAAAgtggtg                                                     >  1:21583/1‑71 (MQ=255)                 ccTCGCCCAGGTTGAAAAGCTGGAAGAAGTCACCGCCGAACTGCACCAGATGTGCCTGAAAGTGGTGGaaa                                                 <  1:386474/71‑1 (MQ=255)                  cTCGCCCAGGTTGAAAAGCTGGAAGAAGTCACCGCCGAACTGCACCAGATGTGCCTGAAAGTGGTGGaaaa                                                <  1:989621/71‑1 (MQ=255)                  cTCGCCCAGGTTGAAAAGCTGGAAGAAGTCACCGCCGAACTGCACCAGATGTGCCTGAAAGTGGTGGaaa                                                 >  1:3076688/1‑70 (MQ=255)                    cGCCCAGGTTGAAAAGCTGGAAGAAGTCACCGCCGAACTGCACCAGATGTGCCTGAAAGTGGTGGAAAAAg                                              >  1:2201115/1‑71 (MQ=255)                     gCCCAGGTTGAATAGCTGGAAGAAGTCACCGCCGAACTGCACCAGATGTGCCTGAAAGTGGTGGaaaaa                                               <  1:1118196/69‑1 (MQ=255)                      cccAGGTTGAAAAGCTGGAAGAAGTCACCGCCGAACTGCACCAGATGTGCCTGAAAGTGGTGGAAAAAGTg                                            <  1:1855799/71‑1 (MQ=255)                       ccAGGTTGAAAAGCTGGAAGAAGTCACCGCCGAACTGCACCAGATGTGCCTGAAAGTGGTGGAAAAAGTg                                            <  1:536568/70‑1 (MQ=255)                         aGGTTGAAAAGCTGGAAGAAGTCACCGCCGAACTGCACCAGATGTGCCTGAAAGTGCTGGaaaa                                                <  1:2371526/64‑1 (MQ=255)                          ggTTGAAAAGCTGGAAGAAGTCACCGCCGAACTGCACCAGATGTGCCTGAAAGTGGTGGAAAAAGTGATCg                                        >  1:2722511/1‑71 (MQ=255)                          ggTTGAAAAGCTGGAAGAAGTCACCGCCGAACTGCACCAGATGTGCCTGAAAGTGGTGGAAAAAGTGATCg                                        >  1:1501350/1‑71 (MQ=255)                          ggTTGAAAAGCTGGAAGAAGTCACCGCCGAACTGCACCAGATGTGCCTGAAAGTGGTGGAAAAAGTGATCg                                        <  1:2789628/71‑1 (MQ=255)                          ggTTGAAAAGCTGGAAGAAGTCACCGCCGAACTGCACCAGATGTGCCTGAAAGTGGTGGAAAAAGTGATCg                                        >  1:2009476/1‑71 (MQ=255)                           gTTGAAAAGCTGGAAGAAGTCACCGCCGAACTGCACCAGATGTGCCTGAAAGTGGTGGAAAAAGt                                             <  1:1087524/65‑1 (MQ=255)                           gTTGAAAAGCTGGAAGAAGTCACCGCCGAACTGCACCAGATGTGCCTGAAAGTGGTGGAAAAAGTGATCGc                                       >  1:915132/1‑71 (MQ=255)                            ttGAAAAGCTGGAAGAAGTCACCGCCGAACTGCACCAGATGTGCCTGaaa                                                           >  1:608976/1‑50 (MQ=255)                             tGAAAAGCTGGAAGAAGTCACCGCCGAACTGCACCAGATGTGCCTGAAAGTGGTGGAAAAAGTGATCGc                                       >  1:2333513/1‑69 (MQ=255)                              gAAAAGCTGGAAGAAGTCACCGCCGAACTGCACCAGATGTGTCTg                                                              >  1:660252/1‑45 (MQ=255)                              gAAAAGCTGGAAGAAGTCACCGCCGAACTGCACCAGATGTGTCTg                                                              >  1:454774/1‑45 (MQ=255)                               aaaaGCTGGAAGAAGTCACCGCCGAACTGCACCAGATGTGCCTGAAAGTGGTGGAAAAAGTGATCGCCAGc                                   <  1:297210/71‑1 (MQ=255)                               aaaaGCTGGAAGAAGTCACCGCCGAACTGCACCAGATGTGCCTGAAAGTGGTGGAAAAAGTGATCGCCAGc                                   <  1:2794959/71‑1 (MQ=255)                               aaaaGCTGGAAGAAGTCACCGCCGAACTGCACCAGATGTGCCTGAAAGTGGTGGAAAAAGTGATCGCCAGc                                   <  1:505462/71‑1 (MQ=255)                               aaaaGCTGGAAGAAGTCACCGCCGAACTGCACCAGATGTGCCTGAAAGTGGTGGAAAAAGTGATCGCCAGc                                   <  1:2744066/71‑1 (MQ=255)                                aaaGCTGGAAGAAGTCACCGCCGAACTGCACCAGATGTGCCTGAAAGTGGTGGAAAAAGTGATCGCCAGCg                                  <  1:242703/71‑1 (MQ=255)                                   gCTGGAAGAAGTCACCGCCGAACTGCACCAGATGTGCCTGAAAGTGGTGGAAAAAGTGATCGCCAGCGATg                               <  1:919091/71‑1 (MQ=255)                                       gaagaaGTCACCGCCGAACTGCACCAGATGTGCCTGAAAGTGGTGGAAAAAGTGATCGCCAGCGATGAg                             <  1:2825414/69‑1 (MQ=255)                                       gaagaaGTCACCGCCGAACTGCACCAGATGTGCCTGAAAGTGGTGGAAAAAGTGATCGCCAGCGATGAg                             <  1:193380/69‑1 (MQ=255)                                         agaagTCACCGCCGAACTGCACCAGATGTGCCTGAAAGTGGTGGAAAAAGTGATCGCCAGCGATGAGCtga                         <  1:2119811/71‑1 (MQ=255)                                         agaagTCACCGCCGAACTGCACCAGATGTGCCTGAAAGTGGTGGAAAAAGTGATCGCCAGCGATGAGCtga                         <  1:345840/71‑1 (MQ=255)                                           aagTCACCGCCGAACTGCACCAGATGTGTCTGAAAGTGGTGGAAAAAGTGATCGCCAGCGAt                                <  1:2105380/62‑1 (MQ=255)                                           aagTCACCGCCGAACTGCACCAGATGTGTCTGAAAGTGGTGGAAAAAGTGATCGCCAGCGAt                                <  1:1546590/62‑1 (MQ=255)                                           aagTCACCGCCGAACTGCACCAGATGTGTCTGAAAGTGGTGGAAAAAGTGATCGCCAGCGAt                                <  1:1088287/62‑1 (MQ=255)                                           aagTCACCGCCGAACTGCACCAGATGTGCCTGAAAGTGGTGGaaa                                                 <  1:909396/45‑1 (MQ=255)                                           aagTCACCGCCGAACTGCACCAGATGTGCCTGAAAGTGGTGGAAAAAg                                              >  1:1828205/1‑48 (MQ=255)                                           aagTCACCGCCGAACTGCACCAGATGTGCCTGAAAGTGGTGGAAAAAGTGATCGCCAGc                                   <  1:2825701/59‑1 (MQ=255)                                           aagTCACCGCCGAACTGCACCAGATGTCTCTGAAAGTGGTGGAAAAAGTGATCGCCAGCGAt                                <  1:2987226/62‑1 (MQ=255)                                            agTCACCGCCGAACTGCACCAGATGTGCCTGAAAGTGGTGGAAAAAGTGATCGCCAGCGATGAGCtgatga                      >  1:3175345/1‑71 (MQ=255)                                            agTCACCGCCGAACTGCACCAGATGTGCCTGAAAGTGGTGGAAAAAGTGATCGCCAGCGATGAGCtgatga                      >  1:2983133/1‑71 (MQ=255)                                            agTCACCGCCGAACTGCACCAGATGTGCCTGAAAGTGGTGGAAAAAGTGATCGCCAGCGATGAGCtgatga                      >  1:2839656/1‑71 (MQ=255)                                             gTCACCGCCGAACTGCACCAGATGTGTCTGAAAGTGGTGGAAAAAGTGATCGCCAGCGATg                               >  1:1356519/1‑61 (MQ=255)                                               cACCGCCGAACTGCACCAGATGTGCCTGAAAGTGGTGGAAAAAGTGATCg                                        >  1:2028067/1‑50 (MQ=255)                                                aCCGCCGAACTGCACCAGATGTGCCTGAAAGTGGTGGAAAAAGTGATCGCCAGCGATGAGCTGATGACCa                   >  1:543541/1‑70 (MQ=255)                                                 ccgccgAACTGCACCAGATGTGCCTGAAAGTGGTGGAAAAAGTGATCGcc                                      <  1:3057794/50‑1 (MQ=255)                                                 ccgccgAACTGCACCAGATGTGCCTGAAAGTGGTGGAAAAAGTGATCGcc                                      <  1:593666/50‑1 (MQ=255)                                                   gccgAACTGCACCAGATGTGCCTGAAAGTGGTGGAAAAAGTGATCGCCAGc                                   <  1:2426885/51‑1 (MQ=255)                                                   gccgAACTGCACCAGATGTGCCTGAAAGTGGTGGAAAAAGTGATCGCCAGc                                   <  1:1007033/51‑1 (MQ=255)                                                   gccgAACTGCACCAGATGTGCCTGAAAGTGGTGGAAAAAGTGATCGCCAGCGATGAGCtga                         <  1:1095655/61‑1 (MQ=255)                                                   gccgAACTGCACCAGATGTGCCTGAAAGTGGTGGAAAAAGTGATCGCCAGCGATGAGCtga                         <  1:2853307/61‑1 (MQ=255)                                                      gAACTGCACCAGATGTGCCTGAAAGTGGTGGAAAAAGTGATCGCCAGCGATGAGCTGATGACCAAATTCCg            <  1:672404/71‑1 (MQ=255)                                                       aaCTGCACCAGATGTGCCTGAAAGTGGTGGAAAAAGTGATCGCCAGc                                   <  1:2314134/47‑1 (MQ=255)                                                        aCTGCACCAGATGTGTCTGAAAGTGGTGGAAAAAGTGATCGCCAGCGAt                                <  1:127335/49‑1 (MQ=255)                                                         cTGCACCAGATGTGTCTGAAAGTGGTGGAAAAAGTGATc                                         <  1:1167006/39‑1 (MQ=255)                                                         cTGCACCAGATGTGTCTGAAAGTGGTGGAAAAAGTGATc                                         <  1:1743188/39‑1 (MQ=255)                                                         cTGCACCAGATGTGCCTGAAAGTGGTGGAAAAAGTGATCGCCAGCGATGAGCTGATGACCaaa                 >  1:1029004/1‑63 (MQ=255)                                                         cTGCACCAGATGTGCCTGAAAGTGGTGGAAAAAGTGATCGCCAGCGATGAGCTGATGACCAAATTCCGc           <  1:1101019/69‑1 (MQ=255)                                                          tGCACCAGATGTGCCTGAAAGTGGTGGAAAAAGTg                                            <  1:1025195/35‑1 (MQ=255)                                                          tGCACCAGATGTGCCTGAAAGTGGTGGAAAAAGTg                                            <  1:367641/35‑1 (MQ=255)                                                          tGCACCAGATGTGCCTGAAAGTGGTGGAAAAAGTg                                            <  1:1614280/35‑1 (MQ=255)                                                            cACCAGATGTGCCTGAAAGTGGTGGAAAAAGTGATCGCCAGCGATGAGCtga                         >  1:3108294/1‑52 (MQ=255)                                                            cACCAGATGTGCCTGAAAGTGGTGGAAAAAGTGATCGCCAGCGATGAGCtga                         >  1:1421083/1‑52 (MQ=255)                                                              ccAGATGTGCCTGAAAGTGGTGGAAAAAGTGATCGCCAGCGATGAGCTGATGACCAAATTCCGCATTCCaa    <  1:972774/71‑1 (MQ=255)                                                              ccAGATGTGCCTGAAAGTGGTGGAAAAAGTGATCGCCAGCGATGAGCTGATGACCAAATTCCGCATTCCaa    <  1:3122196/71‑1 (MQ=255)                                                              ccAGATGTGCCTGAAAGTGGTGGAAAAAGTGATCGCCAGCGATGAGCTGATGACCAAATTCCGCATTCCaa    <  1:678553/71‑1 (MQ=255)                                                                  aTGTGCCTGAAAGTGGTGGAAAAAGTGATCGCCAGc                                   <  1:2053101/36‑1 (MQ=255)                                                                   tgtgCCTGAAAGTGGTGGAAAAAGTGATCGCCAGCGATGAGCTGATGACCaaa                 >  1:723090/1‑53 (MQ=255)                                                                    gtgCCTGAAAGTGGTGGAAAAAGTGATCGCCAGCGATGAGCtg                          <  1:633511/43‑1 (MQ=255)                                                                    gtgCCTGAAAGTGGTGGAAAAAGTGATCGCCAGCGATGAGCtg                          <  1:1398067/43‑1 (MQ=255)                                                                    gtgCCTGAAAGTGGTGGAAAAAGTGATCGCCAGCGATGAGCtg                          <  1:1544782/43‑1 (MQ=255)                                                                    gtgCCTGAAAGTGGTGGAAAAAGTGATCGCCAGCGATGAGCTGATGACCAAATTCCGCATTCCaaaa  >  1:1228375/1‑67 (MQ=255)                                                                    gtgCCTGAAAGTGGTGGAAAAAGTGATCGCCAGCGATGAGCTGATGACCAAATTCCGCATTCCaaaa  >  1:1932661/1‑67 (MQ=255)                                                                       |                                                                CTTACTACAAGTTGACCCTCGCCCAGGTTGAAAAGCTGGAAGAAGTCACCGCCGAACTGCACCAGATGTGCCTGAAAGTGGTGGAAAAAGTGATCGCCAGCGATGAGCTGATGACCAAATTCCGCATTCCAAAA  >  minE/2035975‑2036108

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