Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	1,279,189	G→A	8.3%	W323* (TGG→TAG)	yeeJ →	adhesin

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	1,279,189	0	G	A	8.3%	147.1 / 3.6	48	W323* (TGG→TAG)	yeeJ	adhesin
Reads supporting (aligned to +/- strand):  ref base G (20/24);  new base A (2/2);  total (22/26)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.85e-01

AGCGCACCTGAACTGGACAACGATTATGAAGCACGCCCGGCCAATGGCTGGGATGTACGCGCAGAAAGCTGGCTACCCGCCTGGCCGCACCTTGGCGGTAAACTGGTCTATGAACAGTATTATGGCGATGAAGTGGCCC  >  minE/1279119‑1279257                                                                       |                                                                     agcgcaCCTGAACTGGACAACGATTATGAAGCACGCCCGGCCAATGGCTGGGATGTACGCGCAGAAAGCTg                                                                      <  1:3966774/71‑1 (MQ=255)         tGAACTGGACAACGATTATGAAGCACGCCCGGCCAATGGCTGGGATGTACGCGCAGAAAGCTGGCTACCCg                                                              >  1:1757504/1‑71 (MQ=255)          gAACTGGACAACGATTATGAAGCACGCCCGGCCAATGGCTGGGATGTACGCGCAGAAAGCTGGCTACCCGc                                                             >  1:2566779/1‑71 (MQ=255)           aaCTGGACAACGATTATGAAGCACGCCCGGCCAATGGCTGGGATGTACGCGCAGAAAGCTGGCTACCCg                                                              <  1:4109024/69‑1 (MQ=255)             cTGGACAACGATTATGAAGCACGCCCGGCCAATGGCTGGGATGTACGCGCAGAAAGCTGGCTACCCGCCTg                                                          <  1:1985323/71‑1 (MQ=255)                 aCAACGATTATGAAGCACGCCCGGCCAATGGCTGGGATGTACGCGCAGAAAGCTGGCTACCCGCCTGGCCg                                                      <  1:2761269/71‑1 (MQ=255)                 aCAACGATTATGAAGCACGCCCGGCCAATGGCTGGGATGTACGCGCAGAAAGCTGGCTACCCGCCTGGCCg                                                      <  1:2104505/71‑1 (MQ=255)                   aaCGATTATGAAGCACGCCCGGCCAATGGCTGGGATGTACGCGCAGAAAGCTGGCTACCCGCCTg                                                          <  1:956121/65‑1 (MQ=255)                   aaCGATTATGAAGCACGACCGGCCAATGGCTGGGATGTACGCGCAGAAAGCTGGCTACCCGCCt                                                           >  1:3196251/1‑64 (MQ=255)                    aCGATTATGAAGCACGCCCGGCCAATGGCTGGGATGTACGCGCAGAAAGCTGGCTACCCGCCTGGCCGCAc                                                   >  1:1625532/1‑71 (MQ=255)                     ggATTATCAGGAACGCCCGGCGAATGGCTGGGATATTCGTGCTGAGGGCTATTTACCCGCCTGGCCGCagc                                                  >  1:3126947/2‑69 (MQ=255)                      gATTATGAAGCACGCCCGGCCAATGGCTGGGATGTACGCGCAGAAAGCTGGCTACCCGcc                                                            >  1:581064/1‑60 (MQ=255)                      gATTATGAAGCACGCCCGGCCAATGGCTGGGATGTACGCGCAGAAAGCTGGCTACCCGCCTGGCCGCACCt                                                 >  1:910732/1‑71 (MQ=255)                      gATTATCAGGAACGCCCGGCGAATGGCTGGGATATTCGTGCTGAGGGCTATTTACCCGCCTGGCCGCAGCt                                                 <  1:2662854/71‑1 (MQ=255)                       aTTATCAGGAACGCCCGGCGAATGGCTGGGATATTCGTGCTGAGGGCTATTTACCCGCCTGGCCGc                                                     <  1:1419477/66‑1 (MQ=255)                              aGCACGCCCGGCCAATGGCTGGGATGTACGCGCAGAAAGCTGGCTACCCGCCTGGCCGCACCTTGGCGGTa                                         <  1:485619/71‑1 (MQ=255)                              aGCACGCCCGGCCAATGGCTGGGATGTACGCGCAGAAAGCTGGCTACCCGCCTGGCCGCACCTTGGCGGTa                                         <  1:2608193/71‑1 (MQ=255)                              aGCACGCCCGGCCAATGGCTGGGATGTACGCGCAGAAAGCTGGCTACCCGCCTGGCCGCACCTTGGCGGTa                                         <  1:2513920/71‑1 (MQ=255)                               gCACGCCCGGCCAATGGCTGGGATGTACGCGCAGAAAGCTGGCTACCCGCCTGGCCGCACCTTGGCGGTa                                         <  1:2702910/70‑1 (MQ=255)                                        gCCAATGGCTGGGATGTACGCGCAGAAAGCTGGCTACCCGCCTGGCCGCACCTTGGCGGTaaa                                       >  1:2985413/1‑63 (MQ=255)                                         ccAATGGCTGGGATGTACGCGCAGAAAGCTGGCTACCCGCCTGGCCGCACCTTGGCGGTa                                         <  1:2165892/60‑1 (MQ=255)                                          cAATGGCTGGGATGTACGCGCAGAAAGCTGGCTACCCGCCTGGCCGCACCTTGGCGGTAAACTGGTCTAt                              >  1:2391338/1‑70 (MQ=255)                                          cAATGGCTGGGATGTACGCGCAGAAAGCTGGCTACCCGCCTGGCCGCACCTTGGCGGTAAACTGGTCTATg                             <  1:3805810/71‑1 (MQ=255)                                            aTGGCTGGGATGTACGCGCAGAAAGCTGGCTACCCGCCTGGCCGCACCTTGGCGGTa                                         <  1:639173/57‑1 (MQ=255)                                             tGGCTGGGATATTCGTGCTGAGGGCTATTTACCCGCCTGGCCGCagc                                                  >  1:3804651/1‑45 (MQ=255)                                              ggCTGGGATGTACGCGCAGAAAGCTGGCTACCCGCCTGGCCGCACCTTGGCGGTaaa                                       >  1:2665002/1‑57 (MQ=255)                                              ggCTGGGATGTACGCGCAGAAAGCTGGCTACCCGCCTGGCCGCACCTTGGCGGTAAACTGGTCTa                               >  1:3931117/1‑65 (MQ=255)                                               gCTGGGATGTACGCGCAGAAAGCTGGCTACCCGCCTGGCCGCACCTTGGCGGTAAACTGGTCTATGAACAg                        <  1:2589805/71‑1 (MQ=255)                                               gCTGGGATGTACGCGCAGAAAGCTGGCTACCCGCCTGGCCGCACCTTGGCGGTAAACTGGTCTATGAACAg                        <  1:2032599/71‑1 (MQ=255)                                                cTGGGATGTACGCGCAGAAAGCTGGCTACCCGCCTGGCCGCACCTTGGCGGTAAACTGGTCTATGAACAGt                       <  1:792457/71‑1 (MQ=255)                                                    gATGTACGCGCAGAAAGCTGGCTACCCGCCTGGCCGCACCTTGGCGGTAAACTGGt                                  >  1:4045751/1‑56 (MQ=255)                                                    gATGTACGCGCAGAAAGCTGGCTACCCGCCTGGCCGCACCTTGGCGGTAAACTGGTCTATGAACAGtatta                   >  1:3454334/1‑71 (MQ=255)                                                     aTGTACGCGCAGAAAGCTGGCTACCCGCCTGGCCGCACCTTGGCGGTAAACTGGTCTATGAAc                          >  1:539588/1‑63 (MQ=255)                                                      tGTACGCGCAGAAAGCTGGCTACCCGCCTGGCCGCACCTTGGCGGTAAACTGGTCTATGAACAGta                      <  1:1465955/66‑1 (MQ=255)                                                         aCGCGCAGAAAGCTGGCTACCCGCCTGGCCGCACCTTGGc                                             >  1:2928907/1‑40 (MQ=255)                                                         aCGCGCAGAAAGCTGGCTACCCGCCTGGCCGCACCTTGGCGGTAAACTGGt                                  >  1:1348917/1‑51 (MQ=255)                                                          cgcgCAGAAAGCTGGCTACCCGCCTGGCCGCACCTTGGCGGTAAACTGGTCTATGAACAGTATTATGGCGa             >  1:438868/1‑71 (MQ=255)                                                           gcgcAGAAAGCTGGCTACCCGCCTGGCCGCACCTTGGCGGTAAACTGGTCTATGAACAGTATTATGGCGAt            >  1:2203185/1‑71 (MQ=255)                                                           gcgcAGAAAGCTGGCTACCCGCCTGGCCGCACCTTGGCGGTAAACTGGTCTATGAACAGTATTATGGCGAt            <  1:3917097/71‑1 (MQ=255)                                                            cgcAGAAAGCTGGCTACCCGCCTGGCCGCACCTTGGCGGTAAACTGGTc                                 <  1:1461261/49‑1 (MQ=255)                                                            cgcAGAAAGCTGGCTACCCGCCTGGCCGCACCTTGGCGGTAAACTGGTCTATGAACAGTATTATGGCGATg           <  1:3566021/71‑1 (MQ=255)                                                             gcAGAAAGCTGGCTACCCGCCTGGCCGCACCTTGGCGGTAAACTGGTCTATGAACAGTATTATGGCGATg           <  1:1852795/70‑1 (MQ=255)                                                             gcAGAAAGCTGGCTACCCGCCTGGCCGCACCTTGGCGGTAAACTGGTCTATGAACAGTATTATGGCGATg           <  1:1401545/70‑1 (MQ=255)                                                                gAAAGCTGGCTACCCGCCTGGCCGCACCTTGGCGGTAAACTGGTCTATGAAc                          <  1:1865654/52‑1 (MQ=255)                                                                  aaGCTGGCTACCCGCCTGGCCGCACCTTGGCGGTAAACTGGTCTATGAACAGTATTATGGc               >  1:1910476/1‑61 (MQ=255)                                                                  aaGCTGGCTACCCGCCTGGCCGCACCTTGGCGGTAAACTGGTCTATGAACAGTATTATGGCGATGAAGTgg     <  1:1878218/71‑1 (MQ=255)                                                                  aaGCTGGCTACCCGCCTGGCCGCACCTTGGCGGTAAACTGGTCTATGAACAGTATTATGGCGATGAAGTgg     >  1:2354644/1‑71 (MQ=255)                                                                     cTGGCTACCCGCCTGGCCGCACCTTGGCGGTAAACTGGTCTATGAACAGTATTATGGCGATGAAGTGGccc  >  1:921553/1‑71 (MQ=255)                                                                      tGGCTACCCGCCTGGCCGCACCTTGGCGGTAAACTGGTCTATGAACAGTATTATGGCGAt            <  1:247926/60‑1 (MQ=255)                                                                       |                                                                     AGCGCACCTGAACTGGACAACGATTATGAAGCACGCCCGGCCAATGGCTGGGATGTACGCGCAGAAAGCTGGCTACCCGCCTGGCCGCACCTTGGCGGTAAACTGGTCTATGAACAGTATTATGGCGATGAAGTGGCCC  >  minE/1279119‑1279257

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