Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	1,665,156	C→T	35.6%	L22L (CTG→CTA)	guaA ←	GMP synthetase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	1,665,156	0	C	T	35.6%	58.0 / 48.0	84	L22L (CTG→CTA)	guaA	GMP synthetase
Reads supporting (aligned to +/- strand):  ref base C (29/25);  new base T (13/17);  total (42/42)
Fisher's exact test for biased strand distribution p-value = 4.95e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.88e-01

ACATCCCACGCCCACAGTTCGCAGTAAACACCCAGCTCACGCACGCGGCGCGCAACCAGTTGAGTGTACTGAGAACCGAAGTCCAGAATGAGGATGCGATGCTTATGAATGTTTTCCGT  >  minE/1665100‑1665218                                                         |                                                               acaTCCCACGCCCACAGTTCGCAGTAAACACCCAGCTCACGCACGCGGCGCGCAACCAGTTGAgtgt                                                      <  1:2018267/67‑1 (MQ=255) acaTCCCACGCCCACAGTTCGCAGTAAACACCCAGCTCACGCACGCGGCGCGCAACCAGTTGAgtgt                                                      <  1:2032011/67‑1 (MQ=255) acaTCCCACGCCCACAGTTCGCAGTAAACACCCAGCTCACGCACGCGGCGCGCAACCAGTTGAgtgt                                                      <  1:3237200/67‑1 (MQ=255) acaTCCCACGCCCACAGTTCGCAGTAAACACCCAGCTCACGCACGCGGCGCGCAACCAGTTGAgtgt                                                      <  1:2686900/67‑1 (MQ=255) acaTCCCACGCCCACAGTTCGCAGTAAACACCCAGCTCACGCACGCGGCGCGCAACCAGTTGAgtgt                                                      <  1:2683034/67‑1 (MQ=255) acaTCCCACGCCCACAGTTCGCAGTAAACACCCAGCTCACGCACGCGGCGCGCAACCAGTTGAgtgt                                                      <  1:2397642/67‑1 (MQ=255)   aTCCCACGCCCACAGTTCGCAGTAAACACCCAGCTCACGCACGCGGCGCGCAACCAGTTGAGTGTACt                                                   <  1:216836/68‑1 (MQ=255)                     gCAGTAAACACCCAGCTCACGCACGCGGCGCGCAACCAGTTGAGTGTACTGAGAACCGAAGTCCAGa                                  >  1:2503961/1‑67 (MQ=255)                     gCAGTAAACACCCAGCTCACGCACGCGGCGCGCAACCAGTTGAGTGTACTGAGAACCGAAGTCCAGa                                  >  1:76511/1‑67 (MQ=255)                     gCAGTAAACACCCAGCTCACGCACGCGGCGCGCAACCAGTTGAGTGTACTGAGAACCGAAGTCCAGa                                  >  1:2262748/1‑67 (MQ=255)                     gCAGTAAACACCCAGCTCACGCACGCGGCGCGCAACCAGTTGAGTGTACTGAGAACCGAAGTCCAGa                                  >  1:743808/1‑67 (MQ=255)                     gCAGTAAACACCCAGCTCACGCACGCGGCGCGCAACCAGTTGAGTGTACTGAGAACCGAAGTCCAGa                                  >  1:2397617/1‑67 (MQ=255)                     gCAGTAAACACCCAGCTCACGCACGCGGCGCGCAACCAGTTGAGTGTACTGAGAACCGAAGTCCAGa                                  >  1:366368/1‑67 (MQ=255)                     gCAGTAAACACCCAGCTCACGCACGCGGCGCGCAACCAGTTGAGTGTACTGAGAACCGAAGTCCAGa                                  >  1:105853/1‑67 (MQ=255)                     gCAGTAAACACCCAGCTCACGCACGCGGCGCGCAACCAGTTGAGTGTACTGAGAACCGAAGTCCAGa                                  >  1:1005779/1‑67 (MQ=255)                        gTAAACACCCAGCTCACGCACGCGGCGCGCAACCAGTTGAGTGTACTGAGAACCGAAGTCCAGAATGAg                             <  1:457374/69‑1 (MQ=255)                        gTAAACACCCAGCTCACGCACGCGGCGCGCAACCAGTTGAGTGTACTGAGAACCGAAGTCCAGAATGAg                             <  1:2721566/69‑1 (MQ=255)                        gTAAACACCCAGCTCACGCACGCGGCGCGCAACCAGTTGAGTGTACTGAGAACCGAAGTCCAGAATGAg                             <  1:410725/69‑1 (MQ=255)                           aaCACCCAGCTCACGCACGCGGCGCGCAACCAGTTGAGTGTACTGAGAACCGAAGTCCAg                                   <  1:1102684/60‑1 (MQ=255)                           aaCACCCAGCTCACGCACGCGGCGCGCAACCAGTTGAGTGTACTGAGAACCGAAGTCCAg                                   <  1:3414669/60‑1 (MQ=255)                           aaCACCCAGCTCACGCACGCGGCGCGCAACCAGTTGAGTGTACTGAGAACCGAAGTCCAg                                   <  1:1604918/60‑1 (MQ=255)                           aaCACCCAGCTCACGCACGCGGCGCGCAACCAGTTGAGTGTACTGAGAACCGAAGTCCAg                                   <  1:588244/60‑1 (MQ=255)                           aaCACCCAGCTCACGCACGCGGCGCGCAACCAGTTGAGTGTACTGAGAACCGAAGTCCAGAATGAgg                            >  1:3330019/1‑67 (MQ=255)                           aaCACCCAGCTCACGCACGCGGCGCGCAACCAGTTGAGTGTACTGAGAACCGAAGTCCAGAATGAGgat                          >  1:488133/1‑69 (MQ=255)                           aaCACCCAGCTCACGCACGCGGCGCGCAACCAGTTGAGTGTACTGAGAACCGAAGTCCAGAATGAGgat                          >  1:2981347/1‑69 (MQ=255)                           aaCACCCAGCTCACGCACGCGGCGCGCAACCAGTTGAGTGTACTGAGAACCGAAGTCCAGAATGAGgat                          >  1:538095/1‑69 (MQ=255)                           aaCACCCAGCTCACGCACGCGGCGCGCAACCAGTTGAGTGTACTGAGAACCGAAGTCCAGAATGAGgat                          >  1:584910/1‑69 (MQ=255)                           aaCACCCAGCTCACGCACGCGGCGCGCAACCAGTTGAGTGTACTGAGAACCGAAGTCCAGAATGAGgat                          >  1:3040966/1‑69 (MQ=255)                           aaCACCCAGCTCACGCACGCGGCGCGCAACCAGTTGAGTGTACTGAGAACCGAAGTCCAGAATGAGgat                          >  1:1343444/1‑69 (MQ=255)                           aaCACCCAGCTCACGCACGCGGCGCGCAACCAGTTGAGTGTACTGAGAACCGAAGTCCAGAATGAGgat                          >  1:1927012/1‑69 (MQ=255)                           aaCACCCAGCTCACGCACGCGGCGCGCAACCAGTTGAGTGTACTGAGAACCGAAGTCCAGAATGAGgat                          >  1:1377857/1‑69 (MQ=255)                           aaCACCCAGCTCACGCACGCGGCGCGCAACCAGTTGAGTGTACTGAGAACCGAAGTCCAGAATGAGgat                          >  1:2667623/1‑69 (MQ=255)                           aaCACCCAGCTCACGCACGCGGCGCGCAACCAGTTGAGTGTACTGAGAACCGAAGTCCAGAATGAGgat                          >  1:2532999/1‑69 (MQ=255)                           aaCACCCAGCTCACGCACGCGGCGCGCAACCAGTTGAGTGTACTGAGAACCGAAGTCCAGAATGAGgat                          >  1:532694/1‑69 (MQ=255)                           aaCACCCAGCTCACGCACGCGGCGCGCAACCAGTTGAGTGTACTGAGAACCGAAGTCCAGAATGAGgat                          >  1:2377242/1‑69 (MQ=255)                           aaCACCCAGCTCACGCACGCGGCGCGCAACCAGTTGAGTGTACTGAGAACCGAAGTCCAGAATGAGgat                          >  1:1706666/1‑69 (MQ=255)                           aaCACCCAGCTCACGCACGCGGCGCGCAACCAGTTGAGTGTACTGAGAACCGAAGTCCAGAATGAGgat                          >  1:3356027/1‑69 (MQ=255)                           aaCACCCAGCTCACGCACGCGGCGCGCAACCAGTTGAGTGTACTGAGAACCGAAGTCCAGAATGAGgat                          >  1:1950085/1‑69 (MQ=255)                           aaCACCCAGCTCACGCACGCGGCGCGCAACCAGTTGAGTGTACTGAGAACCGAAGTCCAGAATGAGga                           >  1:1583735/1‑68 (MQ=255)                           aaCACCCAGCTCACGCACGCGGCGCGCAACCAGTTGAGTGTACTGAGAACCGAAGTCCAGAATGAGga                           >  1:386522/1‑68 (MQ=255)                             cacCCAGCTCACGCACGCGGCGCGCAACCAGTTGAGTGTACTGAGAACCGAAGTCCAGAATGAGgatgc                        <  1:672468/69‑1 (MQ=255)                             cacCCAGCTCACGCACGCGGCGCGCAACCAGTTGAGTGTACTGAGAACCGAAGTCCAGAATGAGgatgc                        <  1:2792444/69‑1 (MQ=255)                                 cAGCTCACGCACGCGGCGCGCAACTAGTTGAGTGTACTGAGAACCGAAGTCCAGAATGAGgatgcgatg                    <  1:2690223/69‑1 (MQ=255)                                 cAGCTCACGCACGCGGCGCGCAACTAGTTGAGTGTACTGAGAACCGAAGTCCAGAATGAGgatgcgatg                    <  1:3143988/69‑1 (MQ=255)                                 cAGCTCACGCACGCGGCGCGCAACTAGTTGAGTGTACTGAGAACCGAAGTCCAGAATGAGgatgcgatg                    <  1:4765/69‑1 (MQ=255)                                 cAGCTCACGCACGCGGCGCGCAACTAGTTGAGTGTACTGAGAACCGAAGTCCAGAATGAGgatgcgatg                    <  1:3044237/69‑1 (MQ=255)                                 cAGCTCACGCACGCGGCGCGCAACTAGTTGAGTGTACTGAGAACCGAAGTCCAGAATGAGgatgcgatg                    <  1:3233527/69‑1 (MQ=255)                                 cAGCTCACGCACGCGGCGCGCAACTAGTTGAGTGTACTGAGAACCGAAGTCCAGAATGAGgatgcgatg                    <  1:833578/69‑1 (MQ=255)                                 cAGCTCACGCACGCGGCGCGCAACTAGTTGAGTGTACTGAGAACCGAAGTCCAGAATGAGgatgcgatg                    <  1:3037746/69‑1 (MQ=255)                                 cAGCTCACGCACGCGGCGCGCAACTAGTTGAGTGTACTGAGAACCGAAGTCCAGAATGAGgatgcgatg                    <  1:288452/69‑1 (MQ=255)                                 cAGCTCACGCACGCGGCGCGCAACTAGTTGAGTGTACTGAGAACCGAAGTCCAGAATGAGgatgcgatg                    <  1:2829381/69‑1 (MQ=255)                                 cAGCTCACGCACGCGGCGCGCAACTAGTTGAGTGTACTGAGAACCGAAGTCCAGAATGAGgatgcgatg                    <  1:1500708/69‑1 (MQ=255)                                 cAGCTCACGCACGCGGCGCGCAACTAGTTGAGTGTACTGAGAACCGAAGTCCAGAATGAGgatgcgatg                    <  1:1774044/69‑1 (MQ=255)                                 cAGCTCACGCACGCGGCGCGCAACTAGTTGAGTGTACTGAGAACCGAAGTCCAGAATGAGgatgcgatg                    <  1:787446/69‑1 (MQ=255)                                 cAGCTCACGCACGCGGCGCGCAACTAGTTGAGTGTACTGAGAACCGAAGTCCAGAATGAGgatgcgatg                    <  1:2043146/69‑1 (MQ=255)                                 cAGCTCACGCACGCGGCGCGCAACTAGTTGAGTGTACTGAGAACCGAAGTCCAGAATGAGgatgcgatg                    <  1:1806216/69‑1 (MQ=255)                                 cAGCTCACGCACGCGGCGCGCAACTAGTTGAGTGTACTGAGAACCGAAGTCCAGAATGAGgatgcgatg                    <  1:827316/69‑1 (MQ=255)                                 cAGCTCACGCACGCGGCGCGCAACTAGTTGAGTGAACTGAGAACCGAAGTCCAGAATGAGgatgcgatg                    <  1:1445739/69‑1 (MQ=255)                                 cAGCTCACGAACGCGGCGCGCAACTAGTTGAGTGTACTGAGAACCGAAGTCCAGAATGAGgatgcgatg                    <  1:1684348/69‑1 (MQ=255)                                   gCTCACGCACGCGGCGCGCAACCAGTTGAGTGTACTGAGAACCGAAGTCCAGAATGAGGATGCGATGCt                  <  1:714147/69‑1 (MQ=255)                                   gCTCACGCACGCGGCGCGCAACCAGTTGAGTGTACTGAGAACCGAAGTCCAGAATGAGGATGCGATGCt                  <  1:2423740/69‑1 (MQ=255)                                   gCTCACGCACGCGGCGCGCAACCAGTTGAGTGTACTGAGAACCGAAGTCCAGAATGAGGATGCGATGCt                  <  1:156378/69‑1 (MQ=255)                                   gCTCACGCACGCGGCGCGCAACCAGTTGAGTGTACTGAGAACCGAAGTCCAGAATGAGGATGCGATGCt                  <  1:2371283/69‑1 (MQ=255)                                   gCTCACGCACGCGGCGCGCAACCAGTTGAGTGTACTGAGAACCGAAGTCCAGAATGAGGATGCGATGCt                  <  1:1099465/69‑1 (MQ=255)                                      cacgcacgCGGCGCGCACCCAGTTGAGTGTACTGAGAACCGAAGTCCAGAATGAGGATGCGATGCTTa                >  1:99609/1‑68 (MQ=255)                                      cacgcacgCGGCGCGCAACCAGTTGAGTGTACTGAGAACCGAAGTCCAGAATGAGGATGCGATCCTTa                >  1:3124491/1‑68 (MQ=255)                                          cacgcGGCGCGCAACCAGTTGAGTCTACTGAGAACCGAAGTCCAGAATGAGgat                          >  1:316716/1‑54 (MQ=255)                                            cgcgGCGCGCAACCAGTTGAGTGTACTGAGAACCGAAGTCCAGAATGAGGATGCGATGCt                  <  1:2337314/60‑1 (MQ=255)                                               ggcgCGCAACTAGTTGAGTGTACTGAGAACCGAAGTCCAGAATGAGGATGCGATGCTTATGAATGtttt      >  1:677822/1‑69 (MQ=255)                                               ggcgCGCAACTAGTTGAGTGTACTGAGAACCGAAGTCCAGAATGAGGATGCGATGCTTATGAATGtttt      >  1:906186/1‑69 (MQ=255)                                               ggcgCGCAACTAGTTGAGTGTACTGAGAACCGAAGTCCAGAATGAGGATGCGATGCTTATGAATGtttt      >  1:1898058/1‑69 (MQ=255)                                               ggcgCGCAACTAGTTGAGTGTACTGAGAACCGAAGTCCAGAATGAGGATGCGATGCTTATGAATGtttt      >  1:2144317/1‑69 (MQ=255)                                               ggcgCGCAACTAGTTGAGTGTACTGAGAACCGAAGTCCAGAATGAGGATGCGATGCTTATGAATGtttt      >  1:2193786/1‑69 (MQ=255)                                               ggcgCGCAACTAGTTGAGTGTACTGAGAACCGAAGTCCAGAATGAGGATGCGATGCTTATGAATGtttt      >  1:2446642/1‑69 (MQ=255)                                               ggcgCGCAACTAGTTGAGTGTACTGAGAACCGAAGTCCAGAATGAGGATGCGATGCTTATGAATGtttt      >  1:694744/1‑69 (MQ=255)                                               ggcgCGCAACTAGTTGAGTGTACTGAGAACCGAAGTCCAGAATGAGGATGCGATGCTTATGAATGtttt      >  1:650100/1‑69 (MQ=255)                                               ggcgCGCAACTAGTTGAGTGTACTGAGAACCGAAGTCCAGAATGAGGATGCGATGCTTATGAATGtttt      >  1:2464617/1‑69 (MQ=255)                                               ggcgCGCAACTAGTTGAGTGTACTGAGAACCGAAGTCCAGAATGAGGATGCGATGCTTATGAATGtttt      >  1:2868168/1‑69 (MQ=255)                                               ggcgCGCAACTAGTTGAGTGTACTGAGAACCGAAGTCCAGAATGAGGATGCGATGCTTATGAATGtttt      >  1:541253/1‑69 (MQ=255)                                               ggcgCGCAACTAGTTGAGTGTACTGAGAACCGAAGTCCAGAATGAGGATGCGATGCTTATGAATGtttt      >  1:28847/1‑69 (MQ=255)                                               ggcgCGCAACTAGTTGAGTGTACTGAGAACCGAAGTCCAGAATGAGGATGCGATGCTTATGAATGtttt      >  1:3056186/1‑69 (MQ=255)                                                   cgcAACCAGTTGAGTGTACTGAGAACCGAAGTCCAGAATGAGGATGCGATGCTTATGAATGTTTTCCGt  <  1:2452413/69‑1 (MQ=255)                                                   cgcAACCAGTTGAGTGTACTGAGAACCGAAGTCCAGAATGAGGATGCGATGCTTATGAATGTTTTCCGt  <  1:1270946/69‑1 (MQ=255)                                                   cgcAACCAGTTGAGTGTACTGAGAACCGAAGTCCAGAATGAGGATGCGATGCTTATGAATGTTTTCCGt  <  1:1053844/69‑1 (MQ=255)                                                         |                                                               ACATCCCACGCCCACAGTTCGCAGTAAACACCCAGCTCACGCACGCGGCGCGCAACCAGTTGAGTGTACTGAGAACCGAAGTCCAGAATGAGGATGCGATGCTTATGAATGTTTTCCGT  >  minE/1665100‑1665218

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