Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	W3110S.gb	330,092	A→G	N469S (AAC→AGC)	betT →	choline transporter of high affinity

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	W3110S.gb	330,092	0	A	G	100.0%	230.4 / NA	62	N469S (AAC→AGC)	betT	choline transporter of high affinity
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (53/9);  total (53/9)

TGCGCGTCTTCTGGTCGGTGGCGATTGGCCTGCTGACGCTCGGCATGCTGATGACTAACGGGATATCCGCGCTGCAAAACACCACGGTGATTATGGGGCTGCCGTTCAGCTTTGTG  >  W3110S.gb/330035‑330150                                                          |                                                           tgCGCGTCTTCTGGTCGGTGGCGATTGGCCTGCTGACGCTCGGCATGCTGATGACTAGCggg                                                        <  1:326813/62‑1 (MQ=255) tgCGCGTCTTCTGGTCGGTGGCGATTGGCCTGCTGACGCTCGGCATGCTGATGACTAGCggg                                                        <  1:3714982/62‑1 (MQ=255)  gcgcGTCTTCTGGTCGGTGGCGATTGGCCTGCTGACGCTCGGCATGCTGATGACTAGCggg                                                        <  1:3612360/61‑1 (MQ=255)     cgTCTTCTGGTCGGTGGCGATTGGCCTGCTGACGCTCGGCATGCTGATGACTAGCGGGata                                                     <  1:2245226/61‑1 (MQ=255)     cgTCTTCTGGTCGGTGGCGATTGGCCTGCTGACGCTCGGCATGCTGATGACTAGCGGGata                                                     <  1:2532300/61‑1 (MQ=255)     cgTCTTCTGGTCGGTGGCGATTGGCCTGCTGACGCTCGGCATGCTGATGACTAGCGGGata                                                     <  1:3048610/61‑1 (MQ=255)     cgTCTTCTGGTCGGTGGCGATTGGCCTGCTGACGCTCGGCATGCTGATGACTAGCGGGata                                                     <  1:1885319/61‑1 (MQ=255)                cGGTGGCGATTGGCCTGCTGACGCTCGGCATGCTGATGACTAGCGGGATATCCGCGCTGCaa                                         >  1:1737839/1‑62 (MQ=255)                cGGTGGCGATTGGCCTGCTGACGCTCGGCATGCTGATGACTAGCGGGATATCCGCGCTGCaa                                         >  1:2319696/1‑62 (MQ=255)                cGGTGGCGATTGGCCTGCTGACGCTCGGCATGCTGATGACTAGCGGGATATCCGCGCTGCaa                                         >  1:2303407/1‑62 (MQ=255)                cGGTGGCGATTGGCCTGCTGACGCTCGGCATGCTGATGACTAGCGGGATATCCGCGCTGCaa                                         >  1:256015/1‑62 (MQ=255)                cGGTGGCGATTGGCCTGCTGACGCTCGGCATGCTGATGACTAGCGGGATATCCGCGCTGCaa                                         >  1:322617/1‑62 (MQ=255)                cGGTGGCGATTGGCCTGCTGACGCTCGGCATGCTGATGACTAGCGGGATATCCGCGCTGCaa                                         >  1:3833707/1‑62 (MQ=255)                cGGTGGCGATTGGCCTGCTGACGCTCGGCATGCTGATGACTAGCGGGATATCCGCGCTGCaa                                         >  1:1636686/1‑62 (MQ=255)                cGGTGGCGATTGGCCTGCTGACGCTCGGCATGCTGATGACTAGCGGGATATCCGCGCTGCaa                                         >  1:140637/1‑62 (MQ=255)                cGGTGGCGATTGGCCTGCTGACGCTCGGCATGCTGATGACTAGCGGGATATCCGCGCTGCaa                                         >  1:1343759/1‑62 (MQ=255)                cGGTGGCGATTGGCCTGCTGACGCTCGGCATGCTGATGACTAGCGGGATATCCGCGCTGCaa                                         >  1:1082800/1‑62 (MQ=255)                cGGTGGCGATTGGCCTGCTGACGCTCGGCATGCTGATGACTAGCGGGATATCCGCGCTGCa                                          >  1:513263/1‑61 (MQ=255)                       gATTGGCCTGCTGACGCTCGGCATGCTGATGACTAGCGGGATATCCGCGCTGCAAAAcacca                                  <  1:466669/62‑1 (MQ=255)                        aTTGGCCTGCTGACGCTCGGCATGCTGATGACTAGCGGGATATCCGCGCTGCaa                                         >  1:1709775/1‑54 (MQ=255)                        aTTGGCCTGCTGACGCTCGGCATGCTGATGACTAGCGGGATATCCGCGCTGCaa                                         >  1:2985096/1‑54 (MQ=255)                          tGGCCTGCTGACGCTCGGCATGCTGATGACTAGCGGGATATCCGCGCTGCAAAACACCACg                                >  1:812189/1‑61 (MQ=255)                          tGGCCTGCTGACGCTCGGCATGCTGATGACTAGCGGGATATCCGCGCTGCAAAACACCACg                                >  1:3712045/1‑61 (MQ=255)                          tGGCCTGCTGACGCTCGGCATGCTGATGACTAGCGGGATATCCGCGCTGCAAAACACCACg                                >  1:4560754/1‑61 (MQ=255)                          tGGCCTGCTGACGCTCGGCATGCTGATGACTAGCGGGATATCCGCGCTGCAAAACACCACg                                >  1:2601216/1‑61 (MQ=255)                          tGGCCTGCTGACGCTCGGCATGCTGATGACTAGCGGGATATCCGCGCTGCAAAACACCACg                                >  1:694167/1‑61 (MQ=255)                          tGGCCTGCTGACGCTCGGCATGCTGATGACTAGCGGGATATCCGCGCTGCAAAACACCACg                                >  1:703696/1‑61 (MQ=255)                          tGGCCTGCTGACGCTCGGCATGCTGATGACTAGCGGGATATCCGCGCTGCAAAACACCACg                                >  1:2200514/1‑61 (MQ=255)                          tGGCCTGCTGACGCTCGGCATGCTGATGACTAGCGGGATATCCGCGCTGCAAAACACCACg                                >  1:2140803/1‑61 (MQ=255)                          tGGCCTGCTGACGCTCGGCATGCTGATGACTAGCGGGATATCCGCGCTGCAAAACACCACg                                >  1:1994457/1‑61 (MQ=255)                           ggCCTGCTGACGCTCGGCATGCTGATGACTAGCGGGATATCCGCGCTGCaa                                         >  1:1727976/1‑51 (MQ=255)                           ggCCTGCTGACGCTCGGCATGCTGATGACTAGCGGGATATCCGCGCTGCaa                                         >  1:1317997/1‑51 (MQ=255)                           ggCCTGCTGACGCTCGGCATGCTGATGACTAGCGGGATATCCGCGCTGCaa                                         >  1:423190/1‑51 (MQ=255)                           ggCCTGCTGACGCTCGGCATGCTGATGACTAGCGGGATATCCGCGCTGCaa                                         >  1:2366381/1‑51 (MQ=255)                           ggCCTGCTGACGCTCGGCATGCTGATGACTAGCGGGATATCCGCGCTGCaa                                         >  1:1851530/1‑51 (MQ=255)                           ggCCTGCTGACGCTCGGCATGCTGATGACTAGCGGGATATCCGCGCTGCAAAACACCACGGt                              >  1:1996032/1‑62 (MQ=255)                           ggCCTGCTGACGCTCGGCATGCTGATGACTAGCGGGATATCCGCGCTGCAAAACACCACGGt                              >  1:1537158/1‑62 (MQ=255)                           ggCCTGCTGACGCTCGGCATGCTGATGACTAGCGGGATATCCGCGCTGCAAAACACCACGGt                              >  1:2131686/1‑62 (MQ=255)                           ggCCTGCTGACGCTCGGCATGCTGATGACTAGCGGGATATCCGCGCTGCAAAACACCACGGt                              >  1:3983280/1‑62 (MQ=255)                           ggCCTGCTGACGCTCGGCATGCTGATGACTAGCGGGATATCCGCGCTGCAAAACACCACGGt                              >  1:1822645/1‑62 (MQ=255)                                  tgACGCTCGGCATGCTGATGACTAGCGGGATATCCGCGCTGCaaa                                        >  1:3803017/1‑45 (MQ=255)                                  tgACGCTCGGCATGCTGATGACTAGCGGGATATCCGCGCTGCAAAACACCACGGTGATTATg                       >  1:593777/1‑62 (MQ=255)                                  tgACGCTCGGCATGCTGATGACTAGCGGGATATCCGCGCTGCAAAACACCACGGTGATTATg                       >  1:342812/1‑62 (MQ=255)                                  tgACGCTCGGCATGCTGATGACTAGCGGGATATCCGCGCTGCAAAACACCACGGTGATTATg                       >  1:2757988/1‑62 (MQ=255)                                  tgACGCTCGGCATGCTGATGACTAGCGGGATATCCGCGCTGCAAAACACCACGGTGATTATg                       >  1:3798945/1‑62 (MQ=255)                                        tCGGCATGCTGATGACTAGCGGGATATCCGCGCTGCAAAACACCACGGTGATTATGGGGc                   >  1:3971348/1‑60 (MQ=255)                                        tCGGCATGCTGATGACTAGCGGGATATCCGCGCTGCAAAACACCACGGTGATTATGGGGCTg                 >  1:1689937/1‑62 (MQ=255)                                        tCGGCATGCTGATGACTAGCGGGATATCCGCGCTGCAAAACACCACGGTGATTATGGGGCTg                 >  1:4304078/1‑62 (MQ=255)                                        tCGGCATGCTGATGACTAGCGGGATATCCGCGCTGCAAAACACCACGGTGATTATGGGGCTg                 >  1:3184980/1‑62 (MQ=255)                                                cTGATGACTAGCGGGATATCCGCGCTGCAAAAcacca                                  <  1:2555837/37‑1 (MQ=255)                                                       cTAGCGGGATATCCGCGCTGCAAAACACCACGGTGa                            >  1:4631780/1‑36 (MQ=255)                                                       cTAGCGGGATATCCGCGCTGCAAAACACCACGGTGAt                           >  1:2862063/1‑37 (MQ=255)                                                       cTAGCGGGATATCCGCGCTGCAAAACACCACGGTGAt                           >  1:2101449/1‑37 (MQ=255)                                                       cTAGCGGGATATCCGCGCTGCAAAACACCACGGTGAt                           >  1:3571521/1‑37 (MQ=255)                                                       cTAGCGGGATATCCGCGCTGCAAAACACCACGGTGAt                           >  1:4076914/1‑37 (MQ=255)                                                       cTAGCGGGATATCCGCGCTGCAAAACACCACGGTGAt                           >  1:4291552/1‑37 (MQ=255)                                                       cTAGCGGGATATCCGCGCTGCAAAACACCACGGTGAt                           >  1:2890080/1‑37 (MQ=255)                                                       cTAGCGGGATATCCGCGCTGCAAAACACCACGGTGAt                           >  1:1525588/1‑37 (MQ=255)                                                       cTAGCGGGATATCCGCGCTGCAAAACACCACGGTGAt                           >  1:1437202/1‑37 (MQ=255)                                                       cTAGCGGGATATCCGCGCTGCAAAACACCACGGTGAt                           >  1:1167099/1‑37 (MQ=255)                                                       cTAGCGGGATATCCGCGCTGCAAAACACCACGGTGAt                           >  1:1045783/1‑37 (MQ=255)                                                       cTAGCGGGATATCCGCGCTGCAAAACACCACGGTGATTATGGGGCTGCCGTTCAGCTTtgtg  >  1:1876132/1‑62 (MQ=255)                                                          |                                                           TGCGCGTCTTCTGGTCGGTGGCGATTGGCCTGCTGACGCTCGGCATGCTGATGACTAACGGGATATCCGCGCTGCAAAACACCACGGTGATTATGGGGCTGCCGTTCAGCTTTGTG  >  W3110S.gb/330035‑330150

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