Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	1,469,857	C→T	6.3%	G622D (GGC→GAC)	yfaL ←	adhesin

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	1,469,857	0	C	T	6.3%	205.5 / 3.1	63	G622D (GGC→GAC)	yfaL	adhesin
Reads supporting (aligned to +/- strand):  ref base C (23/36);  new base T (2/2);  total (25/38)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CACTGCGTGTCTACCCCCGCATCAACTGCCACTTCACCGTCGGCGCGCATTTCAATATCACGTCCGTGGCCTGCCGCATCTGAACGCGTGTCAACGTCGGAGGTAATCTCCATCACGCTTTCTTGTTGCTTATCGGTAAA  >  minE/1469788‑1469927                                                                      |                                                                       cACTGCGTGTCTACCCCCGCATCAACTGCCACTTCACCGTCGGCGCGCATTTCAATATCACGTCCGTGGcc                                                                       >  1:2399449/1‑71 (MQ=255)   cTGCGTGTCTACCCCCGCATCAACTGCCACTTCACCGTCGGCGCGCATTTCAATATCACGTCCGTGGCCTg                                                                     >  1:2298794/1‑71 (MQ=255)     gCGTGTCTACCCCCGCATCAACTGCCACTTCACCGTCGGCGCGCATTTCAATATCACGTCCGTGGCCTGc                                                                    >  1:4098735/1‑70 (MQ=255)        tgtCTACCCCCGCATCAACTGCCACTTCACCGTCGGCGCGCATTTCAATATCACGTCCGTGGCCTGCCGCa                                                                <  1:3452737/71‑1 (MQ=255)              cccccGCATCAACTGCCACTTCACCGTCGGCGCGCATTTCAATATCACGTCCGTGGCCTGCCGCATCTg                                                            <  1:3903102/69‑1 (MQ=255)              cccccGCATCAACTGCCACTTCACCGTCGGCGCGCATTTCAATATCACGTCCGTGGCCTGCCGCATCTGa                                                           <  1:1887824/70‑1 (MQ=255)               ccccGCATCAACTGCCACTTCACCGTCGGCGCGCATTTCAATATCACGTCCGTGGCCTGCCGCATCTGAAc                                                         <  1:475798/71‑1 (MQ=255)                cccGCATCAACTGCCACTTCACCGTCGGCGCGCATTTCAATATCACGTCCGTGGCCTGCCGCATCTGAAcg                                                        <  1:2446237/71‑1 (MQ=255)                cccGCATCAACTGCCACTTCACCGTCGGCGCGCATTTCAATATCACGTCCGTGGCCTGCCGCATCTGAAcg                                                        <  1:1794377/71‑1 (MQ=255)                 ccGCATCAACTGCCACTTCACCGTCGGCGCGCATTTCAATATCACGTCCGTGGCCTGCCGCATCTGAAcgc                                                       <  1:3787043/71‑1 (MQ=255)                  cGCATCAACTGCCACTTCACCGTCGGCGCGCATTTCAATATCACGTCCGTGGCCTGCCGCATCTGAAcgcg                                                      >  1:1245636/1‑71 (MQ=255)                   gCATCAACTGCCACTTCACCGTCGGCGCGCATTTCAATATCACGTCCGTGGCCTGCCGCATCTGAACGCgt                                                     <  1:3310235/71‑1 (MQ=255)                    cATCAACTGCCACTTCACCGTCGGCGCGCATTTCAATATCACGTCCGTGGCCTGcc                                                                   >  1:241997/1‑56 (MQ=255)                     aTCAACTGCCACTTCACCGTCGGCGCGCATTTCAATATCACGTCCGTGGCCTGCCGCATCTGaa                                                          >  1:3831383/1‑64 (MQ=255)                      tCAACTGCCACTTCACCGTCGGCGCGCATTTCAATATCACGTCCGTGGCCt                                                                      <  1:1962571/51‑1 (MQ=255)                      tCAACTGCCACTTCACCGTCGGCGCGCATTTCAATATCACGTCCGTGGCCTGCCGCATCTGAACGCGTGTc                                                  <  1:317619/71‑1 (MQ=255)                      tCAACTGCCACTTCACCGTCGGCGCGCATTTCAATATCACGTCCGTGGCCTGCCGCATCTGAACGCGTGTc                                                  <  1:204359/71‑1 (MQ=255)                      tCAACTGCCACTTCACCGTCGGCGCGCATTTCAATATCACGTCCGTGGCCTGCCGCATCTGAACGCGTGTc                                                  <  1:2282426/71‑1 (MQ=255)                      tCAACTGCCACTTCACCGTCGGCGCGCATTTCAATATCACGTCCGTGGCCTGCCGCATCTGAACGCGTGTc                                                  <  1:3880022/71‑1 (MQ=255)                      tCAACTGCCACTTCACCGTCGGCGCGCATTTCAATATCACGTCCGTGGCCTGCCGCATCTGAACGCGTGTc                                                  <  1:2592087/71‑1 (MQ=255)                      tCAACTGCCACTTCACCGTCGGCGCGCATTTCAATATCACGTCCGTGGCCTGCCGCATCTGAACGCGTGTc                                                  <  1:2541297/71‑1 (MQ=255)                       cAACTGCCACTTCACCGTCGGCGCGCATTTCAATATCACGTCCGTGGCCTGCCGCATCTGAACGCGTGTc                                                  >  1:1095679/1‑70 (MQ=255)                        aaCTGCCACTTCACCGTCGGCGCGCATTTCAATATCACGTCCGTGGCCTGCCGCATCTGAACGCGTGTCaa                                                >  1:631704/1‑71 (MQ=255)                        aaCTGCCACTTCACCGTCGGCGCGCATTTCAATATCACGTCCGTGGCCTGCCGCATCTGAACGCGTGTCaa                                                >  1:2565649/1‑71 (MQ=255)                        aaCTGCCACTTCACCGTCGGCGCGCATTTCAATATCACGTCCGTGGCCTGCCGCATCTGAACGCGTGTCaa                                                >  1:3631781/1‑71 (MQ=255)                         aCTGCCACTTCACCGTCGGCGCGCATTTCAATATGACGTCCGTGGCCTGCCGCATCTGAAGGCGTGTCAAc                                               <  1:1791990/71‑1 (MQ=255)                            gCCACTTCACCGTCGGCGCGCATTTCAATATCACGTCCGTGGCCTGCCGCATCTGAACGCGTGTc                                                  <  1:2657126/65‑1 (MQ=255)                             ccACTTCACCGTCGGCGCGCATTTCAATATCACGTCCGTGGCCTGCCGCATCTGAACGCGTGTCAACGt                                             <  1:968305/69‑1 (MQ=255)                              cACTTCACCGTCGGCGCGCATTTCAATATCACGTCCGTGGCCTGCCGCATCTGAACGCGTGTCAACGTCgg                                          <  1:2047683/71‑1 (MQ=255)                              cACTTCACCGTCGGCGCGCATTTCAATATCACGTCCGTGGCCTGCCGCATCTGAACGCGTGTCAACGTCgg                                          <  1:3260543/71‑1 (MQ=255)                              cACTTCACCGTCGGCGCGCATTTCAATATCACGTCCGTGGCCTGCCGCATCTGAACGCGTGGCAACGTCgg                                          <  1:2492040/71‑1 (MQ=255)                                cTTCACCGTCGGCGCGCATTTCAATATCACGTCCGTGGCCTGCCGCATCTGAACGCGTGTCAACGTCgg                                          <  1:2990380/69‑1 (MQ=255)                                cTTCACCGTCGGCGCGCATTTCAATATCACGTCCGTGGCCTGCCGCATCTGAACGCGTGTCAACGTCGGAg                                        <  1:351539/71‑1 (MQ=255)                                  tCACCGTCGGCGCGCATTTCAATATCACGTCCGTGGTCTGCCGCATCTGAACGCGTGTCAACGTCGGAgg                                       >  1:3371031/1‑70 (MQ=255)                                  tCACCGTCGGCGCGCATTTCAATATCACGTCCGTGGCCTGCCGCATCTGAACGCGTGTCAACGTCGGAgg                                       >  1:728895/1‑70 (MQ=255)                                   cACCGTCGGCGCGCATTTCAATATCACGTCCGTGGCCTGCCGCATCTg                                                            <  1:3922731/48‑1 (MQ=255)                                   cACCGTCGGCGCGCATTTCAATATCACGTCCGTGGCCTGCCGCATCTGAACGCGTGTc                                                  >  1:412844/1‑58 (MQ=255)                                    aCCGTCGGCGCGCATTTCAATATCACGTCCGTGGCCTGCCGCATCTGAACGCGTGTCAACGTCGGAGGTaa                                    >  1:2507273/1‑71 (MQ=255)                                     ccGTCGGCGCGCATTTCAATATCACGTCCGTGGTCTGCCGCATCTGAACGCGTGTCAACGTCGGAg                                        <  1:1239445/66‑1 (MQ=255)                                     ccGTCGGCGCGCATTTCAATATCACGTCCGTGGCCTGCCGCATCTGAACGCGTGTc                                                  <  1:3493939/56‑1 (MQ=255)                                       gTCGGCGCGCATTTCAATATCACGTCCGTGGCCTGCCGCATCTGaa                                                          >  1:3756679/1‑46 (MQ=255)                                       gTCGGCGCGCATTTCAATATCACGTCCGTGGCCTGCCGCATCTGaa                                                          >  1:1155231/1‑46 (MQ=255)                                        tCGGCGCGCATTTCAATATCACGTCCGTGGCCTGCCGCATCTGAACGCGTGTCAACGTCGGAGGTAAtctc                                <  1:365225/71‑1 (MQ=255)                                          ggCGCGCATTTCAATATCACGTCCGTGGCCTGCCGCATCTGAACGCGTGTCAACGTCGGAGGTAATCTcc                               <  1:3965800/70‑1 (MQ=255)                                             gcgcATTTCAATATCACGTCCGTGGCCTGCCGCATCTg                                                            <  1:3160362/38‑1 (MQ=255)                                             gcgcATTTCAATATCACGTCCGTGGCCTGCCGCATCTGAACGCGTGTCAACGTCGGAGGTAATCTCCATCa                           >  1:2963905/1‑71 (MQ=255)                                                cATTTCAATATCACGTCCGTGGCCTGCCGCATCTGAACGCGTGTCAACGTCGGAGGTAATCTCCATCACGc                        >  1:1852278/1‑71 (MQ=255)                                                cATTTCAATATCACGTCCGTGGCCTGCCGCATCTGAACGCGTGTCAACGTCGGAGGTAATCTCCATCACGc                        <  1:2765424/71‑1 (MQ=255)                                                  tttCAATATCACGTCCGTGGCCTGCCGCATCTGAACGCGTGTCaa                                                >  1:3471387/1‑45 (MQ=255)                                                  tttCAATATCACGTCCGTGGCCTGCCGCATCTGAACGCGTGTCAACGTCGGAgg                                       <  1:2786449/54‑1 (MQ=255)                                                    tCAATATCACGTCCGTGGCCTGCCGCATCTGAACGCGTGTCAACGTCGGAGGTAATCTCCATCACGCTTTc                    >  1:1229851/1‑71 (MQ=255)                                                      aaTATCACGTCCGTGGCCTGCCGCATCTGAACGCTTGTCAACGTCGGAGGTAATCTCCATCACGCTTTCtt                  <  1:3612127/71‑1 (MQ=255)                                                      aaTATCACGTCCGTGGCCTGCCGCATCTGAACGCGTGTCAACGTCGGAGGTAATCTCCATCACGCTTTCtt                  >  1:2714263/1‑71 (MQ=255)                                                       atatCACGTCCGTGGTCTGCCGCATCTGAACGCgtg                                                    <  1:1980083/36‑1 (MQ=255)                                                       atatCACGTCCGTGGCCTGCCGCATCTGAACGCGTGTCAACGTCGGAGGTAATCTCCATCACGCTTTCttg                 <  1:2295196/71‑1 (MQ=255)                                                       atatCACGTCCGTGGCCTGCCGCATCTGAACGCGTGTCAACGTCGGAGGTAATCTCCATCACGCTTTCttg                 <  1:824433/71‑1 (MQ=255)                                                        tatCACGTCCGTGGCCTGCCGCATCTGAACGCGTGTCAACGTCGGAGGTAAt                                   >  1:3942588/1‑52 (MQ=255)                                                         atCACGTCCGTGGCCTGCCGCATCTGAACGCGTGTCa                                                 <  1:1945910/37‑1 (MQ=255)                                                          tCACGTCCGTGGCCTGCCGCATCTGAACGCGTGTCAACGTCGGAGGTAATCTCCATCACGCTTTCttgttg              >  1:1989671/1‑71 (MQ=255)                                                            aCGTCCGTGGTCTGCCGCATCTGAACGCGTGTCAACGTCGGAGGTAATCTCCATCACGCTTTCTTGTTGCt            >  1:3125039/1‑71 (MQ=255)                                                             cGTCCGTGGCCTGCCGCATCTGAACGCGTGTCAACGTCGGAGGTAATCTCCATCACGCTTTCt                   >  1:506207/1‑63 (MQ=255)                                                                 cGTGGCCTGCCGCATCTGAACGCGTGTCAACGTCGGAGGTAATCTCCATCACGCTTTCTTGTTGCTTATc        <  1:1446767/70‑1 (MQ=255)                                                                 cGTGGCCTGCCGCATCTGAACGCGTGTCAACGTCGGAGGTAATCTCCATCACGCTTTCTTGTTGCTTATCg       >  1:3047071/1‑71 (MQ=255)                                                                  gTGGCCTGCCGCATCTGAACGCGTGTCAACGTCGGAGGTAATCTCCAt                             <  1:2808053/48‑1 (MQ=255)                                                                      ccTGCCGCATCTGAACGCGTGTCAACGTCGGAGGTAATCTCCATCACGCTTTCTTGTTGCTTATCGGTaaa  <  1:3843044/71‑1 (MQ=255)                                                                      |                                                                       CACTGCGTGTCTACCCCCGCATCAACTGCCACTTCACCGTCGGCGCGCATTTCAATATCACGTCCGTGGCCTGCCGCATCTGAACGCGTGTCAACGTCGGAGGTAATCTCCATCACGCTTTCTTGTTGCTTATCGGTAAA  >  minE/1469788‑1469927

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