Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	W3110S.gb	4,580,803	C→T	pseudogene (1987/2721 nt)	yjiV →	ECK4334:JW5954:b4486; conserved hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	W3110S.gb	4,580,803	0	C	T	100.0%	108.9 / NA	35	pseudogene (1987/2721 nt)	yjiV	ECK4334:JW5954:b4486; conserved hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (22/13);  total (22/13)

TGACGCTGGCCTCTCGCAGCCAGCAGGCAACGATTCCGGGGCCGGAGTGGCACCTGAATGACGAACTGGTGGTGCGTAGTCTCGGTTATAAAACGGTAGA  >  W3110S.gb/4580762‑4580861                                          |                                                           tgaCGCTGGCCTCTCGCAGCCAGCAGGCAACGATTCCGGGGTCGGAGTGGCACCTGAATGa                                         >  1:3091922/1‑61 (MQ=255)         gCCTCTCGCAGCCAGCAGGTAACGATTCCGGGGTCgg                                                         >  1:1233190/1‑37 (MQ=37)          ccTCTCGCAGCCAGCAGGCAACGATTCCGGGGTCGGAGTGGCTCCTGAATGACGAACtggtg                               >  1:2563845/1‑62 (MQ=255)          ccTCTCGCAGCCAGCAGGCAACGATTCCGGGGTCGGAGTGGCACCTGAATGACGAACtggtg                               >  1:598219/1‑62 (MQ=255)          ccTCTCGCAGCCAGCAGGCAACGATTCCGGGGTCGGAGTGGCACCTGAATGACGAACtggtg                               >  1:592786/1‑62 (MQ=255)          ccTCTCGCAGCCAGCAGGCAACGATTCCGGGGTCGGAGTGGCACCTGAATGACGAACtggtg                               >  1:1595269/1‑62 (MQ=255)          ccTCTCGCAGCCAGCAGGCAACGATTCCGGGGTCGGAGTGGCACCTGAATGACGAACtggtg                               >  1:470092/1‑62 (MQ=255)          ccTCTCGCAGCCAGCAGGCAACGATTCCGGGGTCGGAGTGGCACCTGAATGACGAACtggtg                               >  1:4422168/1‑62 (MQ=255)          ccTCTCGCAGCCAGCAGGCAACGATTCCGGGGTCGGAGTGGCACCTGAATGACGAACtggtg                               >  1:4173618/1‑62 (MQ=255)                          ggCAACGATTCCGGGGTCGGAGTGGCACCTGAATGACGAACTGGTGGTGCGTAGTCTCg                  <  1:7951/59‑1 (MQ=255)                                gATTCCGGGGTCGGAGTGGCACCTGAATGACGAACTGGTGGTGCGTAGTCTCGGTTATaaaa         <  1:3009275/62‑1 (MQ=255)                                gATTCCGGGGTCGGAGTGGCACCTGAATGACGAACTGGTGGTGCGTAGTCTCGGTTATaaaa         <  1:544352/62‑1 (MQ=255)                                gATTCCGGGGTCGGAGTGGCACCTGAATGACGAACTGGTGGTGCGTAGTCTCGGTTATaaaa         <  1:4627495/62‑1 (MQ=255)                                gATTCCGGGGTCGGAGTGGCACCTGAATGACGAACTGGTGGTGCGTAGTCTCGGTTATaaaa         <  1:2325093/62‑1 (MQ=255)                                gATTCCGGGGTCGGAGTGGCACCTGAATGACGAACTGGTGGTGCGTAGTCTCGGTTATaaaa         <  1:2643910/62‑1 (MQ=255)                                  ttCCGGGGTCGGAGTGGCACCTGAATGACGAACTGGTGGTGCGTAGTCTCGGTTATAAAAc        >  1:2836911/1‑61 (MQ=255)                                    ccGGGGTCGGAGTGGCACCTGAATGACGAACTGGTGGTGCGTAGTCTCGGTTATAAAACGGt     <  1:1406429/62‑1 (MQ=255)                                    ccGGGGTCGGAGTGGCACCTGAATGACGAACTGGTGGTGCGTAGTCTCGGTTATAAAACGGt     <  1:1516490/62‑1 (MQ=255)                                    ccGGGGTCGGAGTGGCACCTGAATGACGAACTGGTGGTGCGTAGTCTCGGTTATAAAACGGt     <  1:1530761/62‑1 (MQ=255)                                    ccGGGGTCGGAGTGGCACCTGAATGACGAACTGGTGGTGCGTAGTCTCGGTTATAAAACGGt     <  1:1710757/62‑1 (MQ=255)                                    ccGGGGTCGGAGTGGCACCTGAATGACGAACTGGTGGTGCGTAGTCTCGGTTATAAAACGGt     <  1:3029125/62‑1 (MQ=255)                                    ccGGGGTCGGAGTGGCACCTGAATGACGAACTGGTGGTGCGTAGTCTCGGTTATAAAACGGt     <  1:2984656/62‑1 (MQ=255)                                     cGGGGTCGGAGTGGCACCTGAATGACGAACTGGTGGTGCGTAGTCTCGGTTATAAAACGGt     >  1:2490437/1‑61 (MQ=255)                                     cGGGGTCGGAGTGGCACCTGAATGACGAACTGGTGGTGCGTAGTCTCGGTTATAAAACGGt     >  1:2171768/1‑61 (MQ=255)                                     cGGGGTCGGAGTGGCACCTGAATGACGAACTGGTGGTGCGTAGTCTCGGTTATAAAACGGt     <  1:2825246/61‑1 (MQ=255)                                       gggTCGGAGTGGCACCTGAATGACGAACTGGTGGTGCGTAGTCTCGGTTATAAAACGGTAGa  >  1:3212808/1‑62 (MQ=255)                                       gggTCGGAGTGGCACCTGAATGACGAACTGGTGGTGCGTAGTCTCGGTTATAAAACGGTAGa  >  1:3905753/1‑62 (MQ=255)                                       gggTCGGAGTGGCACCTGAATGACGAACTGGTGGTGCGTAGTCTCGGTTATAAAACGGTAGa  >  1:403492/1‑62 (MQ=255)                                       gggTCGGAGTGGCACCTGAATGACGAACTGGTGGTGCGTAGTCTCGGTTATAAAACGGTAGa  >  1:4219313/1‑62 (MQ=255)                                       gggTCGGAGTGGCACCTGAATGACGAACTGGTGGTGCGTAGTCTCGGTTATAAAACGGTAGa  >  1:4257702/1‑62 (MQ=255)                                       gggTCGGAGTGGCACCTGAATGACGAACTGGTGGTGCGTAGTCTCGGTTATAAAACGGTAGa  >  1:439995/1‑62 (MQ=255)                                       gggTCGGAGTGGCACCTGAATGACGAACTGGTGGTGCGTAGTCTCGGTTATAAAACGGTAGa  >  1:2432636/1‑62 (MQ=255)                                       gggTCGGAGTGGCACCTGAATGACGAACTGGTGGTGCGTAGTCTCGGTTATAAAACGGTAGa  >  1:2684337/1‑62 (MQ=255)                                       gggTCGGAGTGGCACCTGAATGACGAACTGGTGGTGCGTAGTCTCGGTTATAAAACGGTAGa  >  1:2882815/1‑62 (MQ=255)                                       gggTCGGAGTGGCACCTGAATGACGAACTGGTGGTGCGTAGTCTCGGTTATAAAACGGTAGa  >  1:665082/1‑62 (MQ=255)                                          |                                                           TGACGCTGGCCTCTCGCAGCCAGCAGGCAACGATTCCGGGGCCGGAGTGGCACCTGAATGACGAACTGGTGGTGCGTAGTCTCGGTTATAAAACGGTAGA  >  W3110S.gb/4580762‑4580861

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