Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	1,266,325	A→C	37.0%	I154L (ATT→CTT)	tyrP →	tyrosine transporter

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	1,266,325	0	A	C	37.0%	55.3 / 70.2	81	I154L (ATT→CTT)	tyrP	tyrosine transporter
Reads supporting (aligned to +/- strand):  ref base A (13/38);  new base C (10/20);  total (23/58)
Fisher's exact test for biased strand distribution p-value = 4.57e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GAAATGGTGGCCTCTCCGCATGTTGAGCTGGCAGTGCATTTATTTGCTGATTTAGCCCTCGCCACGTCATTTCTCGGCGTTGCGTTAGGCTTATTTGATTATCTGGCTGATTTATTTCAGCGTTCAAA  >  minE/1266257‑1266384                                                                     |                                                            gAAATGGTGGCCTCTCCGCATGTTGAGCTGGCAGTGCATTTATTTGCTGATTTAGCCCTCGCCACTTCa                                                             <  1:636233/69‑1 (MQ=255) gAAATGGTGGCCTCTCCGCATGTTGAGCTGGCAGTGCATTTATTTGCTGATTTAGCCCTCGCCACGTCa                                                             <  1:214438/69‑1 (MQ=255) gAAATGGTGGCCTCTCCGCATGTTGAGCTGGCAGTGCATTTATTTGCTGATTTAGCCCTCGCCACGTCa                                                             <  1:2799503/69‑1 (MQ=255)     tggtggCCTCTCCGCATGTTGAGCTGGCAGTGCATTTATTTGCTGATTTAGCCCTCGCCACGTCATTtc                                                         >  1:2345942/1‑69 (MQ=255)                 cGCATGTTGAGCTGGCAGTGCATTTATTTGCTGATTTAGCCCTCGCCACGTCATTTCTCGgcgttgc                                               <  1:2310668/67‑1 (MQ=255)                 cGCATGTTGAGCTGGCAGTGCATTTATTTGCTGATTTAGCCCTCGCCACGTCATTTCTCGgcgttgc                                               <  1:2857440/67‑1 (MQ=255)                 cGCATGTTGAGCTGGCAGTGCATTTATTTGCTGATTTAGCCCTCGCCACGTCATTTCTCGgcgttgc                                               <  1:2029651/67‑1 (MQ=255)                 cGCATGTTGAGCTGGCAGTGCATTTATTTGCTGATTTAGCCCTCGCCACGTCATTTCTCGgcgttgc                                               <  1:915419/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:2255958/68‑1 (MQ=255)                                 aGTGCATTTATTTGCTGATTTAGCCCTCGCCACGTCATTTCTCGGCGTTGCGTTAGGCTTATTTGAtt                              <  1:28347/68‑1 (MQ=255)                                 aGTGCATTTATTTGCTGATTTAGCCCTCGCCACGTCATTTCTCGGCGTTGCGTTAGGCTTATTTGAtt                              <  1:1931726/68‑1 (MQ=255)                                 aGTGCATTTATTTGCTGATTTAGCCCTCGCCACGTCATTTCTCGGCGTTGCGTTAGGCTTATTTGAtt                              <  1:1846031/68‑1 (MQ=255)                                 aGTGCATTTATTTGCTGATTTAGCCCTCGCCACGTCATTTCTCGGCGTTGCGTTAGGCTTATTTGAtt                              <  1:2500210/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:921447/61‑1 (MQ=255)                                     cATTTATTTGCTGATTTAGCCCTCGCCACGTCATTTCTCGGCGTTGCGTTAGGCTTATTTg                                 <  1:2533200/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:629184/1‑69 (MQ=255)                                     cATTTATTTGCTGATTTAGCCCTCGCCACGTCATTTCTCGGCGTTGCGTTAGGCTTATTTGATTATCTg                         >  1:2384551/1‑69 (MQ=255)                                     cATTTATTTGCTGATTTAGCCCTCGCCACGTCATTTCTCGGCGTTGCGTTAGGCTTATTTGATTATCTg                         >  1:2881429/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:2193034/1‑66 (MQ=255)                                                tGATTTAGCCCTCGCCACGTCCTTCCTCGGCGTAGCGTTAGGCTTATTTGATTATCTGGCAGatttgt               >  1:1516468/1‑66 (MQ=255)                                                tGATTTAGCCCTCGCCACGTCCTTCCTCGGCGTAGCGTTAGGCTTATTTGATTATCTGGCAGatttgt               >  1:2662378/1‑66 (MQ=255)                                                tGATTTAGCCCTCGCCACGTCCTTCCTCGGCGTAGCGTTAGGCTTATTTGATTATCTGGCAGatttgt               >  1:2630148/1‑66 (MQ=255)                                                tGATTTAGCCCTCGCCACGTCCTTCCTCGGCGTAGCGTTAGGCTTATTTGATTATCTGGCAGatttgt               >  1:3218342/1‑66 (MQ=255)                                                tGATTTAGCCCTCGCCACGTCCTTCCTCGGCGTAGCGTTAGGCTTATTTGATTATCTGGCAGatttgt               >  1:2411078/1‑66 (MQ=255)                                                tGATTTAGCCCTCGCCACGTCCTTCCTCGGCGTAGCGTTAGGCTTATTTGATTATCTGGCAGatttgt               >  1:43571/1‑66 (MQ=255)                                                tGATTTAGCCCTCGCCACGTCCTTCCTCGGCGTAGCGTTAGGCTTATTTGATTATCTGGCAGatttgt               >  1:3274818/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:3144134/69‑1 (MQ=255)                                                         ccTCGCCACGTCATTTCTCGGCGTTGCGTTAGGCTTATTTGATTATCTGGCTGATTTATTTCAGCGTTc     <  1:623533/69‑1 (MQ=255)                                                         ccTCGCCACGTCATTTCTCGGCGTTGCGTTAGGCTTATTTGATTATCTGGCTGATTTATTTCAGCGTTc     <  1:3239450/69‑1 (MQ=255)                                                         ccTCGCCACGTCATTTCTCGGCGTTGCGTTAGGCTTATTTGATTATCTGGCTGATTTATTTCAGCGTTc     <  1:1150597/69‑1 (MQ=255)                                                         ccTCGCCACGTCATTTCTCGGCGTTGCGTTAGGCTTATTTGATTATCTGGCTGATTTATTTCAGCGTTc     <  1:188318/69‑1 (MQ=255)                                                         ccTCGCCACGTCATTTCTCGGCGTTGCGTTAGGCTTATTTGATTATCTGGCTGATTTATTTCAGCGTTc     <  1:2699817/69‑1 (MQ=255)                                                         ccTCGCCACGTCATTTCTCGGCGTTGCGTTAGGCTTATTTGATTATCTGGCTGATTTATTTCAGCGTTc     <  1:1080130/69‑1 (MQ=255)                                                          cTCGCCACGTCATTTCTCGGCGTTGCGTTAGGCTTATTTGATTATCTGGCTGATTTATTTCAGCGTTCa    <  1:820825/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:3136750/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:944699/65‑1 (MQ=25)                                                                aCGTCCTTCCTCGGCGTAGCGTTAGGCTTATTTGATTATCTGGCAGATTTGTTTCAGCGTTCaaa  <  1:713655/65‑1 (MQ=25)                                                                aCGTCCTTCCTCGGCGTAGCGTTAGGCTTATTTGATTATCTGGCAGATTTGTTTCAGCGTTCaaa  <  1:1354473/65‑1 (MQ=25)                                                                aCGTCCTTCCTCGGCGTAGCGTTAGGCTTATTTGATTATCTGGCAGATTTGTTTCAGCGTTCaaa  <  1:2762721/65‑1 (MQ=25)                                                                aCGTCCTTCCTCGGCGTAGCGTTAGGCTTATTTGATTATCTGGCAGATTTGTTTCAGCGTTCaaa  <  1:1828988/65‑1 (MQ=25)                                                                aCGTCCTTCCTCGGCGTAGCGTTAGGCTTATTTGATTATCTGGCAGATTTGTTTCAGCGTTCaaa  <  1:191684/65‑1 (MQ=25)                                                                aCGTCCTTCCTCGGCGTAGCGTTAGGCTTATTTGATTATCTGGCAGATTTGTTTCAGCGTTCaaa  <  1:1945261/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:2320293/65‑1 (MQ=25)                                                                aCGTCCTTCCTCGGCGTAGCGTTAGGCTTATTTGATTATCTGGCAGATTTGTTTCAGCGTTCaaa  <  1:3230574/65‑1 (MQ=25)                                                                aCGTCCTTCCTCGGCGTAGCGTTAGGCTTATTTGATTATCTGGCAGATTTGTTTCAGCGTTCaaa  <  1:2485707/65‑1 (MQ=25)                                                                aCGTCCTTCCTCGGCGTAGCGTTAGGCTTATTTGATTATCTGGCAGATTTGTTTCAGCGTTCaaa  <  1:2526049/65‑1 (MQ=25)                                                                aCGTCCTTCCTCGGCGTAGCGTTAGGCTTATTTGATTATCTGGCAGATTTGTTTCAGCGTTCaaa  <  1:3082391/65‑1 (MQ=25)                                                                aCGTCCTTCCTCGGCGTAGCGTTAGGCTTATTTGATTATCTGGCAGATTTGTTTCAGCGTTCaaa  <  1:3021034/65‑1 (MQ=25)                                                                aCGTCCTTCCTCGGCGTAGCGTTAGGCTTATTTGATTATCTGGCAGATTTGTTTCAGCGTTCaaa  <  1:2693473/65‑1 (MQ=25)                                                                aCGTCCTTCCTCGGCGTAGCGTTAGGCTTATTTGATTATCTGGCAGATTTGTTTCAGCCTTCaaa  <  1:2674440/65‑1 (MQ=21)                                                                     |                                                            GAAATGGTGGCCTCTCCGCATGTTGAGCTGGCAGTGCATTTATTTGCTGATTTAGCCCTCGCCACGTCATTTCTCGGCGTTGCGTTAGGCTTATTTGATTATCTGGCTGATTTATTTCAGCGTTCAAA  >  minE/1266257‑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: ‑