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

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

New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913	= 3815858	0 (0.000)	41 (0.990)	33/264	0.3	100%	pseudogene (24/48 nt)	rph	ribonuclease PH (defective);enzyme; Degradation of RNA; RNase PH
	?	NC_000913	3815941 =	0 (0.000)					pseudogene (609/669 nt)	rph	ribonuclease PH (defective);enzyme; Degradation of RNA; RNase PH

AAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATTCCCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCAGATCGCAAACCGCTTCGCCGTTCACAATTCCGACAGAAACTGCGGCTACCATCCCTT  >  NC_000913/3815802‑3816083                                                                                                                                            |                                                                                                                                               cctgcttgctacgcgcaaattcaataacctggcgctgatatggtttcatgccttcgctcctcatcttacttttctacagacaaaaaaaaggcgactcatcagtcgccttaaaaatcagtttgccagcgccgccttctgcCGTc                                                                                                                                             <  1:271203‑M2/4‑1 (MQ=255)               gcaaattcaataaactggcgctgatatggtttcatgccttcgctcctcatcttacttttctacagacaaaaaaaaggcgactcatcagtcgccttaaaaatcagtttgccagcgccgccttctgcCGTCCCCTGCACTTCAat                                                                                                                               >  1:602725‑M2/126‑143 (MQ=255)               gcaaattcaataaactggcgctgatatggtttcatgccttcgctcctcatcttacttttctacagacaaaaaaaaggcgactcatcagtcgccttaaaaatcagtttgccagcgccgccttctgcCGTCCCCTGCACTTCAat                                                                                                                               >  1:25696‑M2/126‑143 (MQ=255)                      caagtgtctggcgcgggtatgggttcgtgccttcgctcctcatcttattgtactacagacaaaaaaaaggcgactcatcagtcgccttaaagatcagtttgccagcgccgccttctgcCGTCCCCTGCACTTCAATGATgcgc                                                                                                                        <  2:324198‑M2/25‑1 (MQ=255)                               ggcgctgttctgttttcgcgccgtcgccccccttctcccctttctacagacacaaaaaaggcgactcatcagtcgccttaaaaatcagtttgccagcgccgccttctgcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCg                                                                                                               <  2:692631‑M2/34‑1 (MQ=255)                               ggcgctgatatggtgtcatgcctgccctcctcttgtttgttctagccaaacaaaaaaaaggcgactcatcagtcgccttaaaaatcagtttgccagcgccgccttctgcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCg                                                                                                               <  1:795163‑M2/34‑1 (MQ=255)                                              gcgtgccttcgctcctgatcttccttttctacagacaaaaaaaaggcgactcatcagtcgccttaaaaatcagtttgccagcgccgccttctgcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTc                                                                                                <  1:734998‑M2/49‑1 (MQ=255)                                                                 cttacttttctacagacaaaaaaaaggcgactcatcagtcgccttaaaaatcagtttgccagcgccgccctctgaCGGTCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGCCTCTTCGGCag                                                                             >  1:371129‑M2/76‑143 (MQ=255)                                                                  ttacttttctacagacaaaaaaaaggcgactcatcagtcgccttaaaaatcagtttgccagcgccgacttctgcCGTCCCCTGCACTTCAATGATGCGCCCGGCTTCGGTCACCACTCCGTTCATGTCGGTCTCTGCGGcagg                                                                            >  2:377277‑M2/75‑142 (MQ=255)                                                                      ttttctacagacaaaaaaaaggcgactcatcagtcgccttaaacatcagttcgccagcgccgccttctgcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCt                                                                        <  1:454879‑M2/73‑1 (MQ=255)                                                                      ttttctacagacaaaaaaaaggcgactcatcagtcgccttaaaaatcagtttgccagcgccgccttctgcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCGTGTCGGTCGCTGCggcagagagc                                                                        >  1:522012‑M2/71‑140 (MQ=255)                                                                      ttttctacagacaaaaaaaaggcgactcatcagtcgccttaaaaatcagtttgccagcgccgccttctgcCGTCCCATGCACTTCAATGATGCGCCAGTCGTCGGTCATCCCTACGCTCATGTCGGTCTCTACGGCAGAGtcc                                                                        >  1:449494‑M2/71‑142 (MQ=255)                                                                          ccacagaccaagaaaaggcgtctcatcagtcgccttaaaaatcagtttgccagcgccgccttctgcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCaa                                                                    <  2:469081‑M2/77‑1 (MQ=255)                                                                           tacatacaaataaaaggcgtctcatcagtcgccttaaaaatcagtttgccagcgccgccttctgcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAAc                                                                   <  1:640371‑M2/78‑1 (MQ=255)                                                                           tacagacaaaaaaaaggcgactcatcagtcgccttaaaaatcagtttgccagcgccgccttctgcCGTCCCCTGCACTTCAATGATGCGCCCGTCCTCGGTCATCACTACGTTCATGTCGGTCCCTCCGGCACCGTCTccacc                                                                   >  2:495833‑M2/66‑141 (MQ=255)                                                                                acaaaaaaaaggcgactcatcagtcgccttaaaaatcagtttgccagcgccgccttctgcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTAtt                                                              <  1:93934‑M2/83‑1 (MQ=255)                                                                                acaaaaaaaaggcgactcatcagtcgccttaaaaatcagtttgccagcgccgccttctgcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACCTAtt                                                              >  1:773103‑M2/61‑143 (MQ=255)                                                                                 caaaaaaaaggcgactcatcagtcgccttaaaaatcagtttgccagcgccgccttctgcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTc                                                             <  2:25696‑M2/84‑1 (MQ=255)                                                                                      agcaggcgactcatcaggcgccttaaagatcagtttgccagcgccgccttctgcCGTCCCCTGCACTTCAATTATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCAGAt                                                        <  2:253149‑M2/89‑1 (MQ=255)                                                                                            cgactcatcagtcgccttaaaaatcagtttgccagcgccgccttctgcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCAGATCGCaaa                                                  >  2:67427‑M2/49‑143 (MQ=255)                                                                                            cgactcatcagtcgccttaaaaatcagtttgccagcgccgccttctgcCGTCCCCTGCACTTAAATGAGTCACCCGTCTTCCGTCATCACTACGTTCATGTCGGTATCTGCGGCAGCCTCCTCAACGGCTTCAACAACGCaaa                                                  >  2:58385‑M2/49‑143 (MQ=255)                                                                                             gactcatcagtcgccttaaaaatcagtttgccagcgccgccttctgcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGCGTCTTCCACGAGTGGCAGatcgcgcac                                                 >  2:427504‑M2/48‑139 (MQ=255)                                                                                             gactcatcagtcgccttaaaaatcagtttgccagcgccgccttctgcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGCGTCGACCTCAGAGTCTTCAACGTATTCCAGATCGCAAAc                                                 >  1:757780‑M2/48‑143 (MQ=255)                                                                                                  atcagtcgccttaaaaatcagtttgccagcgccgccttctgcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCAGATCGCAAACCGCtt                                            <  1:141185‑M2/101‑1 (MQ=255)                                                                                                  atcagtcgcctgaaaaatcagtttgccagcgccgccttctgcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCAGGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCAGATCGCAAACCGctc                                            >  2:360472‑M2/43‑142 (MQ=255)                                                                                                   tcagtcgccttaaaaatcagtttgccagcgccgccttctgcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCAGATCGCAAACCGCTTc                                           >  1:494260‑M2/42‑143 (MQ=255)                                                                                                   tcagtcgccttaaaaatcagtttgccagcgccgccttctgcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTACCTCAGTAATGGTCTTCAACGTATTCCGGGTCGCAAGCCGCTTc                                           >  1:544249‑M2/42‑143 (MQ=255)                                                                                                   gcagtcgccttaaaaatcagcttgccagcgccgcgtcctgcCGTCGCCTGCACGTCAATGATGCGCCCGTCTTCGGTCCTCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCAGATCGCAAACCGCTTc                                           <  1:334390‑M2/102‑1 (MQ=255)                                                                                                    cagtcgcattaaaaatcagtttgccagcgccgccttctgcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCAGATCGCAAACCGCTTCg                                          <  2:609989‑M2/103‑1 (MQ=255)                                                                                                      gtcgccttaaaaatcagtttgccagcgccgccttctgcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCAGATCGCAAACTGCTTCGcc                                        >  1:651760‑M2/39‑143 (MQ=255)                                                                                                          ccttaaaaatcagtttgccagcgccgccttctgcCGTCCCCTGCACTTCAATGCTGCGCCCGTCTTCGGTCATCACGACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACTTATTCCAGATCGCAAACCGCTTCGCCGTTc                                    <  2:107215‑M2/109‑1 (MQ=255)                                                                                                            ttaaaaatcagtttgccagcgccgccttctgcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCCCTGCGGCATAGTCGTCAACCGATTCCAGACCGCAAACCGCTTCGCCGTTcac                                  >  2:262592‑M2/33‑143 (MQ=255)                                                                                                             taaaaatcagtttgccagcgccgccttctgcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCAGATCGCAAACCGCTTCGCCGTTcaca                                 >  1:637117‑M2/32‑143 (MQ=255)                                                                                                              aaaaatcagtttgccagcgccgccttctgcCGTCCCCTGCACTTCAATGATGCGCGCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCAGATCACACACCCCTTCGACGTTCACaa                                >  2:301450‑M2/31‑143 (MQ=255)                                                                                                              aaaaatcagtttgccagcgccgccttctgcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGGCTCTGCGGCAGAGTCTTCAACGTATTCAAGATCGCAAACCGCTTCGCCGCTCACaa                                >  2:298261‑M2/31‑143 (MQ=255)                                                                                                                          ggcagcgccgccttctgcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCAGATCGCAAACCGCTTCGCCGTTCACAATTCCGACAGaaa                    <  2:221674‑M2/125‑1 (MQ=255)                                                                                                                          accagcgccgccttctgcCGTCCCCTGCACTTCAATGATGCGTCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCAGATCGCAAACCGCTTCGCCGTTCACAATTCCGACAGaaa                    <  2:103816‑M2/125‑1 (MQ=255)                                                                                                                            cagcgccgccttctgcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCAGATCGCAAACCGCTTCGCCGTTCACAATTCCGACAGAAACt                  >  1:526437‑M2/17‑143 (MQ=255)                                                                                                                             agcgccgccttctgcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCAGATCGCAAACCGCTTCGCCGTTCACAATTCCGACAGAAACTg                 >  1:686821‑M2/16‑143 (MQ=255)                                                                                                                              gcgccgccttctgcCTTCTCCTTCACTTCAATGATTCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTGCAACGTATTCCAGATCGCAAACCGCTTCGCCGTTCACAATTCCGACAGAAACTGc                <  1:241326‑M2/129‑1 (MQ=255)                                                                                                                                 ccgccttctgcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCAGATCGCAAACCGCTTCGCCGTTCACAATTCCGACAGAAACTGCGGc             <  1:550109‑M2/132‑1 (MQ=255)                                                                                                                                      ttctgcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCAGATCGCAAACCGCTTCGCCGTTCCCAATTCCGGCAGACACTGCGGCTACCa        >  1:135667‑M2/7‑143 (MQ=255)                                                                                                                                            cGTCCCCTGCACTTCAATGTTGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGATTTTTCACCGTATTCCGGTTCGCAAACCGCTTCGCCGTTCACAATTCCGACAGAAACTGCGGCTACCATCCCtt  <  1:410302/143‑1 (MQ=255)                                                                                                                                            cGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCAGATCGCAAACCGCTTCGCCGTTCACACTCCCGCCAGAAACTGCGCCTACCAGCCCtt  >  2:27838/1‑143 (MQ=255)                                                                                                                                            |                                                                                                                                               AAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATTCCCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCAGATCGCAAACCGCTTCGCCGTTCACAATTCCGACAGAAACTGCGGCTACCATCCCTT  >  NC_000913/3815802‑3816083

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