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	32 [0]	[0] 33	[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)	29 (0.600)	25/262	1.2	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

ACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATTCCCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCAGATCGCAAACCGCTTCGCCGTTCACAATTCCGACAGAAACTGCGGCTACCATCCCTT  >  NC_000913/3815799‑3816083                                                                                                                                               |                                                                                                                                               aaacctgcgtgctaagcgcaaattcaataatctggcgctgatatggtttcatgccttcgctcctcatcttacttttctacagacaaaaaaaaggcgactcatcagtcgccttaaaaatcagtttgccagcgccgccttctgcc                                                                                                                                                <  2:167364‑M2/1‑1 (MQ=255)     ctgcttgctaagcgcaaattcaataaactggcgctgatatggtttcatgccttcgctcctcatcttacttttctacagacaaaaaaaaggcgactcatcagtcgccttaaaaatcagtttgccagcgccgccttctgcCGTcc                                                                                                                                            >  2:928633‑M2/139‑143 (MQ=255)             taagcgcaaattcaataaactggcgctgatatggtttcatgccttcgctcctcatcttccttttctacagacaaaaaaaaggcgactcatcagtcgccttaaaaatcaatttaccagcgccgccttctgcCGTCCACTGCACt                                                                                                                                    >  1:480827‑M2/131‑143 (MQ=255)                 cgaaaattcaataaactggcgctgatatggtttcatgccttcgctcctcatcttacttttctacagacaaaaaaaaggcgactcatcagtcgccttaaaaatcagtttgccagcgccgccttctgcCGTCCCCTGCACTTCaa                                                                                                                                <  2:291457‑M2/17‑1 (MQ=255)                  gcaaattcaataagctggcgctgatatggtttcatgccttcgctcctcatcttacttttctacagacaaaaaaaaggcgactcatcagtcgccttaaaaatcagtttgccagcgccgccttctgcCGTCCCCTGCACTTCAat                                                                                                                               <  1:560531‑M2/18‑1 (MQ=255)                     aattcaataaactggcgctgatatggtttcatgccttcgctcctcatcttacttttctacagacaaaaaaaaggcgactcatcagtcgccttaaaaatcagttagccagcgccgccttctgcCGTCCCCTGCACTTCAatgat                                                                                                                            >  2:911256‑M2/123‑143 (MQ=255)                          aataaactggcgctgatatggtttcatgccttcgctcctcatcttacttttctacagacaaaaaaaaggcgactcatcagtcgccttaaaaatcagtttgccagcgccgccttctgcCGTCCCCTGCACTTCAATGATGCGcc                                                                                                                       >  1:438777‑M2/118‑143 (MQ=255)                          aataaactggcgctgatatggtttcatgccttcgctcctcaacttacttttctacagacaaaaaaaaggcgactcatcagtcgccttaaaaatcagtttgccagcgccgccttctgcCGTCCCCTGCACTTCAATGATGCGcc                                                                                                                       <  1:202467‑M2/26‑1 (MQ=255)                                ctggcgctgatatggtttcatgccttcgctcctcatcttacttttctacagacaaaaaaaaggcgactcatcagtcgccttaaaaatcagtttgccagcgccgccttctgcCGTCCCCTGCACTTCAATGATGCGCCCGTCtt                                                                                                                 <  2:530269‑M2/32‑1 (MQ=255)                                 tggcgctgatatggtttcatgccttcgctcctcatcttacttttctacagacaaaaaaaaggcgactcatcagtcgccttaaaaatcagtttgccagcgccgccttctgcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTc                                                                                                                <  1:334998‑M2/33‑1 (MQ=255)                                  ggcgctgatatggtttcatgccttcgctcctcatcttacttttctacagacaaacaaaaggcgactcatcagtcgccttaaaaatcagtttgccagcgccgccttctgcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCg                                                                                                               <  2:215078‑M2/34‑1 (MQ=255)                                       tgatatggtttcatgccttcgctcctcatcttacttttctacagacaaaaaaaaggcgactcatcagtcgccttaaaaatcagtttgccagcgccgccttctgcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGtcat                                                                                                          <  1:48481‑M2/39‑1 (MQ=255)                                       tgatatggtttcatgccttcgctcctcatcttacttttctacagacaaaaaaaaggcgactcatcagtcgccttaaaaatcagtttgccagcgccgccttctgcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGtcat                                                                                                          <  2:239581‑M2/39‑1 (MQ=255)                                                      ccttcgctcctcagcttacttttctgccgacaaaaaaaagccgacgcgtctgtcgccttaaaaatcagtttgccagcgccgccttctgcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTc                                                                                           <  1:353887‑M2/54‑1 (MQ=255)                                                        ttcgctcctcatcttacttttctacagacaaaaaaaaggcgactcatcagtcgccttaaaaatcagtttgccagcgccgccttctgcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGCTCATGTCgg                                                                                         >  2:545328‑M2/88‑143 (MQ=255)                                                         tcgctcctcatcttacttttctacagacaaaaaaaaggcgactcatcagtcgccttaaaaatcagtttgccagcgccgccttctgcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCCTCACTACGTTCATGTCGGt                                                                                        >  2:889611‑M2/87‑143 (MQ=255)                                                                     tttcttttctacagacaaaaaaaaggcgactcatcagtcgccttaaaaatcagtttgccagcgccgccttctgcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCaga                                                                            <  2:361134‑M2/69‑1 (MQ=255)                                                                     ttacgtttctacagaccataaaaaggcgactcatcagtcgccttaaaactcagtttgccagcgccgcctgctgcCGTCCCCTGCACTTCAATGATGCGCCCGCCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCaga                                                                            <  2:253979‑M2/69‑1 (MQ=255)                                                                                    caaaaaaaaggcgactcatcagtcgccttaaaaatcagtttgccagcgccgccttctgcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTc                                                             <  1:748837‑M2/84‑1 (MQ=255)                                                                                                 actcatcagtcgccttaaaaatcagtttgccagcgccgccttctgcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTCAAACGTATTCCAGATCTcacacg                                                >  1:668173‑M2/47‑142 (MQ=255)                                                                                                     atcagtcgccttaaaaatcagtttgccagcgccgccttctgcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGCGCGTTCAACGTATTCCAGATCGCAAACCGCtt                                            >  2:462391‑M2/43‑143 (MQ=255)                                                                                                     atcagtcgccttaaaaatcagtttgccagcgccgccttctgcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCAGATCGCAAACCGCtt                                            >  2:796998‑M2/43‑143 (MQ=255)                                                                                                               ttaaaagtcagtttgccagcgccgccttctgcCGCCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCAGATCGCAAACCGCTTCGCCGTTcac                                  <  1:754112‑M2/111‑1 (MQ=255)                                                                                                                    aatcagtttgccagcgccgccttctgcCGCCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACGACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCAGATCGCAAACCGCTTCGCCGTTCACAATTc                             <  2:438777‑M2/116‑1 (MQ=255)                                                                                                                     atcagtttgccagcgccgccttctgcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGCGTCTTCAACGTATTCCAGCTCGCACACCGCTTCGCCGTTCACAATTcc                            >  1:198203‑M2/27‑143 (MQ=255)                                                                                                                          tttgccagcgccgccttctgcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCAGATCGCAAACCGCTTCGCCATTCACAATTCCGACAg                       >  1:119676‑M2/22‑143 (MQ=255)                                                                                                                                agcgccgccttctgcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCAGATCGCAAACCGCTTCGCCGTTCACAATTCCGACAGAAACTg                 >  2:345301‑M2/16‑143 (MQ=255)                                                                                                                                agcgccgccttctgcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCAGATCGCAAACCGCTTCGCCGTTCACAATTCCGACAGAAACTg                 >  1:566461‑M2/16‑143 (MQ=255)                                                                                                                                 gcgccgccttctgcCGTCCCCTGCACTGCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCAGATCGCAAACCGCTTCGCCGTTCACAATTCCGACAGAAACTGc                <  2:480827‑M2/129‑1 (MQ=255)                                                                                                                                       ccttctgcCGTCCCCTGCACTTCAAGGATGCGCCCGTCTGCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCAGATCGCAAACCGCTTCGCCGTTCACAATTCCGACAGAAACTGCGGCTAc          <  1:789711‑M2/135‑1 (MQ=255)                                                                                                                                         ttctgcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCAGATCGCAAACCGCTTCGCCGTTCACAATTCCGACAGAAACTGCGGCTACCa        <  1:796998‑M2/137‑1 (MQ=255)                                                                                                                                           ctgcCGTCCCCTGCTCCTCCCCGATGCGCCCGTCTTCGGTCAGCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCAGATCGCAAACCGCTTCGCCGTTCACAATTCCGACAGAAACTGCGGCTACCATc      <  2:452243‑M2/139‑1 (MQ=255)                                                                                                                                               cGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCACAGTCTTCAACGTATTCCAGATCGCAAACCGCTTCGCCGTTCACAATTCCGACAGAAACTGCGGCTACCCTCCCtt  >  1:304906/1‑143 (MQ=255)                                                                                                                                               |                                                                                                                                               ACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATTCCCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCAGATCGCAAACCGCTTCGCCGTTCACAATTCCGACAGAAACTGCGGCTACCATCCCTT  >  NC_000913/3815799‑3816083

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