Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	2,505,784	T→G	45.1%	R80R (CGT→CGG)	yicJ →	predicted transporter

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	2,505,784	0	T	G	45.1%	‑3.1 / 91.9	80	R80R (CGT→CGG)	yicJ	predicted transporter
Reads supporting (aligned to +/- strand):  ref base T (9/35);  new base G (17/19);  total (26/54)
Fisher's exact test for biased strand distribution p-value = 1.62e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.00e-01

TCCGATCCTTGCATGGGGTTGTTGGCCGATCGAACGCGCTCTCGCTGGGGTAAATTTCGTCCGTGGGTACTGTTTGGCGCACTGCCATTCGGGATCGTCTGTGTACTGGCCTATAGCACGC  >  minE/2505725‑2505845                                                            |                                                              tCCGATCCTTGCATGGGGTTGTTGGCCGATCGAACGCGCTCTCGCTGGGGTAAATTTCGTCCGTGGGTa                                                      >  1:1169619/1‑69 (MQ=255)   cGATCCTTGCATGGGGTTGTTGGCCGATCGAACGCGCTCTCGCTGGGGTAAATTTCGGCCGTGGGTGCt                                                    >  1:2695105/1‑69 (MQ=255)   cGATCCTTGCATGGGGTTGTTGGCCGATCGAACGCGCTCTCGCTGGGGTAAATTTCGGCCGTGGGTGCt                                                    >  1:3195905/1‑69 (MQ=255)     aTCCTTGCATGGGGTTGTTGGCCGATCGAACGCGCTCTCGCTGGGGTAAATTTCGTCCGTGGGTACTg                                                   <  1:1643904/68‑1 (MQ=255)     aTCCTTGCATGGGGTTGTTGGCCGATCGAACGCGCTCTCGCTGGGGTAAATTTCGTCCGTGGGTACTg                                                   <  1:3162068/68‑1 (MQ=255)     aTCCTTGCATGGGGTTGTTGGCCGATCGAACGCGCTCTCGCTGGGGTAAATTTCGTCCGTGGGTACTg                                                   <  1:573449/68‑1 (MQ=255)     aTCCTTGCATGGGGTTGTTGGCCGATCGAACGCGCTCTCGCTGGGGTAAATTTCGTCCGTGGGTACTg                                                   <  1:1436174/68‑1 (MQ=255)     aTCCTTGCATGGGGTTGTTGGCCGATCGAACGCGCTCTCGCTGGGGTAAATTTCGTCCGTGGGTACTg                                                   <  1:1389874/68‑1 (MQ=255)     aTCCTTGCATGGGGTTGTTGGCCGATCGAACGCGCTCTCGCTGGGGTAAATTTCGTCCGTGGGTACTg                                                   <  1:1147256/68‑1 (MQ=255)     aTCCTTGCATGGGGTTGTTGGCCGATCGAACGCGCTCTCGCTGGGGTAAATTTCGTCCGTGGGTACTg                                                   <  1:1132447/68‑1 (MQ=255)       ccTTGCATGGGGTTGTTGGCCGATCGAACGCGCTCTCGCTGGGGTAAATTTCGTCCGTGGGTACTg                                                   <  1:379648/66‑1 (MQ=255)              tGGGGTTGTTGGCCGATCGAACGCGCTCTCGCTGGGGTAAATTTCGTCCGTGGGTACTGTTTGGCGCa                                          <  1:283863/68‑1 (MQ=255)              tGGGGTTGTTGGCCGATCGAACGCGCTCTCGCTGGGGTAAATTTCGTCCGTGGGTACTGTTTGGCGCa                                          <  1:112376/68‑1 (MQ=255)              tGGGGTTGTTGGCCGATCGAACGCGCTCTCGCTGGGGTAAATTTCGTCCGTGGGTACTGTTTGGCGCa                                          <  1:2398177/68‑1 (MQ=255)              tGGGGTTGTTGGCCGATCGAACGCGCTCTCGCTGGGGTAAATTTCGTCCGTGGGTACTGTTTGGCGCa                                          <  1:749989/68‑1 (MQ=255)                   ttgttgGCCGATCGAACGCGCTCTCGCTGGGGTAAATTTCGTCCGTGGGTACTGTTTGGCGCACTGCCa                                    >  1:1536908/1‑69 (MQ=255)                          ccGATCGAACGCGCTCTCGCTGGGGTAAATTTCGTCCGTGGGTACTGTTTGGCGCACTGCCATTCggg                              <  1:1359924/68‑1 (MQ=255)                          ccGATCGAACGCGCTCTCGCTGGGGTAAATTTCGTCCGTGGGTACTGTTTGGCGCACTGCCATTCggg                              <  1:2419671/68‑1 (MQ=255)                          ccGATCGAACGCGCTCTCGCTGGGGTAAATTTCGTCCGTGGGTACTGTTTGGCGCACTGCCATTCggg                              <  1:891462/68‑1 (MQ=255)                              tCGAACGCGCTCTCGCTGGGGTAAATTTCGGCCGTGGGTGCTGttt                                                <  1:2823034/46‑1 (MQ=255)                              tCGAACGCGCTCTCGCTGGGGTAAATTTCGGCCGTGGGTGCTGttt                                                <  1:2416915/46‑1 (MQ=255)                              tCGAACGCGCTCTCGCTGGGGTAAATTTCGGCCGTGGGTGCTGttt                                                <  1:2534627/46‑1 (MQ=255)                              tCGAACGCGCTCTCGCTGGGGTAAATTTCGGCCGTGGGTGCTGttt                                                <  1:906794/46‑1 (MQ=255)                              tCGAACGCGCTCTCGCTGGGGTAAATTTCGGCCGTGGGTGCTGttt                                                <  1:1037932/46‑1 (MQ=255)                              tCGAACGCGCTCTCGCTGGGGTAAATTTCGGCCGTGGGTGCTGttt                                                <  1:3156977/46‑1 (MQ=255)                              tCGAACGCGCTCTCGCTGGGGTAAATTTCGGCCGTGGGTGCTGttt                                                <  1:1067643/46‑1 (MQ=255)                              tCGAACGCGCTCTCGCTGGGGTAAATTTCGGCCGTGGGTGCTGttt                                                <  1:94988/46‑1 (MQ=255)                              tCGAACGCGCTCTCGCTGGGGTAAATTTCGGCCGTGGGTGCTGttt                                                <  1:2120763/46‑1 (MQ=255)                              tCGAACGCGCTCTCGCTGGGGTAAATTTCGGCCGTGGGTGCTGttt                                                <  1:1885993/46‑1 (MQ=255)                              tCGAACGCGCTCTCGCTGGGGTAAATTTCGGCCGTGGGTGCTGttt                                                <  1:3371449/46‑1 (MQ=255)                              tCGAACGCGCTCTCGCTGGGGTAAATTTCGGCCGTGGGTGCTGttt                                                <  1:1844615/46‑1 (MQ=255)                              tCGAACGCGCTCTCGCTGGGGTAAATTTCGGCCGTGGGTGCTGttt                                                <  1:1841407/46‑1 (MQ=255)                              tCGAACGCGCTCTCGCTGGGGTAAATTTCGGCCGTGGGTGCTGttt                                                <  1:594125/46‑1 (MQ=255)                              tCGAACGCGCTCTCGCTGGGGTAAATTTCGGCCGTGGGTGCTGttt                                                <  1:600748/46‑1 (MQ=255)                              tCGAACGCGCTCTCGCTGGGGTAAATTTCGGCCGTGGGTGCTGttt                                                <  1:639460/46‑1 (MQ=255)                              tCGAACGCGCTCTCGCTGGGGTAAATTTCGGCCGTGGGTGCTGttt                                                <  1:713797/46‑1 (MQ=255)                              tCGAACGCGCTCTCGCTGGGGTAAATTTCGGCCGTGGGTGCTGttt                                                <  1:782176/46‑1 (MQ=255)                              tCGAACACGCTCTCGCTGGGGTAAATTTCGGCCGTGGGTGCTGttt                                                <  1:1871797/46‑1 (MQ=255)                               cGAACGCGCTCTCGCTGGGGTAAATTTCGTCCGTGGGTACTGTTTGGCGCACTGCCATTCGGGATCGtg                        >  1:1384219/1‑68 (MQ=255)                               cGAACGCGCTCTCGCTGGGGTAAATTTCGTCCGTGGGTACTGTTTGGCGCACTGCCATTCGGGATCGTc                        >  1:1863231/1‑69 (MQ=255)                               cGAACGCGCTCTCGCTGGGGTAAATTTCGTCCGTGGGTACTGTTTGGCGCACTGCCATTCGGGATCGTc                        >  1:252201/1‑69 (MQ=255)                               cGAACGCGCTCTCGCTGGGGTAAATTTCGTCCGTGGGTACTGTTTGGCGCACTGCCATTCGGGATCGTc                        >  1:2554788/1‑69 (MQ=255)                                   cgcgCTCTCGCTGGGGTAAATTTCGTCCGTGGGTACTGTTTGGCGCACTGCCATTCGGGATCGTCtgtg                    <  1:3240022/69‑1 (MQ=255)                                   cgcgCTCTCGCTGGGGTAAATTTCGTCCGTGGGTACTGTTTGGCGCACTGCCATTCGGGATCGTCtgtg                    <  1:1047890/69‑1 (MQ=255)                                   cgcgCTCTCGCTGGGGTAAATTTCGTCCGTGGGTACTGTTTGGCGCACTGCCATTCGGGATCGTCtgtg                    <  1:2776921/69‑1 (MQ=255)                                   cgcgCTCTCGCTGGGGTAAATTTCGTCCGTGGGTACTGTTTGGCGCACTGCCATTCGGGATCGTCtgtg                    <  1:1432349/69‑1 (MQ=255)                                   cgcgCTCTCGCTGGGGTAAATTTCGTCCGTGGGTACTGTTTGGCGCACTGCCATTCGGGATCGTCtgtg                    <  1:1489618/69‑1 (MQ=255)                                            gCTGGGGTAAATTTCGTCCGTGGGTACTGTTTGGCGCACTGCCATTCTGGATCGTCTGTGTa                  <  1:1623383/62‑1 (MQ=255)                                            gCTGGGGTAAATTTCGTCCGTGGGTACTGTTTGGCGCACTGCCATTCGGGATCGTCTGTGTa                  <  1:2466288/62‑1 (MQ=255)                                            gCTGGGGTAAATTTCGTCCGTGGGTACTGTTTGGCGCACTGCCATTCGGGATCGTCTGTGTa                  <  1:2394496/62‑1 (MQ=255)                                            gCTGGGGTAAATTTCGTCCGTGGGTACTGTTTGGCGCACTGCCATTCGGGATCGTCTGTGTa                  <  1:2882624/62‑1 (MQ=255)                                            gCTGGGGTAAATTTCGTCCGTGGGTACTGTTTGGCGCACTGCCATTCGGGATCGTCTGTGTa                  <  1:1457860/62‑1 (MQ=255)                                            gCTGGGGTAAATTTCGTCCGTGGGTACTGTTTGGCGCACTGCCATTCGGGATCGTCTGTGTa                  <  1:139163/62‑1 (MQ=255)                                                     aaTTTCGTCCGTGGGTACTGTTTGGCGCACTGCCATTCGGGATCGTCTGTGTACTGGCCtat         >  1:1748819/1‑62 (MQ=255)                                                     aaTTTCGTCCGTGGGTACTGTTTGGCGCACTGCCATTCGGGATCGTCTGTGTACTGGCCTATAGCAc    >  1:1339423/1‑67 (MQ=255)                                                     aaTTTCGTCCGTGGGTACTGTTTGGCGCACTGCCATTCGGGATCGTCTGTGTACTGGCCTATAGCAc    >  1:2298224/1‑67 (MQ=255)                                                     aaTTTCGGCCGTGGGTGCTGTTTGGCGCACTGCCATTCGGGa                             >  1:293600/1‑42 (MQ=38)                                                     aaTTTCGGCCGTGGGTGCTGTTTGGCGCACTGCCATTCGGGATCGTCTGTGTACt                >  1:1548129/1‑55 (MQ=255)                                                     aaTTTCGGCCGTGGGTGCTGTTTGGCGCACTGCCATTCGGGATCGTCTGTGTACTGGCCTATAGCAc    >  1:1092328/1‑67 (MQ=255)                                                     aaTTTCGGCCGTGGGTGCTGTTTGGCGCACTGCCATTCGGGATCGTCTGTGTACTGGCCTATAGCAc    >  1:2796414/1‑67 (MQ=255)                                                     aaTTTCGGCCGTGGGTGCTGTTTGGCGCACTGCCATTCGGGATCGTCTGTGTACTGGCCTATAGCAc    >  1:749/1‑67 (MQ=255)                                                     aaTTTCGGCCGTGGGTGCTGTTTGGCGCACTGCCATTCGGGATCGTCTGTGTACTGGCCTATAGCAc    >  1:2890025/1‑67 (MQ=255)                                                     aaTTTCGGCCGTGGGTGCTGTTTGGCGCACTGCCATTCGGGATCGTCTGTGTACTGGCCTATAGCAc    >  1:645302/1‑67 (MQ=255)                                                     aaTTTCGGCCGTGGGTGCTGTTTGGCGCACTGCCATTCGGGATCGTCTGTGTACTGGCCTATAGCAc    >  1:1588477/1‑67 (MQ=255)                                                     aaTTTCGGCCGTGGGTGCTGTTTGGCGCACTGCCATTCGGGATCGTCTGTGTACTGGCCTATAGCAc    >  1:2663519/1‑67 (MQ=255)                                                     aaTTTCGGCCGTGGGTGCTGTTTGGCGCACTGCCATTCGGGATCGTCTGTGTACTGGCCTATAGCAc    >  1:2301332/1‑67 (MQ=255)                                                     aaTTTCGGCCGTGGGTGCTGTTTGGCGCACTGCCATTCGGGATCGTCTGTGTACTGGCCTATAGCAc    >  1:2319701/1‑67 (MQ=255)                                                     aaTTTCGGCCGTGGGTGCTGTTTGGCGCACTGCCATTCGGGATCGTCTGTGTACTGGCCTATAGCAc    >  1:3383660/1‑67 (MQ=255)                                                     aaTTTCGGCCGTGGGTGCTGTTTGGCGCACTGCCATTCGGGATCGTCTGTGTACTGGCCTATAGCAc    >  1:1387816/1‑67 (MQ=255)                                                     aaTTTCGGCCGTGGGTGCTGTTTGGCGCACTGCCATTCGGGATCGTCTGTGTACTGGCCTATAGCAc    >  1:1969684/1‑67 (MQ=255)                                                     aaTTTCGGCCGTGGGTGCTGTTTGGCGCACTGCCATTCGGGATCGTCTGTGTACTGGCCTATAGCAc    >  1:32089/1‑67 (MQ=255)                                                       tttCGTCCGTGGGTACTGTTTGGCGCACTGCCATTCGGGATCGTCTGTGTACTGGCCTATAGCACGc  <  1:804154/67‑1 (MQ=255)                                                       tttCGTCCGTGGGTACTGTTTGGCGCACTGCCATTCGGGATCGTCTGTGTACTGGCCTATAGCACGc  <  1:986886/67‑1 (MQ=255)                                                       tttCGTCCGTGGGTACTGTTTGGCGCACTGCCATTCGGGATCGTCTGTGTACTGGCCTATAGCACGc  <  1:781555/67‑1 (MQ=255)                                                       tttCGTCCGTGGGTACTGTTTGGCGCACTGCCATTCGGGATCGTCTGTGTACTGGCCTATAGCACGc  <  1:370518/67‑1 (MQ=255)                                                       tttCGTCCGTGGGTACTGTTTGGCGCACTGCCATTCGGGATCGTCTGTGTACTGGCCTATAGCACGc  <  1:367914/67‑1 (MQ=255)                                                       tttCGTCCGTGGGTACTGTTTGGCGCACTGCCATTCGGGATCGTCTGTGTACTGGCCTATAGCACGc  <  1:337526/67‑1 (MQ=255)                                                       tttCGTCCGTGGGTACTGTTTGGCGCACTGCCATTCGGGATCGTCTGTGTACTGGCCTATAGCACGc  <  1:1881210/67‑1 (MQ=255)                                                       tttCGTCCGTGGGTACTGTTTGGCGCACTGCCATTCGGGATCGTCTGTGTACTGGCCTATAGCACGc  <  1:3117307/67‑1 (MQ=255)                                                       tttCGTCCGTGGGTACTGTTTGGCGCACTGCCATTCGGGATCGTCTGTGTACTGGCCTATAGCACGc  <  1:2972322/67‑1 (MQ=255)                                                            |                                                              TCCGATCCTTGCATGGGGTTGTTGGCCGATCGAACGCGCTCTCGCTGGGGTAAATTTCGTCCGTGGGTACTGTTTGGCGCACTGCCATTCGGGATCGTCTGTGTACTGGCCTATAGCACGC  >  minE/2505725‑2505845

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