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%	119.2 / NA	38	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 (19/19);  total (19/19)

ATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATTCCCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACG  >  NC_000913/3815743‑3815986                                                                                                                         |                                                                                                                            aTAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCGGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGTTACAATGgatt                                                                                                                  >  2:1087347/1‑132 (MQ=255) aTAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCGGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGTTACAATGgatt                                                                                                                  <  1:1087347/132‑1 (MQ=255)         gcgcTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCGGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGTTACAATGGATTCGATTCCCCTCGGGc                                                                                                   <  2:714498/139‑1 (MQ=255)         gcgcTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCGGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGTTACAATGGATTCGATTCCCCTCGGGc                                                                                                   >  2:871349/1‑139 (MQ=255)             tGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCGGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGTTACAATGGATTCGATTCCCCTCGGGCCaga                                                                                               <  2:945088/139‑1 (MQ=255)                tatGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCGGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGTTACAATGGATTCGATTCCCCTCGGGCCAGAGcc                                                                                            <  2:945851/139‑1 (MQ=255)                      ttCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCGGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGTTACAATGGATTCGATTCCCCTCGGGCCAGAGCCAACAAg                                                                                      <  2:555050/139‑1 (MQ=255)                         aTGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCGGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGTTACAATGgattcga                                                                                                               <  2:436950/111‑1 (MQ=255)                         aTGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCGGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGTTACAATGgattcga                                                                                                               >  1:436950/1‑111 (MQ=255)                         aTGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCGGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGTTACAATGGATTCGATTCCCCTCGGGCCAGAGCCAACAAGATg                                                                                   >  1:936674/1‑139 (MQ=255)                             cTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCGGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGTTACAATGGATTCGATTCCCCTCGGGCCAGAGCCAACAAGATGAGTa                                                                               >  1:97913/1‑139 (MQ=255)                                cGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCAACTCATCGGTCGCATTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGTTACAATGGATTCGATTCCCCTCGGGCCAGAGCCAACAAGATGAGTAGct                                                                            >  1:908061/1‑139 (MQ=255)                                             aCTTTTCTACAGACAAAAAAAAGGCGACTCATCGGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGTTACAATGGATTCGATTCCCCt                                                                                                        >  1:620316/1‑98 (MQ=255)                                             aCTTTTCTACAGACAAAAAAAAGGCGACTCATCGGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGTTACAATGGATTCGATTCCCCt                                                                                                        <  2:620316/98‑1 (MQ=255)                                                 ttCTACAGACAAAAAAAAGGCGACTCATCGGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGTTACAATGGATTCGATTCCCCTCGGGc                                                                                                   >  2:126/1‑99 (MQ=255)                                                 ttCTACAGACAAAAAAAAGGCGACTCATCGGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGTTACAATGGATTCGATTCCCCTCGGGc                                                                                                   <  1:126/99‑1 (MQ=255)                                                    tACAGACAAAAAAAAGGCGACTCATCGGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGTTACAATGGATTCGATTCCCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCg                                                        >  2:727095/1‑139 (MQ=255)                                                      cAGACAAAAAAAAGGCGACTCATCGGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGTTACAATGGATTCGATTCCCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGcc                                                      <  2:936674/139‑1 (MQ=255)                                                                   ggCGACTCATCGGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGTTACAATGGATTCGATTCCCCt                                                                                                        >  1:442420/1‑76 (MQ=255)                                                                   ggCGACTCATCGGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGTTACAATGGATTCGATTCCCCt                                                                                                        <  2:442420/76‑1 (MQ=255)                                                                   ggCGACTCATCGGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGTTACAATGGATTCGATTCCCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGc                                                           <  1:1303855/121‑1 (MQ=255)                                                                   ggCGACTCATCGGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGTTACAATGGATTCGATTCCCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGc                                                           >  2:1303855/1‑121 (MQ=255)                                                                   ggCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGTTACAATGGATTCGATTCCCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTcccc                                         <  1:871349/139‑1 (MQ=255)                                                                                   ccTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGTTACAATGGATTCGATTCCCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCAt                                                                      <  1:221118/94‑1 (MQ=255)                                                                                   ccTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGTTACAATGGATTCGATTCCCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCAt                                                                      >  2:221118/1‑94 (MQ=255)                                                                                   ccTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGTTACAATGGATTCGATTCCCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATgc                         >  1:1084555/1‑139 (MQ=255)                                                                                     ttAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGTTACAATGGATTCGATTCCCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATgcgc                       >  1:1046890/1‑139 (MQ=255)                                                                                           aTCAGTTTGCCAGCGCCGCCTTCTGCGTCGTTACAATGGATTCGATTCCCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCt                 <  2:97913/139‑1 (MQ=255)                                                                                            tCAGTTTGCCAGCGCCGCCTTCTGCGTCGTTACAATGGATTCGATTCCCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCtt                <  1:727095/139‑1 (MQ=255)                                                                                                  tGCCAGCGCCGCCTTCTGCGTCGTTACAATGGATTCGATTCCCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTg                                                 >  1:948272/1‑100 (MQ=255)                                                                                                  tGCCAGCGCCGCCTTCTGCGTCGTTACAATGGATTCGATTCCCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTg                                                 <  2:948272/100‑1 (MQ=255)                                                                                                  tGCCAGCGCCGCCTTCTGCGTCGTTACAATGGATTCGATTCCCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGtca          <  1:1045559/139‑1 (MQ=255)                                                                                                   gCCAGCGCCGCCTTCTGCGTCGTTACAATGGATTCGATTCCCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTccc                                          <  1:733815/106‑1 (MQ=255)                                                                                                   gCCAGCGCCGCCTTCTGCGTCGTTACAATGGATTCGATTCCCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTccc                                          >  2:733815/1‑106 (MQ=255)                                                                                                        cgccgcCTTCTGCGTCGTTACAATGGATTCGATTCCCCTCGGGCCAGAGcc                                                                                            <  1:863359/51‑1 (MQ=255)                                                                                                        cgccgcCTTCTGCGTCGTTACAATGGATTCGATTCCCCTCGGGCCAGAGcc                                                                                            >  2:863359/1‑51 (MQ=255)                                                                                                          ccgccTTCTGCGTCGTTACAATGGATTCGATTCCCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACg  >  1:273343/1‑139 (MQ=255)                                                                                                          ccgccTTCTGCGTCGTTACAATGGATTCGATTCCCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACg  >  1:1319603/1‑139 (MQ=255)                                                                                                                    cgtcgtTACAATGGATTCGATTCCCCTCGGGCCAGAGCCAACAAGATGAGTAg                                                                              >  2:1148415/1‑53 (MQ=255)                                                                                                                    cgtcgtTACAATGGATTCGATTCCCCTCGGGCCAGAGCCAACAAGATGAGTAg                                                                              <  1:1148415/53‑1 (MQ=255)                                                                                                                         |                                                                                                                            ATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATTCCCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACG  >  NC_000913/3815743‑3815986

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