Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	1,266,328	T→C	36.2%	S155P (TCT→CCT)	tyrP →	tyrosine transporter

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	1,266,328	0	T	C	36.2%	49.6 / 48.0	80	S155P (TCT→CCT)	tyrP	tyrosine transporter
Reads supporting (aligned to +/- strand):  ref base T (12/39);  new base C (10/19);  total (22/58)
Fisher's exact test for biased strand distribution p-value = 3.10e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.99e-01

TGGTGGCCTCTCCGCATGTTGAGCTGGCAGTGCATTTATTTGCTGATTTAGCCCTCGCCACGTCATTTCTCGGCGTTGCGTTAGGCTTATTTGATTATCTGGCTGATTTATTTCAGCGTTCAAA  >  minE/1266261‑1266384                                                                    |                                                         tggtggCCTCTCCGCATGTTGAGCTGGCAGTGCATTTATTTGCTGATTTAGCCCTCGCCACGTCATTtc                                                         >  1:2345942/1‑69 (MQ=255)             cGCATGTTGAGCTGGCAGTGCATTTATTTGCTGATTTAGCCCTCGCCACGTCATTTCTCGgcgttgc                                               <  1:2029651/67‑1 (MQ=255)             cGCATGTTGAGCTGGCAGTGCATTTATTTGCTGATTTAGCCCTCGCCACGTCATTTCTCGgcgttgc                                               <  1:915419/67‑1 (MQ=255)             cGCATGTTGAGCTGGCAGTGCATTTATTTGCTGATTTAGCCCTCGCCACGTCATTTCTCGgcgttgc                                               <  1:2310668/67‑1 (MQ=255)             cGCATGTTGAGCTGGCAGTGCATTTATTTGCTGATTTAGCCCTCGCCACGTCATTTCTCGgcgttgc                                               <  1:2857440/67‑1 (MQ=255)                aTGTTGAGCTGGCAGTGCATTTATTTGCTGATTTAGCCCTCGCCACGTCATTTCTCGGCGTTGCGTTAg                                          <  1:1295268/69‑1 (MQ=255)                aTGTTGAGCTGGCAGTGCATTTATTTGCTGATTTAGCCCTCGCCACGTCATTTCTCGGCGTTGCGTTAg                                          <  1:3424257/69‑1 (MQ=255)                aTGTTGAGCTGGCAGTGCATTTATTTGCTGATTTAGCCCTCGCCACGTCATTTCTCGGCGTTGCGTTAg                                          <  1:2511730/69‑1 (MQ=255)                             aGTGCATTTATTTGCTGATTTAGCCCTCGCCACGTCATTTCTCGGCGTTGCGTTAGGCTTATTTGAtt                              <  1:1931726/68‑1 (MQ=255)                             aGTGCATTTATTTGCTGATTTAGCCCTCGCCACGTCATTTCTCGGCGTTGCGTTAGGCTTATTTGAtt                              <  1:2255958/68‑1 (MQ=255)                             aGTGCATTTATTTGCTGATTTAGCCCTCGCCACGTCATTTCTCGGCGTTGCGTTAGGCTTATTTGAtt                              <  1:1846031/68‑1 (MQ=255)                             aGTGCATTTATTTGCTGATTTAGCCCTCGCCACGTCATTTCTCGGCGTTGCGTTAGGCTTATTTGAtt                              <  1:2500210/68‑1 (MQ=255)                             aGTGCATTTATTTGCTGATTTAGCCCTCGCCACGTCATTTCTCGGCGTTGCGTTAGGCTTATTTGAtt                              <  1:28347/68‑1 (MQ=255)                             aGTGCATTTATTTGCTGATTTAGCCCTCGCCACGTCATTTCTCGGCGTTGCGTTAGGCTTATTTGAt                               <  1:1678200/67‑1 (MQ=255)                             aGTGCATTTATTTGCTGATTTAGCCCTCGCCACGTCATTTCTCGGCGTTGCGTTAGGCTTATTTGAt                               <  1:141483/67‑1 (MQ=255)                             aGTGCATTTATTTGCTGATTTAGCCCTCGCCACGTCATTTCTCGGCGTTGCGTTAGGCTTATTTGAt                               <  1:2719786/67‑1 (MQ=255)                                 cATTTATTTGCTGATTTAGCCCTCGCCACGTCCTTTCTCGGCGTTGCGTTAGGCTTATTTg                                 <  1:660801/61‑1 (MQ=255)                                 cATTTATTTGCTGATTTAGCCCTCGCCACGTCATTTCTCGGCGTTGCGTTAGGCTTATTTg                                 <  1:2533200/61‑1 (MQ=255)                                 cATTTATTTGCTGATTTAGCCCTCGCCACGTCATTTCTCGGCGTTGCGTTAGGCTTATTTg                                 <  1:921447/61‑1 (MQ=255)                                 cATTTATTTGCTGATTTAGCCCTCGCCACGTCATTTCTCGGCGTTGCGTTAGGCTTATTTg                                 <  1:1849410/61‑1 (MQ=255)                                 cATTTATTTGCTGATTTAGCCCTCGCCACGTCATTTCTCGGCGTTGCGTTAGGCTTATTTg                                 <  1:51963/61‑1 (MQ=255)                                 cATTTATTTGCTGATTTAGCCCTCGCCACGTCATTTCTCGGCGTTGCGTTAGGCTTATTTGATTTTCt                          >  1:2608118/1‑68 (MQ=255)                                 cATTTATTTGCTGATTTAGCCCTCGCCACGTCATTTCTCGGCGTTGCGTTAGGCTTATTTGATTATCTg                         >  1:1547878/1‑69 (MQ=255)                                 cATTTATTTGCTGATTTAGCCCTCGCCACGTCATTTCTCGGCGTTGCGTTAGGCTTATTTGATTATCTg                         >  1:2881429/1‑69 (MQ=255)                                 cATTTATTTGCTGATTTAGCCCTCGCCACGTCATTTCTCGGCGTTGCGTTAGGCTTATTTGATTATCTg                         >  1:2384551/1‑69 (MQ=255)                                 cATTTATTTGCTGATTTAGCCCTCGCCACGTCATTTCTCGGCGTTGCGTTAGGCTTATTTGATTATCTg                         >  1:629184/1‑69 (MQ=255)                                 cATTTATTTGCTGATTTAGCCCTCGCCACGTCATTTCTCGGCGTTGCGTTAGGCTTATTTGATTATCTg                         >  1:2493470/1‑69 (MQ=255)                                  atttatttGCTGATTTAGCCCTCGCCACGTCATTTCTCGGCGTTGCGTTAGGCTTATTTg                                 <  1:380717/60‑1 (MQ=255)                                     tatttGCTGATTTAGCCCTCGCCACGTCATTTCTCGGCGTTGCGTTAGGCTTATTTGATTATCTGGCTg                     >  1:3294520/1‑69 (MQ=255)                                            tGATTTAGCCCTCGCCACGTCCTTCCTCGGCGTAGCGTTAGGCTTATTTGATTATCTGGCAGatttgt               >  1:2411078/1‑66 (MQ=255)                                            tGATTTAGCCCTCGCCACGTCCTTCCTCGGCGTAGCGTTAGGCTTATTTGATTATCTGGCAGatttgt               >  1:2630148/1‑66 (MQ=255)                                            tGATTTAGCCCTCGCCACGTCCTTCCTCGGCGTAGCGTTAGGCTTATTTGATTATCTGGCAGatttgt               >  1:1516468/1‑66 (MQ=255)                                            tGATTTAGCCCTCGCCACGTCCTTCCTCGGCGTAGCGTTAGGCTTATTTGATTATCTGGCAGatttgt               >  1:2662378/1‑66 (MQ=255)                                            tGATTTAGCCCTCGCCACGTCCTTCCTCGGCGTAGCGTTAGGCTTATTTGATTATCTGGCAGatttgt               >  1:43571/1‑66 (MQ=255)                                            tGATTTAGCCCTCGCCACGTCCTTCCTCGGCGTAGCGTTAGGCTTATTTGATTATCTGGCAGatttgt               >  1:3218342/1‑66 (MQ=255)                                            tGATTTAGCCCTCGCCACGTCCTTCCTCGGCGTAGCGTTAGGCTTATTTGATTATCTGGCAGatttgt               >  1:3274818/1‑66 (MQ=255)                                            tGATTTAGCCCTCGCCACGTCCTTCCTCGGCGTAGCGTTAGGCTTATTTGATTATCTGGCAGatttgt               >  1:2193034/1‑66 (MQ=255)                                            tGATTTAGCCCTCGCCACGTCCTTCCTCGGCGTAGCGTTAGGCTTATTTGATTATCTGGCAGatttgt               >  1:2160319/1‑66 (MQ=255)                                            tGATTTAGCCCTCGCCACGTCCTTCCTCGGCGTAGCGTTAGGCTTATTTGATTATCTGGCAGatttg                >  1:2570144/1‑66 (MQ=255)                                                     ccTCGCCACGTCATTTCTCGGCGTTGCGTTAGGCTTATTTGATTATCTGGCTTATTTATTTCAGCGTTc     <  1:3399623/69‑1 (MQ=255)                                                     ccTCGCCACGTCATTTCTCGGCGTTGCGTTAGGCTTATTTGATTATCTGGCTGATTTATTTCAGCGTTc     <  1:1150597/69‑1 (MQ=255)                                                     ccTCGCCACGTCATTTCTCGGCGTTGCGTTAGGCTTATTTGATTATCTGGCTGATTTATTTCAGCGTTc     <  1:623533/69‑1 (MQ=255)                                                     ccTCGCCACGTCATTTCTCGGCGTTGCGTTAGGCTTATTTGATTATCTGGCTGATTTATTTCAGCGTTc     <  1:3144134/69‑1 (MQ=255)                                                     ccTCGCCACGTCATTTCTCGGCGTTGCGTTAGGCTTATTTGATTATCTGGCTGATTTATTTCAGCGTTc     <  1:3239450/69‑1 (MQ=255)                                                     ccTCGCCACGTCATTTCTCGGCGTTGCGTTAGGCTTATTTGATTATCTGGCTGATTTATTTCAGCGTTc     <  1:188318/69‑1 (MQ=255)                                                     ccTCGCCACGTCATTTCTCGGCGTTGCGTTAGGCTTATTTGATTATCTGGCTGATTTATTTCAGCGTTc     <  1:1080130/69‑1 (MQ=255)                                                     ccTCGCCACGTCATTTCTCGGCGTTGCGTTAGGCTTATTTGATTATCTGGCTGATTTATTTCAGCGTTc     <  1:2699817/69‑1 (MQ=255)                                                      cTCGCCACGTCATTTCTCGGCGTTGCGTTAGGCTTATTTGATTATCTGGCTGATTTATTTCAGCGTTCa    <  1:3136750/69‑1 (MQ=255)                                                      cTCGCCACGTCATTTCTCGGCGTTGCGTTAGGCTTATTTGATTATCTGGCTGATTTATTTCAGCGTTCa    <  1:1896990/69‑1 (MQ=255)                                                      cTCGCCACGTCATTTCTCGGCGTTGCGTTAGGCTTATTTGATTATCTGGCTGATTTATTTCAGCGTTCa    <  1:1601620/69‑1 (MQ=255)                                                      cTCGCCACGTCATTTCTCGGCGTTGCGTTAGGCTTATTTGATTATCTGGCTGATTTATTTCAGCGTTCa    <  1:3214387/69‑1 (MQ=255)                                                      cTCGCCACGTCATTTCTCGGCGTTGCGTTAGGCTTATTTGATTATCTGGCTGATTTATTTCAGCGTTCa    <  1:820825/69‑1 (MQ=255)                                                      cTCGCCACGTCATTTCTCGGCGTTGCGTTAGGCTTATTTGATTATCTGGCTGATTTATTTCAGCGTTCa    <  1:856281/69‑1 (MQ=255)                                                      cTCGCCACGTCATTTCTCGGCGTTGCGTTAGGCTTATTTGATTATCTGGCTGATTTATTTCAGCGTTCa    <  1:2541142/69‑1 (MQ=255)                                                      cTCGCCACGTCATTTCTCGGCGTTGCGTTAGGCTTATTTGATTATCTGGCTGATTTATTTCAGCGTTCa    <  1:2791876/69‑1 (MQ=255)                                                      cTCGCCACGTCATTTCTCGGCGTTGCGTTAGGCTTATTTGATTATCTGGCTGATTTATTTCAGCGTTCa    <  1:2790353/69‑1 (MQ=255)                                                       tCGCCACGTCATTTCTCGGCGTTGCGTTAGGCTTATTTGATTATCTGGCTGATTTATTTCAGCGTTCa    <  1:1220904/68‑1 (MQ=255)                                                         gCCACGTCATTTCTCGGCGTTGCGTTAGGCTTATTTGAtt                              >  1:2955289/1‑40 (MQ=255)                                                         gCCACGTCATTTCTCGGCGTTGCGTTAGGCTTATTTGAtt                              >  1:811621/1‑40 (MQ=255)                                                         gCCACGTCATTTCTCGGCGTTGCGTTAGGCTTATTTGAtt                              >  1:672591/1‑40 (MQ=255)                                                         gCCACGTCATTTCTCGGCGTTGCGTTAGGCTTATTTGAtt                              >  1:1428495/1‑40 (MQ=255)                                                         gCCACGTCATTTCTCGGCGTTGCGTTAGGCTTATTTGAtt                              >  1:519850/1‑40 (MQ=255)                                                            aCGTCCTTCCTCGTCGTAGCGTTAGGCTTATTTGATTATCTGGCAGATTTGTTTCAGCGTTCaaa  <  1:1574527/65‑1 (MQ=25)                                                            aCGTCCTTCCTCGGCGTAGCGTTAGGCTTATTTGATTATCTGGCAGATTTGTTTCAGCGTTCaaa  <  1:1945261/65‑1 (MQ=25)                                                            aCGTCCTTCCTCGGCGTAGCGTTAGGCTTATTTGATTATCTGGCAGATTTGTTTCAGCGTTCaaa  <  1:713655/65‑1 (MQ=25)                                                            aCGTCCTTCCTCGGCGTAGCGTTAGGCTTATTTGATTATCTGGCAGATTTGTTTCAGCGTTCaaa  <  1:1354473/65‑1 (MQ=25)                                                            aCGTCCTTCCTCGGCGTAGCGTTAGGCTTATTTGATTATCTGGCAGATTTGTTTCAGCGTTCaaa  <  1:944699/65‑1 (MQ=25)                                                            aCGTCCTTCCTCGGCGTAGCGTTAGGCTTATTTGATTATCTGGCAGATTTGTTTCAGCGTTCaaa  <  1:1828988/65‑1 (MQ=25)                                                            aCGTCCTTCCTCGGCGTAGCGTTAGGCTTATTTGATTATCTGGCAGATTTGTTTCAGCGTTCaaa  <  1:191684/65‑1 (MQ=25)                                                            aCGTCCTTCCTCGGCGTAGCGTTAGGCTTATTTGATTATCTGGCAGATTTGTTTCAGCGTTCaaa  <  1:3334992/65‑1 (MQ=25)                                                            aCGTCCTTCCTCGGCGTAGCGTTAGGCTTATTTGATTATCTGGCAGATTTGTTTCAGCGTTCaaa  <  1:2257788/65‑1 (MQ=25)                                                            aCGTCCTTCCTCGGCGTAGCGTTAGGCTTATTTGATTATCTGGCAGATTTGTTTCAGCGTTCaaa  <  1:2312632/65‑1 (MQ=25)                                                            aCGTCCTTCCTCGGCGTAGCGTTAGGCTTATTTGATTATCTGGCAGATTTGTTTCAGCGTTCaaa  <  1:3230574/65‑1 (MQ=25)                                                            aCGTCCTTCCTCGGCGTAGCGTTAGGCTTATTTGATTATCTGGCAGATTTGTTTCAGCGTTCaaa  <  1:2320293/65‑1 (MQ=25)                                                            aCGTCCTTCCTCGGCGTAGCGTTAGGCTTATTTGATTATCTGGCAGATTTGTTTCAGCGTTCaaa  <  1:2485707/65‑1 (MQ=25)                                                            aCGTCCTTCCTCGGCGTAGCGTTAGGCTTATTTGATTATCTGGCAGATTTGTTTCAGCGTTCaaa  <  1:3082391/65‑1 (MQ=25)                                                            aCGTCCTTCCTCGGCGTAGCGTTAGGCTTATTTGATTATCTGGCAGATTTGTTTCAGCGTTCaaa  <  1:3021034/65‑1 (MQ=25)                                                            aCGTCCTTCCTCGGCGTAGCGTTAGGCTTATTTGATTATCTGGCAGATTTGTTTCAGCGTTCaaa  <  1:2526049/65‑1 (MQ=25)                                                            aCGTCCTTCCTCGGCGTAGCGTTAGGCTTATTTGATTATCTGGCAGATTTGTTTCAGCGTTCaaa  <  1:2762721/65‑1 (MQ=25)                                                            aCGTCCTTCCTCGGCGTAGCGTTAGGCTTATTTGATTATCTGGCAGATTTGTTTCAGCGTTCaaa  <  1:2693473/65‑1 (MQ=25)                                                            aCGTCCTTCCTCGGCGTAGCGTTAGGCTTATTTGATTATCTGGCAGATTTGTTTCAGCCTTCaaa  <  1:2674440/65‑1 (MQ=21)                                                                    |                                                         TGGTGGCCTCTCCGCATGTTGAGCTGGCAGTGCATTTATTTGCTGATTTAGCCCTCGCCACGTCATTTCTCGGCGTTGCGTTAGGCTTATTTGATTATCTGGCTGATTTATTTCAGCGTTCAAA  >  minE/1266261‑1266384

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