Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	2,086,490	C→A	9.2%	P89T (CCG→ACG)	yqjG →	predicted S‑transferase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	2,086,490	0	C	A	9.2%	183.3 / 4.3	54	P89T (CCG→ACG)	yqjG	predicted S‑transferase
Reads supporting (aligned to +/- strand):  ref base C (20/29);  new base A (3/2);  total (23/31)
Fisher's exact test for biased strand distribution p-value = 6.40e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.36e-01

CACCGCACGCTGATCATGCGCAAACTCAAAGGACTGGAACCGTTTATTTCCGTTTCCGTAGTGAACCCGCTGATGCTGGAAAACGGTTGGACCTTTGATGACAGTTTTCCGGGAGCAACCGGAGACACGCTCTATCA  >  minE/2086424‑2086560                                                                   |                                                                       caccgcacGCTGATCATGCGCAAACTCAAAGGACTGGAACCGTTTATTTCCGTTTCCGTAGTGAACCCGCt                                                                    <  1:2370768/71‑1 (MQ=255)    cgcacGCTGATCATGCGCAAACTCAAAGGACTGGAACCGTTTATTTCCGTTTCCGTAGTGAACACGCTGAt                                                                 >  1:1522650/1‑71 (MQ=255)      cacGCTGATCATGCGCAAACTCAAAGGACTGGAACCGTTTATTTCCGTTTCCGTAGTGAACCCg                                                                      >  1:2845204/1‑64 (MQ=255)      cacGCTGATCATGCGCAAACTCAAAGGACTGGAACCGTTTATTTCCGTTTCCGTAGTGAACCCGCTGATGc                                                               >  1:1479561/1‑71 (MQ=255)      cacGCTGATCATGCGCAAACTCAAAGGACTGGAACCGTTTATTTCCGTTTCCGTAGTGAACCCGCTGATGc                                                               >  1:3151165/1‑71 (MQ=255)      cacGCTGATCATGCGCAAACTCAAAGGACTGGAACCGTTTATTTCCGTTTCCGTAGTGAACCCGCTGATGc                                                               >  1:1543364/1‑71 (MQ=255)      cacGCTGATCATGCGCAAACTCAAAGGACTGGAACCGTTTATTTCCGTTTCCGTAGTGAACACGCTGATGc                                                               >  1:3265435/1‑71 (MQ=255)             aTCATGCGCAAACTCAAAGGACTGGAACCGTTTATTTCCGTTTCCGTAGTGAACCCGCTGATGCTGGaaaa                                                        >  1:3109157/1‑71 (MQ=255)             aTCATGCGCAAACTCAAAGGACTGGAACCGTTTATTTCCGTTTCCGTAGTGAACCCGCTGATGCTGGaaaa                                                        >  1:1815035/1‑71 (MQ=255)                aTGCGCAAACTCAAAGGACTGGAACCGTTTATTTCCGTTTCCGTAGTGAACCCGCt                                                                    <  1:120371/56‑1 (MQ=255)                aTGCGCAAACTCAAAGGACTGGAACCGTTTATTTCCGTTTCCGTAGTGAACCCGCTGATGCTGGAAAACg                                                      <  1:2933196/70‑1 (MQ=255)                aTGCGCAAACTCAAAGGACTGGAACCGTTTATTTCCGTTTCCGTAGTGAACCCGCTGATGCTGGAAAACg                                                      <  1:288568/70‑1 (MQ=255)                aTGCGCAAACTCAAAGGACTGGAACCGTTTATTTCCGTTTCCGTAGTGAACCCGCTGATGCTGGAAAACg                                                      <  1:1291459/70‑1 (MQ=255)                aTGCGCAAACTCAAAGGACTGGAACCGTTTATTTCCGTTTCCGTAGTGAACACGCTGATGCTGGAAAACgg                                                     <  1:1846126/71‑1 (MQ=255)                 tGCGCAAACTCAAAGGACTGGAACCGTTTATTTCCGTTTCCGTAGTGAACCCGCTGATGCTGGAAAACGGt                                                    <  1:2983661/71‑1 (MQ=255)                  gcgcAAACTCAAAGGACTGGAACCGTTTATTTCCGTTTCCGTAGTGAACCCGCTGAt                                                                 <  1:32853/57‑1 (MQ=255)                  gcgcAAACTCAAAGGACTGGAACCGTTTATTTCCGTTTCCGTAGTGAACCCGCTGAt                                                                 <  1:3157304/57‑1 (MQ=255)                  gcgcAAACTCAAAGGACTGGAACCGTTTATTTCCGTTTCCGTAGTGAACCCGCTGAt                                                                 <  1:675006/57‑1 (MQ=255)                   cgcAAACTCAAAGGACTGGAACCGTTTATTTCCGTTTCCGTAGTGAACCCGCTGATGCTGGAAAACGGTTg                                                  <  1:681915/71‑1 (MQ=255)                    gcAAACTCAAAGGACTGGAACCGTTTATTTCCGTTTCCGTAGTGAACCCGCTGATGc                                                               >  1:3112408/1‑57 (MQ=255)                       aaCTCAAAGGACTGGAACCGTTTATTTCCGTTTCCGTAGTGAACCCGCTGATGCTGGAAAACGGTTGGAcc                                              <  1:849654/71‑1 (MQ=255)                         cTCAAAGGACTGGAACCGTTTATTTCCGTTTCCGTAGTGAACCAGCTGATGCTGGAAAACGGTTGGACCtt                                            >  1:568762/1‑71 (MQ=255)                         cTCAAAGGACTGGAACCGTTTATTTCCGTTTCCGTAGTGAACACGCTGATGCTGGAAAACGGTTGGACCtt                                            >  1:1793173/1‑71 (MQ=255)                          tCAAAGGACTGGAACCGTTTATTTCCGTTTCCGTAGTGAACCCGCTGATGCTGGAAAACg                                                      <  1:562204/60‑1 (MQ=255)                            aaaGGACTGGAACCGTTTATTTCCGTTTCCGTAGTGAACCCGCTGATGCTGGAAAACGGTTGGACCTTtga                                         <  1:2446023/71‑1 (MQ=255)                             aaGGACTGGAACCGTTTATTTCCGTTTCCGTAGTGAACCCGCTGATGCTGGAAAACGGTTGGACCTTtgat                                        <  1:1848581/71‑1 (MQ=255)                               ggACTGGAACCGTTTATTTCCGTTTCCGTAGTGAACACGCTGATGCTGGAAAACGGTTGGACCTTtgatga                                      <  1:2978072/71‑1 (MQ=255)                                gACTGGAACCGTTTATTTCCGTTTCCGTAGTGAACCGGCTGATGCTGGAAAACGGTTGGACCTTTGATGAc                                     <  1:473993/71‑1 (MQ=255)                                     gAACCGTTTATTTCCGTTTCCGTAGTGAACCCGCTGATGCTGGAAAACGGTTGGACCTTTGATGACAGttt                                >  1:1635582/1‑71 (MQ=255)                                        ccGTTTATTTCCGTTTCCGTAGTGAACCCGCTGATGCTGGAAAACGGTTGGACCTTTGATGACAGTTTTcc                             <  1:785258/71‑1 (MQ=255)                                        ccGTTTATTTCCGTTTCCGTAGTGAACCCGCTGATGCTGGAAAACGGTTGGACCTTTGATGACAGTTTTcc                             <  1:1186043/71‑1 (MQ=255)                                         cGTTTATTTCCGTTTCCGTAGTGAACCCGCTGATGCTGGAAAACGGTTGGACCTTTGATGACAGTTTTCCg                            <  1:296968/71‑1 (MQ=255)                                          gTTTATTTCCGTTTCCGTAGTGAACCCGCTGATGCTGGAAAACGGTTGGACCTTTGATGACAGTTTTCCgg                           <  1:1984626/71‑1 (MQ=255)                                          gTTTATTTCCGTTTACGTAGTGAACCCGCTGATGCTGGAAAACGGTTGGACCTTtg                                          >  1:1296458/1‑56 (MQ=255)                                             tATTTCCGTTTCCGTAGTGAACCCGCTGATGCTGGaaaa                                                        >  1:867820/1‑39 (MQ=255)                                             tATTTCCGTTTCCGTAGTGAACCCGCTGATGCTGGAAAACGGTTGGACCTTTGATGACAGTTTTCCGGGa                         <  1:2476354/70‑1 (MQ=255)                                               tttccgtttccgTAGTGAACCCGCTGATGCTGGAAAACGGTTGGACCTTTGATGACAGTTTTCCGGGAGc                       <  1:1863328/70‑1 (MQ=255)                                               tttccgtttccgTAGTGAACCCGCTGATGCTGGAAAACGGTTGGACCTTTGATGACAGTTTTCCGGGAGc                       <  1:838721/70‑1 (MQ=255)                                               tttccgtttccgTAGTGAACCCGCTGATGCTGGAAAACGGTTGGACCTTTGATGACAGTTTTCCGGGAGc                       <  1:840861/70‑1 (MQ=255)                                               tttccgtttccgTAGTGAACCCGCTGATGCTGGAAAACGGTTGGACCTTTGATGACAGTTTTCCGGGAGc                       <  1:3036281/70‑1 (MQ=255)                                                ttccgtttccgtAGTGAACCCGCTGATGCTGGAAAACGGTTGGACCTTTGATGACAGTTTTCCGGGAGCaa                     >  1:3274486/1‑71 (MQ=255)                                                ttccgtttccgtAGTGAACCCGCTGATGCTGGAAAACGGTTGGACCTTTGATGACAGTTTTCCGGGAGCaa                     >  1:858499/1‑71 (MQ=255)                                                 tccgtttccgtAGTGAACCCGCTGATGCTGGAAAACGGTTGGACCTTTGATGACAGTTTTCCGGGAGCaa                     >  1:1980608/1‑70 (MQ=255)                                                 tccgtttccgtAGTGAACCCGCTGATGCTGGAAAACGGTTGGACCTTTGATGACAGTTTTCCGGGAGCaa                     >  1:475698/1‑70 (MQ=255)                                                     tttccgtAGTGAACCCGCTGATGCTGGAAAACGGTTgg                                                 <  1:1006749/38‑1 (MQ=255)                                                        ccgtAGTGAACCCGCTGATGCTGGAAAACGGTTGGACCTTTGATGACAGTTTTCCGGGAGCAACCGgaga              >  1:1951435/1‑70 (MQ=255)                                                        ccgtAGTGAACCCGCTGATGCTGGAAAACGGTTGGACCTTTGATGACAGTTTTCCGGGAGCAACCGgaga              >  1:45707/1‑70 (MQ=255)                                                          gtAGTGAACCCGCTGATGCTGGAAAACGGTTGGACCTTTGATGACAGTTTTCCGGGAGCAACCGGAGacac           >  1:2485364/1‑71 (MQ=255)                                                            aGTGAACCCGCTGATGCTGGAAAACGGTTGGACCTTTGATGACAGTTTTCCGGGAGCAACCGGAGACACGc         >  1:1479228/1‑71 (MQ=255)                                                              tGAACCCGCTGATGCTGGAAAACGGTTGGACCTTTGATGACAGTTTTCCGGGAGCAACCg                  <  1:1722262/60‑1 (MQ=255)                                                              tGAACCCGCTGATGCTGGAAAACGGTTGGACCTTTGATGACAGTTTTCCGGGAGCAACCg                  <  1:879832/60‑1 (MQ=255)                                                              tGAACCCGCTGATGCTGGAAAACGGTTGGACCTTTGATGACAGTTTTCCGGGAGCAACCg                  <  1:1514365/60‑1 (MQ=255)                                                              tGAACCCGCTGATGCTGGAAAACGGTTGGACCTTTGATGACAGTTTTCCGGGAGCAACCg                  <  1:263930/60‑1 (MQ=255)                                                                aaCCCGCTGATGCTGGAAAACGGTTGGACCTTTGATGACAGtttt                               >  1:2690770/1‑45 (MQ=255)                                                                   ccGCTGATGCTGGAAAACGGTTGGACCTTTGATGACAGTTTTCCGGGAGCAACCGga                <  1:1558629/57‑1 (MQ=255)                                                                   ccGCTGATGCTGGAAAACGGTTGGACCTTTGATGACAGTTTTCCGGGAGCAACCGGAGACACGCTCTATCa  <  1:1576298/71‑1 (MQ=255)                                                                   |                                                                       CACCGCACGCTGATCATGCGCAAACTCAAAGGACTGGAACCGTTTATTTCCGTTTCCGTAGTGAACCCGCTGATGCTGGAAAACGGTTGGACCTTTGATGACAGTTTTCCGGGAGCAACCGGAGACACGCTCTATCA  >  minE/2086424‑2086560

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