Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,815,863	C→T	pseudogene (19/48 nt)	rph ←	ribonuclease PH (defective),enzyme, Degradation of RNA, RNase PH

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,815,863	0	C	T	100.0%	126.8 / NA	40	pseudogene (19/48 nt)	rph	ribonuclease PH (defective),enzyme, Degradation of RNA, RNase PH
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (18/22);  total (18/22)

CAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATTCCCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCG  >  NC_000913/3815741‑3815995                                                                                                                           |                                                                                                                                     cAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCGGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGTTACAATGgattcgatt                                                                                                                      <  2:870462/139‑1 (MQ=255) cAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCGGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGTTACAATGgattcgatt                                                                                                                      <  2:168030/139‑1 (MQ=255)                 atatGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCGGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGTTACAATGGATTCGATTCCCCTCGGGCCAGAGc                                                                                                      >  1:139532/1‑139 (MQ=255)                   atGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCGGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGTTACAATGGATTCGATTCCCCTCGGGCCAGAGCca                                                                                                    <  1:560832/139‑1 (MQ=255)                           aTGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCGGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGTTACAATGGATTCGATTCCCCTCGGGc                                                                                                            <  2:207279/123‑1 (MQ=255)                           aTGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCGGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGTTACAATGGATTCGATTCCCCTCGGGc                                                                                                            >  1:207279/1‑123 (MQ=255)                                ttCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCGGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGTTACAATGGATTCGATTCCCCTCGGGCCAGAGCCAACAAGATGAGTAg                                                                                       <  2:766824/139‑1 (MQ=255)                                ttCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCGGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGTTACAATGGATTCGATTCCCCTCGGGCCAGAGCCAACAAGATGAGTAg                                                                                       <  2:1172168/139‑1 (MQ=255)                                   gCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCGGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGCCGTTACAATGGATTCGATTCCCCTCGGGCCAGAGCCAACAAGATGAGTAGctc                                                                                    >  2:1050368/1‑139 (MQ=255)                                       ctcATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCGGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGTTACAATGGATTCGATTCCCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCa                                                                                >  2:948121/1‑139 (MQ=255)                                                 ttttCTACAGACAAAAAAAAGGCGACTCATCGGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGTTACAATGGATTCGATTCCCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACg                                                                      >  2:619102/1‑139 (MQ=255)                                                   ttCTACAGACAAAGAAAATGCGTCCCATCGGTCGGCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGTTACAATGGATTCGATTCCCCCCGGGCCAGAGCCAACAAGATGAGTAGCTCTGCATGGGTGAACGGc                                                                    <  2:1145957/139‑1 (MQ=255)                                                    tCTACAGACAAAAAAAAGGCGACTCATCGGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGTTACAATGgattcgatt                                                                                                                      <  2:323797/88‑1 (MQ=255)                                                    tCTACAGACAAAAAAAAGGCGACTCATCGGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGTTACAATGgattcgatt                                                                                                                      >  1:323797/1‑88 (MQ=255)                                                     cTACAGACAAAAAAAAGGCGACTCATCGGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGTTACa                                                                                                                                  <  1:534776/75‑1 (MQ=255)                                                     cTACAGACAAAAAAAAGGCGACTCATCGGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGTTACa                                                                                                                                  >  2:534776/1‑75 (MQ=255)                                                        cAGACAAAAAAAAGGCGACTCATCGGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGTTACAATGGATTCGATTCCCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGcc                                                               <  1:1246071/139‑1 (MQ=255)                                                           aCAAAAAAAAGGCGACTCATCGGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGTTACAATGGATTCGATTCCCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTc                                                            >  1:102848/1‑139 (MQ=255)                                                           aCAAAAAAAAGGCGACTCATCGGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGTTACAATGGATTCGATTCCCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTc                                                            >  2:1191104/1‑139 (MQ=255)                                                                     ggCGACTCATCGGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGTTACAATGGATTCGATTCCCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTc                                                     <  1:1053243/136‑1 (MQ=255)                                                                     ggCGACTCATCGGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGTTACAATGGATTCGATTCCCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTc                                                     >  2:1053243/1‑136 (MQ=255)                                                                       cGACTCATCGGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGTTACAATGGATTCGATTCCCCt                                                                                                                 <  1:1048590/74‑1 (MQ=255)                                                                       cGACTCATCGGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGTTACAATGGATTCGATTCCCCt                                                                                                                 >  2:1048590/1‑74 (MQ=255)                                                                       cGACTCATCGGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGTTACAATGGATTCGATTCCCCTCGGGCCAGAGCca                                                                                                    <  1:484220/87‑1 (MQ=255)                                                                       cGACTCATCGGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGTTACAATGGATTCGATTCCCCTCGGGCCAGAGCca                                                                                                    >  2:484220/1‑87 (MQ=255)                                                                         aCTCATCGGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGTTACAATGGATTCGATTCCCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCtt                                                                                  >  2:1106083/1‑103 (MQ=255)                                                                         aCTCATCGGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGTTACAATGGATTCGATTCCCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCtt                                                                                  <  1:1106083/103‑1 (MQ=255)                                                                                     ccTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGTTACAATGGATTCGATTCCCCTCGGGCCAGAGCcaaca                                                                                                 <  1:59016/76‑1 (MQ=255)                                                                                     ccTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGTTACAATGGATTCGATTCCCCTCGGGCCAGAGCcaaca                                                                                                 >  2:59016/1‑76 (MQ=255)                                                                                      cTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGTTACAATGGATTCGATTCCCCt                                                                                                                 <  2:1187211/59‑1 (MQ=255)                                                                                      cTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGTTACAATGGATTCGATTCCCCt                                                                                                                 >  1:1187211/1‑59 (MQ=255)                                                                                                aGTTTGCCAGCGCCGCCTTCTGCGTCGTTACAATGGATTCGATTCCCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCg                       <  1:619102/139‑1 (MQ=255)                                                                                                  tttGCCAGCGCCGCCTTCTGCGTCGTTACAATGGATTCGATTCCCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTg                                                >  1:429257/1‑112 (MQ=255)                                                                                                  tttGCCAGCGCCGCCTTCTGCGTCGTTACAATGGATTCGATTCCCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTg                                                <  2:429257/112‑1 (MQ=255)                                                                                                    tGCCAGCGCCGCCTTCTGCGTCGTTACAATGGATTCGATTCCCCTCGGGCCAGAGcc                                                                                                     >  1:812464/1‑57 (MQ=255)                                                                                                    tGCCAGCGCCGCCTTCTGCGTCGTTACAATGGATTCGATTCCCCTCGGGCCAGAGcc                                                                                                     <  2:812464/57‑1 (MQ=255)                                                                                                    tGCCAGCGCCGCCTTCTGCGTCGTTACAATGGATTCGATTCCCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGtca                   <  2:86831/139‑1 (MQ=255)                                                                                                     gCCAGCGCCGCCTTCTGCGTCGTTACAATGGATTCGATTCCCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGtcat                  <  2:341131/139‑1 (MQ=255)                                                                                                          cgccgcCTTCTGCGTCGTTACAATGGATTCGATTCCCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTa             >  2:664293/1‑139 (MQ=255)                                                                                                                     gCGTCGTTACAATGGATTCGATTCCCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCg  <  2:139532/139‑1 (MQ=255)                                                                                                                           |                                                                                                                                     CAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATTCCCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCG  >  NC_000913/3815741‑3815995

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