Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	W3110S.gb	4,591,445	G→A	T18I (ACA→ATA)	hsdR ←	endonuclease R

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	W3110S.gb	4,591,445	0	G	A	100.0%	127.1 / NA	40	T18I (ACA→ATA)	hsdR	endonuclease R
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (35/5);  total (35/5)

ATGGCATAAGTGAAGTCGTTGACGCCCTTCAGGAATTCAAAATTGGATTTATTCATCATTGTTATTAATCCATTGCTGTGCGGGCCTGTCCAAATATTTAAGGCCCATAACATCTCAT  >  W3110S.gb/4591385‑4591502                                                             |                                                          aTGGCATAAGTGAAGTCGTT‑ACGCCCTTCAGGAATTCAAAATTGGATTTATTCATCattatt                                                         <  1:4176526/62‑1 (MQ=255)  tGGCATAAGTGAAGTCGTTGACGCCCTTCAGGAATTCAAAATTGGATTTATTCATCattatt                                                         <  1:3957021/62‑1 (MQ=255)  tGGCATAAGTGAAGTCGTTGACGCCCTTCAGGAATTCAAAATTGGATTTATTCATCattatt                                                         <  1:3739371/62‑1 (MQ=255)  tGGCATAAGTGAAGTCGTTGACGCCCTTCAGGAATTCAAAATTGGATTTATTCATCattatt                                                         <  1:3982433/62‑1 (MQ=255)    gCATAAGTGAAGTCGTTGACGCCCTTCAGGAATTCAAAATTGGATTTATTCATCAttatta                                                        >  1:266799/1‑61 (MQ=255)    gCATAAGTGAAGTCGTTGACGCCCTTCAGGAATTCAAAATTGGATTTATTCATCATtattat                                                       <  1:3272973/62‑1 (MQ=255)        aaGTGAAGTCGTTGACGCCCTTCAGGAATTCAAAATTGGATTTATTCATCATTATTATTAAt                                                   <  1:3093892/62‑1 (MQ=255)                tCGTTGACGCCCTTCAGGAATTCAAAATTGGATTTATTCATCATTATTATTAATCCATTGCt                                           >  1:4636364/1‑62 (MQ=255)                tCGTTGACGCCCTTCAGGAATTCAAAATTGGATTTATTCATCATTATTATTAATCCATTGCt                                           >  1:3738675/1‑62 (MQ=255)                tCGTTGACGCCCTTCAGGAATTCAAAATTGGATTTATTCATCATTATTATTAATCCATTGCt                                           >  1:3138806/1‑62 (MQ=255)                tCGTTGACGCCCTTCAGGAATTCAAAATTGGATTTATTCATCATTATTATTAATCCATTGCt                                           >  1:4320449/1‑62 (MQ=255)                tCGTTGACGCCCTTCAGGAATTCAAAATTGGATTTATTCATCATTATTATTAATCCATTGCt                                           >  1:4364451/1‑62 (MQ=255)                tCGTTGACGCCCTTCAGGAATTCAAAATTGGATTTATTCATCATTATTATTAATCCATTGCt                                           >  1:1358852/1‑62 (MQ=255)                    tGACGCCCTTCAGGAATTCAAAATTGGATTTATTCATCattatt                                                         <  1:4132973/44‑1 (MQ=255)                        gCCCTTCAGGAATTCAAAATTGGATTTATTCATCattatt                                                         <  1:4224490/40‑1 (MQ=255)                        gCCCTTCAGGAATTCAAAATTGGATTTATTCATCATTATTATTAATCCATTGCTGTGCGGGc                                   >  1:3974059/1‑62 (MQ=255)                         cccTTCAGGAATTCAAAATTGGATTTATTCATCATTATTATTAATCCATTGCTGTGCGGGc                                   >  1:3934748/1‑61 (MQ=255)                         cccTTCAGGAATTCAAAATTGGATTTATTCATCATTATTATTAATCCATTGCTGTGCGGGc                                   >  1:2045498/1‑61 (MQ=255)                         cccTTCAGGAATTCAAAATTGGATTTATTCATCATTATTATTAATCCATTGCTGTGCGGGc                                   >  1:261385/1‑61 (MQ=255)                         cccTTCAGGAATTCAAAATTGGATTTATTCATCATTATTATTAATCCATTGCTGTGCGGGc                                   >  1:1290038/1‑61 (MQ=255)                         cccTTCAGGAATTCAAAATTGGATTTATTCATCATTATTATTAATCCATTGCTGTGCGGGc                                   >  1:2188214/1‑61 (MQ=255)                         cccTTCAGGAATTCAAAATTGGATTTATTCATCATTATTATTAATCCATTGCTGTGCGGGc                                   >  1:850069/1‑61 (MQ=255)                         cccTTCAGGAATTCAAAATTGGATTTATTCATCATTATTATTAATCCATTGCTGTGCGGGc                                   >  1:660760/1‑61 (MQ=255)                         cccTTCAGGAATTCAAAATTGGATTTATTCATCATTATTATTAATCCATTGCTGTGCGGGc                                   >  1:3761698/1‑61 (MQ=255)                                     tCAAAATTGGATTTATTCATCATTATTATTAATCCATTGCTGTGCGGGCCTGTCCAAatat                       >  1:1544865/1‑61 (MQ=255)                                     tCAAAATTGGATTTATTCATCATTATTATTAATCCATTGCTGTGCGGGCCTGTCCAAatat                       >  1:734934/1‑61 (MQ=255)                                     tCAAAATTGGATTTATTCATCATTATTATTAATCCATTGCTGTGCGGGCCTGTCCAAatat                       >  1:2101079/1‑61 (MQ=255)                                     tCAAAATTGGATTTATTCATCATTATTATTAATCCATTGCTGTGCGGGCCTGTCCAAATAtt                      <  1:1036539/62‑1 (MQ=255)                                     tCAAAATTGGATTTATTCATCATTATTATTAATCCATTGCTGTGCGGGCCTGTCCAAATAtt                      <  1:4500016/62‑1 (MQ=255)                                           ttGGATTTATTCATCATTATTATTAATCCATTGCTGTGCGGGCCTGTCCAAATATTTAAGGc                <  1:2746680/62‑1 (MQ=255)                                           ttGGATTTATTCATCATTATTATTAATCCATTGCTGTGCGGGCCTGTCCAAATATTTAAGGc                <  1:4585178/62‑1 (MQ=255)                                                  tATTCATCATTATTATTAATCCATTGCTGTGCGGGCCTg                                >  1:4490314/1‑39 (MQ=255)                                                  tATTCATCATTATTATTAATCCATTGCTGTGCGGGCCTg                                >  1:621932/1‑39 (MQ=255)                                                  tATTCATCATTATTATTAATCCATTGCTGTGCGGGCCTg                                >  1:3057690/1‑39 (MQ=255)                                                  tATTCATCATTATTATTAATCCATTGCTGTGCGGGCCTg                                >  1:3009842/1‑39 (MQ=255)                                                  tATTCATCATTATTATTAATCCATTGCTGTGCGGGCCTg                                >  1:821665/1‑39 (MQ=255)                                                  tATTCATCATTATTATTAATCCATTGCTGTGCGGGCCTg                                >  1:2397557/1‑39 (MQ=255)                                                  tATTCATCATTATTATTAATCCATTGCTGTGCGGGCCTg                                >  1:2386647/1‑39 (MQ=255)                                                  tATTCATCATTATTATTAATCCATTGCTGTGCGGGCCTg                                >  1:2355382/1‑39 (MQ=255)                                                  tATTCATCATTATTATTAATCCATTGCTGTGCGGGCCTg                                >  1:1604672/1‑39 (MQ=255)                                                  tATTCATCATTATTATTAATCCATTGCTGTGCGGGCCTg                                >  1:983543/1‑39 (MQ=255)                                                    ttCATCATTATTATTAATCCATTGCTGTGCGGGCCTGTCCAAATATTTAAGGCCCATAAcat       >  1:769987/1‑62 (MQ=255)                                                    ttCATCATTATTATTAATCCATTGCTGTGCGGGCCTGTCCAAATATTTAAGGCCCATAAcat       >  1:1184092/1‑62 (MQ=255)                                                    ttCATCATTATTATTAATCCATTGCTGTGCGGGCCTGTCCAAATATTTAAGGCCCATAAcat       >  1:3063871/1‑62 (MQ=255)                                                    ttCATCATTATTATTAATCCATTGCTGTGCGGGCCTGTCCAAATATTTAAGGCCCATAAcat       >  1:1395305/1‑62 (MQ=255)                                                    ttCATCATTATTATTAATCCATTGCTGTGCGGGCCTGTCCAAATATTTAAGGCCCATAAca        >  1:2415407/1‑61 (MQ=255)                                                    ttCATCAGTATTATTAATCCATTGCTGTGCGGGCCTGTCCAAATATTTAAGGCCCATAAcat       >  1:1423890/1‑62 (MQ=255)                                                        tcatTATTATTAATCCATTGCTGTTCGGGCCTGTCCAAATATTTAAGGCCCATAAcatctca   >  1:3823377/1‑62 (MQ=255)                                                          attattattAATCCATTGCTGTGCGGGCCTGTCCAAATATTTAAGGCCCATAAcatctcat  >  1:4387804/1‑61 (MQ=255)                                                          attattattAATCCATTGCTGTGCGGGCCTGTCCAAATATTTAAGGCCCATAAcatctcat  >  1:3890233/1‑61 (MQ=255)                                                          attattattAATCCATTGCTGTGCGGGCCTGTCCAAATATTTAAGGCCCATAAcatctcat  >  1:355911/1‑61 (MQ=255)                                                          attattattAATCCATTGCTGTGCGGGCCTGTCCAAATATTTAAGGCCCATAAcatctcat  >  1:2742300/1‑61 (MQ=255)                                                          attattattAATCCATTGCTGTGCGGGCCTGTCCAAATATTTAAGGCCCATAAcatctcat  >  1:2426917/1‑61 (MQ=255)                                                             |                                                          ATGGCATAAGTGAAGTCGTTGACGCCCTTCAGGAATTCAAAATTGGATTTATTCATCATTGTTATTAATCCATTGCTGTGCGGGCCTGTCCAAATATTTAAGGCCCATAACATCTCAT  >  W3110S.gb/4591385‑4591502

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