Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	1,191,560	T→C	18.2%	*116W (TGA→TGG)	rnd ←	ribonuclease D

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	1,191,560	0	T	C	18.2%	147.9 / 27.7	66	*116W (TGA→TGG)	rnd	ribonuclease D
Reads supporting (aligned to +/- strand):  ref base T (29/25);  new base C (6/6);  total (35/31)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.80e-01

CCAGTCGGTGCGCGATTCACTCTTGTCCAGCGTAACGCCGGAATACTCTTCCACCATGGAAGCGAAACCTCATGACATCGGGCGTCCGCAGAAGGCAGCAAGGATTTGCGTGTCAATCAAGGGTTGTGGTAAT  >  minE/1191491‑1191623                                                                      |                                                                ccagTCGGTGCGCGATTCACTCTTGTCCAGCGTAACGCCGGAATACTCTTCCACCATGGAAGCGAAACCTc                                                                <  1:3868063/71‑1 (MQ=255) ccagTCGGTGCGCGATTCACTCTTGTCCAGCGTAACGCCGGAATACTCTTCCACCATGGAAGCGAAACCTc                                                                <  1:2508753/71‑1 (MQ=255)      cGGTGCGCGATTCACTCTTGTCCAGCGTAACGCCGGAATACTCTTCCACCATGGAAGCGAAACCTCATGa                                                            >  1:3252754/1‑70 (MQ=255)      cGGTGCGCGATTCACTCTTGTCCAGCGTAACGCCGGAATACTCTTCCACCATGGAAGCGAAACCTCATGa                                                            >  1:2250582/1‑70 (MQ=255)       ggTGCGCGATTCACTCTTGTCCAGCGTAACGCCGGAATACTCTTCCACCATGGAAGCGAAACCTCATGACa                                                          >  1:571447/1‑71 (MQ=255)        gTGCGCGATTCACTCTTGTCCAGCGTAACGCCGGAATACTCTTCCACCATGGAAGCGAAACCCCATGACa                                                          <  1:4031445/70‑1 (MQ=255)               aTTCACTCTTGTCCAGCGTAACGCCGGAATACTCTTCCACCATGGAAGCGAAACCTCAt                                                              <  1:1400561/59‑1 (MQ=255)               aTTCACTCTTGTCCAGCGTAACGCCGGAATACTCTTCCACCATGGAAGCGAAACCTCATg                                                             <  1:387054/60‑1 (MQ=255)               aTTCACTCTTGTCCAGCGTAACGCCGGAATACTCTTCCACCATGGAAGCGAAACCTCATGACATc                                                        <  1:3593070/65‑1 (MQ=255)                ttCACTCTTGTCCAGCGTAACGCCGGAATACTCTTCCACCATGGAAGCGAAACCTCATGACATCGGGCGTc                                                 <  1:1435529/71‑1 (MQ=255)                 tCACTCTTGTCCAGCGTAACGCCGGAATACTCTTCCACCATGGAAGCGAAACCTCATGa                                                            >  1:87317/1‑59 (MQ=255)                 tCACTCTTGTCCAGCGTAACGCCGGAATACTCTTCCACCATGGAAGCGAAACCTCATGACATCg                                                       >  1:723806/1‑64 (MQ=255)                  cACTCTTGTCCAGCGTAACGCCGGAATACTCTTCCACCATGGAAGCGAAACCCCATGACATCGGGCGTcc                                                >  1:3854360/1‑70 (MQ=255)                    ctctTGTCCAGCGTAACGCCGGAATACTCTTCCACCATGGAAGCGAAACCTCATg                                                             <  1:1590755/55‑1 (MQ=255)                    ctctTGTCCAGCGTAACGCCGGAATACTCTTCCACCATGGAAGCGAAACCTCATGACATCgg                                                      >  1:1956589/1‑62 (MQ=255)                     tctTGTCCAGCGTAACGCCGGAATACTCTTCCACCATGGAAGCGAAACCTc                                                                <  1:1656106/51‑1 (MQ=255)                     tctTGTCCAGCGTAACGCCGGAATACTCTTCCACCATGGAAGCGAAACCCCAt                                                              <  1:2558576/53‑1 (MQ=255)                      ctTGTCCAGCGTAACGCCGGAATACTCTTCCACCATGGAAGCGAAAGCTCATGACATCGGGCGTCCGCAg                                            >  1:836885/1‑70 (MQ=255)                      ctTGTCCAGCGTAACGCCGGAATACTCTTCCACCATGGAAGCGAAACCt                                                                 >  1:3987960/1‑49 (MQ=255)                      ctTGTCCAGCGTAACGCCGGAATACTCTTCCACCATGGAAGCGAAACCTCATGACATCGGGCGTCCGCa                                             >  1:4120970/1‑69 (MQ=255)                      ctTGTCCAGCGTAACGCCGGAATACTCTTCCACCATGGAAGCGAAACCTCATGACATCGGGCGTCCGCAg                                            >  1:860418/1‑70 (MQ=255)                            cAGCGTAACGCCGGAATACTCTTCCACCATGGAAGCGAAACCCCATGACATCGGGCGTCCGCAGAAGgca                                      <  1:2968132/70‑1 (MQ=255)                                 tAACGCCGGAATACTCTTCCACCATGGAAGCGAAACCTCATGACATCGGGCGTCCGCAGAAGGCAGCAAgg                                <  1:161125/71‑1 (MQ=255)                                 tAACGCCGGAATACTCTTCCACCATGGAAGCGAAACCTCATGACATCGGGCGTCCGCAGAAGGCAGCAAgg                                <  1:153139/71‑1 (MQ=255)                                  aaCGCCGGAATACTCTTCCACCATGGAAGCGAAACCTc                                                                <  1:250998/38‑1 (MQ=255)                                   aCGCCGGAATACTCTTCCACCATGGAAGCGAAACCTCATGACATCGGGCGTCCGCAGAAGGCAGCAAGGa                               <  1:1549529/70‑1 (MQ=255)                                   aCGCCGGAATACTCTTCCACCATGGAAGCGAAACCTCATGACATCGGGCGTCCGCAGAAGGCAGCAAGGAt                              >  1:1291643/1‑71 (MQ=255)                                       cGGAATACTCTTCCACCATGGAAGCGAAACCTCATGAc                                                           >  1:2382649/1‑38 (MQ=255)                                         gAATACTCTTCCACCATGGAAGCGAAACCTCATGACATCGGGCGTCCGCAGaa                                          >  1:3504209/1‑53 (MQ=255)                                          aaTACTCTTCCACCATGGAAGCGAAACCTCATGACATCGGGCGTCCGCAGAAGGCAGCAAGGATTTGCgtg                       >  1:1460451/1‑71 (MQ=255)                                            tACTCTTCCACCATGGAAGCGAAACCTCATGACATCgg                                                      <  1:3574815/38‑1 (MQ=255)                                            tACTCTTCCACCATGGAAGCGAAACCCCATGACATCGGGCGTCCGCAGAAGGCAGCAAGGATTTGCGTGtc                     <  1:930391/71‑1 (MQ=255)                                            tACTCTTCCACCATGGAAGCGAAACCCCATGACATCGGGCGTCCGCAGAAGGCAGCAAGGATTTGCGTGtc                     <  1:450531/71‑1 (MQ=255)                                              ctctTCCACCATGGAAGCGAAACCTCATGACATCGGGCGTCCGCAGaa                                          >  1:1850848/1‑48 (MQ=255)                                              ctctTCCACCATGGAAGCGAAACCTCATGACATCGGGCGTCCGCAGaa                                          >  1:3054172/1‑48 (MQ=255)                                              ctctTCCACCATGGAAGCGAAACCTCATGACATCGGGCGTCCGCAGAAGGCAGCAAGGATTTGCGTGtcaa                   >  1:3981045/1‑71 (MQ=255)                                                ctTCCACCATGGAAGCGAAACCCCATGACATCGGGCGTCCGCAGAAGGCAGCAAGGATTTGCGTGtcaat                  >  1:936052/1‑70 (MQ=255)                                                ctTCCACCATGGAAGCGAAACCCCATGACATCGGGCGTCCGCAGAAGGCAGCAAGGATTTGCGTGtcaat                  >  1:2964648/1‑70 (MQ=255)                                                 ttCCACCATGGAAGCGAAACCTCATGACATCGGGCGTCCGCAGAAGGCAGCAAGGATTTGCGTGtcaat                  <  1:352709/69‑1 (MQ=255)                                                 ttCCACCATGGAAGCGAAACCTCATGACATCGGGCGTCCGCAGAAGGCAGCAAGGATTTGCGTGtcaat                  >  1:2415579/1‑69 (MQ=255)                                                   ccaccaTGGAAGCGAAACCTCATGACATCGGGCGTCCGCAGAAGgcagc                                    <  1:654118/49‑1 (MQ=255)                                                   ccaccaTGGAAGCGAAACCTCATGACATCGGGCGTCCGCAGAAGgca                                      >  1:150613/1‑47 (MQ=255)                                                   ccaccaTGGAAGCGAAACCTCATGACATCGGGCGTCCGCAGAAGGCAGCAAGGATTTGCGTGtcaatc                 >  1:3409601/1‑68 (MQ=255)                                                   ccaccaTGGAAGCGAAACCTCATGACATCGGGCGTCCGCAGAAGGCAGCAAGGATTTGCGTGtcaatc                 >  1:277296/1‑68 (MQ=255)                                                   ccaccaTGGAAGCGAAACCTCATGACATCGGGCGTCCGCAGAAGGCAGCAAGGATTTGCGTGTCAATCAAg              >  1:759438/1‑71 (MQ=255)                                                   ccaccaTGGAAGCGAAACCCCATGACATCGGGCGTCCGCAGaa                                          >  1:1942010/1‑43 (MQ=255)                                                   ccaccaTGGAAGCGAAACCCCATGACATCGGGCGTCCGCAGAAGGCAGCAAGGATTTGCGTGTCAATCAAg              >  1:2282588/1‑71 (MQ=255)                                                    caccaTGGAAGCGAAACCTCATGACATCGGGCGTCCGCAGAAGGCAGCAAgg                                <  1:404210/52‑1 (MQ=255)                                                    caccaTGGAAGCGAAACCCCATGACATCGGGCGTCCGCAGAAGGCAGCAAGGATTTGCGTGTCAATCAAg              <  1:1961885/70‑1 (MQ=255)                                                      ccaTGGAAGCGAAACCTCATGACATCGGGCGTCCGCAg                                            <  1:2261676/38‑1 (MQ=255)                                                      ccaTGGAAGCGAAACCTCATGACATCGGGCGTCCGCAGAAGGCAGCAAGGATTTGCGTGtc                     >  1:3457526/1‑61 (MQ=255)                                                       caTGGAAGCGAAACCTCATGACATCGGGCGTCCGCAGAAGGCAGCAAGGATTTGCGTGTCAATCAAggg            <  1:3624921/69‑1 (MQ=255)                                                       caTGGAAGCGAAACCTCATGACATCGGGCGTCCGCAGAAGGCAGCAAGGATTTGCGTGTCAATCAAggg            <  1:2339977/69‑1 (MQ=255)                                                       caTGGAAGCGAAACCTCATGACATCGGGCGTCCGCAGAAGGCAGCAAGGATTTGCGTGTCAATCAAGGGtt          >  1:1720908/1‑71 (MQ=255)                                                        aTGGAAGCGAAACCTCATGACATCTGGCGTCCGCAGaa                                          >  1:2449337/1‑38 (MQ=255)                                                        aTGGAAGCGAAACCTCATGACATCGGGCGTCCGCAGAAGGCAGCAAGGATTTGCGTGTCAATCAAGGGTtg         >  1:1507822/1‑71 (MQ=255)                                                         tGGAAGCGAAACCTCATGACATCGGGCGTCCGCAGAAGGCAGCAAGGATTTGCGTGtcaatc                 <  1:1826449/62‑1 (MQ=255)                                                           gAAGCGAAACCTCATGACATCGGGCGTCCGCAGAAGGCAGCAAGGATTTGCGTGtcaatc                 >  1:499270/1‑60 (MQ=255)                                                           gAAGCGAAACCTCATGACATCGGGCGTCCGCAGAAGGCAGCAAGGATTTGCGTGTCAATCAAGGGTTGTgg      <  1:3675252/71‑1 (MQ=255)                                                           gAAGCGAAACCTCATGACATCGGGCGTCCGCAGAAGGCAGCAAGGATTTGCGTGTCAATCAAGGGTTGTgg      >  1:3317724/1‑71 (MQ=255)                                                               cGAAACCTCATGACATCGGGCGTCCGCAGAAGGCAGCAAGGATTTGCGTGTCAATCAAGGGTTGTGGTAAt  <  1:4004933/71‑1 (MQ=255)                                                                gAAACCTCATGACATCGGGCGTCCGCAGAAGGCAGCAAGGATTTGCGTGtcaa                   >  1:484417/1‑53 (MQ=255)                                                                 aaaCCCCATGACATCGGGCGTCCGCAGAAGGCAGCAAGGATTTGCGTGTCAATCAAGGGTTGTGGTaa   >  1:2347118/1‑68 (MQ=255)                                                                    ccTCATGACATCGGGCGTCCGCAGAAGGCAGCAAGGAttt                            <  1:1032928/40‑1 (MQ=255)                                                                    ccTCATGACATCGGGCGTCCGCAGAAGGCAGCAAGGATTTGCg                         >  1:1218305/1‑43 (MQ=255)                                                                    ccTCATGACATCGGGCGTCCGCAGAAGGCAGCAAGGATTTGCGTGTCAATCAAGGGTTGTGGTaa   <  1:696459/65‑1 (MQ=255)                                                                     cTCATGACATCGGGCGTCCGCAGAAGGCAGCAAGGATTTGCGTGtcaa                   <  1:1853458/48‑1 (MQ=255)                                                                      tCATGACATCGGGCGTCCGCAGAAGGCAGCAAGGATTTGCGTGtca                    <  1:1421783/46‑1 (MQ=255)                                                                      |                                                                CCAGTCGGTGCGCGATTCACTCTTGTCCAGCGTAACGCCGGAATACTCTTCCACCATGGAAGCGAAACCTCATGACATCGGGCGTCCGCAGAAGGCAGCAAGGATTTGCGTGTCAATCAAGGGTTGTGGTAAT  >  minE/1191491‑1191623

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