Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	602,686	C→T	52.2%	S39S (TCC→TCT)	ycaD →	predicted transporter

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	602,686	0	C	T	52.2%	‑0.3 / 60.5	69	S39S (TCC→TCT)	ycaD	predicted transporter
Reads supporting (aligned to +/- strand):  ref base C (3/30);  new base T (17/19);  total (20/49)
Fisher's exact test for biased strand distribution p-value = 5.28e-04
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GTTGACTCTGGCGATTGCGGTGTTAAATACACTCGTGCCGCTTTGGCTCGCCCAGGAACACATGTCCACATGGCAGGTAGGCGTTGTCAGCTCATCCTATTTTACCGGCAACCTTGTCGGTACATTGCTGAC  >  minE/602620‑602751                                                                   |                                                                  gttgACTCTGGCGATTGCGGTGTTAAATACACTCGTGCCGCTTTGGCTCGCCCAGGAACACATGTCca                                                                  <  1:819065/68‑1 (MQ=255) gttgACTCTGGCGATTGCGGTGTTAAATACACTCGTGCCGCTTTGGCTCGCCCAGGAACACATGTCca                                                                  <  1:1700311/68‑1 (MQ=255) gttgACTCTGGCGATTGCGGTGTTAAATACACTCGTGCCGCTTTGGCTCGCCCAGGAACACATGTCca                                                                  <  1:551492/68‑1 (MQ=255) gttgACTCTGGCGATTGCGGTGTTAAATACACTCGTGCCGCTTTGGCTCGCCCAGGAACACATGTCca                                                                  <  1:1597583/68‑1 (MQ=255) gttgACTCTGGCGATTGCGGTGCTAAATACACTCGTGCCGCTTTGGCTCGCCCAGGAACACATGTCca                                                                  <  1:60939/68‑1 (MQ=255)           gCGATTGCGGTGTTAAATACACTCGTGCCGCTTTGGCTCGCCCAGGAACACATGTCCACATGGCAgg                                                         <  1:190209/67‑1 (MQ=255)           gCGATTGCGGTGTTAAATACACTCGTGCCGCTTTGGCTCGCCCAGGAACACATGTCCACATGGCAgg                                                         <  1:3025909/67‑1 (MQ=255)            cGATTGCGGTGTTAAATACACTCGTGCCGCTTTGGCTCGCCCAGGAACACATGTCTACATGGCAGGTAg                                                      >  1:1897125/1‑69 (MQ=255)            cGATTGCGGTGTTAAATACACTCGTGCCGCTTTGGCTCGCCCAGGAACACATGTCTACATGGCAGGTAg                                                      >  1:1780813/1‑69 (MQ=255)            cGATTGCGGTGTTAAATACACTCGTGCCGCTTTGGCTCGCCCAGGAACACATGTCTACATGGCAGGTAg                                                      >  1:137314/1‑69 (MQ=255)            cGATTGCGGTGTTAAATACACTCGTGCCGCTTTGGCTCGCCCAGGAACACATGTCTACATGGCAGGTAg                                                      >  1:2576949/1‑69 (MQ=255)            cGATTGCGGTGTTAAATACACTCGTGCCGCTTTGGCTCGCCCAGGAACACATGTCTACATGGCAGGTAg                                                      >  1:2394228/1‑69 (MQ=255)            cGATTGCGGTGTTAAATACACTCGTGCCGCTTTGGCTCGCCCAGGAACACATGTCCACATGGCAgg                                                         <  1:3155191/66‑1 (MQ=255)               ttGCGGTGTTAAATACACTCGTGCCGCTTTGGCTCGCCCAGGAACACATGTCTACATGGCATGTAGGCg                                                   <  1:1168070/69‑1 (MQ=255)               ttGCGGTGTTAAATACACTCGTGCCGCTTTGGCTCGCCCAGGAACACATGTCTACATGGCAGGTAGGCg                                                   <  1:1454816/69‑1 (MQ=255)               ttGCGGTGTTAAATACACTCGTGCCGCTTTGGCTCGCCCAGGAACACATGTCTACATGGCAGGTAGGCg                                                   <  1:2609432/69‑1 (MQ=255)               ttGCGGTGTTAAATACACTCGTGCCGCTTTGGCTCGCCCAGGAACACATGTCTACATGGCAGGTAGGCg                                                   <  1:2722527/69‑1 (MQ=255)                tGCGGTGTTAAATACACTCGTGCCGCTTTGGCTCGCCCAGGAACACATGTCTACATGGCAGGTAGGCg                                                   <  1:760859/68‑1 (MQ=255)                    gtgtTAAATACACTCGTGCCGCTTTGGCTCGCCCAGGAACACATGTCTACATGGCAGGTAGGCGTTGTc                                              <  1:2552897/69‑1 (MQ=255)                    gtgtTAAATACACTCGTGCCGCTTTGGCTCGCCCAGGAACACATGTCTACATGGCAGGTAGGCGTTGTc                                              <  1:3013623/69‑1 (MQ=255)                    gtgtTAAATACACTCGTGCCGCTTTGGCTCGCCCAGGAACACATGTCTACATGGCAGGTAGGCGTTGTc                                              <  1:2948499/69‑1 (MQ=255)                    gtgtTAAATACACTCGTGCCGCTTTGGCTCGCCCAGGAACACATGTCTACATGGCAGGTAGGCGTTGTc                                              <  1:3014087/69‑1 (MQ=255)                    gtgtTAAATACACTCGTGCCGCTTTGGCTCGCCCAGGAACACATGTCTACATGGCAGGTAGGCGTTGTc                                              <  1:2571926/69‑1 (MQ=255)                    gtgtTAAATACACTCGTGCCGCTTTGGCTCGCCCAGGAACACATGTCTACATGGCAGGTAGGCGTTGTc                                              <  1:2760935/69‑1 (MQ=255)                    gtgtTAAATACACTCGTGCCGCTTTGGCTCGCCCAGGAACACATGTCTACATGGCAGGTAGGCGTTGTc                                              <  1:310347/69‑1 (MQ=255)                    gtgtTAAATACACTCGTGCCGCTTTGGCTCGCCCAGGAACACATGTCTACATGGCAGGTAGGCGTTGTc                                              <  1:1807842/69‑1 (MQ=255)                    gtgtTAAATACACTCGTGCCGCTTTGGCTCGCCCAGGAACACATGTCTACATGGCAGGTAGGCGTTGTc                                              <  1:1759511/69‑1 (MQ=255)                    gtgtTAAATACACTCGTGCCGCTTTGGCTCGCCCAGGAACACATGTCTACATGGCAGGTAGGCGTTGTc                                              <  1:124087/69‑1 (MQ=255)                    gtgtTAAATACACTCGTGCCGCTTTGGCTCGCCCAGGAACACATGTCTACATGGCAGGTAGGCGTTGTc                                              <  1:635476/69‑1 (MQ=255)                    gtgtTAAATACACTCGTGCCGCTTTGGCTCGCCCAGGAACACATGTCTACATGGCAGGTAGGCGTTGTc                                              <  1:1603494/69‑1 (MQ=255)                    gtgtTAAATACACTCGTGCCGCTTTGGCTCGCCCAGGAACACATGTCTACATGGCAGGTAGGCGTTGTc                                              <  1:1448261/69‑1 (MQ=255)                    gtgtTAAATACACTCGTGCCGCTTTGGCTCGCCCAGGAACACATGTCTACATGGCAGGTAGGCGTTGTc                                              <  1:1477355/69‑1 (MQ=255)                           aTACACTCGTGCCGCTTTGGCTCGCCCAGGAACATATGTCTACATgg                                                             >  1:2502728/1‑47 (MQ=255)                           aTACACTCGTGCCGCTTTGGCTCGCCCAGGAACACATGTCTACATgg                                                             >  1:1810881/1‑47 (MQ=255)                           aTACACTCGTGCCGCTTTGGCTCGCCCAGGAACACATGTCTACATgg                                                             >  1:287165/1‑47 (MQ=255)                           aTACACTCGTGCCGCTTTGGCTCGCCCAGGAACACATGTCTACATgg                                                             >  1:290306/1‑47 (MQ=255)                           aTACACTCGTGCCGCTTTGGCTCGCCCAGGAACACATGTCTACATgg                                                             >  1:1404711/1‑47 (MQ=255)                           aTACACTCGTGCCGCTTTGGCTCGCCCAGGAACACATGTCTACATgg                                                             >  1:1596747/1‑47 (MQ=255)                           aTACACTCGTGCCGCTTTGGCTCGCCCAGGAACACATGTCTACATgg                                                             >  1:2534738/1‑47 (MQ=255)                           aTACACTCGTGCCGCTTTGGCTCGCCCAGGAACACATGTCTACATgg                                                             >  1:554598/1‑47 (MQ=255)                           aTACACTCGTGCCGCTTTGGCTCGCCCAGGAACACATGTCTACATgg                                                             >  1:362381/1‑47 (MQ=255)                           aTACACTCGTGCCGCTTTGGCTCGCCCAGGAACACATGTCTACATgg                                                             >  1:2055150/1‑47 (MQ=255)                           aTACACTCGTGCCGCTTTGGCTCGCCCAGGAACACATGTCTACATgg                                                             >  1:2053069/1‑47 (MQ=255)                           aTACACTCGTGCCGCTTTGGCTCGCCCAGGAACACATGTCTACATg                                                              >  1:118598/1‑46 (MQ=255)                            tACACTCGTGCCGCTTTGGCTCGCCCAGGAACACATGTCCACATGGCAGGTAGGCGTTGTCAGCTCATc                                      <  1:862047/69‑1 (MQ=255)                            tACACTCGTGCCGCTTTGGCTCGCCCAGGAACACATGTCCACATGGCAGGTAGGCGTTGTCAGCTCATc                                      <  1:1821421/69‑1 (MQ=255)                             acacTCGTGCCGCTTTGGCTCGCCCAGGAACACATGTCCACATGGCAGGTAGGCGTTGTCAGCTCATcc                                     <  1:3353072/69‑1 (MQ=255)                             acacTCGTGCCGCTTTGGCTCGCCCAGGAACACATGTCCACATGGCAGGTAGGCGTTGTCAGCTCATcc                                     <  1:678920/69‑1 (MQ=255)                             acacTCGTGCCGCTTTGGCTCGCCCAGGAACACATGTCCACATGGCAGGTAGGCGTTGTCAGCTCATcc                                     <  1:326535/69‑1 (MQ=255)                             acacTCGTGCCGCTTTGGCTCGCCCAGGAACACATGTCCACATGGCAGGTAGGCGTTGTCAGCTCATcc                                     <  1:3255362/69‑1 (MQ=255)                             acacTCGTGCCGCTTTGGCTCGCCCAGGAACACATGTCCACATGGCAGGTAGGCGTTGTCAGCTCATcc                                     <  1:3116120/69‑1 (MQ=255)                             acacTCGTGCCGCTTTGGCTCGCCCAGGAACACATGTCCACATGGCAGGTAGGCGTTGTCAGCTCATcc                                     <  1:255836/69‑1 (MQ=255)                             acacTCGTGCCGCTTTGGCTCGCCCAGGAACACATGTCCACATGGCAGGTAGGCGTTGTCAGCTCATcc                                     <  1:2653621/69‑1 (MQ=255)                                cTCGTGCCGCTTTGGCTCGCCCAGGAACACATGTCCACATGGCAGGTAGGCGTTGTCAGCTCATCCTa                                   <  1:424965/68‑1 (MQ=255)                                cTCGTGCCGCTTTGGCTCGCCCAGGAACACATGTCCACATGGCAGGTAGGCGTTGTCAGCTCATCCTa                                   <  1:1434836/68‑1 (MQ=255)                                  cGTGCCGCTTTGGCTCGCCCAGGAACACATGTCCACATGGCAGGTAGGCGTTGTCAGCTCATCCTAttt                                <  1:848157/69‑1 (MQ=255)                                            tGGCTCGCCCAGGATCACATGTCCACATGGCAGGTAGGCGTTGTCAGCTCATc                                      <  1:559428/53‑1 (MQ=255)                                                    ccATGAACACATGTCCACATGGCAGGTAGGCGTTGTCAGCTCATCCTATTTTACCGGCAACCTTGTCgg              <  1:1756026/69‑1 (MQ=255)                                                    ccAGGAACACATGTCCACATGGCAGGTAGGCGTTGTCAGCTCATCCTATTTTACCGGCAACCTTGTCgg              <  1:2166486/69‑1 (MQ=255)                                                    ccAGGAACACATGTCCACATGGCAGGTAGGCGTTGTCAGCTCATCCTATTTTACCGGCAACCTTGTCgg              <  1:1200046/69‑1 (MQ=255)                                                    ccAGGAACACATGTCCACATGGCAGGTAGGCGTTGTCAGCTCATCCTATTTTACCGGCAACCTTGTCgg              <  1:1200260/69‑1 (MQ=255)                                                    ccAGGAACACATGTCCACATGGCAGGTAGGCGTTGTCAGCTCATCCTATTTTACCGGCAACCTTGTCgg              <  1:2937887/69‑1 (MQ=255)                                                    ccAGGAACACATGTCCACATGGCAGGTAGGCGTTGTCAGCTCATCCTATTTTACCGGCAACCTTGTCgg              <  1:2629312/69‑1 (MQ=255)                                                    ccAGGAACACATGTCCACATGGCAGGTAGGCGTTGTCAGCTCATCCTATTTTACCGGCAACCTTGTCgg              <  1:1653057/69‑1 (MQ=255)                                                    ccAGGAACACATGTCCACATGGCAGGTAGGCGTTGTCAGCTCATCCTATTTTACCGGCAACCTTGTCgg              <  1:1339968/69‑1 (MQ=255)                                                     cAGGAACACATGTCCACATGGCTGGTAGGCGTTGTCAGCTCATcc                                     <  1:674842/45‑1 (MQ=255)                                                               tGTCCACATGGCAGGTAGGCGTTGTCAGCTCATCCTa                                   <  1:2781375/37‑1 (MQ=255)                                                               tGTCCACATGGCAGGTAGGCGTTGTCAGCTCATCCTa                                   <  1:1643644/37‑1 (MQ=255)                                                               tGTCCACATGGCAGGTAGGCGTTGTCAGCTCATCCTa                                   <  1:297423/37‑1 (MQ=255)                                                               tGTCCACATGGCAGGTAGGCGTTGTCAGCTCATCCTa                                   <  1:2013775/37‑1 (MQ=255)                                                               tGTCCACATGGCAGGTAGGCGTTGTCAGCTCATCCTa                                   <  1:2145882/37‑1 (MQ=255)                                                                gTCCACATGGCAGGTAGGCGTTGTCAGCTCATCCTATTTTACCGGCAACCTTGTCGGTACATTGCTGAc  >  1:374292/1‑69 (MQ=255)                                                                gTCCACATGGCAGGTAGGCGTTGTCAGCTCATCCTATTTTACCGGCAACCTTGTCGGTACATTGCTGAc  >  1:2775003/1‑69 (MQ=255)                                                                gTCCACATGGCAGGTAGGCGTTGTCAGCTCATCCTATTTTACCGGCAACCTTGTCGGTACATTGCTGAc  >  1:1775879/1‑69 (MQ=255)                                                                   |                                                                  GTTGACTCTGGCGATTGCGGTGTTAAATACACTCGTGCCGCTTTGGCTCGCCCAGGAACACATGTCCACATGGCAGGTAGGCGTTGTCAGCTCATCCTATTTTACCGGCAACCTTGTCGGTACATTGCTGAC  >  minE/602620‑602751

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