Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	2,557,631	G→A	18.6%	A95A (GCG→GCA)	secB →	protein export chaperone

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	2,557,631	0	G	A	18.6%	133.4 / 18.2	85	A95A (GCG→GCA)	secB	protein export chaperone
Reads supporting (aligned to +/- strand):  ref base G (36/33);  new base A (5/11);  total (41/44)
Fisher's exact test for biased strand distribution p-value = 1.69e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CGTTCCTGTGTGAAGTTCAGCAGGGCGGTATTTTCTCCATCGCGGGTATCGAAGGCACCCAGATGGCGCATTGCCTGGGAGCATACTGCCCGAACATTCTGTTCCCGTATGCTCGTGAGTGCATCACCAGCATGGT  >  minE/2557564‑2557699                                                                    |                                                                     cgTTCCTGTGTGAAGTTCAGCAGGGCGGTATTTTCTCCATCGCGGGTATCGAAGGCACCCAGATGgcgc                                                                     <  1:3375026/69‑1 (MQ=255)     ccTGTGTGAAGTTCAGCAGGGCGGTATTTTCTCCATCGCGGGTATCGAAGGCACCCAGATGGCGCATTg                                                                 >  1:1854638/1‑69 (MQ=255)     ccTGTGTGAAGTTCAGCAGGGCGGTATTTTCTCCATCGCGGGTATCGAAGGCACCCAGATGGCGCATTg                                                                 >  1:688756/1‑69 (MQ=255)        gtgtgAAGTTCAGCAGGGCGGTATTTTCTCCATCGCGGGTATCGAAGGCACCCAGATGGCGCATTGCCt                                                              <  1:2378852/69‑1 (MQ=255)        gtgtgAAGTTCAGCAGGGCGGTATTTTCTCCATCGCGGGTATCGAAGGCACCCAGATGGCGCATTGCCt                                                              <  1:1333100/69‑1 (MQ=255)           tgAAGTTCAGCAGGGCGGTATTTTCTCCATCGCGGGTATCGAAGGCACCCAGATGGCGCATTGCCTggg                                                           <  1:1007924/69‑1 (MQ=255)             aaGTTCAGCAGGGCGGTATTTTCTCCATCGCGGGTATCGAAGGCACCCAGATGGCGCATTGCCTGGGAg                                                         <  1:1394194/69‑1 (MQ=255)               gTTCAGCAGGGCGGTATTTTCTCCATCGCGGGTATCGAAGGCACCCAGATGGCGCATTGCCTGGGAGCa                                                       <  1:1817181/69‑1 (MQ=255)               gTTCAGCAGGGCGGTATTTTCTCCATCGCGGGTATCGAAGGCACCCAGATGGCGCATTGCCTGGGAGCa                                                       <  1:3278901/69‑1 (MQ=255)               gTTCAGCAGGGCGGTATTTTCTCCATCGCGGGTATCGAAGGCACCCAGATGGCGCATTGCCTGGGAGCa                                                       <  1:487776/69‑1 (MQ=255)               gTTCAGCAGGGCGGTATTTTCTCCATCGCGGGTATCGAAGGCACCCAGATGGCGCATTGCCTGGGAGCa                                                       <  1:133077/69‑1 (MQ=255)                ttCAGCAGGGCGGTATTTTCTCCATCGCGGGTATCGAAGGCACCCAGATGGCGCATTGCCTGGGAGCa                                                       <  1:2270888/68‑1 (MQ=255)                               ttttCTCCATCGCGGGTATCGAAGGCACCGAGATGGCGCATTGCCTGGGAGCATACTg                                                  >  1:2966808/1‑58 (MQ=255)                               ttttCTCCATCGCGGGTATCGAAGGCACCCAGATGGCGCATTGCCTGGGAGCATACTg                                                  >  1:442836/1‑58 (MQ=255)                               ttttCTCCATCGCGGGTATCGAAGGCACCCAGATGGCGCATTGCCTGGGAGCATACTg                                                  >  1:1521224/1‑58 (MQ=255)                               ttttCTCCATCGCGGGTATCGAAGGCACCCAGATGGCGCATTGCCTGGGAGCATACTg                                                  >  1:152720/1‑58 (MQ=255)                               ttttCTCCATCGCGGGTATCGAAGGCACCCAGATGGCGCATTGCCTGGGAGCATACTg                                                  >  1:1559660/1‑58 (MQ=255)                               ttttCTCCATCGCGGGTATCGAAGGCACCCAGATGGCGCATTGCCTGGGAGCATACTg                                                  >  1:2657825/1‑58 (MQ=255)                               ttttCTCCATCGCGGGTATCGAAGGCACCCAGATGGCGCATTGCCTGGGAGCATACTg                                                  >  1:1116192/1‑58 (MQ=255)                               ttttCTCCATCGCGGGTATCGAAGGCACCCAGATGGCGCATTGCCTGGGAGCATACTg                                                  >  1:2833966/1‑58 (MQ=255)                               ttttCTCCATCGCGGGTATCGAAGGCACCCAGATGGCGCATTGCCTGGGAGCATACTg                                                  >  1:1703506/1‑58 (MQ=255)                                 ttCTCCATCGCGGGTATCTAAGGCACCCAGATGGCGCATTGCCTGGGAg                                                         <  1:2499124/49‑1 (MQ=255)                                 ttCTCCATCGCGGGTATCGAAGGCACCCAGATGGCGCATTGCCt                                                              >  1:2754795/1‑44 (MQ=255)                                 ttCTCCATCGCGGGTATCGAAGGCACCCAGATGGCGCATTGCCTTGGAg                                                         <  1:3070856/49‑1 (MQ=255)                                 ttCTCCATCGCGGGTATCGAAGGCACCCAGATGGCGCATTGCCTGGGAGCa                                                       >  1:2418243/1‑51 (MQ=255)                                 ttCTCCATCGCGGGTATCGAAGGCACCCAGATGGCGCATTGCCTGGGAGCa                                                       >  1:927479/1‑51 (MQ=255)                                 ttCTCCATCGCGGGTATCGAAGGCACCCAGATGGCGCATTGCCTGGGAGCa                                                       >  1:1661213/1‑51 (MQ=255)                                 ttCTCCATCGCGGGTATCGAAGGCACCCAGATGGCGCATTGCCTGGGAGCa                                                       >  1:2697246/1‑51 (MQ=255)                                 ttCTCCATCGCGGGTATCGAAGGCACCCAGATGGCGCATTGCCTGGGAGCATACTGCCCGAACATTCTg                                     >  1:3211663/1‑69 (MQ=255)                                 ttCTCCATCGCGGGTATCGAAGGCACCCAGATGGCGCATTGCCTGGGAGCATACTGCCCGAACATTCTg                                     >  1:2288861/1‑69 (MQ=255)                                 ttCTCCATCGCGGGTATCGAAGGCACCCAGATGGCGCATTGCCTGGGAGCATACTGCCCGAACATTCTg                                     >  1:2217764/1‑69 (MQ=255)                                         cgcgGGTATCGAAGGCACCCAGATGGCGCATTGCCTGGGAGCATACTGCCCGAACATTCTGTTccc                                >  1:2041626/1‑66 (MQ=255)                                         cgcgGGTATCGAAGGCACCCAGATGGCGCATTGCCTGGGAGCATACTGCCCGAACATTCTGTTCCCGt                              >  1:681244/1‑68 (MQ=255)                                         cgcgGGTATCGAAGGCACCCAGATGGCGCATTGCCTGGGAGCATACTGCCCGAACATTCTGTTCCCGTa                             >  1:1950354/1‑69 (MQ=255)                                         cgcgGGTATCGAAGGCACCCAGATGGCGCATTGCCTGGGAGCATACTGCCCGAACATTCTGTTCCCGTa                             >  1:872905/1‑69 (MQ=255)                                         cgcgGGTATCGAAGGCACCCAGATGGCGCATTGCCTGGGAGCATACTGCCCGAACATTCTGTTCCCGTa                             >  1:1608834/1‑69 (MQ=255)                                         cgcgGGTATCGAAGGCACCCAGATGGCGCATTGCCTGGGAGCATACTGCCCGAACATTCTGTTCCCGTa                             >  1:3301975/1‑69 (MQ=255)                                         cgcgGGTATCGAAGGCACCCAGATGGCGCATTGCCTGGGAGCATACTGCCCGAACATTCTGTTCCCGTa                             >  1:3275528/1‑69 (MQ=255)                                            gggTATCGAAGGCACCCAGATGGCACATTGCCTGGGCGCATACTGCCCGAACATTCTGTTCCCGTATGc                          >  1:2733894/1‑69 (MQ=255)                                            gggTATCGAAGGCACCCAGATGGCACATTGCCTGGGCGCATACTGCCCGAACATTCTGTTCCCGTATGc                          >  1:1608694/1‑69 (MQ=255)                                            gggTATCGAAGGCACCCAGATGGCACATTGCCTGGGCGCATACTGCCCGAACATTCTGTTCCCGTATGc                          >  1:2768343/1‑69 (MQ=255)                                            gggTATCGAAGGCACCCAGATGGCACATTGCCTGGGCGCATACTGCCCGAACATTCTGTTCCCGTATGc                          >  1:250646/1‑69 (MQ=255)                                            gggTATCGAAGGCACCCAGATGGCACATTGCCTGGGCGCATACTGCCCGAACATTCTGTTCCCGTATGc                          >  1:1065191/1‑69 (MQ=255)                                               tATCGAAGGCACCCAGATGGCGCATTGCCTGGGAGCATACTGCCCGAACATTCTGTTCCCGTATGCTc                        >  1:2725410/1‑68 (MQ=255)                                               tATCGAAGGCACCCAGATGGCGCATTGCCTGGGAGCATACTGCCCGAACATTCTGTTCCCGTATGCTc                        >  1:215258/1‑68 (MQ=255)                                                aTCGAAGGCACCCAGATGGCGCATTGCCTGGGAGCATACTGCCCGAACATTCTGTTCCCGTATGCTc                        >  1:993058/1‑67 (MQ=255)                                                 tCGAAGGCACCCAGATGGCGCATTGCCTGGGAGCATACTGCCCGAACATTCTGTTCCCGTATGCTCGTg                     <  1:2057310/69‑1 (MQ=255)                                                 tCGAAGGCACCCAGATGGCGCATTGCCTGGGAGCATACTGCCCGAACATTCTGTTCCCGTATGCTCGTg                     <  1:2038342/69‑1 (MQ=255)                                                 tCGAAGGCACCCAGATGGCGCATTGCCTGGGAGCATACTGCCCGAACATTCTGTTCCCGTATGCTCGTg                     <  1:771044/69‑1 (MQ=255)                                                 tCGAAGGCACCCAGATGGCGCATTGCCTGGGAGCATACTGCCCGAACATTCCGTTCCCGTATGCTCGTg                     <  1:1057837/69‑1 (MQ=255)                                                   gAAGGCACCCAGATGGCGCATTGCCTGGGAGCATACTGCCCGAACa                                          <  1:726917/46‑1 (MQ=255)                                                   gAAGGCACCCAGATGGCGCATTGCCTGGGAGCATACTGCCCGAACATTCTGTTc                                  >  1:1656154/1‑54 (MQ=255)                                                   gAAGGCACCCAGATGGCGCATTGCCTGGGAGCATACTGCCCGAACATTCTGTTCCCg                               <  1:254193/57‑1 (MQ=255)                                                   gAAGGCACCCAGATGGCGCATTGCCTGGGAGCATACTGCCCGAACATTCTGTTCCCg                               <  1:1972026/57‑1 (MQ=255)                                                   gAAGGCACCCAGATGGCGCATTGCCTGGGAGCATACTGCCCGAACATTCTGTTCCCg                               <  1:896017/57‑1 (MQ=255)                                                   gAAGGCACCCAGATGGCGCATTGCCTGGGAGCATACTGCCCGAACATTCTGTTCCCg                               <  1:2066819/57‑1 (MQ=255)                                                   gAAGGCACCCAGATGGCGCATTGCCTGGGAGCATACTGCCCGAACATTCTGTTCCCg                               <  1:1333543/57‑1 (MQ=255)                                                   gAAGGCACCCAGATGGCGCATTGCCTGGGAGCATACTGCCCGAACATTCTGTTCCCg                               <  1:862784/57‑1 (MQ=255)                                                   gAAGGCACCCAGATGGCGCATTGCCTGGGAGCATACTGCCCGAACATTCTGTTCCCg                               <  1:1649430/57‑1 (MQ=255)                                                   gAAGGCACCCAGATGGCGCATTGCCTGGGAGCATACTGCCCGAACATTCTGTTCCCGTATGCTCGTGAg                   <  1:668640/69‑1 (MQ=255)                                                   gAAGGCACCCAGATGGCGCATTGCCTGGGAGCATACTGCCCGAACATTCTGTTCCCGTATGCTCGTGAg                   >  1:1652118/1‑69 (MQ=255)                                                   gAAGGCACCCAGATGGCGCATTGCCTGGGAGCATACTGCCCGAACATTCTGTTCCCGTATGCTCGTGAg                   <  1:1365716/69‑1 (MQ=255)                                                   gAAGGCACCCAGATGGCGCATTGCCTGGGAGCATACTGCCCGAACATTCTGTTCCCGTATGCTCGTGAg                   <  1:2020557/69‑1 (MQ=255)                                                   gAAGGCACCCAGATGGCGCATTGCCTGGGAGCATACTGCCCGAACATTCTGTTCCCGTATGCTCGTGAg                   <  1:842861/69‑1 (MQ=255)                                                   gAAGGCACCCAGATGGCGCATTGCCTGGGAGCATACTGCCCGAACATTCTGTTCCCGTATGCTCGTGAg                   <  1:2170587/69‑1 (MQ=255)                                                   gAAGGCACCCAGATGGCGCATTGCCTGGGAGCATACTGCCCGAACATTCTGTTCCCGTATGCTCGTGAg                   <  1:222134/69‑1 (MQ=255)                                                   gAAGGCACCCAGATGGCGCATTGCCTGGGAGCATACTGCCCGAACATTCTGTTCCCGTATGCTCGTGAg                   >  1:2236176/1‑69 (MQ=255)                                                        cACCCAGATGGCACATTGCCTGGGCGCATACTGCCCGAACATTCTGTTCCCGTATGCTCGTGAGTGCAt              <  1:1048471/69‑1 (MQ=255)                                                        cACCCAGATGGCACATTGCCTGGGCGCATACTGCCCGAACATTCTGTTCCCGTATGCTCGTGAGTGCAt              <  1:478169/69‑1 (MQ=255)                                                        cACCCAGATGGCACATTGCCTGGGCGCATACTGCCCGAACATTCTGTTCCCGTATGCTCGTGAGTGCAt              <  1:826880/69‑1 (MQ=255)                                                        cACCCAGATGGCACATTGCCTGGGCGCATACTGCCCGAACATTCTGTTCCCGTATGCTCGTGAGTGCAt              <  1:97286/69‑1 (MQ=255)                                                        cACCCAGATGGCACATTGCCTGGGCGCATACTGCCCGAACATTCTGTTCCCGTATGCTCGTGAGTGCAt              <  1:1516002/69‑1 (MQ=255)                                                        cACCCAGATGGCACATTGCCTGGGCGCATACTGCCCGAACATTCTGTTCCCGTATGCTCGTGAGTGCAt              <  1:3199278/69‑1 (MQ=255)                                                        cACCCAGATGGCACATTGCCTGGGCGCATACTGCCCGAACATTCTGTTCCCGTATGCTCGTGAGTGCAt              <  1:3049050/69‑1 (MQ=255)                                                        cACCCAGATGGCACATTGCCTGGGCGCATACTGCCCGAACATTCTGTTCCCGTATGCTCGTGAGTGCAt              <  1:3009568/69‑1 (MQ=255)                                                        cACCCAGATGGCACATTGCCTGGGCGCATACTGCCCGAACATTCTGTTCCCGTATGCTCGTGAGTGCAt              <  1:2976618/69‑1 (MQ=255)                                                        cACCCAGATGGCACATTGCCTGGGCGCATACTGCCCGAACATTCTGTTCCCGTATGCTCGTGAGTGCAt              <  1:2050772/69‑1 (MQ=255)                                                        cACCCAGATGGCACATTGCCTGGGCGCATACTGCCCGAACATTCTGTTCCCGTATGCTCGTGAGTGCAt              <  1:2461186/69‑1 (MQ=255)                                                          cccAGATGGCGGATTGCCTGGGAGCATACTGCCCGAACATTCTGTTCCCGTATGCTCGTGAGTGCATc             <  1:1990563/68‑1 (MQ=255)                                                                 ggCGCATTGCCTGGGAGCATACTGCCCGAACATTCTGTTCCCGTATGCTCGTGAGTGCATCAc           <  1:424132/63‑1 (MQ=255)                                                                 ggCGCATTGCCTGGGAGCATACTGCCCGAACATTCTGTTCCCGTATGCTCGTGAGTGCATCAc           <  1:1462535/63‑1 (MQ=255)                                                                 ggCGCATTGCCTGGGAGCATACTGCCCGAACATTCTGTTCCCGTATGCTCGTGAGTGCATCAc           <  1:2427240/63‑1 (MQ=255)                                                                 ggCGCATTGCCTGGGAGCATACTGCCCGAACATTCTGTTCCCGTATGCTCGTGAGTGCATCACCAGCAt     >  1:3324936/1‑69 (MQ=255)                                                                 ggCGCATTGCCTGGGAGCATACTGCCCGAACATTCTGTTCCCGTATGCTCGTGAGTGCATCACCAGCAt     >  1:1361268/1‑69 (MQ=255)                                                                 ggCGCATTGCCTGGGAGCATACTGCCCGAACATTCTGTTCCCGTATGCTCGTGAGTGCATCACCAGCAt     >  1:2494988/1‑69 (MQ=255)                                                                 ggCGCATTGCCTGGGAGCATACTGCCCGAACATTCTGTTCCCGTATGCCCGTGAGTGCATCACCAGCAt     >  1:2985941/1‑69 (MQ=255)                                                                    gcATTGCCTGGGAGCATACTGCCCGAACATTCTGTTCCCGTATGCTCGTGAGTGCATCACCAGCATGGt  >  1:2135788/1‑69 (MQ=255)                                                                    |                                                                     CGTTCCTGTGTGAAGTTCAGCAGGGCGGTATTTTCTCCATCGCGGGTATCGAAGGCACCCAGATGGCGCATTGCCTGGGAGCATACTGCCCGAACATTCTGTTCCCGTATGCTCGTGAGTGCATCACCAGCATGGT  >  minE/2557564‑2557699

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