Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,997,779	T→C	pseudogene (210/300 nt)	ygeP ←	protein YgeP

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,997,779	0	T	C	100.0%	66.4 / NA	20	pseudogene (210/300 nt)	ygeP	protein YgeP
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (12/8);  total (12/8)

TTCTTCTGGTGTTTATCTTTGAAAGATTTGTAGGATAAATAACATTAAGTAAAGCCGCATCGGCTTCTCTGTAATGTTTATCTTTTTCCAGATAACCTGGGGATTCTCCCCCTAAATCAAGCTGTAACGTTATTCTTTCATAACCTTTATCTA  >  NC_000913/2997703‑2997855                                                                             |                                                                             ttcttcTGGTGTTTATCTTTGAAAGATTTGTAGGATAAATAACATTAAGTAAAGCCGCATCGGCTTCTCTGTAATGctt                                                                            <  1:547008/79‑1 (MQ=255)    ttctGGTGTTTATCTTTGAAAGATTTGTAGGATAAATAACATTAAGTAAAGCCGCATCGGCTTCTCTGTAATGCTTATCtt                                                                       >  1:547028/1‑81 (MQ=255)     tctGGTGTTTATCTTTGAAAGATTTGTAGGATAAATAACATTAAGTAAAGCCGCATCGGCTTCTCTGTAATGCTTATCttt                                                                      >  2:547043/1‑81 (MQ=255)     tctGGTGTTTATCTTTGAAAGATTTGTAGGATAAATAACATTAAGTAAAGCCGCATCGGCTTCTCTGTAATGCTTATCttt                                                                      >  2:547042/1‑81 (MQ=255)      ctGGTGTTTATCTTTGAAAGATTTGTAGGATAAATAACATTAAGTAAAGCCGCATCGGCTTCTCTGTAATGCTTATCttt                                                                      <  1:547009/80‑1 (MQ=255)         gtgtTTATCTTTGAAAGATTTGTAGGATAAATAACATTAAGTAAAGCCGCATCGGCTTCTCTGTAATGCTTATCTTTTTcc                                                                  <  2:547010/81‑1 (MQ=255)                 cTTTGAAAGATTTGTAGGATAAATAACATTAAGTAAAGCCGCATCGGCTTCTCTGTAATGCTTATCTTTTTCCAGATAAcc                                                          >  1:547049/1‑81 (MQ=255)                           tttGTAGGATAAATAACATTAAGTAAAGCCGCATCGGCTTCTCTGTAATGCTTATCTTTTTCCAGATAACCTGGGGATtct                                                >  2:547040/1‑81 (MQ=255)                             tGTAGGATAAATAACATTAAGTAAAGCCGCATCGGCTTCTCTGTAATGCTTATCTTTTTCCAGATAACCTGGGGATTCTcc                                              >  2:547039/1‑81 (MQ=255)                                 ggATAAATAACATTAAGTAAAGCCGCATCGGCTTCTCTGTAATGCTTATCTTTTTCCAGATAACCTGGGGATTCTccccc                                           >  1:547031/1‑80 (MQ=255)                                                  tAAAGCCGCATCGGCTTCTCTGTAATGCTTATCTTTTTCCAGATAACCTGGGGATTCTCCCCCTAAATCAAGCTGTAACGt                         <  2:547011/81‑1 (MQ=255)                                                  tAAAGCCGCATCGGCTTCTCTGTAATGCTTATCTTTTTCCAGATAACCTGGGGATTCTCCCCCTAAATCAAGCTGTAACGt                         <  2:547012/81‑1 (MQ=255)                                                     aGCCGCATCGGCTTCTCTGTAATGCTTATCTTTTTCCAGATAACCTGGGGATTCTCCCCCTAAATCAAGCTGTAACGTTAt                      >  1:547034/1‑81 (MQ=255)                                                     aGCCGCATCGGCTTCTCTGTAATGCTTATCTTTTTCCAGATAACCTGGGGATTCTCCCCCTAAATCAAGCTGTAACGTTAt                      >  2:547048/1‑81 (MQ=255)                                                           aTCGGCTTCTCTGTAATGCTTATCTTTTTCCAGATAACCTGGGGATTCTCCCCCTAAATCAAGCTGTAACGTTATTCTTTc                <  2:547013/81‑1 (MQ=255)                                                               gCTTCTCTGTAATGCTTATCTTTTTCCAGATAACCTGGGGATTCTCCCCCTAAATCAAGCTGTAACGTTATTCTTTCATaa            >  1:547052/1‑81 (MQ=255)                                                                   ctctGTAATGCTTATCTTTTTCCAGATAACCTGGGGATTCTCCCCCTAAATCAAGCTGTAACGTTATTCTTTCATAACCt         <  2:547014/80‑1 (MQ=255)                                                                     ctGTAATGCTTATCTTTTTCCAGATAACCTGGGGATTCTCCCCCTAAATCAAGCTGTAACGTTATTCTTTCATAACCTTTa      >  1:547050/1‑81 (MQ=255)                                                                      tGTAATGCTTATCTTTTTCCAGATAACCTGGGGATTCTCCCCCTAAATCAAGCTGTAACGTTATTCTTTCATAACCttt       <  2:547015/79‑1 (MQ=255)                                                                        tAATGCTTATCTTTTTCCAGATAACCTGGGGATTCTCCCCCTAAATCAAGCTGTAACGTTATTCTTTCATAACCTTTa      <  2:547016/78‑1 (MQ=255)                                                                         aaTGCTTATCTTTTTCCAGATAACCTGGGGATTCTCCCCCTAAATCAAGCTGTAACGTTATTCTTTCATAACCTTTATCTa  >  1:547044/1‑81 (MQ=255)                                                                             |                                                                             TTCTTCTGGTGTTTATCTTTGAAAGATTTGTAGGATAAATAACATTAAGTAAAGCCGCATCGGCTTCTCTGTAATGTTTATCTTTTTCCAGATAACCTGGGGATTCTCCCCCTAAATCAAGCTGTAACGTTATTCTTTCATAACCTTTATCTA  >  NC_000913/2997703‑2997855

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