Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	1,266,337	T→A	37.6%	C158S (TGC→AGC)	tyrP →	tyrosine transporter

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	1,266,337	0	T	A	37.6%	62.7 / 65.8	77	C158S (TGC→AGC)	tyrP	tyrosine transporter
Reads supporting (aligned to +/- strand):  ref base T (12/36);  new base A (10/19);  total (22/55)
Fisher's exact test for biased strand distribution p-value = 4.39e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.93e-01

CGCATGTTGAGCTGGCAGTGCATTTATTTGCTGATTTAGCCCTCGCCACGTCATTTCTCGGCGTTGCGTTAGGCTTATTTGATTATCTGGCTGATTTATTTCAGCGTTCAAATACCGTTGGTGGACGG  >  minE/1266273‑1266400                                                                 |                                                                cGCATGTTGAGCTGGCAGTGCATTTATTTGCTGATTTAGCCCTCGCCACGTCATTTCTCGgcgttgc                                                               <  1:2857440/67‑1 (MQ=255) cGCATGTTGAGCTGGCAGTGCATTTATTTGCTGATTTAGCCCTCGCCACGTCATTTCTCGgcgttgc                                                               <  1:915419/67‑1 (MQ=255) cGCATGTTGAGCTGGCAGTGCATTTATTTGCTGATTTAGCCCTCGCCACGTCATTTCTCGgcgttgc                                                               <  1:2310668/67‑1 (MQ=255) cGCATGTTGAGCTGGCAGTGCATTTATTTGCTGATTTAGCCCTCGCCACGTCATTTCTCGgcgttgc                                                               <  1:2029651/67‑1 (MQ=255)    aTGTTGAGCTGGCAGTGCATTTATTTGCTGATTTAGCCCTCGCCACGTCATTTCTCGGCGTTGCGTTAg                                                          <  1:2511730/69‑1 (MQ=255)    aTGTTGAGCTGGCAGTGCATTTATTTGCTGATTTAGCCCTCGCCACGTCATTTCTCGGCGTTGCGTTAg                                                          <  1:3424257/69‑1 (MQ=255)    aTGTTGAGCTGGCAGTGCATTTATTTGCTGATTTAGCCCTCGCCACGTCATTTCTCGGCGTTGCGTTAg                                                          <  1:1295268/69‑1 (MQ=255)                 aGTGCATTTATTTGCTGATTTAGCCCTCGCCACGTCATTTCTCGGCGTTGCGTTAGGCTTATTTGAtt                                              <  1:2500210/68‑1 (MQ=255)                 aGTGCATTTATTTGCTGATTTAGCCCTCGCCACGTCATTTCTCGGCGTTGCGTTAGGCTTATTTGAtt                                              <  1:1846031/68‑1 (MQ=255)                 aGTGCATTTATTTGCTGATTTAGCCCTCGCCACGTCATTTCTCGGCGTTGCGTTAGGCTTATTTGAtt                                              <  1:1931726/68‑1 (MQ=255)                 aGTGCATTTATTTGCTGATTTAGCCCTCGCCACGTCATTTCTCGGCGTTGCGTTAGGCTTATTTGAtt                                              <  1:2255958/68‑1 (MQ=255)                 aGTGCATTTATTTGCTGATTTAGCCCTCGCCACGTCATTTCTCGGCGTTGCGTTAGGCTTATTTGAtt                                              <  1:28347/68‑1 (MQ=255)                 aGTGCATTTATTTGCTGATTTAGCCCTCGCCACGTCATTTCTCGGCGTTGCGTTAGGCTTATTTGAt                                               <  1:2719786/67‑1 (MQ=255)                 aGTGCATTTATTTGCTGATTTAGCCCTCGCCACGTCATTTCTCGGCGTTGCGTTAGGCTTATTTGAt                                               <  1:1678200/67‑1 (MQ=255)                 aGTGCATTTATTTGCTGATTTAGCCCTCGCCACGTCATTTCTCGGCGTTGCGTTAGGCTTATTTGAt                                               <  1:141483/67‑1 (MQ=255)                     cATTTATTTGCTGATTTAGCCCTCGCCACGTCCTTTCTCGGCGTTGCGTTAGGCTTATTTg                                                 <  1:660801/61‑1 (MQ=255)                     cATTTATTTGCTGATTTAGCCCTCGCCACGTCATTTCTCGGCGTTGCGTTAGGCTTATTTg                                                 <  1:51963/61‑1 (MQ=255)                     cATTTATTTGCTGATTTAGCCCTCGCCACGTCATTTCTCGGCGTTGCGTTAGGCTTATTTg                                                 <  1:2533200/61‑1 (MQ=255)                     cATTTATTTGCTGATTTAGCCCTCGCCACGTCATTTCTCGGCGTTGCGTTAGGCTTATTTg                                                 <  1:1849410/61‑1 (MQ=255)                     cATTTATTTGCTGATTTAGCCCTCGCCACGTCATTTCTCGGCGTTGCGTTAGGCTTATTTg                                                 <  1:921447/61‑1 (MQ=255)                     cATTTATTTGCTGATTTAGCCCTCGCCACGTCATTTCTCGGCGTTGCGTTAGGCTTATTTGATTTTCt                                          >  1:2608118/1‑68 (MQ=255)                     cATTTATTTGCTGATTTAGCCCTCGCCACGTCATTTCTCGGCGTTGCGTTAGGCTTATTTGATTATCTg                                         >  1:2493470/1‑69 (MQ=255)                     cATTTATTTGCTGATTTAGCCCTCGCCACGTCATTTCTCGGCGTTGCGTTAGGCTTATTTGATTATCTg                                         >  1:629184/1‑69 (MQ=255)                     cATTTATTTGCTGATTTAGCCCTCGCCACGTCATTTCTCGGCGTTGCGTTAGGCTTATTTGATTATCTg                                         >  1:2384551/1‑69 (MQ=255)                     cATTTATTTGCTGATTTAGCCCTCGCCACGTCATTTCTCGGCGTTGCGTTAGGCTTATTTGATTATCTg                                         >  1:2881429/1‑69 (MQ=255)                     cATTTATTTGCTGATTTAGCCCTCGCCACGTCATTTCTCGGCGTTGCGTTAGGCTTATTTGATTATCTg                                         >  1:1547878/1‑69 (MQ=255)                      atttatttGCTGATTTAGCCCTCGCCACGTCATTTCTCGGCGTTGCGTTAGGCTTATTTg                                                 <  1:380717/60‑1 (MQ=255)                         tatttGCTGATTTAGCCCTCGCCACGTCATTTCTCGGCGTTGCGTTAGGCTTATTTGATTATCTGGCTg                                     >  1:3294520/1‑69 (MQ=255)                                tGATTTAGCCCTCGCCACGTCCTTCCTCGGCGTAGCGTTAGGCTTATTTGATTATCTGGCAGatttgt                               >  1:2662378/1‑66 (MQ=255)                                tGATTTAGCCCTCGCCACGTCCTTCCTCGGCGTAGCGTTAGGCTTATTTGATTATCTGGCAGatttgt                               >  1:2411078/1‑66 (MQ=255)                                tGATTTAGCCCTCGCCACGTCCTTCCTCGGCGTAGCGTTAGGCTTATTTGATTATCTGGCAGatttgt                               >  1:3218342/1‑66 (MQ=255)                                tGATTTAGCCCTCGCCACGTCCTTCCTCGGCGTAGCGTTAGGCTTATTTGATTATCTGGCAGatttgt                               >  1:1516468/1‑66 (MQ=255)                                tGATTTAGCCCTCGCCACGTCCTTCCTCGGCGTAGCGTTAGGCTTATTTGATTATCTGGCAGatttgt                               >  1:2193034/1‑66 (MQ=255)                                tGATTTAGCCCTCGCCACGTCCTTCCTCGGCGTAGCGTTAGGCTTATTTGATTATCTGGCAGatttgt                               >  1:43571/1‑66 (MQ=255)                                tGATTTAGCCCTCGCCACGTCCTTCCTCGGCGTAGCGTTAGGCTTATTTGATTATCTGGCAGatttgt                               >  1:2160319/1‑66 (MQ=255)                                tGATTTAGCCCTCGCCACGTCCTTCCTCGGCGTAGCGTTAGGCTTATTTGATTATCTGGCAGatttgt                               >  1:3274818/1‑66 (MQ=255)                                tGATTTAGCCCTCGCCACGTCCTTCCTCGGCGTAGCGTTAGGCTTATTTGATTATCTGGCAGatttgt                               >  1:2630148/1‑66 (MQ=255)                                tGATTTAGCCCTCGCCACGTCCTTCCTCGGCGTAGCGTTAGGCTTATTTGATTATCTGGCAGatttg                                >  1:2570144/1‑66 (MQ=255)                                         ccTCGCCACGTCATTTCTCGGCGTTGCGTTAGGCTTATTTGATTATCTGGCTTATTTATTTCAGCGTTc                     <  1:3399623/69‑1 (MQ=255)                                         ccTCGCCACGTCATTTCTCGGCGTTGCGTTAGGCTTATTTGATTATCTGGCTGATTTATTTCAGCGTTc                     <  1:2699817/69‑1 (MQ=255)                                         ccTCGCCACGTCATTTCTCGGCGTTGCGTTAGGCTTATTTGATTATCTGGCTGATTTATTTCAGCGTTc                     <  1:3239450/69‑1 (MQ=255)                                         ccTCGCCACGTCATTTCTCGGCGTTGCGTTAGGCTTATTTGATTATCTGGCTGATTTATTTCAGCGTTc                     <  1:1080130/69‑1 (MQ=255)                                         ccTCGCCACGTCATTTCTCGGCGTTGCGTTAGGCTTATTTGATTATCTGGCTGATTTATTTCAGCGTTc                     <  1:3144134/69‑1 (MQ=255)                                         ccTCGCCACGTCATTTCTCGGCGTTGCGTTAGGCTTATTTGATTATCTGGCTGATTTATTTCAGCGTTc                     <  1:623533/69‑1 (MQ=255)                                         ccTCGCCACGTCATTTCTCGGCGTTGCGTTAGGCTTATTTGATTATCTGGCTGATTTATTTCAGCGTTc                     <  1:188318/69‑1 (MQ=255)                                         ccTCGCCACGTCATTTCTCGGCGTTGCGTTAGGCTTATTTGATTATCTGGCTGATTTATTTCAGCGTTc                     <  1:1150597/69‑1 (MQ=255)                                          cTCGCCACGTCATTTCTCGGCGTTGCGTTAGGCTTATTTGATTATCTGGCTGATTTATTTCAGCGTTCa                    <  1:3136750/69‑1 (MQ=255)                                          cTCGCCACGTCATTTCTCGGCGTTGCGTTAGGCTTATTTGATTATCTGGCTGATTTATTTCAGCGTTCa                    <  1:2541142/69‑1 (MQ=255)                                          cTCGCCACGTCATTTCTCGGCGTTGCGTTAGGCTTATTTGATTATCTGGCTGATTTATTTCAGCGTTCa                    <  1:1896990/69‑1 (MQ=255)                                          cTCGCCACGTCATTTCTCGGCGTTGCGTTAGGCTTATTTGATTATCTGGCTGATTTATTTCAGCGTTCa                    <  1:3214387/69‑1 (MQ=255)                                          cTCGCCACGTCATTTCTCGGCGTTGCGTTAGGCTTATTTGATTATCTGGCTGATTTATTTCAGCGTTCa                    <  1:820825/69‑1 (MQ=255)                                          cTCGCCACGTCATTTCTCGGCGTTGCGTTAGGCTTATTTGATTATCTGGCTGATTTATTTCAGCGTTCa                    <  1:1601620/69‑1 (MQ=255)                                          cTCGCCACGTCATTTCTCGGCGTTGCGTTAGGCTTATTTGATTATCTGGCTGATTTATTTCAGCGTTCa                    <  1:2790353/69‑1 (MQ=255)                                          cTCGCCACGTCATTTCTCGGCGTTGCGTTAGGCTTATTTGATTATCTGGCTGATTTATTTCAGCGTTCa                    <  1:2791876/69‑1 (MQ=255)                                          cTCGCCACGTCATTTCTCGGCGTTGCGTTAGGCTTATTTGATTATCTGGCTGATTTATTTCAGCGTTCa                    <  1:856281/69‑1 (MQ=255)                                           tCGCCACGTCATTTCTCGGCGTTGCGTTAGGCTTATTTGATTATCTGGCTGATTTATTTCAGCGTTCa                    <  1:1220904/68‑1 (MQ=255)                                             gCCACGTCATTTCTCGGCGTTGCGTTAGGCTTATTTGAtt                                              >  1:519850/1‑40 (MQ=255)                                             gCCACGTCATTTCTCGGCGTTGCGTTAGGCTTATTTGAtt                                              >  1:2955289/1‑40 (MQ=255)                                             gCCACGTCATTTCTCGGCGTTGCGTTAGGCTTATTTGAtt                                              >  1:1428495/1‑40 (MQ=255)                                             gCCACGTCATTTCTCGGCGTTGCGTTAGGCTTATTTGAtt                                              >  1:811621/1‑40 (MQ=255)                                             gCCACGTCATTTCTCGGCGTTGCGTTAGGCTTATTTGAtt                                              >  1:672591/1‑40 (MQ=255)                                                aCGTCCTTCCTCGTCGTAGCGTTAGGCTTATTTGATTATCTGGCAGATTTGTTTCAGCGTTCaaa                  <  1:1574527/65‑1 (MQ=25)                                                aCGTCCTTCCTCGGCGTAGCGTTAGGCTTATTTGATTATCTGGCAGATTTGTTTCAGCGTTCaaa                  <  1:3230574/65‑1 (MQ=25)                                                aCGTCCTTCCTCGGCGTAGCGTTAGGCTTATTTGATTATCTGGCAGATTTGTTTCAGCGTTCaaa                  <  1:713655/65‑1 (MQ=25)                                                aCGTCCTTCCTCGGCGTAGCGTTAGGCTTATTTGATTATCTGGCAGATTTGTTTCAGCGTTCaaa                  <  1:944699/65‑1 (MQ=25)                                                aCGTCCTTCCTCGGCGTAGCGTTAGGCTTATTTGATTATCTGGCAGATTTGTTTCAGCGTTCaaa                  <  1:3334992/65‑1 (MQ=25)                                                aCGTCCTTCCTCGGCGTAGCGTTAGGCTTATTTGATTATCTGGCAGATTTGTTTCAGCGTTCaaa                  <  1:2485707/65‑1 (MQ=25)                                                aCGTCCTTCCTCGGCGTAGCGTTAGGCTTATTTGATTATCTGGCAGATTTGTTTCAGCGTTCaaa                  <  1:1945261/65‑1 (MQ=25)                                                aCGTCCTTCCTCGGCGTAGCGTTAGGCTTATTTGATTATCTGGCAGATTTGTTTCAGCGTTCaaa                  <  1:191684/65‑1 (MQ=25)                                                aCGTCCTTCCTCGGCGTAGCGTTAGGCTTATTTGATTATCTGGCAGATTTGTTTCAGCGTTCaaa                  <  1:1828988/65‑1 (MQ=25)                                                aCGTCCTTCCTCGGCGTAGCGTTAGGCTTATTTGATTATCTGGCAGATTTGTTTCAGCGTTCaaa                  <  1:2257788/65‑1 (MQ=25)                                                aCGTCCTTCCTCGGCGTAGCGTTAGGCTTATTTGATTATCTGGCAGATTTGTTTCAGCGTTCaaa                  <  1:2312632/65‑1 (MQ=25)                                                aCGTCCTTCCTCGGCGTAGCGTTAGGCTTATTTGATTATCTGGCAGATTTGTTTCAGCGTTCaaa                  <  1:2320293/65‑1 (MQ=25)                                                aCGTCCTTCCTCGGCGTAGCGTTAGGCTTATTTGATTATCTGGCAGATTTGTTTCAGCGTTCaaa                  <  1:2526049/65‑1 (MQ=25)                                                aCGTCCTTCCTCGGCGTAGCGTTAGGCTTATTTGATTATCTGGCAGATTTGTTTCAGCGTTCaaa                  <  1:1354473/65‑1 (MQ=25)                                                aCGTCCTTCCTCGGCGTAGCGTTAGGCTTATTTGATTATCTGGCAGATTTGTTTCAGCGTTCaaa                  <  1:2693473/65‑1 (MQ=25)                                                aCGTCCTTCCTCGGCGTAGCGTTAGGCTTATTTGATTATCTGGCAGATTTGTTTCAGCGTTCaaa                  <  1:2762721/65‑1 (MQ=25)                                                aCGTCCTTCCTCGGCGTAGCGTTAGGCTTATTTGATTATCTGGCAGATTTGTTTCAGCGTTCaaa                  <  1:3082391/65‑1 (MQ=25)                                                aCGTCCTTCCTCGGCGTAGCGTTAGGCTTATTTGATTATCTGGCAGATTTGTTTCAGCGTTCaaa                  <  1:3021034/65‑1 (MQ=25)                                                aCGTCCTTCCTCGGCGTAGCGTTAGGCTTATTTGATTATCTGGCAGATTTGTTTCAGCCTTCaaa                  <  1:2674440/65‑1 (MQ=21)                                                         ctcGGCGTTGCGTTAGGCTTATTTGTTTATCTGGCTGATTTATTTCTGCGTTc                     <  1:3289956/53‑1 (MQ=255)                                                              cgtagcgttAGGCTTATTTGATTATCTGGCAGATTTGTTTCAGCGTTCAAATACCGTTGGTGGACga  <  1:1513930/67‑2 (MQ=38)                                                              cgtagcgttAGGCTTATTTGATTATCTGGCAGATTTGTTTCAGCGTTCAAATACCGTTGGTGGACga  <  1:1087926/67‑2 (MQ=38)                                                              cgtagcgttAGGCTTATTTGATTATCTGGCAGATTTGTTTCAGCGTTCAAATACCGTTGGTGGACga  <  1:1613412/67‑2 (MQ=38)                                                              cgtagcgttAGGCTTATTTGATTATCTGGCAGATTTGTTTCAGCGTTCAAATACCGTTGGTGGACga  <  1:1717096/67‑2 (MQ=38)                                                              cgtagcgttAGGCTTATTTGATTATCTGGCAGATTTGTTTCAGCGTTCAAATACCGTTGGTGGACga  <  1:1138731/67‑2 (MQ=38)                                                              cgtagcgttAGGCTTATTTGATTATCTGGCAGATTTGTTTCAGCGTTCAAATACCGTTGGTGGACga  <  1:1774797/67‑2 (MQ=38)                                                              cgtagcgttAGGCTTATTTGATTATCTGGCAGATTTGTTTCAGCGTTCAAATACCGTTGGTGGACga  <  1:1807307/67‑2 (MQ=38)                                                              cgtagcgttAGGCTTATTTGATTATCTGGCAGATTTGTTTCAGCGTTCAAATACCGTTGGTGGACga  <  1:640652/67‑2 (MQ=38)                                                              cgtagcgttAGGCTTATTTGATTATCTGGCAGATTTGTTTCAGCGTTCAAATACCGTTGGTGGACga  <  1:619102/67‑2 (MQ=38)                                                              cgtagcgttAGGCTTATTTGATTATCTGGCAGATTTGTTTCAGCGTTCAAATACCGTTGGTGGACga  <  1:613463/67‑2 (MQ=38)                                                              cgtagcgttAGGCTTATTTGATTATCTGGCAGATTTGTTTCAGCGTTCAAATACCGTTGGTGGACga  <  1:2031850/67‑2 (MQ=38)                                                              cgtagcgttAGGCTTATTTGATTATCTGGCAGATTTGTTTCAGCGTTCAAATACCGTTGGTGGACga  <  1:208607/67‑2 (MQ=38)                                                              cgtagcgttAGGCTTATTTGATTATCTGGCAGATTTGTTTCAGCGTTCAAATACCGTTGGTGGACga  <  1:2167024/67‑2 (MQ=38)                                                              cgtagcgttAGGCTTATTTGATTATCTGGCAGATTTGTTTCAGCGTTCAAATACCGTTGGTGGACga  <  1:3428499/67‑2 (MQ=38)                                                              cgtagcgttAGGCTTATTTGATTATCTGGCAGATTTGTTTCAGCGTTCAAATACCGTTGGTGGACga  <  1:3383131/67‑2 (MQ=38)                                                              cgtagcgttAGGCTTATTTGATTATCTGGCAGATTTGTTTCAGCGTTCAAATACCGTTGGTGGACga  <  1:3302032/67‑2 (MQ=38)                                                              cgtagcgttAGGCTTATTTGATTATCTGGCAGATTTGTTTCAGCGTTCAAATACCGTTGGTGGACga  <  1:3253801/67‑2 (MQ=38)                                                              cgtagcgttAGGCTTATTTGATTATCTGGCAGATTTGTTTCAGCGTTCAAATACCGTTGGTGGACga  <  1:3238596/67‑2 (MQ=38)                                                              cgtagcgttAGGCTTATTTGATTATCTGGCAGATTTGTTTCAGCGTTCAAATACCGTTGGTGGACga  <  1:2856762/67‑2 (MQ=38)                                                                 |                                                                CGCATGTTGAGCTGGCAGTGCATTTATTTGCTGATTTAGCCCTCGCCACGTCATTTCTCGGCGTTGCGTTAGGCTTATTTGATTATCTGGCTGATTTATTTCAGCGTTCAAATACCGTTGGTGGACGG  >  minE/1266273‑1266400

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