Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	W3110S.gb	631,602	A→G	*66* (TAA→TAG)	ybdD →	conserved hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	W3110S.gb	631,602	0	A	G	100.0%	227.6 / NA	61	*66* (TAA→TAG)	ybdD	conserved hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (55/6);  total (55/6)

TTTTTCCGTGAGCGGCAGGACGCGCGCTACGGTGGAAAAGGCGGCGCGCGCTGCTGCTAAATTTGTACGTCAGGCTTTAAACGATTCCACTTTTTTGAACGCTGCACGCAAC  >  W3110S.gb/631543‑631654                                                            |                                                     tttttCCGTGAGCGGCAGGACGCGCGCTACGGTGGAAAAGGCGGCGCGCGCTGCTGCTAGAt                                                    >  1:3314187/1‑62 (MQ=255) tttttCCGTGAGCGGCAGGACGCGCGCTACGGTGGAAAAGGCGGCGCGCGCTGCTGCTAGAt                                                    >  1:2519284/1‑62 (MQ=255) tttttCCGTGAGCGGCAGGACGCGCGCTACGGTGGAAAAGGCGGCGCGCGCTGCTGCTAGAt                                                    >  1:1428974/1‑62 (MQ=255)   tttCCGTGAGCGGCAGGACGCGCGCTACGGTGGAAAAGGCGGCGCGCGCTGCTGCTAGAtt                                                   >  1:3239500/1‑61 (MQ=255)   tttCCGTGAGCGGCAGGACGCGCGCTACGGTGGAAAAGGCGGCGCGCGCTGCTGCTAGAtt                                                   >  1:952770/1‑61 (MQ=255)   tttCCGTGAGCGGCAGGACGCGCGCTACGGTGGAAAAGGCGGCGCGCGCTGCTGCTAGAtt                                                   >  1:746928/1‑61 (MQ=255)   tttCCGTGAGCGGCAGGACGCGCGCTACGGTGGAAAAGGCGGCGCGCGCTGCTGCTAGAtt                                                   >  1:73348/1‑61 (MQ=255)   tttCCGTGAGCGGCAGGACGCGCGCTACGGTGGAAAAGGCGGCGCGCGCTGCTGCTAGAtt                                                   >  1:3649278/1‑61 (MQ=255)   tttCCGTGAGCGGCAGGACGCGCGCTACGGTGGAAAAGGCGGCGCGCGCTGCTGCTAGAtt                                                   >  1:1886044/1‑61 (MQ=255)   tttCCGTGAGCGGCAGGACGCGCGCTACGGTGGAAAAGGCGGCGCGCGCTGCTGCTAGAtt                                                   >  1:1559685/1‑61 (MQ=255)   tttCCGTGAGCGGCAGGACGCGCGCTACGGTGGAAAAGGCGGCGCGCGCTGCTGCTAGAtt                                                   >  1:1171758/1‑61 (MQ=255)         tGAGCGGCAGGACGCGCGCTACGGTGGAAAAGGCGGCGCGCGCTGCTGCTAGATTTGTAc                                              >  1:3752689/1‑60 (MQ=255)         tGAGCGGCAGGACGCGCGCTACGGTGGAAAAGGCGGCGCGCGCTGCTGCTAGATTTGTAc                                              >  1:3448637/1‑60 (MQ=255)         tGAGCGGCAGGACGCGCGCTACGGTGGAAAAGGCGGCGCGCGCTGCTGCTAGATTTGTAc                                              >  1:2371727/1‑60 (MQ=255)         tGAGCGGCAGGACGCGCGCTACGGTGGAAAAGGCGGCGCGCGCTGCTGCTAGATTTGTAc                                              >  1:2280677/1‑60 (MQ=255)          gAGCGGCAGGACGCGCGCTACGGTGGAAAAGGCGGCGCGCGCTGCTGCTAGATTTGTACGTc                                           >  1:1605616/1‑62 (MQ=255)          gAGCGGCAGGACGCGCGCTACGGTGGAAAAGGCGGCGCGCGCTGCTGCTAGATTTGTACGTc                                           >  1:3841362/1‑62 (MQ=255)          gAGCGGCAGGACGCGCGCTACGGTGGAAAAGGCGGCGCGCGCTGCTGCTAGATTTGTACGTc                                           >  1:3982483/1‑62 (MQ=255)          gAGCGGCAGGACGCGCGCTACGGTGGAAAAGGCGGCGCGCGCTGCTGCTAGATTTGTACGTc                                           >  1:4350786/1‑62 (MQ=255)                         cgcTACGGTGGAAAAGGCGGCGCGCGCTGCTGCTAGATTTGTa                                               >  1:1934752/1‑43 (MQ=255)                         cgcTACGGTGGAAAAGGCGGCGCGCGCTGCTGCTAGATTTGTACGTCAGGCTTTAAACGAt                             >  1:1239033/1‑61 (MQ=255)                         cgcTACGGTGGAAAAGGCGGCGCGCGCTGCTGCTAGATTTGTACGTCAGGCTTTAAACGAt                             >  1:1167994/1‑61 (MQ=255)                         cgcTACGGTGGAAAAGGCGGCGCGCGCTGCTGCTAGATTTGTACGTCAGGCTTTAAACGAt                             >  1:1110044/1‑61 (MQ=255)                         cgcTACGGTGGAAAAGGCGGCGCGCGCTGCTGCTAGATTTGTACGTCAGGCTTTAAACGAt                             >  1:597534/1‑61 (MQ=255)                         cgcTACGGTGGAAAAGGCGGCGCGCGCTGCTGCTAGATTTGTACGTCAGGCTTTAAACGAt                             >  1:305725/1‑61 (MQ=255)                         cgcTACGATGGAAAAGGCGGCGCGCGCTGCTGCTAGATTTGTACGTCAGGCTTTAAACGAt                             >  1:199076/1‑61 (MQ=255)                           cTACGGTGGAAAAGGCGGCGCGCGCTGCTGCTAGATTTGTACGTCAGGCTTTAAACGATTcc                          >  1:4529353/1‑62 (MQ=255)                           cTACGGTGGAAAAGGCGGCGCGCGCTGCTGCTAGATTTGTACGTCAGGCTTTAAACGATTcc                          >  1:3261103/1‑62 (MQ=255)                           cTACGGTGGAAAAGGCGGCGCGCGCTGCTGCTAGATTTGTACGTCAGGCTTTAAACGATTcc                          >  1:454843/1‑62 (MQ=255)                           cTACGGTGGAAAAGGCGGCGCGCGCTGCTGCTAGATTTGTACGTCAGGCTTTAAACGATTcc                          >  1:3144788/1‑62 (MQ=255)                            tACGGTGGAAAAGGCGGCGCGCGCTGCTGCTAGATTTGTACGTCAGGCTTTAAACGa                              >  1:3644875/1‑57 (MQ=255)                            tACGGTGGAAAAGGCGGCGCGCGCTGCTGCTAGATTTGTACGTCAGGCTTTAAACGATTCCa                         >  1:2943130/1‑62 (MQ=255)                            tACGGTGGAAAAGGCGGCGCGCGCTGCTGCTAGATTTGTACGTCAGGCTTTAAACGATTCCa                         >  1:1222636/1‑62 (MQ=255)                            tACGGTGGAAAAGGCGGCGCGCGCTGCTGCTAGATTTGTACGTCAGGCTTTAAACGATTCCa                         >  1:3396976/1‑62 (MQ=255)                            tACGGTGGAAAAGGCGGCGCGCGCTGCTGCTAGATTTGTACGTCAGGCTTTAAACGATTCCa                         >  1:1153206/1‑62 (MQ=255)                            tACGGTGGAAAAGGCGGCGCGCGCTGCTGCTAGATTTGTACGTCAGGCTTTAAACGATTCCa                         >  1:3584179/1‑62 (MQ=255)                                gTGGAAAAGGCGGCGCGCGCTGCTGCTAGATTTGTACGTCAGGCTTTAAACGATTCCACttt                     <  1:4595092/62‑1 (MQ=255)                                gTGGAAAAGGCGGCGCGCGCTGCTGCTAGATTTGTACGTCAGGCTTTAAACGATTCCACttt                     <  1:2079016/62‑1 (MQ=255)                                            gcgctcgcTGCTGCTAGATTTGTACGTCAGGCTTTAAACGATTCCACTTTTTTGAACGCTg          >  1:904686/1‑61 (MQ=255)                                            gcgcgcgcTGCTGCTAGATTTGTACGTCAGGCTTTAAACGATTCCACTTTTTTGAACGCTGc         >  1:406867/1‑62 (MQ=255)                                            gcgcgcgcTGCTGCTAGATTTGTACGTCAGGCTTTAAACGATTCCACTTTTTTGAACGCTGc         >  1:4219850/1‑62 (MQ=255)                                            gcgcgcgcTGCTGCTAGATTTGTACGTCAGGCTTTAAACGATTCCACTTTTTTGAACGCTGc         >  1:3195626/1‑62 (MQ=255)                                            gcgcgcgcTGCTGCTAGATTTGTACGTCAGGCTTTAAACGATTCCACTTTTTTGAACGCTGc         >  1:4616317/1‑62 (MQ=255)                                            gcgcgcgcTGCTGCTAGATTTGTACGTCAGGCTTTAAACGATTCCACTTTTTTGAACGCTGc         >  1:69122/1‑62 (MQ=255)                                            gcgcgcgcTGCTGCTAGATTTGTACGTCAGGCTTTAAACGATTCCACTTTTTTGAACGCTGc         >  1:2172957/1‑62 (MQ=255)                                            gcgcgcgcTGCTGCTAGATTTGTACGTCAGGCTTTAAACGATTCCACTTTTTTGAACGCTGc         >  1:737550/1‑62 (MQ=255)                                            gcgcgcgcTGCTGCTAGATATGTACGTCAGGCTTTAAACGATTCAACTTTTTTGAACGCTGc         >  1:1697721/1‑62 (MQ=255)                                               cgcgcTGCTGCTAGATTTGTACGTCAGGCTTTAAACGATTCCACTTTTTTGAACGCTGCAc       >  1:111503/1‑61 (MQ=255)                                               cgcgcTGCTGCTAGATTTGTACGTCAGGCTTTAAACGATTCCACTTTTTTGAACGCTGCAc       >  1:452814/1‑61 (MQ=255)                                               cgcgcTGCTGCTAGATTTGTACGTCAGGCTTTAAACGATTCCACTTTTTTGAACGCTGCAc       >  1:354721/1‑61 (MQ=255)                                               cgcgcTGCTGCTAGATTTGTACGTCAGGCTTTAAACGATTCCACTTTTTTGAACGCTGCAc       >  1:3463248/1‑61 (MQ=255)                                               cgcgcTGCTGCTAGATTTGTACGTCAGGCTTTAAACGATTCCACTTTTTTGAACGCTGCAc       >  1:3041708/1‑61 (MQ=255)                                               cgcgcTGCTGCTAGATTTGTACGTCAGGCTTTAAACGATTCCACTTTTTTGAACGCTGCAc       >  1:2709789/1‑61 (MQ=255)                                               cgcgcTGCTGCTAGATTTGTACGTCAGGCTTTAAACGATTCCACTTTTTTGAACGCTGCAc       >  1:1247610/1‑61 (MQ=255)                                               cgcgcTGCTGCTAGATTTGTACGTCAGGCTTTAAACGATTCCACTTTTTTGAACGCTGCAc       >  1:171502/1‑61 (MQ=255)                                               cgcgcTGCTGCTAGATTTGTACGTCAGGCTTTAAACGAGTCCACTTTTTTGAACGCTGCAc       >  1:4464377/1‑61 (MQ=255)                                               cgcgcTGCCGCTAGATTTGTACGTCAGGCTTTAAACGATTCCACTTTTTTGAACGCTGCAc       >  1:182408/1‑61 (MQ=255)                                                 cgcTGCTGCTAGATTTGTACGTCAGGCTTTAAACGATTCCACTTTTTTGAACGCTGCACGCa    <  1:411850/62‑1 (MQ=255)                                                 cgcTGCTGCTAGATTTGTACGTCAGGCTTTAAACGATTCCACTTTTTTGAACGCTGCACGCa    <  1:2547962/62‑1 (MQ=255)                                                   ctgctgctAGATTTGTACGTCAGGCTTTAAACGATTCCACTTTTTTGAACGCTGCACGCAAc  <  1:4036635/62‑1 (MQ=255)                                                   ctgctgctAGATTTGTACGTCAGGCTTTAAACGATTCCACTTTTTTGAACGCTGCACGCAAc  <  1:151844/62‑1 (MQ=255)                                                            |                                                     TTTTTCCGTGAGCGGCAGGACGCGCGCTACGGTGGAAAAGGCGGCGCGCGCTGCTGCTAAATTTGTACGTCAGGCTTTAAACGATTCCACTTTTTTGAACGCTGCACGCAAC  >  W3110S.gb/631543‑631654

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑