Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	435,077	C→T	28.7%	Q41Q (CAG→CAA)	ybfE ←	LexA regulated protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	435,077	0	C	T	28.7%	102.3 / 40.6	97	Q41Q (CAG→CAA)	ybfE	LexA regulated protein
Reads supporting (aligned to +/- strand):  ref base C (48/21);  new base T (7/21);  total (55/42)
Fisher's exact test for biased strand distribution p-value = 8.92e-05
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.97e-01

GCCTCGACAGCTTCCGCGTTCAGCTTCAGTTCGACACGCTTAAGGCCACGTACTTTGTCGCGTTTTAGCTGGTTGCGTTTATTAATACGCAGCTGTTCATCGCGCGAAAGCGGATTAGTTTTCGGTCG  >  minE/435009‑435136                                                                     |                                                            gcCTCGACAGCTTCCGCGTTCAGCTTCAGTTCGACACGCTTAAGGCCACGTACTTTGTCGCGTTTTAGc                                                             >  1:1975247/1‑69 (MQ=255)          gCTTCCGCGTTCAGCTTCAGTTCGACACGCTTAAGGCCACGTACTTTGTCGCGTTTTAGCTGGTTGCg                                                     >  1:3206089/1‑68 (MQ=255)          gCTTCCGCGTTCAGCTTCAGTTCGACACGCTTAAGGCCACGTACTTTGTCGCGTTTTAGCTGGTTGCg                                                     >  1:1686821/1‑68 (MQ=255)          gCTTCCGCGTTCAGCTTCAGTTCGACACGCTTAAGGCCACGTACTTTGTCGCGTTTTAGCTGGTTGCg                                                     >  1:1652708/1‑68 (MQ=255)          gCTTCCGCGTTCAGCTTCAGTTCGACACGCTTAAGGCCACGTACTTTGTCGCGTTTTAGCTGGTTGCg                                                     >  1:1639631/1‑68 (MQ=255)          gCTTCCGCGTTCAGCTTCAGTTCGACACGCTTAAGGCCACGTACTTTGTCGCGTTTTAGCTGGTTGCg                                                     >  1:1783145/1‑68 (MQ=255)          gCTTCCGCGTTCAGCTTCAGTTCGACACGCTTAAGGCCACGTACTTTGTCGCGTTTTAGCTGGTTGCg                                                     >  1:1812997/1‑68 (MQ=255)          gCTTCCGCGTTCAGCTTCAGTTCGACACGCTTAAGGCCACGTACTTTGTCGCGTTTTAGCTGGTTGCg                                                     >  1:3159563/1‑68 (MQ=255)          gCTTCCGCGTTCAGCTTCAGTTCGACACGCTTAAGGCCACGTACTTTGTCGCGTTTTAGCTGGTTGCg                                                     >  1:2729601/1‑68 (MQ=255)          gCTTCCGCGTTCAGCTTCAGTTCGACACGCTTAAGGCCACGTACTTTGTCGCGTTTTAGCTGGTTGCg                                                     >  1:2845137/1‑68 (MQ=255)          gCTTCCGCGTTCAGCTTCAGTTCGACACGCTTAAGGCCACGTACTTTGTCGCGTTTTAGCTGGTTGCg                                                     >  1:266814/1‑68 (MQ=255)          gCTTCCGCGTTCAGCTTCAGTTCGACACGCTTAAGGCCACGTACTTTGTCGCGTTTTAGCTGGTTGCg                                                     >  1:2277608/1‑68 (MQ=255)          gCTTCCGCGTTCAGCTTCAGTTCGACACGCTTAAGGCCACGTACTTTGTCGCGTTTTAGCTGGTTGCg                                                     >  1:501504/1‑68 (MQ=255)          gCTTCCGCGTTCAGCTTCAGTTCGACACGCTTAAGGCCACGTACTTTGTCGCGTTTTAGCTGGTTGCg                                                     >  1:2075669/1‑68 (MQ=255)           cTTCCGCGTTCAGCTTCAGTTCGACACGCTTAAGGCCACGTACTTTGTCGCGTTTTAGCTGGTTGCGtt                                                   <  1:1419211/69‑1 (MQ=255)           cTTCCGCGTTCAGCTTCAGTTCGACACGCTTAAGGCCACGTACTTTGTCGCGTTTTAGCTGGTTGCGtt                                                   <  1:3021713/69‑1 (MQ=255)           cTTCCGCGTTCAGCTTCAGTTCGACACGCTTAAGGCCACGTACTTTGTCGCGTTTTAGCTGGTTGCGtt                                                   <  1:539191/69‑1 (MQ=255)           cTTCCGCGTTCAGCTTCAGTTCGACACGCTTAAGGCCACGTACTTTGTCGCGTTTTAGCTGGTTGCGtt                                                   <  1:1213197/69‑1 (MQ=255)             tCCGCGTTCAGCTTCAGTTCGACACGCTTAAGGCCACGTACTTTGTCGCGTTTTAGCTGGTTGCGttt                                                  >  1:2967914/1‑68 (MQ=255)             tCCGCGTTCAGCTTCAGTTCGACACGCTTAAGGCCACGTACTTTGTCGCGTTTTAGCTGGTTGCGttt                                                  >  1:3018616/1‑68 (MQ=255)             tCCGCGTTCAGCTTCAGTTCGACACGCTTAAGGCCACGTACTTTGTCGCGTTTTAGCTGGTTGCGttt                                                  >  1:2246183/1‑68 (MQ=255)             tCCGCGTTCAGCTTCAGTTCGACACGCTTAAGGCCACGTACTTTGTCGCGTTTTAGCTGGTTGCGttt                                                  >  1:2960844/1‑68 (MQ=255)             tCCGCGTTCAGCTTCAGTTCGACACGCTTAAGGCCACGTACTTTGTCGCGTTTTAGCTGGTTGCGttt                                                  >  1:641455/1‑68 (MQ=255)             tCCGCGTTCAGCTTCAGTTCGACACGCTTAAGGCCACGTACTTTGTCGCGTTTTAGCTGGTTGCGttt                                                  >  1:2747386/1‑68 (MQ=255)             tCCGCGTTCAGCTTCAGTTCGACACGCTTAAGGCCACGTACTTTGTCGCGTTTTAGCTGGTTGCGttt                                                  >  1:2751414/1‑68 (MQ=255)             tCCGCGTTCAGCTTCAGTTCGACACGCTTAAGGCCACGTACTTTGTCGCGTTTTAGCTGGTTGCGttt                                                  >  1:1023905/1‑68 (MQ=255)             tCCGCGTTCAGCTTCAGTTCGACACGCTTAAGGCCACGTACTTTGTCGCGTTTTAGCTGGTTGCGttt                                                  >  1:1020077/1‑68 (MQ=255)             tCCGCGTTCAGCTTCAGTTCGACACGCTTAAGGCCACGTACTTTGTCGCGTTTTAGCTGGTTGCGttt                                                  >  1:1712695/1‑68 (MQ=255)             tCCGCGTTCAGCTTCAGTTCGACACGCTTAAGGCCACGTACTTTGTCGCGTTTTAGCTGGTTGCGttt                                                  >  1:896960/1‑68 (MQ=255)             tCCGCGTTCAGCTTCAGTTCGACACGCTTAACGCCACGTACTTTGTCGCGTTTTAGCTGGTTGCGttt                                                  >  1:3386079/1‑68 (MQ=255)                gcgTTCAGCTTCAGTTCGACACGCTTTAGGCCACGTACTTTGTCGCGTTTTAGCTGGTTGCGtttttt                                               <  1:1904216/68‑1 (MQ=255)                gcgTTCAGCTTCAGTTCGACACGCTTAAGGCCACGTACTTTGTCGCGTTTTAGCTGGTTGCGTttatt                                               <  1:2225860/68‑1 (MQ=255)                gcgTTCAGCTTCAGTTCGACACGCTTAAGGCCACGTACTTTGTCGCGTTTTAGCTGGTTGCGTttatt                                               <  1:2797047/68‑1 (MQ=255)                gcgTTCAGCTTCAGTTCGACACGCTTAAGGCCACGTACTTTGTCGCGTTTTAGCTGGTTGCGTttatt                                               <  1:2812381/68‑1 (MQ=255)                gcgTTCAGCTTCAGTTCGACACGCTTAAGGCCACGTACTTTGTCGCGTTTTAGCTGGTTGCGTttatt                                               <  1:595753/68‑1 (MQ=255)                gcgTTCAGCTTCAGTTCGACACGCTTAAGGCCACGTACTTTGTCGCGTTTTAGCTGGTTGCGTttatt                                               <  1:2927948/68‑1 (MQ=255)                gcgTTCAGCTTCAGTTCGACACGCTTAAGGCCACGTACTTTGTCGCGTTTTAGCTGGTTGCGTttatt                                               <  1:573861/68‑1 (MQ=255)                gcgTTCAGCTTCAGTTCGACACGCTTAAGGCCACGTACTTTGTCGCGTTTTAGCTGGTTGCGTttatt                                               <  1:1342289/68‑1 (MQ=255)                gcgTTCAGCTTCAGTTCGACACGCTTAAGGCCACGTACTTTGTCGCGTTTTAGCTGGTTGCGTttatt                                               <  1:1360464/68‑1 (MQ=255)                gcgTTCAGCTTCAGTTCGACACGCTTAAGGCCACGTACTTTGTCGCGTTTTAGCTGGTTGCGTttatt                                               <  1:219784/68‑1 (MQ=255)                gcgTTCAGCTTCAGTTCGACACGCTTAAGGCCACGTACTTTGTCGCGTTTTAGCTGGTTGCGTttatt                                               <  1:141948/68‑1 (MQ=255)                gcgTTCAGCTTCAGTTCGACACGCTTAAGGCCACGTACTTTGTCGCGTTTTAGCTGGTTGCGTttatt                                               <  1:2145150/68‑1 (MQ=255)                gcgTTCAGCTTCAGTTCGACACGCTTAAGGCCACGTACTTTGTCGCGTTTTAGCTGGTTGCGTttatt                                               <  1:1484565/68‑1 (MQ=255)                gcgTTCAGCTTCAGTTCGACACGCTTAAGGCCACGTACTTTGTCGCGTTTTAGCTGGTTGCGTttatt                                               <  1:1494184/68‑1 (MQ=255)                gcgTTCAGCTTCAGTTCGACACGCTTAAGGCCACGTACTTTGTCGCGTTTTAGCTGGTTGCGTttatt                                               <  1:1618677/68‑1 (MQ=255)                     cAGCTTCAGTTCGACACGCTTAAGGCCACGTACTTTGTCGCGTTTTAGCTGGTTGCGTTTATTAATa                                           >  1:2644695/1‑67 (MQ=255)                     cAGCTTCAGTTCGACACGCTTAAGGCCACGTACTTTGTCGCGTTTTAGCTGGTTGCGTTTATTAATa                                           >  1:761535/1‑67 (MQ=255)                     cAGCTTCAGTTCGACACGCTTAAGGCCACGTACTTTGTCGCGTTTTAGCTGGTTGCGTTTATTAATa                                           >  1:391601/1‑67 (MQ=255)                     cAGCTTCAGTTCGACACGCTTAAGGCCACGTACTTTGTCGCGTTTTAGCTGGTTGCGTTTATTAATa                                           >  1:2606082/1‑67 (MQ=255)                     cAGCTTCAGTTCGACACGCTTAAGGCCACGTACTTTGTCGCGTTTTAGCTGGTTGCGTTTATTAATa                                           >  1:2939875/1‑67 (MQ=255)                     cAGCTTCAGTTCGACACGCTTAAGGCCACGTACTTTGTCGCGTTTTAGCTGGTTGCGTTTATTAATa                                           >  1:1399054/1‑67 (MQ=255)                     cAGCTTCAGTTCGACACGCTTAAGGCCACGTACTTTGTCGCGTTTTAGCTGGTTGCGTTTATTAATa                                           >  1:424348/1‑67 (MQ=255)                                cGACACGCTTAAGGCCACGTACTTTGTCGCGTTTTAGCTGGTTGCGTTTATTAATACGc                                        <  1:1608973/59‑1 (MQ=255)                                     cGCTTAAGGCCACGTACTTTGTCGCGTTTTAGCTGGTTGCGTTTATTAATACGCAGCTGTTCATcgcg                          >  1:389273/1‑68 (MQ=255)                                     cGCTTAAGGCCACGTACTTTGTCGCGTTTTAGCTGGTTGCGTTTATTAATACGCAGCTGTTCATcgcg                          >  1:2719175/1‑68 (MQ=255)                                     cGCTTAAGGCCACGTACTTTGTCGCGTTTTAGCTGGTTGCGTTTATTAATACGCAGCTGTTCATcgcg                          >  1:3286782/1‑68 (MQ=255)                                             gCCACGTACTTTGTCGCGTTTTAGCTGGTTGCGTTTATTAATACGCAGCTGTTCATcgc                           >  1:1856115/1‑59 (MQ=255)                                             gCCACGTACTTTGTCGCGTTTTAGCTGGTTGCGTTTATTAATACGCAGCTGTTCATCGCGCGAAAGCgg                 >  1:397077/1‑69 (MQ=255)                                             gCCACGTACTTTGTCGCGTTTTAGCTGGTTGCGTTTATTAATACGCAGCTGTTCATCGCGCGAAAGCgg                 >  1:513178/1‑69 (MQ=255)                                             gCCACGTACTTTGTCGCGTTTTAGCTGGTTGCGTTTATTAATACGCAGCTGTTCATCGCGCGAAAGCgg                 >  1:631374/1‑69 (MQ=255)                                             gCCACGTACTTTGTCGCGTTTTAGCTGGTTGCGTTTATTAATACGCAGCTGTTCATCGCGCGAAAGCgg                 >  1:641471/1‑69 (MQ=255)                                             gCCACGTACTTTGTCGCGTTTTAGCTGGTTGCGTTTATTAATACGCAGCTGTTCATCGCGCGAAAGCgg                 >  1:2521940/1‑69 (MQ=255)                                             gCCACGTACTTTGTCGCGTTTTAGCTGGTTGCGTTTATTAATACGCAGCTGTTCATCGCGCGAAAGCgg                 >  1:2576798/1‑69 (MQ=255)                                             gCCACGTACTTTGTCGCGTTTTAGCTGGTTGCGTTTATTAATACGCAGCTGTTCATCGCGCGAAAGCgg                 >  1:1009621/1‑69 (MQ=255)                                             gCCACGTACTTTGTCGCGTTTTAGCTGGTTGCGTTTATTAATACGCAGCTGTTCATCGCGCGAAAGCgg                 >  1:2145172/1‑69 (MQ=255)                                             gCCACGTACTTTGTCGCGTTTTAGCTGGTTGCGTTTATTAATACGCAGCTGTTCATCGCGCGAAAGCgg                 >  1:1532467/1‑69 (MQ=255)                                             gCCACGTACTTTGTCGCGTTTTAGCTGGTTGCGTTTATTAATACGCAGCTGTTCATCGCGCGAAAGCgg                 >  1:174182/1‑69 (MQ=255)                                             gCCACGTACTTTGTCGCGTTTTAGCTGGTTGCGTTTATTAATACGCAGCTGTTCATCGCGCGAAAGCgg                 >  1:1755769/1‑69 (MQ=255)                                             gCCACGTACTTTGTCGCGTTTTAGCTGGTTGCGTTTATTAATACGCAGCTGTTCATCGCGCGAAAGCgc                 >  1:2143018/1‑68 (MQ=255)                                                 cGTACTTTGTCGCGTTTTAGCTGGTTGCGTttatt                                               <  1:2682665/35‑1 (MQ=255)                                                       ttGTCGCGTTTTAGTTGGTTGCGTTTATTAATACGCAGCTGCTCATCGCGCGAAAg                    >  1:1934968/1‑56 (MQ=255)                                                       ttGTCGCGTTTTAGTTGGTTGCGTTTATTAATACGCAGCTGCTCATCGCGCGAAAg                    >  1:447420/1‑56 (MQ=255)                                                       ttGTCGCGTTTTAGTTGGTTGCGTTTATTAATACGCAGCTGCTCATCGCGCGAAAg                    >  1:1398586/1‑56 (MQ=255)                                                       ttGTCGCGTTTTAGTTGGTTGCGTTTATTAATACGCAGCTGCTCATCGCGCGAAAg                    >  1:12876/1‑56 (MQ=255)                                                       ttGTCGCGTTTTAGTTGGTTGCGTTTATTAATACGCAGCTGCTCATCGCGCGAAAg                    >  1:172908/1‑56 (MQ=255)                                                       ttGTCGCGTTTTAGTTGGTTGCGTTTATTAATACGCAGCTGCTCATCGCGCGAAAg                    >  1:3009849/1‑56 (MQ=255)                                                       ttGTCGCGTTTTAGTTGGTTGCGTTTATTAATACGCAGCTGCTCATCGCGCGAAAg                    >  1:1010168/1‑56 (MQ=255)                                                       ttGTCGCGTTTTAGTTGGTTGCGTTTATTAATACGCAGCTGCTCATCGCGCGAAAGCGGATTAGtttt        <  1:1219648/68‑1 (MQ=255)                                                       ttGTCGCGTTTTAGTTGGTTGCGTTTATTAATACGCAGCTGCTCATCGCGCGAAAGCGGATTAGtttt        <  1:1725334/68‑1 (MQ=255)                                                       ttGTCGCGTTTTAGTTGGTTGCGTTTATTAATACGCAGCTGCTCATCGCGCGAAAGCGGATTAGtttt        <  1:1272485/68‑1 (MQ=255)                                                       ttGTCGCGTTTTAGTTGGTTGCGTTTATTAATACGCAGCTGCTCATCGCGCGAAAGCGGATTAGtttt        <  1:115207/68‑1 (MQ=255)                                                       ttGTCGCGTTTTAGTTGGTTGCGTTTATTAATACGCAGCTGCTCATCGCGCGAAAGCGGATTAGtttt        <  1:1140768/68‑1 (MQ=255)                                                       ttGTCGCGTTTTAGTTGGTTGCGTTTATTAATACGCAGCTGCTCATCGCGCGAAAGCGGATTAGtttt        <  1:1480214/68‑1 (MQ=255)                                                       ttGTCGCGTTTTAGTTGGTTGCGTTTATTAATACGCAGCTGCTCATCGCGCGAAAGCGGATTAGtttt        <  1:790954/68‑1 (MQ=255)                                                       ttGTCGCGTTTTAGTTGGTTGCGTTTATTAATACGCAGCTGCTCATCGCGCGAAAGCGGATTAGtttt        <  1:443191/68‑1 (MQ=255)                                                       ttGTCGCGTTTTAGTTGGTTGCGTTTATTAATACGCAGCTGCTCATCGCGCGAAAGCGGATTAGtttt        <  1:1537600/68‑1 (MQ=255)                                                       ttGTCGCGTTTTAGTTGGTTGCGTTTATTAATACGCAGCTGCTCATCGCGCGAAAGCGGATTAGtttt        <  1:158296/68‑1 (MQ=255)                                                       ttGTCGCGTTTTAGTTGGTTGCGTTTATTAATACGCAGCTGCTCATCGCGCGAAAGCGGATTAGtttt        <  1:3363229/68‑1 (MQ=255)                                                       ttGTCGCGTTTTAGTTGGTTGCGTTTATTAATACGCAGCTGCTCATCGCGCGAAAGCGGATTAGtttt        <  1:3174842/68‑1 (MQ=255)                                                       ttGTCGCGTTTTAGTTGGTTGCGTTTATTAATACGCAGCTGCTCATCGCGCGAAAGCGGATTAGtttt        <  1:310764/68‑1 (MQ=255)                                                       ttGTCGCGTTTTAGTTGGTTGCGTTTATTAATACGCAGCTGCTCATCGCGCGAAAGCGGATTAGtttt        <  1:2034053/68‑1 (MQ=255)                                                       ttGTCGCGTTTTAGTTGGTTGCGTTTATTAATACGCAGCTGCTCATCGCGCGAAAGCGGATTAGtttt        <  1:2701606/68‑1 (MQ=255)                                                            gcgTTTTAGTTGGTTGCGTTTATTAATACGCAGCTGCTCATCGCGCGAAAGCGGATTAGTTTTCGgtcg  <  1:1713886/69‑1 (MQ=255)                                                            gcgTTTTAGTTGGTTGCGTTTATTAATACGCAGCTGCTCATCGCGCGAAAGCGGATTAGTTTTCGgtcg  <  1:2783943/69‑1 (MQ=255)                                                            gcgTTTTAGTTGGTTGCGTTTATTAATACGCAGCTGCTCATCGCGCGAAAGCGGATTAGTTTTCGgtcg  <  1:349054/69‑1 (MQ=255)                                                            gcgTTTTAGTTGGTTGCGTTTATTAATACGCAGCTGCTCATCGCGCGAAAGCGGATTAGTTTTCGgtcg  <  1:1060152/69‑1 (MQ=255)                                                            gcgTTTTAGTTGGTTGCGTTTATTAATACGCAGCTGCTCATCGCGCGAAAGCGGATTAGTTTTCGgtcg  <  1:2989943/69‑1 (MQ=255)                                                               ttttAGTTGGTTGCGTTTATTAATACGCAGCTGCTCATCGCGCGAAAGCGGATTAGtttt        <  1:452646/60‑1 (MQ=255)                                                                     |                                                            GCCTCGACAGCTTCCGCGTTCAGCTTCAGTTCGACACGCTTAAGGCCACGTACTTTGTCGCGTTTTAGCTGGTTGCGTTTATTAATACGCAGCTGTTCATCGCGCGAAAGCGGATTAGTTTTCGGTCG  >  minE/435009‑435136

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