Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	1,266,364	T→A	61.0%	*167R (TGA→AGA)	tyrP →	tyrosine transporter

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	1,266,364	0	T	A	61.0%	46.1 / 74.4	82	*167R (TGA→AGA)	tyrP	tyrosine transporter
Reads supporting (aligned to +/- strand):  ref base T (4/28);  new base A (10/40);  total (14/68)
Fisher's exact test for biased strand distribution p-value = 5.49e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TATTTGCTGATTTAGCCCTCGCCACGTCATTTCTCGGCGTTGCGTTAGGCTTATTTGATTATCTGGCTGATTTATTTCAGCGTTCAAATACCGTTGGTGGACGGTTGCAAACTGGTGCAATTACGTTT  >  minE/1266297‑1266424                                                                    |                                                             tatttGCTGATTTAGCCCTCGCCACGTCATTTCTCGGCGTTGCGTTAGGCTTATTTGATTATCTGGCTg                                                             >  1:3294520/1‑69 (MQ=255)        tGATTTAGCCCTCGCCACGTCCTTCCTCGGCGTAGCGTTAGGCTTATTTGATTATCTGGCAGatttgt                                                       >  1:3274818/1‑66 (MQ=255)        tGATTTAGCCCTCGCCACGTCCTTCCTCGGCGTAGCGTTAGGCTTATTTGATTATCTGGCAGatttgt                                                       >  1:2411078/1‑66 (MQ=255)        tGATTTAGCCCTCGCCACGTCCTTCCTCGGCGTAGCGTTAGGCTTATTTGATTATCTGGCAGatttgt                                                       >  1:2193034/1‑66 (MQ=255)        tGATTTAGCCCTCGCCACGTCCTTCCTCGGCGTAGCGTTAGGCTTATTTGATTATCTGGCAGatttgt                                                       >  1:2160319/1‑66 (MQ=255)        tGATTTAGCCCTCGCCACGTCCTTCCTCGGCGTAGCGTTAGGCTTATTTGATTATCTGGCAGatttgt                                                       >  1:3218342/1‑66 (MQ=255)        tGATTTAGCCCTCGCCACGTCCTTCCTCGGCGTAGCGTTAGGCTTATTTGATTATCTGGCAGatttgt                                                       >  1:2662378/1‑66 (MQ=255)        tGATTTAGCCCTCGCCACGTCCTTCCTCGGCGTAGCGTTAGGCTTATTTGATTATCTGGCAGatttgt                                                       >  1:43571/1‑66 (MQ=255)        tGATTTAGCCCTCGCCACGTCCTTCCTCGGCGTAGCGTTAGGCTTATTTGATTATCTGGCAGatttgt                                                       >  1:2630148/1‑66 (MQ=255)        tGATTTAGCCCTCGCCACGTCCTTCCTCGGCGTAGCGTTAGGCTTATTTGATTATCTGGCAGatttgt                                                       >  1:1516468/1‑66 (MQ=255)        tGATTTAGCCCTCGCCACGTCCTTCCTCGGCGTAGCGTTAGGCTTATTTGATTATCTGGCAGatttg                                                        >  1:2570144/1‑66 (MQ=255)                 ccTCGCCACGTCATTTCTCGGCGTTGCGTTAGGCTTATTTGATTATCTGGCTTATTTATTTCAGCGTTc                                             <  1:3399623/69‑1 (MQ=255)                 ccTCGCCACGTCATTTCTCGGCGTTGCGTTAGGCTTATTTGATTATCTGGCTGATTTATTTCAGCGTTc                                             <  1:623533/69‑1 (MQ=255)                 ccTCGCCACGTCATTTCTCGGCGTTGCGTTAGGCTTATTTGATTATCTGGCTGATTTATTTCAGCGTTc                                             <  1:3239450/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)                 ccTCGCCACGTCATTTCTCGGCGTTGCGTTAGGCTTATTTGATTATCTGGCTGATTTATTTCAGCGTTc                                             <  1:1150597/69‑1 (MQ=255)                 ccTCGCCACGTCATTTCTCGGCGTTGCGTTAGGCTTATTTGATTATCTGGCTGATTTATTTCAGCGTTc                                             <  1:3144134/69‑1 (MQ=255)                  cTCGCCACGTCATTTCTCGGCGTTGCGTTAGGCTTATTTGATTATCTGGCTGATTTATTTCAGCGTTCa                                            <  1:820825/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)                  cTCGCCACGTCATTTCTCGGCGTTGCGTTAGGCTTATTTGATTATCTGGCTGATTTATTTCAGCGTTCa                                            <  1:3136750/69‑1 (MQ=255)                  cTCGCCACGTCATTTCTCGGCGTTGCGTTAGGCTTATTTGATTATCTGGCTGATTTATTTCAGCGTTCa                                            <  1:3214387/69‑1 (MQ=255)                  cTCGCCACGTCATTTCTCGGCGTTGCGTTAGGCTTATTTGATTATCTGGCTGATTTATTTCAGCGTTCa                                            <  1:1896990/69‑1 (MQ=255)                  cTCGCCACGTCATTTCTCGGCGTTGCGTTAGGCTTATTTGATTATCTGGCTGATTTATTTCAGCGTTCa                                            <  1:2541142/69‑1 (MQ=255)                  cTCGCCACGTCATTTCTCGGCGTTGCGTTAGGCTTATTTGATTATCTGGCTGATTTATTTCAGCGTTCa                                            <  1:1601620/69‑1 (MQ=255)                   tCGCCACGTCATTTCTCGGCGTTGCGTTAGGCTTATTTGATTATCTGGCTGATTTATTTCAGCGTTCa                                            <  1:1220904/68‑1 (MQ=255)                        aCGTCCTTCCTCGTCGTAGCGTTAGGCTTATTTGATTATCTGGCAGATTTGTTTCAGCGTTCaaa                                          <  1:1574527/65‑1 (MQ=25)                        aCGTCCTTCCTCGGCGTAGCGTTAGGCTTATTTGATTATCTGGCAGATTTGTTTCAGCGTTCaaa                                          <  1:3230574/65‑1 (MQ=25)                        aCGTCCTTCCTCGGCGTAGCGTTAGGCTTATTTGATTATCTGGCAGATTTGTTTCAGCGTTCaaa                                          <  1:2762721/65‑1 (MQ=25)                        aCGTCCTTCCTCGGCGTAGCGTTAGGCTTATTTGATTATCTGGCAGATTTGTTTCAGCGTTCaaa                                          <  1:3021034/65‑1 (MQ=25)                        aCGTCCTTCCTCGGCGTAGCGTTAGGCTTATTTGATTATCTGGCAGATTTGTTTCAGCGTTCaaa                                          <  1:3082391/65‑1 (MQ=25)                        aCGTCCTTCCTCGGCGTAGCGTTAGGCTTATTTGATTATCTGGCAGATTTGTTTCAGCGTTCaaa                                          <  1:2693473/65‑1 (MQ=25)                        aCGTCCTTCCTCGGCGTAGCGTTAGGCTTATTTGATTATCTGGCAGATTTGTTTCAGCGTTCaaa                                          <  1:3334992/65‑1 (MQ=25)                        aCGTCCTTCCTCGGCGTAGCGTTAGGCTTATTTGATTATCTGGCAGATTTGTTTCAGCGTTCaaa                                          <  1:713655/65‑1 (MQ=25)                        aCGTCCTTCCTCGGCGTAGCGTTAGGCTTATTTGATTATCTGGCAGATTTGTTTCAGCGTTCaaa                                          <  1:944699/65‑1 (MQ=25)                        aCGTCCTTCCTCGGCGTAGCGTTAGGCTTATTTGATTATCTGGCAGATTTGTTTCAGCGTTCaaa                                          <  1:2526049/65‑1 (MQ=25)                        aCGTCCTTCCTCGGCGTAGCGTTAGGCTTATTTGATTATCTGGCAGATTTGTTTCAGCGTTCaaa                                          <  1:2485707/65‑1 (MQ=25)                        aCGTCCTTCCTCGGCGTAGCGTTAGGCTTATTTGATTATCTGGCAGATTTGTTTCAGCGTTCaaa                                          <  1:1354473/65‑1 (MQ=25)                        aCGTCCTTCCTCGGCGTAGCGTTAGGCTTATTTGATTATCTGGCAGATTTGTTTCAGCGTTCaaa                                          <  1:1828988/65‑1 (MQ=25)                        aCGTCCTTCCTCGGCGTAGCGTTAGGCTTATTTGATTATCTGGCAGATTTGTTTCAGCGTTCaaa                                          <  1:2320293/65‑1 (MQ=25)                        aCGTCCTTCCTCGGCGTAGCGTTAGGCTTATTTGATTATCTGGCAGATTTGTTTCAGCGTTCaaa                                          <  1:2312632/65‑1 (MQ=25)                        aCGTCCTTCCTCGGCGTAGCGTTAGGCTTATTTGATTATCTGGCAGATTTGTTTCAGCGTTCaaa                                          <  1:191684/65‑1 (MQ=25)                        aCGTCCTTCCTCGGCGTAGCGTTAGGCTTATTTGATTATCTGGCAGATTTGTTTCAGCGTTCaaa                                          <  1:2257788/65‑1 (MQ=25)                        aCGTCCTTCCTCGGCGTAGCGTTAGGCTTATTTGATTATCTGGCAGATTTGTTTCAGCGTTCaaa                                          <  1:1945261/65‑1 (MQ=25)                        aCGTCCTTCCTCGGCGTAGCGTTAGGCTTATTTGATTATCTGGCAGATTTGTTTCAGCCTTCaaa                                          <  1:2674440/65‑1 (MQ=21)                                 ctcGGCGTTGCGTTAGGCTTATTTGTTTATCTGGCTGATTTATTTCTGCGTTc                                             <  1:3289956/53‑1 (MQ=255)                                      cgtagcgttAGGCTTATTTGATTATCTGGCAGATTTGTTTCAGCGTTCAAATACCGTTGGTGGACga                          <  1:3383131/67‑2 (MQ=38)                                      cgtagcgttAGGCTTATTTGATTATCTGGCAGATTTGTTTCAGCGTTCAAATACCGTTGGTGGACga                          <  1:1087926/67‑2 (MQ=38)                                      cgtagcgttAGGCTTATTTGATTATCTGGCAGATTTGTTTCAGCGTTCAAATACCGTTGGTGGACga                          <  1:1138731/67‑2 (MQ=38)                                      cgtagcgttAGGCTTATTTGATTATCTGGCAGATTTGTTTCAGCGTTCAAATACCGTTGGTGGACga                          <  1:1513930/67‑2 (MQ=38)                                      cgtagcgttAGGCTTATTTGATTATCTGGCAGATTTGTTTCAGCGTTCAAATACCGTTGGTGGACga                          <  1:640652/67‑2 (MQ=38)                                      cgtagcgttAGGCTTATTTGATTATCTGGCAGATTTGTTTCAGCGTTCAAATACCGTTGGTGGACga                          <  1:1613412/67‑2 (MQ=38)                                      cgtagcgttAGGCTTATTTGATTATCTGGCAGATTTGTTTCAGCGTTCAAATACCGTTGGTGGACga                          <  1:619102/67‑2 (MQ=38)                                      cgtagcgttAGGCTTATTTGATTATCTGGCAGATTTGTTTCAGCGTTCAAATACCGTTGGTGGACga                          <  1:613463/67‑2 (MQ=38)                                      cgtagcgttAGGCTTATTTGATTATCTGGCAGATTTGTTTCAGCGTTCAAATACCGTTGGTGGACga                          <  1:3428499/67‑2 (MQ=38)                                      cgtagcgttAGGCTTATTTGATTATCTGGCAGATTTGTTTCAGCGTTCAAATACCGTTGGTGGACga                          <  1:2167024/67‑2 (MQ=38)                                      cgtagcgttAGGCTTATTTGATTATCTGGCAGATTTGTTTCAGCGTTCAAATACCGTTGGTGGACga                          <  1:1717096/67‑2 (MQ=38)                                      cgtagcgttAGGCTTATTTGATTATCTGGCAGATTTGTTTCAGCGTTCAAATACCGTTGGTGGACga                          <  1:3302032/67‑2 (MQ=38)                                      cgtagcgttAGGCTTATTTGATTATCTGGCAGATTTGTTTCAGCGTTCAAATACCGTTGGTGGACga                          <  1:1774797/67‑2 (MQ=38)                                      cgtagcgttAGGCTTATTTGATTATCTGGCAGATTTGTTTCAGCGTTCAAATACCGTTGGTGGACga                          <  1:2031850/67‑2 (MQ=38)                                      cgtagcgttAGGCTTATTTGATTATCTGGCAGATTTGTTTCAGCGTTCAAATACCGTTGGTGGACga                          <  1:1807307/67‑2 (MQ=38)                                      cgtagcgttAGGCTTATTTGATTATCTGGCAGATTTGTTTCAGCGTTCAAATACCGTTGGTGGACga                          <  1:208607/67‑2 (MQ=38)                                      cgtagcgttAGGCTTATTTGATTATCTGGCAGATTTGTTTCAGCGTTCAAATACCGTTGGTGGACga                          <  1:2856762/67‑2 (MQ=38)                                      cgtagcgttAGGCTTATTTGATTATCTGGCAGATTTGTTTCAGCGTTCAAATACCGTTGGTGGACga                          <  1:3253801/67‑2 (MQ=38)                                      cgtagcgttAGGCTTATTTGATTATCTGGCAGATTTGTTTCAGCGTTCAAATACCGTTGGTGGACga                          <  1:3238596/67‑2 (MQ=38)                                          gcgttAGGCTTATTTGATTATCTGGCTGATTTATTTCAGCGTTCAAATACCGTTGGTGGACGGTTGCaa                    <  1:2484360/69‑1 (MQ=255)                                          gcgttAGGCTTATTTGATTATCTGGCTGATTTATTTCAGCGTTCAAATACCGTTGGTGGACGGTTGCaa                    <  1:3075446/69‑1 (MQ=255)                                          gcgttAGGCTTATTTGATTATCTGGCTGATTTATTTCAGCGTTCAAATACCGTTGGTGGACGGTTGCaa                    <  1:459688/69‑1 (MQ=255)                                          gcgttAGGCTTATTTGATTATCTGGCTGATTTATTTCAGCGTTCAAATACCGTTGGTGGACGGTTGCaa                    <  1:3254639/69‑1 (MQ=255)                                          gcgttAGGCTTATTTGATTATCTGGCTGATTTATTTCAGCGTTCAAATACCGTTGGTGGACGGTTGCaa                    <  1:3315507/69‑1 (MQ=255)                                            gttAGGCTTATTTGATTATCTGGCTGATTTATTTCAGCGTTCAAATACCGTTGGTGGACGGTTGCaaa                   <  1:873970/68‑1 (MQ=255)                                            gttAGGCTTATTTGATTATCTGGCTGATTTATTTCAGCGTTCAAATACCGTTGGTGGACGGTTGCaaa                   <  1:2395379/68‑1 (MQ=255)                                            gttAGGCTTATTTGATTATCTGGCTGATTTATTTCAGCGTTCAAATACCGTTGGTGGACGGTTGCaaa                   <  1:944445/68‑1 (MQ=255)                                            gttAGGCTTATTTGATTATCTGGCTGATTTATTTCAGCGTTCAAATACCGTTGGTGGACGGTTGCaaa                   <  1:1677091/68‑1 (MQ=255)                                                ggCTTATTTGATTATCTGGCAGATTTGTTTCAGCGTTCaaa                                          <  1:2572609/41‑1 (MQ=37)                                                ggCTTATTTGATTATCTGGCAGATTTGTTTCAGCGTTCAAATACCGTTGGTGGACga                          <  1:570544/57‑2 (MQ=38)                                                  cTTATTTGATTATCTGGCTGATTTATTTCAGCGTTCAAATACCGTTGGTGGACGGTTGCAAACTGGTGc            >  1:3388959/1‑69 (MQ=255)                                                            tATCTGGCTGATTTATTTCAGCGTTCAAATACCGTTGGTGGACGGTTGCAAACTGGTGCAATTACGttt  >  1:20582/1‑69 (MQ=255)                                                            tATCTGGCTGATTTATTTCAGCGTTCAAATACCGTTGGTGGACGGTTGCAAACTGGTGCAATTACGttt  >  1:1656974/1‑69 (MQ=255)                                                                    |                                                             TATTTGCTGATTTAGCCCTCGCCACGTCATTTCTCGGCGTTGCGTTAGGCTTATTTGATTATCTGGCTGATTTATTTCAGCGTTCAAATACCGTTGGTGGACGGTTGCAAACTGGTGCAATTACGTTT  >  minE/1266297‑1266424

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