Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	W3110S.gb	3,824,280	T→C	pseudogene (413/687 nt)	rph →	ECK3633:JW3618:b3643; ribonuclease PH

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	W3110S.gb	3,824,280	0	T	C	100.0%	133.4 / NA	38	pseudogene (413/687 nt)	rph	ECK3633:JW3618:b3643; ribonuclease PH
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (18/20);  total (18/20)

GGCTGATGGTGGCACGCGTACCGCGTCGATTACGGGTGCCTGCGTGGCGCTGGTAGATGCGCTACAGAAGCTGGTGGAAAACGGCAAGCTGAAAACCAATCCGATGAAAGGGATGGTAGC  >  W3110S.gb/3824227‑3824346                                                      |                                                                   ggCTGATGGTGGCACGCGTACCGCGTCGATTACGGGTGCCTGCGTGGCGCTGGCAGATGCGCTACAGAAg                                                    >  1:2663695/1‑70 (MQ=255) ggCTGATGGTGGCACGCGTACCGCGTCGATTACGGGTGCCTGCGTGGCGCTGGCAGATGCGCTACAGAAGc                                                   >  1:1175508/1‑71 (MQ=255) ggCTGATGGTGGCACGCGTACCGCGTCGATTACGGGTGCCTGCGTGGCGCTGGCAGATGCGCTACAGAAGc                                                   >  1:4047367/1‑71 (MQ=255) ggCTGATGGTGGCACGCGTACCGCGTCGATTACGGGTGCCTGCGTGGCGCTGGCAGATGCGCTACAGAAGc                                                   >  1:1846445/1‑71 (MQ=255)   cTGATGGTGGCACGCGTACCGCGTCGATTACGGGTGCCTGCGTGGCGCTGGCAGATGCGCTACAGAAGCtg                                                 >  1:353669/1‑71 (MQ=255)   cTGATGGTGGCACGCGTACCGCGTCGATTACGGGTGCCTGCGTGGCGCTGGCAGATGCGCTACAGAAGCtg                                                 >  1:1299971/1‑71 (MQ=255)   cTGATGGTGGCACGCGTACCGCGTCGATTACGGGTGCCTGCGTGGCGCTGGCAGATGCGCTACAGAAGCtg                                                 >  1:432682/1‑71 (MQ=255)   cTGATGGTGGCACGCGTACCGCGTCGATTACGGGTGCCTGCGTGGCGCTGGCAGATGCGCTACAGAAGCtg                                                 >  1:1606176/1‑71 (MQ=255)   cTGATGGTGGCACGCGTACCGCGTCGATTACGGGTGCCTGCGTGGCGCTGGCAGATGCGCTACAGAAGCtg                                                 >  1:2870387/1‑71 (MQ=255)   cTGATGGTGGCACGCGTACCGCGTCGATTACGGGTGCCTGCGTGGCGCTGGCAGATGCGCTACAGAAGCtg                                                 >  1:44026/1‑71 (MQ=255)   cTGATGGTGGCACGCGTACCGCGTCGATTACGGGTGCCTGCGTGGCGCTGGCAGATGCGCTACAGAAGCtg                                                 >  1:2323041/1‑71 (MQ=255)   cTGATGGTGGCACGCGTACCGCGTCGATTACGGGTGCCTGCGTGGCGCTGGCAGATGCGCTACAGAAGCtg                                                 >  1:263260/1‑71 (MQ=255)      aTGGTGGCACGCGTACCGCGTCGATTACGGGTGCCTTCGTGGCGCTGGCAGATGCGCTACAGAAGCtgg                                                <  1:59376/69‑1 (MQ=255)      aTGGTGGCACGCGTACCGCGTCGATTACGGGTGCCTGCGTGGCGCTGGCAGATGCGCTACAGAAGCtgg                                                <  1:2845737/69‑1 (MQ=255)      aTGGTGGCACGCGTACCGCGTCGATTACGGGTGCCTGCGTGGCGCTGGCAGATGCGCTACAGAAGCtgg                                                <  1:3186425/69‑1 (MQ=255)        ggtggCACGCGTACCGCGTCGATTACGGGTGCCTGCGTGGCGCTGGCAGATGCGCTACAGAAGCtggtgg                                             >  1:3855454/1‑70 (MQ=255)         gtggCACGCGTACCGCGTCGATTACGGGTGCCTGCGTGGCGCTGGCAGATGCGCTACAGAAGCTGTTGGaa                                           <  1:1439840/71‑1 (MQ=255)         gtggCACGCGTACCGCGTCGATTACGGGTGCCTGCGTGGCGCTGGCAGATGCGCTACAGAAGCTGGTGGaa                                           <  1:3224757/71‑1 (MQ=255)         gtggCACGCGTACCGCGTCGATTACGGGTGCCTGCGTGGCGCTGGCAGATGCGCTACAGAAGCTGGTGGaa                                           <  1:3754370/71‑1 (MQ=255)         gtggCACGCGTACCGCGTCGATTACGGGTGCCTGCGTGGCGCTGGCAGATGCGCTACAGAAGCTGGTGGaa                                           <  1:3625324/71‑1 (MQ=255)         gtggCACGCGTACCGCGTCGATTACGGGTGCCTGCGTGGCGCTGGCAGATGCGCTACAGAAGCTGGTGGaa                                           <  1:3311097/71‑1 (MQ=255)         gtggCACGCGTACCGCGTCGATTACGGGTGCCTGCGTGGCGCTGGCAGATGCGCTACAGAAGCTGGTGGaa                                           <  1:1540668/71‑1 (MQ=255)         gtggCACGCGTACCGCGTCGATTACGGGTGCCTGCGTGGCGCTGGCAGATGCGCTACAGAAGCTGGTGGaa                                           <  1:2888450/71‑1 (MQ=255)         gtggCACGCGTACCGCGTCGATTACGGGTGCCTGCGTGGCGCTGGCAGATGCGCTACAGAAGCTGGTGGaa                                           <  1:1932459/71‑1 (MQ=255)         gtggCACGCGTACCGCGTCGATTACGGGTGCCTGCGTGGCGCTGGCAGATGCGCTACAGAAGCTGGTGGa                                            >  1:207501/1‑70 (MQ=255)                 cgTACCGCGTCGATTACGGGTGCCTGCGTGGCGCTGGCAGATGCTCTACAGAAGCtgg                                                <  1:3072431/58‑1 (MQ=39)                                             tGGCGCTGTCAGATGCGCTACAGAAGCTGGTGGAAAACGGCAAGCTGAAAACCAAccc                    <  1:3261315/58‑1 (MQ=255)                                               gctcTGGCAGATGCGCTACAGAAGCTGGTGGAAAACGGCAAGCTGAAAACCAACCCGATGaaa             <  1:3046114/60‑1 (MQ=255)                                               gcgcTGGCAGATGCGCTACAGAAGCTGGTGGAAAACGGCAAGCTGAAAACCAACCCGa                  >  1:3167706/1‑58 (MQ=255)                                               gcgcTGGCAGATGCGCTACAGAAGCTGGTGGAAAACGGCAAGCTGAAAACCAACCCGa                  >  1:1096631/1‑58 (MQ=255)                                               gcgcTGGCAGATGCGCTACAGAAGCTGGTGGAAAACGGCAAGCTGAAAACCAACCCGa                  >  1:1249454/1‑58 (MQ=255)                                               gcgcTGGCAGATGCGCTACAGAAGCTGGTGGAAAACGGCAAGCTGAAAACCAACCCGa                  >  1:3617346/1‑58 (MQ=255)                                               gcgcTGGCAGATGCGCTACAGAAGCTGGTGGAAAACGGCAAGCTGAAAACCAACCCGATGaaa             <  1:938408/63‑1 (MQ=255)                                               gcgcTGGCAGATGCGCTACAGAAGCTGGTGGAAAACGGCAAGCTGAAAACCAACCCGATGaaa             <  1:3429957/63‑1 (MQ=255)                                               gcgcTGGCAGATGCGCTACAGAAGCTGGTGGAAAACGGCAAGCTGAAAACCAACCCGATGaaa             <  1:2920025/63‑1 (MQ=255)                                                  cTGGCAGATGCGCTACAGAAGCTGGTGGAAAACGGCAAGCTGAAAACCAACCCGATGAAAGGGATGGTAGc  <  1:2835262/71‑1 (MQ=255)                                                  cTGGCAGATGCGCTACAGAAGCTGGTGGAAAACGGCAAGCTGAAAACCAACCCGATGAAAGGGATGGTAGc  <  1:1963341/71‑1 (MQ=255)                                                  cTGGCAGATGCGCTACAGAAGCTGGTGGAAAACGGCAAGCTGAAAACCAACCCGATGAAAGGGATGGTAGc  <  1:1427225/71‑1 (MQ=255)                                                      |                                                                   GGCTGATGGTGGCACGCGTACCGCGTCGATTACGGGTGCCTGCGTGGCGCTGGTAGATGCGCTACAGAAGCTGGTGGAAAACGGCAAGCTGAAAACCAATCCGATGAAAGGGATGGTAGC  >  W3110S.gb/3824227‑3824346

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