Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	2,387,702	G→A	14.7%	L238L (CTG→TTG)	yigL ←	predicted hydrolase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	2,387,702	0	G	A	14.7%	257.7 / 30.9	95	L238L (CTG→TTG)	yigL	predicted hydrolase
Reads supporting (aligned to +/- strand):  ref base G (46/35);  new base A (8/6);  total (54/41)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.82e-01

CAGATAATGCGGCACCGCGTCGTCGGCATTAGTACCAATCACTTCCAGCTCGGGATGAAGGTCTTTCAGACGCTGGTGCGCACTGCCCATAATGCAGCCTTTCCCCGCCATCGACAGCATTTCGGCGTCGTTCATCCCG  >  minE/2387634‑2387772                                                                     |                                                                       cAGATAATGCGGCACCGCGTCGTCGGCATTAGTACCAATCACTTCCAGCTCGGGATGAAGGTCTTTCAAAc                                                                      <  1:1655267/71‑1 (MQ=255)  aGATAATGCGGCACCGCGTCGTCGGCATTAGTACCAATCACTTCCAGCTCGGGATGAAGGTCTTTCAGACg                                                                     <  1:2829746/71‑1 (MQ=255)    aTAATGCGGCACCGCGTCGTCGGCATTAGTACCAATCACTTCCAGCTCGGGATGAAGGTCTTTCAGACGCt                                                                   <  1:3564549/71‑1 (MQ=255)     tAATGCGGCACCGCGTCGTCGGCATTAGTACCAATCACTTCCCGCTCGGGATGAAGGTCTTTCAGACGCTg                                                                  >  1:354531/1‑71 (MQ=255)      aaTGCGGCACCGCGTCGTCGGCATTAGTACCAATCACTTCCAGCTCGGGATGAAGGTCTTTCAGACGCTgg                                                                 >  1:1700536/1‑71 (MQ=255)      aaTGCGGCACCGCGTCGTCGGCATTAGTACCAATCACTTCCAGCTCGGGATGAAGGTCTTTCAAACGCTgg                                                                 >  1:3730858/1‑71 (MQ=255)      aaTGCGGCACCGCGTCGTCGGCATTAGTACCAATCACTTCCAGCTCGGGATGAAGGTCTTTCAAACGCTgg                                                                 >  1:27304/1‑71 (MQ=255)           ggCACCGCGTCGTCGGCATTAGTACCAATCACTTCCAGCTCGGGATGAAGGTCTTTCAGACGCTGGTgcgc                                                            <  1:718719/71‑1 (MQ=255)           ggCACCGCGTCGTCGGCATTAGTACCAATCACTTCCAGCTCGGGATGAAGGTCTTTCAGACGCTGGTgcg                                                             >  1:3790815/1‑70 (MQ=255)            gCACCGCGTCGTCGGCATTAGTACCAATCACTTCCAGCTCGGGATGAAGGTCTTTCAGACGCTGGTGCGCa                                                           <  1:2044316/71‑1 (MQ=255)               ccGCGTCGTCGGCATTAGTACCAATCACTTCCAGCTCGGGATGAAGGTCTTTCAAACGCTGGTGCGCACTg                                                        >  1:3165869/1‑71 (MQ=255)                    tcgtcgGCATTAGTACCAATCACTTCCAGCTCGGGATGAAGGTCTTTCAGACg                                                                     <  1:2554966/53‑1 (MQ=255)                    tcgtcgGCATTAGTACCAATCACTTCCAGCTCGGGATGAAGGTCTTTCAGACGCTGGTGCGCACTGCCCAt                                                   <  1:2507613/71‑1 (MQ=255)                    tcgtcgGCATTAGTACCAATCACTTCCAGCTCGGGATGAAGGTCTTTCAGACGCTGGTGCGCACTGCCCAt                                                   <  1:2562827/71‑1 (MQ=255)                        cgGCATTAGTACCAATCACTTCCAGCTCGGGATGAAGGTCTTTCAGACGCTGGTGCGCACTGCCCATAATg                                               >  1:2147197/1‑71 (MQ=255)                         ggCATTAGTACCAATCACTTCCAGCTCGGGATGAAGGTCTTTCAGACGCTGGTGCGCACTGCCCAt                                                   <  1:2765235/66‑1 (MQ=255)                         ggCATTAGTACCAATCACTTCCAGCTCGGGATGAAGGTCTTTCAGACGCTGGTGCGCACTGCCCATAATg                                               <  1:3399928/70‑1 (MQ=255)                         ggCATTAGTACCAATCACTTCCAGCTAGGGATGAAGGTCTTTCAGACGCTGGTGCGCACTg                                                        >  1:2246101/1‑61 (MQ=255)                            aTTAGTACCAATCACTTCCAGCTCGGGATGAAGGTCTTTCAGACGCTGGTGCGCACTGCCCATAATGCAGc                                           >  1:2764146/1‑71 (MQ=255)                             ttAGTACCAATCACTTCCAGCTCGGGATGAAGGTCTTTCAGa                                                                       >  1:3324619/1‑42 (MQ=255)                             ttAGTACCAATCACTTCCAGCTCGGGATGAAGGTCTTTCAGACGCTGGTgcg                                                             <  1:1165181/52‑1 (MQ=255)                             ttAGTACCAATCACTTCCAGCTCGGGATAAAGGTCTTTCAGACGCTGGTGCGCACTGCCCATAATGCAGcc                                          >  1:2229175/1‑71 (MQ=255)                              tAGTACCAATCACTTCCAGCTCGGGATGAAGGTCTTTCAAACGCTGGTGCGCACTGCCCATAATGCAGCCt                                         >  1:2626171/1‑71 (MQ=255)                              tAGTACCAATCACTTCCAGCTCGGGATGAAGGTCTTTCAAACGCTGGTGCGCACTGCCCATAATGCAGCCt                                         >  1:54222/1‑71 (MQ=255)                               aGTACCAATCACTTCCAGCTCGGGATGAAGGTCTTTCAGACGCTGGTGCGCACTGCCCATaa                                                 >  1:912540/1‑62 (MQ=255)                               aGTACCAATCACTTCCAGCTCGGGATGAAGGTCTTTCAGACGCTGGTGCGCACTGCCCATaa                                                 >  1:2819749/1‑62 (MQ=255)                               aGTACCAATCACTTCCAGCTCGGGATGAAGGTCTTTCAGACGCTGGTGCGCACTGCCCATAATGCAGCCtt                                        >  1:1699710/1‑71 (MQ=255)                                gTACCAATCACTTCCAGCTCGGGATGAAGGTCTTTCAGACGCTGGTGCGCACTGCCCATAATGCAGCCtt                                        <  1:2413099/70‑1 (MQ=255)                                 tACCAATCACTTCCAGCTCGGGATGAATGTCTTTCAAACGCt                                                                   <  1:2375132/42‑1 (MQ=255)                                 tACCAATCACTTCCAGCTCGGGATGAAGGTCTTTCAGACGCTGGTGCGCACTGCCCATAATGCAGCCTTTc                                      <  1:2257660/71‑1 (MQ=255)                                 tACCAATCACTTCCAGCTCGGGATGAAGGTCTTTCAAACGCt                                                                   <  1:1979568/42‑1 (MQ=255)                                   ccAATCACTTCCAGCTCGGGATGGAGGTCTTTCAGACGCt                                                                   >  1:3773054/1‑40 (MQ=255)                                   ccAATCACTTCCAGCTCGGGATGAAGGTCTTTCAGACGCTGGTGCTCACTGCCCATAATGCAGCCTTTccc                                    <  1:3828844/71‑1 (MQ=255)                                   ccAATCACTTCCAGCTCGGGATGAAGGTCTTTCAGACGCTGGTGCGCACTGCCCATAATGCa                                             >  1:260662/1‑62 (MQ=255)                                   ccAATCACTTCCAGCTCGGGATGAAGGTCTTTCAGACGCTGGTGCGCACTGCCCATAATGCAGCCTTTccc                                    >  1:935922/1‑71 (MQ=255)                                   ccAATCACTTCCAGCTCGGGATGAAGGTCTTTCAGACGCTGGTGCGCACTGCCCATAATGCAGCCTTTccc                                    >  1:3887524/1‑71 (MQ=255)                                   ccAATCACTTCCAGCTCGGGATGAAGGTCTTTCAGACGCTGGTGCGCACTGCCCATAATGCAGCCTTTccc                                    >  1:3562129/1‑71 (MQ=255)                                   ccAATCACTTCCAGCTCGGGATGAAGGTCTTTCAGACGCTGGTGCGCACTGCCCACAATGCAGCCTTTcc                                     >  1:2360927/1‑70 (MQ=255)                                   ccAATCACTTCCAGCTCGGGATGAAGGTCTTTCAAACGCTgg                                                                 >  1:1481373/1‑42 (MQ=255)                                   ccAATCACTTCCAGCTCGGGATGAAGGTCTTTCAAACGCTgg                                                                 >  1:2732340/1‑42 (MQ=255)                                      aTCACTTCCAGCTCGGGATGAAGGTCTTTCAGACGCTGGTGCGCACTGCCCATAATGCAGCCTTTCCCCGc                                 >  1:4078157/1‑71 (MQ=255)                                      aTCACTTCCAGCTCGAGATGAAGGTCTTTCAGACGCTGGTGCGCACTGCCCATAATGCAGCCTTTCCCCGc                                 >  1:3023229/1‑71 (MQ=255)                                       tCACTTCCAGCTCGGGATGAAGGTCTTTCAGACGCTGGTGCGCACTGCCCAt                                                   <  1:3265339/52‑1 (MQ=255)                                       tCACTTCCAGCTCGGGATGAAGGTCTTTCAGACGCTGGTGCGCACTGCCCATAATGCAGCCTTTCCCCGcc                                <  1:2196150/71‑1 (MQ=255)                                       tCACTTCCAGCTCGGGATGAAGGTCTTTCAGACGCTGGTGCGCACTGCCCATAATGCAGCCTTTCCCCGc                                 >  1:3414385/1‑70 (MQ=255)                                        cACTTCCAGCTCGGGATGAAGGTCTTTCAAACCCTGGTGCGCACTGCCCATAATGc                                              >  1:2826675/1‑56 (MQ=255)                                         aCTTCCAGCTCGGGATGAAGGTCTTTCAGACGCTGGTGCGCACTGCCCATAATGCAGCCTTTc                                      >  1:1923286/1‑63 (MQ=255)                                         aCTTCCAGCTCGGGATGAAGGTCTTTCAGACGCTGGTGCGCACTGCCCATAATGCAGCCTTTCCCCGCCAt                              >  1:1377181/1‑71 (MQ=255)                                          cTTCCAGCTCGGGATGAAGGTCTTTCAAACGCTGGTGCGCACTGCCCATAATGCAGCCTTTCCCCGCCa                               <  1:49747/69‑1 (MQ=255)                                           ttCCAGCTCGGGATGAAGGTCTTTCAGACGCTGGTGCGCACTGCCCATAATGCAGCCTTTCCCCGc                                 <  1:2780468/66‑1 (MQ=255)                                           ttCCAGCTCGGGATGAAGGTCTTTCAGACGCTGGTGCGCACTGCCCATAATGCAGCCTTTCCCCGCCATc                             >  1:2169011/1‑70 (MQ=255)                                            tCCAGCTCGGGATGAAGGTCTTTCAGACGCTGGTGCGCACTGCCCATAATGCa                                             >  1:2844429/1‑53 (MQ=255)                                            tCCAGCTCGGGATGAAGGTCTTTCAGACGCTGGTGCGCACTGCCCATAATGCAGCCTTTCCCCGc                                 <  1:3267713/65‑1 (MQ=255)                                            tCCAGCTCGGGATGAAGGTCTTTCAGACGCTGGTGCGCACTGCCCATAATGCAGCCTTTCCCCGCCATCGa                           >  1:1361599/1‑71 (MQ=255)                                            tCCAGCTCGGGATGAAGGTCTTTCAGACGCTGGTGCGCACTGCCCATAATGCAGCCTTTCCCCGCCATCGa                           >  1:3519437/1‑71 (MQ=255)                                              cAGCTCGGGATGAAGGTCTTTCAGACGCTGGTGCGCACTGCCCATAATGc                                              >  1:1851052/1‑50 (MQ=255)                                              cAGCTCGGGATGAAGGTCTTTCAGACGCTGGTGCGCACTGCCCATAATGc                                              >  1:1690559/1‑50 (MQ=255)                                              cAGCTCGGGATGAAGGTCTTTCAGACGCTGGTGCGCACTGCCCATAATGc                                              >  1:1933299/1‑50 (MQ=255)                                              cAGCTCGGGATGAAGGTCTTTCAGACGCTGGTGCGCACTGCCCATAATGc                                              >  1:3456815/1‑50 (MQ=255)                                              cAGCTCGGGATGAAGGTCTTTCAGACGCTGGTGCGCACTGCCCATAATGCAGCCTTTCCCCGCCATCGa                           >  1:2944695/1‑69 (MQ=255)                                              cAGCTCGGGATGAAGGTCTTTCAGACGCTGGTGCGCACTGCCCATAATGCAGCCTTTCCCCGCCATCGa                           >  1:3588545/1‑69 (MQ=255)                                              cAGCTCGGGATGAAGGTCTTTCAGACGCTGGTGCGCACTGCCCATAATGCAGCCTTTCCCCGCCATCGa                           >  1:3251669/1‑69 (MQ=255)                                              cAGCTCGGGATGAAGGTCTTTCAAACGCTGGTGCGCACTGCCCATAATGCAGCCTTTCCCCGCCATCGAc                          <  1:746538/70‑1 (MQ=255)                                                 cTCGGGATGAAGGTCTTTCAGACGCTGGTGCGCACTGCCCATAATGCAGCCTTTCCCCGCCATCGACAGCa                      >  1:340485/1‑71 (MQ=255)                                                    gggATGAAGGTCTTTCAGACGCTGGTGCGCACTGCCCATaa                                                 >  1:1411314/1‑41 (MQ=255)                                                    gggATGAAGGTCTTTCAGACGCTGGTGCGCACTGCCCATAATGCAGCCTTTCCCCGCCATCGACAGCAttt                   <  1:3985793/71‑1 (MQ=255)                                                    gggATGAAGGTCTTTCAGACGCTGGTGCGCACTGCCCATAATGCAGCCTTTCCCCGCCATCGACAGCAtt                    <  1:2183171/70‑1 (MQ=255)                                                    gggATGAAGGTCTTTCAGACGCTGGTGCGCACTGCCCATAATGCAGCCTTTCCCCGCCATCGACAGCAt                     <  1:475068/69‑1 (MQ=255)                                                    gggATGAAGGTCTTTCAGACGCTGGTGCGCACTGCCCATAATGCAGCCTTTCCCCGCCATCGACAGCAt                     <  1:1379830/69‑1 (MQ=255)                                                    gggATGAAGGGCTTTCAGACGCTGGTGCGCACTGCCCATAATGCAGCCTTTCCCCGCCATCGACAGCAtt                    <  1:1263766/70‑1 (MQ=255)                                                      gATGAAGGTCTTTCAGACGCTGGTGCGCACTGCCCATAATGCAGCCTTTCCCCGCCATCGACAGCATTTCg                 >  1:439963/1‑71 (MQ=255)                                                        tGAAGGTCTTTCAGACGCTGGTGCGCACTGCCCATAATg                                               <  1:3834295/39‑1 (MQ=255)                                                        tGAAGGTCTTTCAGACGCTGGTGCGCACTGCCCATAATGCAGCCTTTCCCCGCCATCGa                           <  1:2751040/59‑1 (MQ=255)                                                        tGAAGGTCTTTCAGACGCTGGTGCGCACTGCCCATAATGCAGCCTTTCCCCGCCATCGACAGCATTTCGGc               <  1:3412566/71‑1 (MQ=255)                                                        tGAAGGTCTTTCAGACGCTGGTGCGCACTGCCCATAATGCAGCCTTTCCCCGCCATCGACAGCATTTCGGc               <  1:2473033/71‑1 (MQ=255)                                                        tGAAGGTCTTTCAGACGCTGGTGCGCACTGCCCATAATGCAGCCTTTCCCCGCCATCGACAGCATTTCGGc               <  1:1972077/71‑1 (MQ=255)                                                         gAAGGTCTTTCAGACGCTGGTGCGCACTGCCCATAATGCAGCCTTTCCCCGCCATCGACAGCATTTCGGcg              >  1:3325807/1‑71 (MQ=255)                                                         gAAGGTCTTTCAGACGCTGGTGCGCACTGCCCATAATGCAGCCTTTCCCCGCCATCGACAGCATTTCGGcg              >  1:1060279/1‑71 (MQ=255)                                                             gTCTTTCAGACGCTGGTGCGCACTGCCCATAATGCAGCCTTTCCCCGCCATCGACAGCAt                     <  1:3652051/60‑1 (MQ=255)                                                             gTCTTTCAGACGCTGGTGCGCACTGCCCATAATGCAGCCTTTCCCCGCCATCGACAGCAt                     <  1:1432326/60‑1 (MQ=255)                                                             gTCTTTCAGACGCTGGTGCGCACTGCCCATAATGCAGCCTTTCCCCGCCATCGACAGCATTTCGGcgtcgt          <  1:1565650/71‑1 (MQ=255)                                                              tCTTTCAGACGCTGGTGCGCACTGCCCATAATGCAGCCTTTCCCCGCCATCGACAGCAt                     <  1:174577/59‑1 (MQ=255)                                                              tCTTTCAGACGCTGGTGCGCACTGCCCATAATGCAGCCTTTCCCCGCCATCGACAGCATTTCg                 >  1:132008/1‑63 (MQ=255)                                                               cTTTCAGACGCTGGTGCGCACTGCCCATAATGCAGCCTTTCCCCGCCATCGACAGCATTTCGGCGTCGtt         >  1:3186899/1‑70 (MQ=255)                                                                tttCAGACGCTGGTGCGCACTGCCCATAATGCAGCCTTTCCCCGCCa                               <  1:3827021/47‑1 (MQ=255)                                                                tttCAGACGCTGGTGCGCACTGCCCATAATGCAGCCTTTCCCCGCCa                               <  1:1858574/47‑1 (MQ=255)                                                                tttCAGACGCTGGTGCGCACTGCCCATAATGCAGCCTTTCCCCGCCATCGACAGCATTTCGGCGTCGTTCa       >  1:4072444/1‑71 (MQ=255)                                                                  tCAGACGCTGGTGCGCACTGCCCATAATGCAGCCTTTCCCCGCCATCGACAGCATTTCGGCGTCGTTc        <  1:1735491/68‑1 (MQ=255)                                                                  tCAGACGCTGGTGCGCACTGCCCATAATGCAGCCTTTCCCCGCCATCGACAGCATTTCGGCGTCGTTCATc     >  1:2033135/1‑71 (MQ=255)                                                                  tCAGACGCTGGTGCGCACTGCCCATAATGCAGCCTTTCCCCGCCATCGACAGCATTTCGGCGTCGTTCATc     >  1:2433847/1‑71 (MQ=255)                                                                  tCAGACGCTGGTGCGCACTGCCCATAATGCAGCCTTTCCCCGCCATCGACAGCATTTCGGCGTCGTTCATc     >  1:509280/1‑71 (MQ=255)                                                                  tCAGACGCTGGTGCGCACTGCCCATAATGCAGCCTTTCCCCGCCATCGACAGCATTTCGGCGTCGTTCATc     >  1:3249798/1‑71 (MQ=255)                                                                  tCAGACGCTGGTGCGCACTGCCCATAATCCAGCCTTTCCCCGCCATCGACAGCATTTCGGCGTCGTTCa       <  1:3892453/69‑1 (MQ=255)                                                                  tCAAACGCTGGTGCGCACTGCCCATAATGCAGCCTTTCCCCGCCATCGACAGCATTTCGGCGTCGTTc        <  1:300256/68‑1 (MQ=255)                                                                    aGACGCTGGTGCGCACTGCCCATAATGCAGCCTTTCCCCGCCATCGACAGCATTTCGGCGTCGTTCATccc   >  1:3281847/1‑71 (MQ=255)                                                                     gACGCTGGTGCGCACTGCCCATAATGCAGCCTTTCCCCGCCATCGACAGCATTTCGGCGTCGTTCATCCCg  <  1:3563095/71‑1 (MQ=255)                                                                     gACGCTGGTGCGCACTGCCCATAATGCAGCCTTTCCCCGCCATCGACAGCATTTCGGCGTCGTTCATCCCg  <  1:240/71‑1 (MQ=255)                                                                     gACGCTGGTGCGCACTGCCCATAATGCAGCCTTTCCCCGCCATCGACAGCATTTCGGCGTCGTTCATCCCg  <  1:770987/71‑1 (MQ=255)                                                                     |                                                                       CAGATAATGCGGCACCGCGTCGTCGGCATTAGTACCAATCACTTCCAGCTCGGGATGAAGGTCTTTCAGACGCTGGTGCGCACTGCCCATAATGCAGCCTTTCCCCGCCATCGACAGCATTTCGGCGTCGTTCATCCCG  >  minE/2387634‑2387772

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