Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	minE	899,984	(T)9→8	intergenic (+41/‑172)	pspE → / → ycjM	thiosulfate:cyanide sulfurtransferase/predicted glucosyltransferase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	899,976	0	T	.	100.0%	247.8 / NA	57	intergenic (+33/‑180)	pspE/ycjM	thiosulfate:cyanide sulfurtransferase/predicted glucosyltransferase
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base . (23/34);  total (23/34)

AAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCGTTTTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAA  >  minE/899912‑900034                                                                 |                                                           aaGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCGttttttt                                                      <  1:2031380/71‑1 (MQ=255) aaGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCGttttttt                                                      <  1:305418/71‑1 (MQ=255) aaGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCGttttttt                                                      <  1:1299518/71‑1 (MQ=255)    aCATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGATCGGCGCTTATGCGCCGttttttttct                                                   <  1:1171558/71‑3 (MQ=255)    aCATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCGttttttttct                                                   <  1:2620105/71‑3 (MQ=255)    aCATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCGttttttttct                                                   >  1:27272/1‑69 (MQ=255)    aCATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCGttttttttct                                                   >  1:3270791/1‑69 (MQ=255)    aCATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCGttttttttct                                                   <  1:3465684/71‑3 (MQ=255)    aCATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCGttttttttct                                                   <  1:3476215/71‑3 (MQ=255)    aCATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCGttttttttct                                                   <  1:3684423/71‑3 (MQ=255)    aCATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCGttttttttct                                                   <  1:4018761/71‑3 (MQ=255)    aCATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCGttttttttct                                                   <  1:4051556/71‑3 (MQ=255)    aCATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCGttttttttct                                                   <  1:4113364/71‑3 (MQ=255)    aCATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCGttttttttct                                                   <  1:4148848/71‑3 (MQ=255)    aCATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCGttttttttct                                                   <  1:732342/71‑3 (MQ=255)    aCATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCGttttttttc                                                    >  1:3471629/1‑69 (MQ=255)    aCATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCGttttttttc                                                    >  1:1355500/1‑69 (MQ=255)    aCATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCGttttttttc                                                    >  1:749665/1‑69 (MQ=255)     cATCGCTATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCGttttttttct                                                   <  1:1098237/70‑3 (MQ=255)     cATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGATCGGCGCTTATGCGCCGttttttttct                                                   <  1:2002141/70‑3 (MQ=255)     cATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCGttttttttct                                                   <  1:942916/70‑3 (MQ=255)     cATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCGttttttttct                                                   <  1:1957416/70‑3 (MQ=255)        cGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCG‑TTTTTTTTCTGTca                                              >  1:3262251/1‑71 (MQ=255)        cGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCG‑TTTTTTTTCTGTca                                              >  1:1717678/1‑71 (MQ=255)        cGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCG‑TTTTTTTTCTGTca                                              >  1:654309/1‑71 (MQ=255)         gCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCGttttttttctg                                                  <  1:3152780/67‑4 (MQ=255)         gCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCGttttttttctg                                                  <  1:2849567/67‑4 (MQ=255)         gCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCGttttttttctg                                                  <  1:347010/67‑4 (MQ=255)         gCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCGttttttttctg                                                  <  1:4215720/67‑4 (MQ=255)         gCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCGttttttttctg                                                  <  1:1451499/67‑4 (MQ=255)         gCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCGttttttttctg                                                  <  1:3983764/67‑4 (MQ=255)         gCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCGttttttttctg                                                  <  1:150511/67‑4 (MQ=255)         gCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCGttttttttctg                                                  <  1:3655115/67‑4 (MQ=255)         gCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCG‑TTTTTTTTCTGTca                                              >  1:2233105/1‑70 (MQ=255)         gCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCG‑TTTTTTTTCTGTca                                              >  1:1945781/1‑70 (MQ=255)         gCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCG‑TTTTTTTTCTGTca                                              >  1:2777999/1‑70 (MQ=255)           aaTGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCG‑TTTTTTTTCTGTCacac                                           <  1:2142269/71‑1 (MQ=255)           aaTGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCG‑TTTTTTTTCTGTCacac                                           <  1:3454610/71‑1 (MQ=255)            aTGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCGtttttt                                                       >  1:76575/1‑59 (MQ=255)            aTGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCG‑TTTTTTTTCTGTcacagc                                          >  1:1153879/1‑69 (MQ=255)            aTGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCG‑TTTTTTTTCTGTCACAcc                                          >  1:3304168/1‑71 (MQ=255)            aTGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCG‑TTTTTTTTCTGTCACAcc                                          >  1:631753/1‑71 (MQ=255)                 gAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCGttttttttc                                                    >  1:3063914/1‑56 (MQ=255)                 gAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCG‑TTTTTTTTCTGTCACACCTTAtt                                     >  1:1715505/1‑71 (MQ=255)                 gAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCG‑TTTTTTTTCTGTCACACCTTAtt                                     >  1:2369709/1‑71 (MQ=255)                 gAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCG‑TTTTTTTTCTGTCACACCTTAtt                                     >  1:2329749/1‑71 (MQ=255)                 gAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCG‑TTTTTTTTCTGTCACACCATAtt                                     >  1:540043/1‑71 (MQ=255)                      tCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCG‑TTTTTTTTCTGTCACACCTTATTTACAt                                <  1:412401/71‑1 (MQ=255)                      tCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCG‑TTTTTTTTCTGTCACACCTTATTTACAt                                <  1:607725/71‑1 (MQ=255)                      tCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCG‑TTTTTTTTCTGTCACACCTTATTTACAt                                <  1:1067768/71‑1 (MQ=255)                      tCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCG‑TTTTTTTTCTGTCACACCTTATTTACAt                                <  1:181406/71‑1 (MQ=255)                         aaGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCG‑TTTTTTTTCTGTCACACCTTATTTACATccc                             <  1:2720422/71‑1 (MQ=255)                         aaGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCG‑TTTTTTTTCTGTCACACCTTATTTACATccc                             <  1:92473/71‑1 (MQ=255)                         aaGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCG‑TTTTTTTTCTGTCACACCTTATTTACATccc                             <  1:191995/71‑1 (MQ=255)                         aaGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCG‑TTTTTTTTCTGTCACACCTTATTTACATccc                             <  1:2329095/71‑1 (MQ=255)                         aaGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCG‑TTTTTTTTCTGTCACACCTTATTTACATccc                             <  1:2298542/71‑1 (MQ=255)                          aGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCG‑TTTTTTTTCTGTcaca                                            >  1:1377802/1‑55 (MQ=255)                           ggTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCG‑TTTTTTTTCTGTCACACCTTATTTACATCCCCa                           <  1:3173457/71‑1 (MQ=255)                           ggTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCG‑TTTTTTTTCTGTCACACCTTATTTACATCCCCa                           <  1:2625975/71‑1 (MQ=255)                           ggTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCG‑TTTTTTTTCTGTCACACCTTATTTACATCCCCa                           <  1:2825054/71‑1 (MQ=255)                           ggTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCG‑TTTTTTTTCTGTCACACCTTATTTACATCCCCa                           <  1:2959309/71‑1 (MQ=255)                              tAAAAGAATTCACCATGAGCGGCGCTTATGCGCCG‑TTTTTTTTCTGTca                                              >  1:2314603/1‑49 (MQ=255)                              tAAAAGAATTCACCATGAGCGGCGCTTATGCGCCG‑TTTTTTTTCTGTca                                              >  1:2267688/1‑49 (MQ=255)                                aaaGAATTCACCATGAGCGGCGCTTATGCGCCG‑TTTTTTTTCTGTCACACCTTATTTACATCCCCATAGat                      <  1:497204/71‑1 (MQ=255)                                aaaGAATTCACCATGAGCGGCGCTTATGCGCCG‑TTTTTTTTCTGTCACACCTTATTTACATCCCCATAGat                      <  1:44576/71‑1 (MQ=255)                                aaaGAATTCACCATGAGCGGCGCTTATGCGCCG‑TTTTTTTTCTGTCACACCTTATTTACATCCCCATAGat                      <  1:2672692/71‑1 (MQ=255)                                aaaGAATTCACCATGAGCGGCGCTTATGCGCCG‑TTTTTTTTCTGTCACACCTTATTTACATCCCCATAGat                      <  1:2219367/71‑1 (MQ=255)                                aaaGAATTCACCATGAGCGGCGCTTATGCGCCG‑TTTTTTTTCTGTCACACCTTATTTACATCCCCATAGat                      <  1:1408715/71‑1 (MQ=255)                                aaaGAATTCACCATGAGCGGCGCTTATGCGCCG‑TTTTTTTTCTGTCACACCTTATTTACATCCCCATAGat                      <  1:3068016/71‑1 (MQ=255)                                aaaGAATTCACCATGAGCGGCGCTTATGCGCCG‑TTTTTTTTCTGTCACACCTTATTTACATCCCCATAGat                      <  1:4012247/71‑1 (MQ=255)                                    aaTTCACCATGAGCGGCGCTTATGCGCCG‑TTTTTTTTCTGTCACACCTTATTTACATCCCCATAGattatt                  >  1:3276802/1‑71 (MQ=255)                                        cACCATGAGCGGCGCTTATGCGCCG‑TTTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGc               <  1:3184112/70‑1 (MQ=255)                                        cACCATGAGCGGCGCTTATGCGCCG‑TTTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGc               <  1:2689765/70‑1 (MQ=255)                                        cACCATGAGCGGCGCTTATGCGCCG‑TTTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGc               <  1:4323159/70‑1 (MQ=255)                                        cACCATGAGCGGCGCTTATGCGCCG‑TTTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGc               <  1:3893923/70‑1 (MQ=255)                                        cACCATGAGCGGCGCTTATGCGCCG‑TTTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGc               <  1:2728269/70‑1 (MQ=255)                                        cACCATGAGCGGCGCTTATGCGCCG‑TTTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGc               <  1:1104260/70‑1 (MQ=255)                                        cACCATGAGCGGCGCTTATGCGCCG‑TTTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGc               <  1:2806180/70‑1 (MQ=255)                                        cACCATGAGCGGCGCTTATGCGCCG‑TTTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGc               <  1:1676693/70‑1 (MQ=255)                                        cACCATGAGCGGCGCTTATGCGCCG‑TTTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCg              >  1:377511/1‑71 (MQ=255)                                        cACCATGAGCGGCGCTTATGCGCCG‑TTTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCg              >  1:74255/1‑71 (MQ=255)                                        cACCATGAGCGGCGCTTATGCGCCG‑TTTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCg              >  1:1654497/1‑71 (MQ=255)                                              gAGCGGAGCTTATGCGCCG‑TTTTTTTTCTGTCACACCTTAtt                                     >  1:559585/1‑42 (MQ=38)                                                 cggcgCTTATGCGCCG‑TTTTTTTTCTGTCACACCTTATTTACATccc                             <  1:4040627/47‑1 (MQ=255)                                                  ggcgCTTCTGCGCCG‑TTTTTTTTCTGTCACACCTTATTTACATCCCCa                           <  1:3997550/48‑1 (MQ=255)                                                  ggcgCTTATGCGCCG‑TTTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTcaca    <  1:2129987/71‑1 (MQ=255)                                                     gcTTATGCGCCG‑TTTTTTTTCTGTCACACCTTATTTACa                                 >  1:1402470/1‑39 (MQ=38)                                                     gcTTATGCGCCG‑TTTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACaaa  >  1:2771192/1‑70 (MQ=255)                                                         aTGCGCCG‑TTTTTTTTCTGTCACACCTTATTTACATCCCCATAGat                      <  1:640612/46‑1 (MQ=255)                                                                 |                                                           AAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCGTTTTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAA  >  minE/899912‑900034

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