Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,624,198	T→C	pseudogene (312/351 nt)	yrhC →	pseudogene fragment

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,624,198	0	T	C	100.0%	88.3 / NA	26	pseudogene (312/351 nt)	yrhC	pseudogene fragment
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (13/13);  total (13/13)

TCGGAGCATTTATACTCGGGCAAATATAAATTTGCAAAAAATGCGAAGATAATGTTTTTAATATTTTGAACAATATGGCATTGAATGATAAGAGTGCTTGCGTTCGAGCTACGGCAATCGAGTCAACGGCCAGGGAAAGATCATCAAAAGTTCACTTTTTGCACTAAATAATTCGCATTTTATGTTTAAAAATTGAGATATTCCTTATTACCTGAAGCTGTTTTTTA  >  NC_000913/3624104‑3624330                                                                                               |                                                                                                                                     tCGGAGCATTTATACTCGGGCAAATAAAAATTTGCAAAAAATGCGAAGATAATGTTTTTAATATTTTGAACAATATGGCATTGAATGATAAGAGCGCTTGCGTTCGAGCTACGGCAATCGAGTCAACGGCCAGGGAAAg                                                                                          >  1:537985/1‑139 (MQ=9)      gCATTTATACTCGGGCAAATAAAAATTTGCAAAAAATGCGAAGATAATGTTTTTAATATTTTGAACAATATGGCATTGAATGATAAGAGCGCTTGCGTTCGAGCTACGGCAATCGAGTCAACGGCCAGGGAAAGatcat                                                                                     <  1:457403/139‑1 (MQ=255)           tataCTCGGGCAAATAAAAATTTGCAAAAAATGCGAAGATAATGTTTTTAATATTTTGAACAATATGGCATTGAATGATAAGAGCGCTTGCGTTCGAGCTACGGCAATCGAGTCAACGGCCAGGGAAAGATCATCaaaa                                                                                >  2:645402/1‑139 (MQ=255)           tataCTCGGGCAAATAAAAATTTGCAAAAAATGCGAAGATAATGTTTTTAATATTTTGAACAATATGGCATTGAATGATAAGAGCGCTTGCGTTCGAGCTACGGCAATCGAGTCAACGGCCAGGGAAAGATCATCaaaa                                                                                >  1:679704/1‑139 (MQ=255)                    gCAAATAAAAATTTGCAAAAAATGCGAAGATAATGTTTTTAATATTTTGAACAATATGGCATTGAATGATAAGAGCGCTTGCGTTCGAGCTACGGCAATCGAGTCAACGGCCAGGGAAAGATCATCAAAAGTTCACttt                                                                       >  1:670525/1‑139 (MQ=255)                    gCAAATAAAAATTTGCAAAAAATGCGAAGATAATGTTTTTAATATTTTGAACAATATGGCATTGAATGATAAGAGCGCTTGCGTTCGAGCTACGGCAATCGAGTCAACGGCCAGGGAAAGATCATCAAAAGTTCACttt                                                                       >  1:853404/1‑139 (MQ=255)                    gCAAATAAAAATTTGCAAAAAATGCGAAGATAATGTTTTTAATATTTTGAACAATATGGCATTGAATGATAAGAGCGCTTGCGTTCGAGCTACGGCAATCGAGTCAACGGCCAGGGAAAGATCATCAAAAGTTAACttt                                                                       >  2:737967/1‑139 (MQ=255)                                                              tatTTTGAACAATATGGCATTGAATGATAAGAGCGCTTGCGTTCGAGCTACGGCAATCGAGTCAACGGCCAGGGAAAGATCATCAAAAGTTCACTTTTTGCACTAAATAATTCGCATTTTATGTTTAAAAATTGAGata                             <  2:853404/139‑1 (MQ=255)                                                               atTTTTAACAATATGCCATTGAATGATAAGATCGCCTGCCTTCGAGCTACGGCAATCGCGTCAACCGCCATTGAAAGATCATCAAAAGTTCACTTTTTGCACTAAATAATTCGCATTTTATGTTTAAAAATTGAGatat                            <  1:737967/139‑1 (MQ=21)                                                               atTTTGAACAATATGGCATTGAATGATAAGAGCGCTTGCGTTCGAGCTACGGCAATCGAGTCAACGGCCAGGGAAAGATCATCAAAAGTTCACTTTTTGCACTAAATAATTCGCATTTTATGTTTAAAAATTGAGatat                            <  2:537985/139‑1 (MQ=255)                                                                 tttGAACAATATGGCATTGAATGATAAGAGCGCTTGCGTTCGAGCTACGGCAATCGAGTCAACGGCCAGGGAAAGATCATCAAAAGTTCACTTTTTGCa                                                                  <  1:1029075/99‑1 (MQ=255)                                                                 tttGAACAATATGGCATTGAATGATAAGAGCGCTTGCGTTCGAGCTACGGCAATCGAGTCAACGGCCAGGGAAAGATCATCAAAAGTTCACTTTTTGCa                                                                  >  2:1029075/1‑99 (MQ=255)                                                                  ttGAACAATATGGCATTGAATGATAAGAGCGCTTGCGTTCGAGCTACGGCAATCGAGTCAACGGCCAGGGAAAGATCATCAAAAGTTCACTTTTTGCACTAAATAATTCGCATTTTATGTTTAAAAATTGAGATATTcc                         <  2:279190/139‑1 (MQ=255)                                                                                  tGAATGATAAGAGCGCTTGCGTTCGAGCTACGGCAATCGAGTCAACGGCCAGGGAAAGATCATCAAAAGTTCACTTTTTGCACTAAATAATTCGCATTTTATGTTTAAAAATTGAGATATTCCTTATTACCTGAAGCTg         >  2:727040/1‑139 (MQ=255)                                                                                    aaTGATAAGAGCGCTTGCGTTCGAGCTACGGCAATCGAGTCAACGGCCAGGGAAAGATCATCAAAAGTTCACTTTTTGCa                                                                  >  1:443274/1‑80 (MQ=255)                                                                                    aaTGATAAGAGCGCTTGCGTTCGAGCTACGGCAATCGAGTCAACGGCCAGGGAAAGATCATCAAAAGTTCACTTTTTGCa                                                                  <  2:443274/80‑1 (MQ=255)                                                                                    aaTGATAAGAGCGCTTGCGTTCGAGCTACGGCAATCGAGTCAACGGCCAGGGAAAGATCATCAAAAGTTCACTTTTTGCACt                                                                >  2:971533/1‑82 (MQ=255)                                                                                    aaTGATAAGAGCGCTTGCGTTCGAGCTACGGCAATCGAGTCAACGGCCAGGGAAAGATCATCAAAAGTTCACTTTTTGCACt                                                                <  1:971533/82‑1 (MQ=255)                                                                                    aaTGATAAGAGCGCTTGCGTTCGAGCTACGGCAATCGAGTCAACGGCCAGGGAAAGATCATCAAAAGTTCACTTTTTGCACTAAATAATTCGCATTTTATGTTTAAAAATTGAGata                             >  2:618975/1‑117 (MQ=255)                                                                                    aaTGATAAGAGCGCTTGCGTTCGAGCTACGGCAATCGAGTCAACGGCCAGGGAAAGATCATCAAAAGTTCACTTTTTGCACTAAATAATTCGCATTTTATGTTTAAAAATTGAGata                             <  1:618975/117‑1 (MQ=255)                                                                                         tAAGAGCGCTTGCGTTCGAGCTACGGCAATCGAGTCAACGGCCAGGGAAAGATCATCAAAAGTTCACTTTTTGCACTAAATAATTCGCATTTTATGTTTAAAAATTGAGata                             >  2:474828/1‑112 (MQ=255)                                                                                         tAAGAGCGCTTGCGTTCGAGCTACGGCAATCGAGTCAACGGCCAGGGAAAGATCATCAAAAGTTCACTTTTTGCACTAAATAATTCGCATTTTATGTTTAAAAATTGAGata                             <  1:474828/112‑1 (MQ=255)                                                                                         tAAGAGCGCTTGCGTTCGAGCTACGGCAATCGAGTCAACGGCCAGGGAAAGATCATCAAAAGTTCACTTTTTGCACTAAATAATTCGCATTTTATGTTTAAAAATTGAGATATTCCTTATTACCTGAAGCTGTTTTTTa  <  2:195299/139‑1 (MQ=255)                                                                                         tAAGAGCGCTTGCGTTCGAGCTACGGCAATCGAGTCAACGGCCAGGGAAAGATCATCAAAAGTTCACTTTTTGCACTAAATAATTCGCATTTTATGTTTAAAAATTGAGATATTCCTTATTACCTGAAGCTGTTTTTTa  <  1:645402/139‑1 (MQ=255)                                                                                           agagCGCTTGCGTTCGAGCTACGGCAATCGAGTCAACGGCCAGGGAAAGATCATCAAAAGTTCACTTTTTGCACTAAATAATTCGCATTTTATGTTTAAAAATTgag                                >  2:124311/1‑107 (MQ=255)                                                                                           agagCGCTTGCGTTCGAGCTACGGCAATCGAGTCAACGGCCAGGGAAAGATCATCAAAAGTTCACTTTTTGCACTAAATAATTCGCATTTTATGTTTAAAAATTgag                                <  1:124311/107‑1 (MQ=255)                                                                                               |                                                                                                                                     TCGGAGCATTTATACTCGGGCAAATATAAATTTGCAAAAAATGCGAAGATAATGTTTTTAATATTTTGAACAATATGGCATTGAATGATAAGAGTGCTTGCGTTCGAGCTACGGCAATCGAGTCAACGGCCAGGGAAAGATCATCAAAAGTTCACTTTTTGCACTAAATAATTCGCATTTTATGTTTAAAAATTGAGATATTCCTTATTACCTGAAGCTGTTTTTTA  >  NC_000913/3624104‑3624330

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