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%	124.6 / 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 (20/20);  total (20/20)

GCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATTCCCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTC  >  NC_000913/3815734‑3815994                                                                                                                                  |                                                                                                                                    gcAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGTGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGTTACAATGga                                                                                                                            >  2:304235/1‑139 (MQ=255)                   cgcTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGTGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGTTACAATGGATTCGATTCCCCt                                                                                                                <  2:826971/133‑1 (MQ=255)                   cgcTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGTGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGTTACAATGGATTCGATTCCCCt                                                                                                                >  1:826971/1‑133 (MQ=255)                                  aTGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGTGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGTTACAATGGATTCGATTCCCCTCGGGCCAGAGCCAACAAGATg                                                                                           >  1:358446/1‑139 (MQ=255)                                  aTGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGTGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGTTACAATGGATTCGATTCCCCTCGGGCCAGAGCCAACAAGATg                                                                                           >  2:246332/1‑139 (MQ=255)                                   tGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGTGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGTTACAATGGATTCGATTCCCCTCGGGCCAGAGCCAACAAg                                                                                              <  1:463178/135‑1 (MQ=255)                                   tGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGTGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGTTACAATGGATTCGATTCCCCTCGGGCCAGAGCCAACAAg                                                                                              >  2:463178/1‑135 (MQ=255)                                       ttCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGTGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGTTACAATGGATTCGATTCCCCTCggg                                                                                                            <  1:189218/117‑1 (MQ=255)                                       ttCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGTGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGTTACAATGGATTCGATTCCCCTCggg                                                                                                            >  2:189218/1‑117 (MQ=255)                                       ttCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGTGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGTTACAATGGATTCGATTCCCCTCGGGCCAGAGcc                                                                                                    <  1:402948/125‑1 (MQ=255)                                       ttCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGTGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGTTACAATGGATTCGATTCCCCTCGGGCCAGAGcc                                                                                                    >  2:402948/1‑125 (MQ=255)                                       ttCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGTGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGTTACAATGGATTCGATTCCCCTCGGGCCAGAGCCAACAAGATGAGTAg                                                                                      <  1:246332/139‑1 (MQ=255)                                       ttCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGTGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGTTACAATGGATTCGATTCCCCTCGGGCCAGAGCCAACAAGATGAGTAg                                                                                      <  2:526732/139‑1 (MQ=255)                                         cGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGTGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGTTACAATGGATTCGATTCCCCTCGGGCCAGAGCCAACAAGATGAGTAGct                                                                                    >  1:198890/1‑139 (MQ=255)                                          gCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGTGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGTTACAATGGATTCGATTCCCCt                                                                                                                >  1:481673/1‑110 (MQ=255)                                          gCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGTGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGTTACAATGGATTCGATTCCCCt                                                                                                                <  2:481673/110‑1 (MQ=255)                                            tcctcATCTTACTTTTCTACAGACAAAAAAAAGGTGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGTTACAATGGATTCGATTCCCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCtt                                                                                 <  1:645089/139‑1 (MQ=255)                                                   cTTACTTTTCTACAGACAAAAAAAAGGTGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGTTACAATGGATTCGATTCCCCTCGGGCCAGAGCCAACAAGATGAGt                                                                                        <  1:580623/125‑1 (MQ=255)                                                   cTTACTTTTCTACAGACAAAAAAAAGGTGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGTTACAATGGATTCGATTCCCCTCGGGCCAGAGCCAACAAGATGAGt                                                                                        >  2:580623/1‑125 (MQ=255)                                                          ttCTACAGACAAAAAAAAGGTGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGTTACAATGGATTCGATTCCCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGc                                                                   >  1:675266/1‑139 (MQ=255)                                                               cAGACAAAAAAAAGGTGACTCATCAGTCGCCTTTAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGTTACAATGGATTCGATTCCCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGcc                                                              <  2:198890/139‑1 (MQ=255)                                                                    aaaaaaaaGGTGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGTTACAATGGATTCGATTCCCCTCGGGCCAGAGCCAACAAGATGAGTAg                                                                                      <  1:513891/110‑1 (MQ=255)                                                                    aaaaaaaaGGTGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGTTACAATGGATTCGATTCCCCTCGGGCCAGAGCCAACAAGATGAGTAg                                                                                      >  2:513891/1‑110 (MQ=255)                                                                    aaaaaaaaGGTGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGTTACAATGGATTCGATTCCCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGcc                                                              <  2:276584/134‑1 (MQ=255)                                                                    aaaaaaaaGGTGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGTTACAATGGATTCGATTCCCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGcc                                                              >  1:276584/1‑134 (MQ=255)                                                                            ggtgACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGTTACAATGGATTCGATTCCCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTcccc                                                 <  2:444618/136‑1 (MQ=255)                                                                              tggttcatcaGTCGCCTTAAAAAACAGTTTGCCAGCGCCGCCTTCTGCGTCGTTACAATGGATTCGATTCCCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTg                                               <  1:741253/135‑1 (MQ=255)                                                                              tgACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGTTACAATGGATTCGATTCCCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTg                                               <  2:358446/138‑1 (MQ=255)                                                                              tgACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGTTACAATGGATTCGATTCCCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTg                                               <  2:788948/138‑1 (MQ=255)                                                                                aCTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGTTACAATGGATTCGATTCCCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTccc                                                  >  1:90402/1‑134 (MQ=255)                                                                                aCTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGTTACAATGGATTCGATTCCCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTccc                                                  <  2:90402/134‑1 (MQ=255)                                                                                           gCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGTTACAATGGATTCGATTCCCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGa                                                                        >  2:638656/1‑101 (MQ=255)                                                                                           gCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGTTACAATGGATTCGATTCCCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGa                                                                        <  1:638656/101‑1 (MQ=255)                                                                                              ttAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGTTACAATGGATTCGATTCCCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATgcgc                               >  1:278085/1‑139 (MQ=255)                                                                                                         tttGCCAGCGCCGCCTTCTGCGTCGTTACAATGGATTCGATTCCCCTCGGGCCAGAGCCAACAAGATGAg                                                                                         >  1:222207/1‑70 (MQ=255)                                                                                                         tttGCCAGCGCCGCCTTCTGCGTCGTTACAATGGATTCGATTCCCCTCGGGCCAGAGCCAACAAGATGAg                                                                                         <  2:222207/70‑1 (MQ=255)                                                                                                         tttGCCAGCGCCGCCTTCTGCGTCGTTACAATGGATTCGATTCCCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGt                    >  2:59654/1‑139 (MQ=255)                                                                                                                  gccgccTTCTGCGTCGTTACAATGGATTCGATTCCCCTCGGGCCAGAGCCAACAAGATGAGTAg                                                                                      <  1:445953/64‑1 (MQ=255)                                                                                                                  gccgccTTCTGCGTCGTTACAATGGATTCAATTCCCCTCGGGCCAGAGCCAACAAGATGAGTAg                                                                                      >  2:445953/1‑64 (MQ=255)                                                                                                                           tGCGTCGTTACAATGGATTCGATTCCCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTc  >  2:679947/1‑139 (MQ=255)                                                                                                                                  |                                                                                                                                    GCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATTCCCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTC  >  NC_000913/3815734‑3815994

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