Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
MC JC	NC_000913_3_hsa_tpiA	3,815,859	Δ82 bp		[rph]–[rph]	[rph], [rph]

Missing coverage evidence...
			seq id	start	end	size	←reads	reads→	gene	description
*	*	÷	NC_000913_3_hsa_tpiA	3815859	3815940	82	50 [0]	[0] 50	[rph]–[rph]	[rph],[rph]

New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913_3_hsa_tpiA	= 3815858	0 (0.000)	47 (0.840)	43/262	0.4	100%	pseudogene (24/48 nt)	rph	ribonuclease PH (defective);enzyme; Degradation of RNA; RNase PH
	?	NC_000913_3_hsa_tpiA	3815941 =	0 (0.000)					pseudogene (609/669 nt)	rph	ribonuclease PH (defective);enzyme; Degradation of RNA; RNase PH

TGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATTCCCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCG  >  NC_000913_3_hsa_tpiA/3815722‑3816004                                                                                                                                         |                                                                                                                                                   tgcttgctAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcac                                                                                                                                        >  1:197353‑M1/1‑137 (MQ=255)   cttgctAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcactt                                                                                                                                      >  2:690032‑M1/1‑135 (MQ=255)      gctAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaa                                                                                                                                   <  1:605883‑M1/149‑18 (MQ=255)       ctAATCGGAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCAATACATTGATCTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaat                                                                                                                                  <  1:244913‑M1/149‑19 (MQ=255)       ctAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCTCTTGAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaat                                                                                                                                  >  1:510992‑M1/1‑131 (MQ=255)        tAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatg                                                                                                                                 >  1:662047‑M1/1‑130 (MQ=255)              cAAATTCAATAAACTGGCACTGATATGGTTTCATGCCTTCGCTCCTCATCTTAATTTTCAACAAAAAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgc                                                                                                                           <  2:612470‑M1/149‑26 (MQ=255)               aaaTTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTCCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcc                                                                                                                          >  1:196399‑M1/1‑123 (MQ=255)                aaTTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgccc                                                                                                                         <  1:344025‑M1/149‑28 (MQ=255)                aaTTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgccc                                                                                                                         <  2:197353‑M1/149‑28 (MQ=255)                     aaTAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgctgcgcccgtctt                                                                                                                    >  1:744245‑M1/1‑117 (MQ=255)                      aTAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttc                                                                                                                   <  2:451228‑M1/149‑34 (MQ=255)                      aTAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttc                                                                                                                   >  1:165468‑M1/1‑116 (MQ=255)                      aTAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgagaccgtattc                                                                                                                   >  1:719597‑M1/1‑116 (MQ=255)                       tAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCAGCTTACTTTTCTACAGACAAAAAAAAAGCGACTCTTCACTCGCCTTAAAAATCCCTTTGCCAGCGCCGCCTTCTGCcctcccctgcacttcattgctgcgcccgccttcg                                                                                                                  >  2:721128‑M1/1‑115 (MQ=255)                       gacactgGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcg                                                                                                                  <  1:360456‑M1/146‑35 (MQ=255)                          aCTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtc                                                                                                               <  1:690032‑M1/149‑38 (MQ=255)                              ccgctGATATGGTTTCCTGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatca                                                                                                           <  1:319111‑M1/148‑42 (MQ=255)                                gcTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcact                                                                                                         <  1:868038‑M1/149‑44 (MQ=255)                                   gATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactacg                                                                                                      >  1:708589‑M1/1‑103 (MQ=255)                                     tatGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactacgtt                                                                                                    <  1:215206‑M1/149‑49 (MQ=255)                                        ggTTTCCTGCCTTCGCCCCTCATCTTACTTTTCCACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactacgttcat                                                                                                 <  1:1623‑M1/149‑52 (MQ=255)                                           ttCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactacgttcatgtc                                                                                              <  1:430570‑M1/149‑55 (MQ=255)                                                   ttCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactacgttcatgtcggtctctg                                                                                      <  2:16114‑M1/149‑63 (MQ=255)                                                       ctcctcATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactacgttcatgtcggtctctgcggc                                                                                  >  2:443900‑M1/1‑83 (MQ=255)                                                          ctcATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcactgcaatgctgcgctcgtcttcggtcatcactacgttcatgtcggtctctgcggcaga                                                                               <  1:941212‑M1/149‑70 (MQ=255)                                                           tcATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatcagcccgtcttcgggcatcactacggtcatgtcggtctatgtggcggcg                                                                              >  1:532410‑M1/1‑79 (MQ=255)                                                                ttACTTATCTACAGAAAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactacgttcatgtcggtctctgcggcagagtcttc                                                                         <  2:593003‑M1/149‑76 (MQ=255)                                                                   tttttcTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactacgttcatgtcggtctctgcggcagagtcttcaac                                                                      <  1:551602‑M1/148‑79 (MQ=255)                                                                          aCAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactacgttcatgtcggtctctgcggcagagtcttcaacgtattcc                                                               <  1:443900‑M1/149‑86 (MQ=255)                                                                            aGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactacgttcatgtcggtctctgcagcagagtcttccacgtaatccag                                                             >  2:363413‑M1/1‑62 (MQ=255)                                                                              aCAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactacgttcatgtcggtctctgcggcagagtcttcaacgtattccagat                                                           >  1:266765‑M1/1‑60 (MQ=255)                                                                                aaaaaaaaGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactacgttcatgtcggtctctgcggcagagtcttcaacgtattccagatcg                                                         >  2:482379‑M1/1‑58 (MQ=255)                                                                                aaaaaaaaGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactacgttcatgtcggtctctgcggcagagtcttcaacgtattccagatcg                                                         >  2:909272‑M1/1‑58 (MQ=255)                                                                                  aaaaaaGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactacgttcatgtcggtctctgcggcagagtcttcaacgtattccagatcgca                                                       >  2:917794‑M1/1‑56 (MQ=255)                                                                                      aaGGCGACTGATGAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactacgttcatgtcggtctctgcggcagagtcttcaacgtattccagagcgcaaacc                                                   >  1:1013409‑M1/1‑52 (MQ=255)                                                                                           gACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactacgttcatgtcggtctctgcggcagagtcttcaacgtattccagatcgcaaaccgcttc                                              <  1:917794‑M1/149‑103 (MQ=255)                                                                                                atcaGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCtgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactacgttcatgtcggtctctgcggcagagtcttcaacgtattccagatcgcaaaccgcttcgccgt                                         >  1:941628‑M1/1‑42 (MQ=255)                                                                                                atcaGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactacgttcatgtcggtctctgcggcagagtcttcaacgtattccagatcgcaaaccgcttcgccgt                                         <  2:624398‑M1/149‑108 (MQ=255)                                                                                                     tCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactacgttcatgtcggtctctgcggcagg                                                                               >  1:542349‑M1/1‑37 (MQ=255)                                                                                                     tCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactacgttcatgtcggtctctgcggcagg                                                                               <  2:542349‑M1/106‑70 (MQ=255)                                                                                                              aaaTCCGTTTGCCAGCTCTCTCTTCTGCcgtcccctgcacttcactgatgcgcccgtcttcggtcatcactacgttcatgtcggtctctgcggcagagtcttcaacgtattccagatcgcaaaccgcttcgccgttcacaattccgaca                           <  2:695050‑M1/149‑122 (MQ=255)                                                                                                                    gTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactacgttcatgtcggtctctgcggcagagtcttcaacgtattccagatcgcaaaccgcttcgccgttcacacttacgacagaacct                     >  2:198112‑M1/1‑22 (MQ=255)                                                                                                                      ttGCCATCGCCGCCTTCTGCtgtcccccgcacttcaatgatgcgcccgtcttcggtcatcactacgttcatgtcggtctctgcggcagagtcttcaacgtattccagatcgcaaaccgcttcgccgttcacaattccgacagaaactgc                   <  1:191745‑M1/149‑130 (MQ=255)                                                                                                                       tGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactacgttcatgtcggtctctgcggcagagtcttcaacgtattccagatcgcaaaccgcttcgccgttcacaattccgacagaaactgcg                  <  1:482379‑M1/149‑131 (MQ=255)                                                                                                                            gcgcCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactacgttcatgtcggtctctgcggcagagtcttcaacgtattccagatctcaaaccgcttcgccgttcacaattccgacagaaacgtaggctac             >  1:334051‑M1/1‑14 (MQ=255)                                                                                                                              gccgccTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactacgttcatgtcggtctctgcggcagagtcttcaacgtattccagatcgcaaaccgcttcgccgttcacaattccgacagaaactgcggctacca           >  1:796248‑M1/1‑12 (MQ=255)                                                                                                                                       tGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactacgttcatgtcggtctctgcggcagagtcttcaacgtattccagatcgcaaaccgcttcgccgttcacaattccgacagaaactgcggctaccatccctttca  >  2:38560‑M1/1‑3 (MQ=255)                                                                                                                                         ccgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactacgttcatgtcggtctctgcggcagagtcttcaacgtattccagatcgcaaaccgcttcgccgttcacaattccgacagaaactgcg                  >  2:710947‑M1/1‑1 (MQ=255)                                                                                                                                         ccgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactacgttcatgtcggtctctgcggcagagtcttcaacgtattccagatcgcaaaccgcttcgccgttcacaattccgacagaaactgcg                  <  1:710947‑M1/131‑131 (MQ=255)                                                                                                                                         |                                                                                                                                                   TGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATTCCCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCG  >  NC_000913_3_hsa_tpiA/3815722‑3816004

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 22 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑