Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	W3110S.gb	3,704,876	C→T	M164I (ATG→ATA)	trkD ←	potassium transporter

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	W3110S.gb	3,704,876	0	C	T	100.0%	150.5 / NA	45	M164I (ATG→ATA)	trkD	potassium transporter
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (24/21);  total (24/21)

CCAGTCCTGCCAGAATCAAAAACCAGGTCAGCATGATCGGCGCAAACAGCTTACCGACCATAGCCGTGCCATGTTTTTGAATCATAAATAATAAAGTGAGAACGATAATTGAGAGG  >  W3110S.gb/3704818‑3704933                                                           |                                                          ccAGTCCTGCCAGAATCAAAAACCAGGTCAGCATGATCGGCGCAAACAGCTTACCGACtata                                                        <  1:3603519/62‑1 (MQ=255) ccAGTCCTGCCAGAATCAAAAACCAGGTCAGCATGATCGGCGCAAACAGCTTACCGACtata                                                        <  1:2821159/62‑1 (MQ=255)  cAGTCCTGCCAGAATCAAAAACCAGGTCAGCATGATCGGCGCAAACAGCTTACCGACtata                                                        <  1:2128051/61‑1 (MQ=255)     tCCTGCCAGAATCAAAAACCAGGTCAGCATGATCGGCGCAAACAGCTTACCGACTATAGCCg                                                    >  1:106118/1‑62 (MQ=255)     tCCTGCCAGAATCAAAAACCAGGTCAGCATGATCGGCGCAAACAGCTTACCGACTATAGCCg                                                    >  1:4695808/1‑62 (MQ=255)     tCCTGCCAGAATCAAAAACCAGGTCAGCATGATCGGCGCAAACAGCTTACCGACTATAGCCg                                                    >  1:3506724/1‑62 (MQ=255)     tCCTGCCAGAATCAAAAACCAGGTCAGCATGATCGGCGCAAACAGCTTACCGACTATAGCCg                                                    >  1:3040420/1‑62 (MQ=255)     tCCTGCCAGAATCAAAAACCAGGTCAGCATGATCGGCGCAAACAGCTTACCGACTATAGCCg                                                    >  1:2836618/1‑62 (MQ=255)     tCCTGCCAGAATCAAAAACCAGGTCAGCATGATCGGCGCAAACAGCTTACCGACTATAGCCg                                                    >  1:2542647/1‑62 (MQ=255)     tCCTGCCAGAATCAAAAACCAGGTCAGCATGATCGGCGCAAACAGCTTACCGACTATAGCCg                                                    >  1:197263/1‑62 (MQ=255)     tCCTGCCAGAATCAAAAACCAGGTCAGCATGATCGGCGCAAACAGCTTACCGACTATAGCCg                                                    >  1:1395076/1‑62 (MQ=255)     tCCTGCCAGAATCAAAAACCAGGTCAGCATGATCGGCGCAAACAGCTTACCGACTATAGCCg                                                    >  1:1407061/1‑62 (MQ=255)             gAATCAAAAACCAGGTCAGCATGATCGGCGCAAACAGCTTACCGACTATAGCCGTGCCATGt                                            >  1:1131102/1‑62 (MQ=255)             gAATCAAAAACCAGGTCAGCATGATCGGCGCAAACAGCTTACCGACTATAGCCGTGCCATGt                                            >  1:3018132/1‑62 (MQ=255)             gAATCAAAAACCAGGTCAGCATGATCGGCGCAAACAGCTTACCGACTATAGCCGTGCCATGt                                            >  1:1303920/1‑62 (MQ=255)                         aGGTCAGCATGATCGGCGCAAACAGCTTACCGACTATAGCCGTGCCATGTTTTTGAATCata                                <  1:894156/62‑1 (MQ=255)                         aGGTCAGCATGATCGGCGCAAACAGCTTACCGACTATAGCCGTGCCATGTTTTTGAATCata                                <  1:1207577/62‑1 (MQ=255)                         aGGTCAGCATGATCGGCGCAAACAGCTTACCGACTATAGCCGTGCCATGTTTTTGAATCata                                <  1:1920743/62‑1 (MQ=255)                         aGGTCAGCATGATCGGCGCAAACAGCTTACCGACTATAGCCGTGCCATGTTTTTGAATCata                                <  1:1852073/62‑1 (MQ=255)                                    aTCGGCGCAAACAGCTTACCGACTATAGCCGTGCCATGTTTTTGAATCATAAATAATaaa                       >  1:75294/1‑60 (MQ=255)                                    aTCGGCGCAAACAGCTTACCGACTATAGCCGTGCCATGTTTTTGAATCATAAATAATaaa                       >  1:2043608/1‑60 (MQ=255)                                    aTCGGCGCAAACAGCTTACCGACTATAGCCGTGCCATGTTTTTGAATCATAAATAATaaa                       >  1:3576430/1‑60 (MQ=255)                                    aTCGGCGCAAACAGCTTACCGACTATAGCCGTGCCATGTTTTTGAATCATAAATAATaaa                       >  1:3539223/1‑60 (MQ=255)                                    aTCGGCGCAAACAGCTTACCGACTATAGCCGTGCCATGTTTTTGAATCATAAATAATaaa                       >  1:3200040/1‑60 (MQ=255)                                    aTCGGCGCAAACAGCTTACCGACTATAGCCGTGCCATGTTTTTGAATCATAAATAATaaa                       >  1:310557/1‑60 (MQ=255)                                    aTCGGCGCAAACAGCTTACCGACTATAGCCGTGCCATGTTTTTGAATCATAAATAATaaa                       >  1:3088063/1‑60 (MQ=255)                                         cgcAAACAGCTTACCGACTATAGCCGTGCCATGTTTTTGAATCATAAATAATAAAGTgaga                 <  1:1077851/61‑1 (MQ=255)                                         cgcAAACAGCTTACCGACTATAGCCGTGCCATGTTTTTGAATCATAAATAATAAAGTgaga                 <  1:1447092/61‑1 (MQ=255)                                         cgcAAACAGCTTACCGACTATAGCCGTGCCATGTTTTTGAATCATAAATAATAAAGTgaga                 <  1:3162552/61‑1 (MQ=255)                                         cgcAAACAGCTTACCGACTATAGCCGTGCCATGTTTTTGAATCATAAATAATAAAGTgaga                 <  1:2970833/61‑1 (MQ=255)                                          gcAAACAGCTTACCGACTATAGCCGTGACATGTTTTTGAATCATAAATAATAAAGTgaga                 <  1:1886640/60‑1 (MQ=255)                                           cAAACAGCTTACCGACTATAGCCGTGCCATGTTTTTGAATCATAAATAATAAAGTGAGAAc               >  1:199346/1‑61 (MQ=255)                                           cAAACAGCTTACCGACTATAGCCGTGCCATGTTTTTGAATCATAAATAATAAAGTGAGAAc               >  1:2048844/1‑61 (MQ=255)                                           cAAACAGCTTACCGACTATAGCCGTGCCATGTTTTTGAATCATAAATAATAAAGTGAGAAc               >  1:1337692/1‑61 (MQ=255)                                           cAAACAGCTTACCGACTATAGCCGTGCCATGTTTTTGAATCATAAATAATAAAGTGAGAACg              <  1:1060202/62‑1 (MQ=255)                                           cAAACAGCTTACCGACTATAGCCGTGCCATGTTTTTGAATCATAAATAATAAAGTGAGAACg              <  1:3912001/62‑1 (MQ=255)                                            aaaCAGCTTACCGACTATAGCCGTGCCATGTTTTTGAATCATAAATAATAAAGTGAGAACg              <  1:1578005/61‑1 (MQ=255)                                              aCAGCTTACCGACTATAGCCGTGCCATGTTTTTGAATCATAAATAATAAAGTGAGAACGATa           >  1:1224947/1‑62 (MQ=255)                                                  cTTACCGACTATAGCCGTGCCATGTTTTTGAATCATAAATAATAAAGTGAGAACGATAATTg       <  1:1676439/62‑1 (MQ=255)                                                  cTTACCGACTATAGCCGTGCCATGTTTTTGAATCATAAATAATAAAGTGAGAACGATAATTg       <  1:4459617/62‑1 (MQ=255)                                                   ttACCGACTATAGCCGTGCCATGTTTTTGAATCATAAATAATAAAGTGAGAACGATAATTg       <  1:2677934/61‑1 (MQ=255)                                                      ccGACTATAGCCGTGCCATGTTTTTGAATCATAAATAATAAAGTGAGAACGATa           >  1:646660/1‑54 (MQ=255)                                                        gACTATAGCCGTGCCATGTTTTTGAATCTTAAATAATAAAGTGAGAACGATAATTGAGAgg  <  1:2204744/61‑1 (MQ=255)                                                        gACTATAGCCGTGCCATGTTTTTGAATCATAAATAATAAAGTTAGAACGATAATTGAGAgg  <  1:6761/61‑1 (MQ=255)                                                        gACTATAGCCGTGCCATGTTTTTGAATCATAAATAATAAAGTGAGAACGATAATTGAGAgg  <  1:4422442/61‑1 (MQ=255)                                                        gACTATAGCCGTGCCATGTTTTTGAATCATAAATAATAAAGTGAGAACGATAATTGAGAgg  <  1:4432663/61‑1 (MQ=255)                                                        gACTATAGCCGTGCCATGTTTTTGAATCATAAATAATAAAGTGAGAACGATAATTGAGAgg  <  1:1540982/61‑1 (MQ=255)                                                        gACTATAGCCGTGCCATGTTTTTGAATCATAAATAATAAAGTGAGAACGATAATTGAGAgg  <  1:1410565/61‑1 (MQ=255)                                                           |                                                          CCAGTCCTGCCAGAATCAAAAACCAGGTCAGCATGATCGGCGCAAACAGCTTACCGACCATAGCCGTGCCATGTTTTTGAATCATAAATAATAAAGTGAGAACGATAATTGAGAGG  >  W3110S.gb/3704818‑3704933

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