Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	165,472	A→G	43.4%	R668R (CGT→CGC)	glnD ←	uridylyltransferase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	165,472	0	A	G	43.4%	19.8 / 64.6	83	R668R (CGT→CGC)	glnD	uridylyltransferase
Reads supporting (aligned to +/- strand):  ref base A (14/33);  new base G (24/12);  total (38/45)
Fisher's exact test for biased strand distribution p-value = 9.90e-04
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.39e-01

AATGGCGGGCATGCCAGGCCAGTTGATTTGGGCTATGGCGGACAAAATAGTTAGCACGACAGCGTGACCAAATTTGGTGCAGCGCCTCTTCGTCGATGTTATCCATGCGCAGTAGTGCC  >  minE/165417‑165535                                                        |                                                                aaTGGCGGGCATGCCAGGCCAGTTGATTTGGGCTATGGCGGACAAAATAGTTAGCACGACAGCGTGAc                                                     >  1:2810924/1‑68 (MQ=255) aaTGGCGGGCATGCCAGGCCAGTTGATTTGGGCTATGGCGGACAAAATAGTTAGCACGACAGCGTGAc                                                     >  1:1515772/1‑68 (MQ=255) aaTGGCGGGCATGCCAGGCCAGTTGATTTGGGCTATGGCGGACAAAATAGTTAGCACGACAGCGTGAc                                                     >  1:3183291/1‑68 (MQ=255) aaTGGCGGGCATGCCAGGCCAGTTGATTTGGGCTATGGCGGACAAAATAGTTAGCACGACAGCGTGAc                                                     >  1:1800546/1‑68 (MQ=255)          cATGCCAGGCCAGTTGATTTGGGCTATGGCGGACAAAATAGTTAGCACGACAGCGTGACCAAATTTgg                                            >  1:363803/1‑68 (MQ=255)          cATGCCAGGCCAGTTGATTTGGGCTATGGCGGACAAAATAGTTAGCACGACAGCGTGACCAAATTTgg                                            >  1:836760/1‑68 (MQ=255)          cATGCCAGGCCAGTTGATTTGGGCTATGGCGGACAAAATAGTTAGCACGACAGCGTGACCAAATTTgg                                            >  1:2592802/1‑68 (MQ=255)          cATGCCAGGCCAGTTGATTTGGGCTATGGCGGACAAAATAGTTAGCACGACAGCGTGACCAAATTTgg                                            >  1:3093036/1‑68 (MQ=255)          cATGCCAGGCCAGTTGATTTGGGCTATGGCGGACAAAATAGTTAGCACGACAGCGTGACCAAATTTgg                                            >  1:6781/1‑68 (MQ=255)          cATGCCAGGCCAGTTGATTTGGGCTATGGCGGACAAAATAGTTAGCACGACAGCGTGACCAAATTTgg                                            >  1:3396130/1‑68 (MQ=255)          cATGCCAGGCCAGTTGATTTGGGCTATGGCGGACAAAATAGTTAGCACGACAGCGTGACCAAATTTgg                                            >  1:1510822/1‑68 (MQ=255)              ccaggccagTTGATTTGGGCTATGGCGGACAAAATAGTTAGCACGACAGCGTGACCAAATTTGGTGCAg                                       >  1:1439493/1‑69 (MQ=255)                                gCTATGGCGGACAAAATAGTTAGCACGACAGCGTGACCAAATTTGGTGCAGCGCCTCTTCGTCGATGt                      <  1:1773663/68‑1 (MQ=255)                                 cTATGGCGGACAAAATAGTTAGCACGACAGCGTGACCAAATTTGGTGCAGCGCCTCTTCGTCGATGt                      <  1:1612571/67‑1 (MQ=255)                                 cTATGGCGGACAAAATAGTTAGCACGACAGCGTGACCAAATTTGGTGCAGCGCCTCTTCGTCGATGTTa                    <  1:20441/69‑1 (MQ=255)                                   aTGGCGGACAAAATAGTTAGCACGACAGCGTGACCAAATTTGGTGCAGCGCCTCTTCGTCGATGTTATc                  <  1:2556386/69‑1 (MQ=255)                                   aTGGCGGACAAAATAGTTAGCACGACAGCGTGACCAAATTTGGTGCAGCGCCTCTTCGTCGATGTTATc                  <  1:1077938/69‑1 (MQ=255)                                   aTGGCGGACAAAATAGTTAGCACGACAGCGTGACCAAATTTGGTGCAGCGCCTCTTCGTCGATGTTATc                  <  1:76125/69‑1 (MQ=255)                                     ggCGGACAAAATAGTTAGCACGACAGCGTGACCAAATTTGGTGCAGCGCCTCTTCGTCGATGTTATc                  <  1:2309861/67‑1 (MQ=255)                                            aaaaTAGTTAGCACGACAGCGTGACCAAATTTGGTGCAGCGCCTCTTCGTCGATGTTATCCATGCGCag         <  1:2336030/69‑1 (MQ=255)                                            aaaaTAGTTAGCACGACAGCGTGACCAAATTTGGTGCAGCGCCTCTTCGTCGATGTTATCCATGCGCag         <  1:2294987/69‑1 (MQ=255)                                            aaaaTAGTTAGCACGACAGCGTGACCAAATTTGGTGCAGCGCCTCTTCGTCGATGTTATCCATGCGCag         <  1:2693984/69‑1 (MQ=255)                                            aaaaTAGTTAGCACGACAGCGTGACCAAATTTGGTGCAGCGCCTCTTCGTCGATGTTATCCATGCGCag         <  1:2953892/69‑1 (MQ=255)                                            aaaaTAGTTAGCACGACAGCGTGACCAAATTTGGTGCAGCGCCTCTTCGTCGATGTTATCCATGCGCag         <  1:2049398/69‑1 (MQ=255)                                            aaaaTAGTTAGCACGACAGCGTGACCAAATTTGGTGCAGCGCCTCTTCGTCGATGTTATCCATGCGCag         <  1:1881957/69‑1 (MQ=255)                                            aaaaTAGTTAGCACGACAGCGTGACCAAATTTGGTGCAGCGCCTCTTCGTCGATGTTATCCATGCGCag         <  1:71362/69‑1 (MQ=255)                                            aaaaTAGTTAGCACGACAGCGTGACCAAATTTGGTGCAGCGCCTCTTCGTCGATGTTATCCATGCGCag         <  1:185670/69‑1 (MQ=255)                                            aaaaTAGTTAGCACGACAGCGTGACCAAATTTGGTGCAGCGCCTCTTCGTCGATGTTATCCATGCGCag         <  1:1805756/69‑1 (MQ=255)                                            aaaaTAGTTAGCACGACAGCGTGACCAAATTTGGTGCAGCGCCTCTTCGTCGATGTTATCCATGCGCag         <  1:1534838/69‑1 (MQ=255)                                             aaaTAGTTAGCACGACAGCGTGACCAAATTTGGTGCAGCGCCTCTTCGTCGATGTTATCCATGCGCagt        <  1:1936526/69‑1 (MQ=255)                                             aaaTAGTTAGCACGACAGCGTGACCAAATTTGGTGCAGCGCCTCTTCGTCGATGTTATCCATGCGCagt        <  1:2646222/69‑1 (MQ=255)                                             aaaTAGTTAGCACGACAGCGTGACCAAATTTGGTGCAGCGCCTCTTCGTCGATGTTATCCATGCGCagt        <  1:2741850/69‑1 (MQ=255)                                             aaaTAGTTAGCACGACAGCGTGACCAAATTTGGTGCAGCGCCTCTTCGTCGATGTTATCCATGCGCagt        <  1:572497/69‑1 (MQ=255)                                               aTAGTTAGCGCGACAGCGTGACCAAATTTGGTGCAGCGCCTCTTCGTCGATGTTATCCATGCGCagtag      <  1:1079614/69‑1 (MQ=255)                                               aTAGTTAGCGCGACAGCGTGACCAAATTTGGTGCAGCGCCTCTTCGTCGATGTTATCCATGCGCagtag      <  1:829557/69‑1 (MQ=255)                                               aTAGTTAGCGCGACAGCGTGACCAAATTTGGTGCAGCGCCTCTTCGTCGATGTTATCCATGCGCagtag      <  1:827608/69‑1 (MQ=255)                                               aTAGTTAGCGCGACAGCGTGACCAAATTTGGTGCAGCGCCTCTTCGTCGATGTTATCCATGCGCagtag      <  1:746184/69‑1 (MQ=255)                                               aTAGTTAGCGCGACAGCGTGACCAAATTTGGTGCAGCGCCTCTTCGTCGATGTTATCCATGCGCagtag      <  1:62055/69‑1 (MQ=255)                                               aTAGTTAGCGCGACAGCGTGACCAAATTTGGTGCAGCGCCTCTTCGTCGATGTTATCCATGCGCagtag      <  1:350183/69‑1 (MQ=255)                                               aTAGTTAGCGCGACAGCGTGACCAAATTTGGTGCAGCGCCTCTTCGTCGATGTTATCCATGCGCagtag      <  1:365779/69‑1 (MQ=255)                                               aTAGTTAGCGCGACAGCGTGACCAAATTTGGTGCAGCGCCTCTTCGTCGATGTTATCCATGCGCagtag      <  1:1486359/69‑1 (MQ=255)                                               aTAGTTAGCGCGACAGCGTGACCAAATTTGGTGCAGCGCCTCTTCGTCGATGTTATCCATGCGCagtag      <  1:1427333/69‑1 (MQ=255)                                               aTAGTTAGCGCGACAGCGTGACCAAATTTGGTGCAGCGCCTCTTCGTCGATGTTATCCATGCGCagtag      <  1:1114898/69‑1 (MQ=255)                                               aTAGTTAGCGCGACAGCGTGACCAAATTTGGTGCAGCGCCTCTTCGTCGATGTTATCCATGCGCagtag      <  1:1236788/69‑1 (MQ=255)                                               aTAGTTAGCGCGACAGCGTGACCAAATTTGGTGCAGCGCCTCTTCGTCGATGTTATCCATGCGCagtag      <  1:1878434/69‑1 (MQ=255)                                               aTAGTTAGCACGACAGCGTGACCAAATTTGGTGCAGCGCCTCTTCGTCGATGTTATCCATGCGCagtag      <  1:2710721/69‑1 (MQ=255)                                               aTAGTTAGCACGACAGCGTGACCAAATTTGGTGCAGCGCCTCTTCGTCGATGTTATCCATGCGCagtag      <  1:869321/69‑1 (MQ=255)                                                tAGTTAGCACGACAGAGTGACCAAATTTGGTGCAGCGCCTCTTCGTCGATGTTATCCATGCGCagtag      <  1:3259098/68‑1 (MQ=255)                                                  gTTAGCGCGACAGCGTGACCAAATTTGGTGCAGCGCCTCTtcg                             >  1:1080238/1‑43 (MQ=255)                                                  gTTAGCGCGACAGCGTGACCAAATTTGGTGCAGCGCCTCTTCGTCGATGTTAt                   >  1:929524/1‑53 (MQ=255)                                                  gTTAGCGCGACAGCGTGACCAAATTTGGTGCAGCGCCTCTTCGTCGATGTTATCCATGCGCagtagt     >  1:780655/1‑67 (MQ=255)                                                  gTTAGCGCGACAGCGTGACCAAATTTGGTGCAGCGCCTCTTCGTCGATGTTATCCATGCGCagtagt     >  1:2485883/1‑67 (MQ=255)                                                  gTTAGCGCGACAGCGTGACCAAATTTGGTGCAGCGCCTCTTCGTCGATGTTATCCATGCGCagtagt     >  1:1327772/1‑67 (MQ=255)                                                  gTTAGCGCGACAGCGTGACCAAATTTGGTGCAGCGCCTCTTCGTCGATGTTATCCATGCGCagtagt     >  1:701768/1‑67 (MQ=255)                                                  gTTAGCGCGACAGCGTGACCAAATTTGGTGCAGCGCCTCTTCGTCGATGTTATCCATGCGCagtagt     >  1:1294509/1‑67 (MQ=255)                                                  gTTAGCGCGACAGCGTGACCAAATTTGGTGCAGCGCCTCTTCGTCGATGTTATCCATGCGCagtagt     >  1:1693460/1‑67 (MQ=255)                                                  gTTAGCGCGACAGCGTGACCAAATTTGGTGCAGCGCCTCTTCGTCGATGTTATCCATGCGCagtagt     >  1:862215/1‑67 (MQ=255)                                                  gTTAGCGCGACAGCGTGACCAAATTTGGTGCAGCGCCTCTTCGTCGATGTTATCCATGCGCagtagt     >  1:867728/1‑67 (MQ=255)                                                  gTTAGCGCGACAGCGTGACCAAATTTGGTGCAGCGCCTCTTCGTCGATGTTATCCATGCGCagtagt     >  1:923842/1‑67 (MQ=255)                                                  gTTAGCGCGACAGCGTGACCAAATTTGGTGCAGCGCCTCTTCGTCGATGTTATCCATGCGCagtagt     >  1:948694/1‑67 (MQ=255)                                                  gTTAGCGCGACAGCGTGACCAAATTTGGTGCAGCGCCTCTTCGTCGATGTTATCCATGCGCagtagt     >  1:1801248/1‑67 (MQ=255)                                                  gTTAGCGCGACAGCGTGACCAAATTTGGTGCAGCGCCTCTTCGTCGATGTTATCCATGCGCagtagt     >  1:2127278/1‑67 (MQ=255)                                                  gTTAGCGCGACAGCGTGACCAAATTTGGTGCAGCGCCTCTTCGTCGATGTTATCCATGCGCagtagt     >  1:2100972/1‑67 (MQ=255)                                                  gTTAGCGCGACAGCGTGACCAAATTTGGTGCAGCGCCTCTTCGTCGATGTTATCCATGCGCagtagt     >  1:2822076/1‑67 (MQ=255)                                                  gTTAGCGCGACAGCGTGACCAAATTTGGTGCAGCGCCTCTTCGTCGATGTTATCCATGCGCagtagt     >  1:2939184/1‑67 (MQ=255)                                                  gTTAGCGCGACAGCGTGACCAAATTTGGTGCAGCGCCTCTTCGTCGATGTTATCCATGCGCagtagt     >  1:2098346/1‑67 (MQ=255)                                                  gTTAGCGCGACAGCGTGACCAAATTTGGTGCAGCGCCTCTTCGTCGATGTTATCCATGCGCagtagt     >  1:2980962/1‑67 (MQ=255)                                                  gTTAGCGCGACAGCGTGACCAAATTTGGTGCAGCGCCTCTTCGTCGATGTTATCCATGCGCagtagt     >  1:1906289/1‑67 (MQ=255)                                                  gTTAGCGCGACAGCGTGACCAAATTTGGTGCAGCGCCTCTTCGTCGATGTTATCCATGCGCagtagt     >  1:3370684/1‑67 (MQ=255)                                                  gTTAGCGCGACAGCGTGACCAAATTTGGTGCAGCGCCTCTTCGTCGATGTTATCCATGCGCagtagt     >  1:1809545/1‑67 (MQ=255)                                                  gTTAGCGCGACAGCGTGACCAAATTTGGTGCAGCGCCTCTTCGTCGATGTTATCCATGCGCagtagt     >  1:358970/1‑67 (MQ=255)                                                  gTTAGCGCGACAGCGTGACCAAATTTGGTGCAGCGCCTCTTCGTCGATGTTATCCATGCGCagtag      >  1:651595/1‑66 (MQ=255)                                                  gTTAGCACGACAGCGTGACCAAATTTGGTGCAGCGCCTCTTCGTCGATGTTATCCATGCGCagtagt     >  1:1195924/1‑67 (MQ=255)                                                  gTTAGCACGACAGCGTGACCAAATTTGGTGCAGCGCCTCTTCGTCGATGTTATCCATGCGCagtagt     >  1:3431557/1‑67 (MQ=255)                                                  gTTAGCACGACAGCGTGACCAAATTTGGTGCAGCGCCTCTTCGTCGATGTTATCCATGCGCAGTagtgc   <  1:3390502/69‑1 (MQ=255)                                                    tAGCACGACAGCGTGACCAAATTTGGTGCAGCGCCTCTTCGTCGATGTTATCCATGCGCAGTagtgcc  <  1:793475/68‑1 (MQ=255)                                                    tAGCACGACAGCGTGACCAAATTTGGTGCAGCGCCTCTTCGTCGATGTTATCCATGCGCAGTagtgcc  <  1:2945895/68‑1 (MQ=255)                                                    tAGCACGACAGCGTGACCAAATTTGGTGCAGCGCCTCTTCGTCGATGTTATCCATGCGCAGTagtgcc  <  1:2945413/68‑1 (MQ=255)                                                    tAGCACGACAGCGTGACCAAATTTGGTGCAGCGCCTCTTCGTCGATGTTATCCATGCGCAGTagtgcc  <  1:2264732/68‑1 (MQ=255)                                                    tAGCACGACAGCGTGACCAAATTTGGTGCAGCGCCTCTTCGTCGATGTTATCCATGCGCAGTagtgcc  <  1:2342205/68‑1 (MQ=255)                                                    tAGCACGACAGCGTGACCAAATTTGGTGCAGCGCCTCTTCGTCGATGTTATCCATGCGCAGTagtccc  <  1:1464825/68‑1 (MQ=255)                                                     aGCACGACAGCGTGACCAAATTTGGTGCAGCGCCTCTTCGTCGATGTTATCCATGCGCAGTagtgcc  <  1:1314507/67‑1 (MQ=255)                                                     aGCACGACAGCGTGACCAAATTTGGTGCAGCGCCTCTTCGTCGATGTTATCCATGCGCAGTagtgcc  <  1:553679/67‑1 (MQ=255)                                                        |                                                                AATGGCGGGCATGCCAGGCCAGTTGATTTGGGCTATGGCGGACAAAATAGTTAGCACGACAGCGTGACCAAATTTGGTGCAGCGCCTCTTCGTCGATGTTATCCATGCGCAGTAGTGCC  >  minE/165417‑165535

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑