Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	909,913:1	+T	35.4%	coding (2090/2268 nt)	ycjT →	predicted hydrolase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	909,913	1	.	T	35.4%	61.5 / 117.5	65	coding (2090/2268 nt)	ycjT	predicted hydrolase
Reads supporting (aligned to +/- strand):  ref base . (26/16);  new base T (8/15);  total (34/31)
Fisher's exact test for biased strand distribution p-value = 4.25e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.23e-01

GTGAATTGCATCTCAATCCGGCGTTACCTGAGCAGTGGCAACAGTTGTCTTTCCCTCTGTTCTGGCA‑GGGCTGCGAATTACAGGTCACTCTTGACGCGCAGCGTATTGCGATTCGAACTTCTGCGCCCGTT  >  minE/909847‑909977                                                                    |                                                                 gTGAATTGCATCTCAATCCGGCGTTACCTGAGCAGTGGCAACAGTTGTCTTTCCCTCTGTTCTGGCA‑tgg                                                               >  1:3429183/1‑70 (MQ=255)   gAATTGCATCTCAATCCGGCGTTACCTGAGCAGTGGCAACAGTTGTCTTTCCCTCTGTTCTGGCA‑GGGCTg                                                            >  1:2395261/1‑71 (MQ=255)    aaTTGCATCTCAATCCGGCGTTACCTGAGCAGTGGCAACAGTTGTCTTTCCCTCTGTTCTGGCATGGGCt                                                             <  1:1828232/70‑1 (MQ=255)    aaTTGCATCTCAATCCGGCGTTACCTGAGCAGTGGCAACAGTTGTCTTTCCCTCTGTTCTGGCA‑GGGc                                                              >  1:1155243/1‑68 (MQ=255)     aTTGCATCTCAATCCGGCGTTACCTGAGCAGTGGCAACAGTTGTCTTTCCCTCTGTTCTGGCA‑GGTCTGCg                                                          <  1:1181408/71‑1 (MQ=255)     aTTGCATCTCAATCCGGCGTTACCTGAGCAGTGGCAACAGTTGTCTTTCCCTCTGTTCTGGCA‑GGGCTGc                                                           >  1:3814872/1‑70 (MQ=255)           tctcAATCCGGCGTTACCTGAGCAGTGGCAACAGTTGTCTTTCCCTCTGTTCTGGCATGGGCTGCGAAtt                                                      <  1:1132580/70‑1 (MQ=255)             tcAATCCGGCGTTACCTGAGCAGTGGCAACAGTTGTCTTTCCCTCTGTTCTGGCA‑GGGCt                                                             >  1:865769/1‑60 (MQ=255)             tcAATCCGGCGTTACCTGAGCAGTGGCAACAGTTGTCTTTCCCTCTGTTCTGGCA‑GGGCt                                                             >  1:631894/1‑60 (MQ=255)             tcAATCCGGCGTTACCTGAGCAGTGGCAACAGTTGTCTTTCCCTCTGTTCTGGCA‑GGGCTGCGAATTa                                                     >  1:1529559/1‑68 (MQ=255)             tcAATCCGGCGTTACCTGAGCAGTGGCAACAGTTGTCTTTCCCTCTGTTCTGGCA‑GGGCTGCGAATTa                                                     >  1:1579705/1‑68 (MQ=255)             tcAATCCGGCGTTACCTGAGCAGTGGCAACAGTTGTCTTTCCCTCTGTTCTGGCA‑GGGCTGCGAATTACAg                                                  >  1:1204861/1‑71 (MQ=255)             tcAATCCGGCGTTACCTGAGCAGTGGCAACAGTTGTCTTTCCCTCTGTTCTGGCA‑GGGCTGCGAATTACAg                                                  >  1:2963869/1‑71 (MQ=255)             tcAATCCGGCGTTACCTGAGCAGTGGCAACAGTTGTCTTTCCCTCTGTTCTGGCA‑GGGCTGCGAATTACAg                                                  >  1:448815/1‑71 (MQ=255)              cAATCCGGCGTTACCTGAGCAGTGGCAACAGTTGTCTTTCCCTCTGTTCTGGCA‑GGGCTGCGAATTACAgg                                                 >  1:1818045/1‑71 (MQ=255)              cAATCCGGCGTTACCTGAGCAGTGGCAACAGTTGTCTTTCCCTCTGTTCTGGCA‑GGGCTGCGAATTACAgg                                                 >  1:3706460/1‑71 (MQ=255)                   cGGCGTTACCTGAGCAGTGGCAACAGTTGTCTTTCCCTCTGTTCTGGCATGGGCTGCGAATTACAGGTCAc                                             <  1:3252686/71‑1 (MQ=255)                   cGGCGTTACCTGAGCAGTGGCAACAGTTGTCTTTCCCTCTGTTCTGGCATGGGCTGCGAATTACAGGTCAc                                             <  1:1456927/71‑1 (MQ=255)                       gTTACCTGAGCAGTGGCAACAGTTGTCTTTCCATCTGTTCTGGCATGGGCTGCGAATTACAGGTCActct                                          >  1:3575545/1‑70 (MQ=255)                         tACCTGAGCAGTGGCAACAGTTGTCTTTCCCTCTGTTCTGGCA‑GGGCTGCGAATTACAGGTCACTCTTGAc                                      >  1:3877207/1‑71 (MQ=255)                          aCCTGAGCAGTGGCAACAGTTGTCTTTCCCTCTGTTCTGGCA‑GGGCTGCGaa                                                        >  1:2044110/1‑52 (MQ=255)                          aCCTGAGCAGTGGCAACAGTTGTCTTTCCCTCTGTTCTGGCA‑GGGCTGCGAATTACAGGTCACTCTTGAc                                      <  1:2551400/70‑1 (MQ=255)                              gAGCAGTGGCAACAGTTGTCTTTCCCTCTGTTCTGGCA‑GGGCTGCGaa                                                        >  1:491171/1‑48 (MQ=255)                              gAGCAGTGGCAACAGTTGTCTTTCACTCTGTTCTGGCA‑GGGCTGCGAATTACAGGTCACTCTTGAcgcg                                   >  1:2205571/1‑69 (MQ=255)                                gCAGTGGCAACAGTTGTCTTTCCCTCTGTTCTGGCA‑GGGCTGCGAATTACAGGTCACTCTTGACGCGCAGc                               >  1:1232723/1‑71 (MQ=255)                                 cAGTGGCAACAGTTGTCTTTCCCTCTGTTCTGGCA‑GGGCTGCGAATTACAGGTCACTCTTGACGCGCAGCg                              <  1:2901939/71‑1 (MQ=255)                                 cAGTGGCAACAGTTGTCTTTCCCTCTGTTCTGGCA‑GGGCTGCGAATTACAGGTCACTCTTGACGCGCAGCg                              <  1:2473123/71‑1 (MQ=255)                                   gTGGCAACAGTTGTCTTTCCCTCTGTTCTGGCA‑GGGCTGCGAATTACAGGTCACTCTTGACGCGCAGCGt                             >  1:1992626/1‑70 (MQ=255)                                    tGGCAACAGTTGTCTTTCCCTCTGTTCTGGCATGGGCTGCGAATTACAGGTCACTCTTGACGCGCAGCg                              <  1:1626985/69‑1 (MQ=255)                                    tGGCAACAGTTGTCTTTCCCTCTGTTCTGGCA‑GGGCTGCGAAt                                                       <  1:1699915/43‑1 (MQ=255)                                      gCAACAGTTGTCTTTCCCTCTGTTCTGGCA‑GGGCTGCGAATTACa                                                   >  1:2382014/1‑45 (MQ=255)                                      gCAACAGTTGTCTTTCCCTCTGTTCTGGCA‑GGGCTGCGAATTACAGGTCACTCTTGACGCGCAGCGTa                            >  1:702006/1‑68 (MQ=255)                                      gCAACAGTTGTCTTTCCCTCTGTTCTGGCA‑GGGCTGCGAATTACAGGTCACTCTTGACGCGCAGCGTATTg                         >  1:3872648/1‑71 (MQ=255)                                       cAACAGTTGTCTTTCCCTCTGTTCTGGCATGGGCTGCGAATTACAGGTCACTCTTGAcgcg                                   <  1:555689/61‑1 (MQ=255)                                        aaCAGTTGTCTTTCCCTCTGTTCTGGCATGGGCTGCGaa                                                        >  1:3405742/1‑39 (MQ=37)                                        aaCAGTTGTCTTTCCCTCTGTTCTGGCA‑GGGCTGCg                                                          <  1:187652/36‑1 (MQ=255)                                        aaCAGTTGTCTTTCCCTCTGTTCTGGCA‑GGGCTGCg                                                          <  1:884413/36‑1 (MQ=255)                                         aCAGTTGTCTTTCCCTCTGTTCTGGCATGGGCTGCGAATTACAGGTCACTCTTGa                                       >  1:1931525/1‑55 (MQ=255)                                          cAGTTGTCTTTCCCTCTGTTCTGGCATGGGCTGCGAATTACAGGTCACTCTTGACGCGCAGCGTATTGCGa                      <  1:897644/71‑1 (MQ=255)                                           aGTTGTCTTTCCCTCTGTTCTGGCATGGGCTGCGAATTACAGGTCACTCTTGACGCGCAGCGTATTGCGAt                     >  1:2871347/1‑71 (MQ=255)                                            gTTGTCTTTCCCTCTGTTCTGGCATGGGCTGCGAATTACAGGTCACTCTTGACGCGCAGCg                              <  1:3945176/61‑1 (MQ=255)                                            gTTGTCTTTCCCTCTGTTCTGGCA‑GGGCTGCGAATTa                                                     <  1:347308/37‑1 (MQ=255)                                            gTTGTCTTTCCCTCTGTTCTGGCA‑GGGCTGCGAATTa                                                     <  1:3553198/37‑1 (MQ=255)                                            gTTGTCTTTCCCTCTGTTCTGGCA‑GGGCTGCGAATTACAGGTCACTCTTGACGCGCAGCGTATTGCGATTc                   <  1:3701076/71‑1 (MQ=255)                                              tGTCTTTCCCTCTGTTCTGGCATGGGCTGCGAATTACAGGTCACTCTTGACGCGCAGCGTATTGCGAttt                   >  1:1329403/1‑69 (MQ=255)                                              tGTCTTTCCCTCTGTTCTGGCA‑GGGCTGCGAATTACAGGTCACTCTTGACGCGCAGCGTATTGc                        >  1:3302091/1‑64 (MQ=255)                                               gTCTTTCCCTCTGTTCTGGCATGGGCTGCGAATTACAGGTCACTCTTGACGCGCAGCGTATTGCGATTCg                  >  1:1211303/1‑70 (MQ=255)                                               gTCTTTCCCTCTGTTCTGGCA‑GGGCTGCGAATTACATGTCACTCTTGa                                       >  1:3713403/1‑48 (MQ=255)                                                 cTTTCCCTCTGTTCTGGCA‑GGGCTGCGAATTACAGGTCACTCTTGACGCGCAGCGTATTGCGATTCGAACt              >  1:1710063/1‑71 (MQ=255)                                                  tttCCCTCTGTTCTGGCATGGGCTGCGAATTACAGGTCACTCTTGACGCGCAGCGTATTGCGATTCGAACt              >  1:3066631/1‑71 (MQ=255)                                                    tCCCTCTGTTCTGGCATGGGCTGCGAATTACAGGTCACTCTTGACGCGCAGCGTATTGCGATTCGAACt              <  1:2121110/69‑1 (MQ=255)                                                     cccTCTGTTCTGGCATGGGCTGCGAATTACAGGTCACTCTTGACGCGCAGCGTATTGCGATTCGAACTTc            <  1:847024/70‑1 (MQ=255)                                                     cccTCTGTTCTGGCA‑GGGCTGCGAATTACAGGTCACTCTTGACGCGCAGct                              <  1:1373954/51‑2 (MQ=255)                                                     cccTCTGTTCTGGCA‑GGGCTGCGAATTACAGGTCACTCTTGACGCGCAGCg                              <  1:3511170/51‑1 (MQ=255)                                                     cccTCTGTTCTGGCA‑GGGCTGCGAATTACAGGTCACTCTTGACGCGCAGCg                              <  1:3981919/51‑1 (MQ=255)                                                      ccTCTGTTCTGGCA‑GGGCTGCGAATTACAGGTCACTCTTGACGCGCAGCGTATTGCGATTCGAACtt             <  1:305919/67‑1 (MQ=255)                                                        tctGTTCTGGCATGGGCTGCGAATTACAGGTCACTCtt                                         <  1:2196674/38‑1 (MQ=38)                                                        tctGTTCTGGCATGGGCTGCGAATTACAGGTCACTCtt                                         <  1:2369888/38‑1 (MQ=38)                                                        tctGTTCTGGCATGGGCTGCGAATTACAGGTCACTCTTGAcgc                                    <  1:1661736/43‑1 (MQ=255)                                                        tctGTTCTGGCATGGGCTGCGAATTACAGGTCACTCTTGACGCGCAGCg                              <  1:2683475/49‑1 (MQ=255)                                                         ctGTTCTGGCATGGGCTGCGAATTACAGGTCACTCTTGACGCGCAGCGTATTGCGATTCGAACtt             >  1:4072661/1‑65 (MQ=255)                                                          tGTTCTGGCATGGGCTGCGAATTACAGGTCACTCTTTACGCGCAGCg                              <  1:543084/47‑1 (MQ=255)                                                           gTTCTGGCA‑GGGCTGCGAATTACAGGTCACTCTTGACGCGCAGCGTATTGCGATTCGAACTTc            <  1:2253085/63‑1 (MQ=255)                                                           gTTCTGGCA‑GGGCTGCGAATTACAGGTCACTCTTGACGCGCAGCGTATTGCGATTCGAACTTc            <  1:3020667/63‑1 (MQ=255)                                                           gTTCTGGCA‑GGGCTGCGAATTACAGGTCACTCTTGACGCGCAGCGTATTGCGATTCGAACTTCTGCGCCCg    >  1:1904032/1‑71 (MQ=255)                                                             tCTGGCA‑GGGCTGCGAATTACAGGTCACTCTTGACGCGCAGCGTATTGCGATTCGAACTTCTGCGCCCGtt  >  1:1160401/1‑71 (MQ=255)                                                                    |                                                                 GTGAATTGCATCTCAATCCGGCGTTACCTGAGCAGTGGCAACAGTTGTCTTTCCCTCTGTTCTGGCA‑GGGCTGCGAATTACAGGTCACTCTTGACGCGCAGCGTATTGCGATTCGAACTTCTGCGCCCGTT  >  minE/909847‑909977

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