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	43 [0]	[0] 43	[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)	39 (1.140)	31/264	0.1	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

ACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATTCCCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCAGATCGCAAACCGCTTCGCCGTTCACAATTCCGACAGAAACTGCGGCTACCATCCCT  >  NC_000913/3815799‑3816082                                                                                                                                               |                                                                                                                                              acacctgcttgctaagcgcaaattcaataaactggcgctgatatggtttcatgccttcgctcctcatcttacttttctacagacaaaaaaaaggcgactcatcagtcgccttaaaaatcagtttgccagcgccgccttctgcc                                                                                                                                               <  1:576090‑M2/1‑1 (MQ=255)   tcctgcttgctaagcgcaaattcaataaactggcgcttatatggtttcatgccttcgctcctcatcttacttttctacagacaaaaaaaaggcgactcatcagtcgccttaaaaatcagtttgccagcgccgccttctgcCGt                                                                                                                                             <  2:125099‑M2/3‑1 (MQ=255)      tgcttgctaagcgcaaattcaataaactggcgctgatatggtttcatgccttcgctcctcatcttacttttctacagacaaaaaaaaggcgactcatcagtcgccttaaaaatcagtttgccagcgccgccttctgcCGTccc                                                                                                                                          <  1:423333‑M2/6‑1 (MQ=255)      tgcttgctaagcgcaaattcaataaactggcgctgatatggtttcatgccctcgctcctcttcttacttttctacagacaaaaaaaaggcgactcatcagtcgccttaaaaatcagtttgccagcgccgccttctgcCGTccc                                                                                                                                          <  1:383488‑M2/6‑1 (MQ=255)       gcttgctaagcgcaaattcaataaactggcgctgatatggtttcatgccttcgctcctcatcttacttttctacagacaaaaaaaaggcgactcatcagtcgccttaaaaatcagtttaccagcgccgccttctgcCGTcccc                                                                                                                                         >  2:421834‑M2/137‑143 (MQ=255)        cttgctaagcgcaaattcaataaactggcgctgatatggtttcatgccttcgctcctcatcttacttttctacagacaaaaaaaaggcgactcatcagtcgccttaaaaatcagtttgccagcgccgccttctgcCGTCCCCt                                                                                                                                        <  2:546495‑M2/8‑1 (MQ=255)           gctaagcgcaaattcaataaactggcgctgatatggtttcatgccttcgctcctcatcttacttttctacagacaaaaaaaaggcgactcatcagtcgccttaaaaatcagtttgccagcgccgccttctgcCGTCCCCTGCa                                                                                                                                     <  2:150470‑M2/11‑1 (MQ=255)             taagcgcaaattcaataaactggcgctgatatggtttcatgccttcgctcctcatcttacttttctacagacaaaaaaaaggcgactcatcagtcgccttaaaaatcagtttgccagcgccgccttctacCGTACCCTGCACt                                                                                                                                   >  2:565552‑M2/131‑143 (MQ=255)                  gcaaattcaataaactgtcgcttatatggtgtcatgccttcgctcctcatcttactattctacagacaaaaaaaaggcgactcatcagtcgccttaaaaatcagtttgccagcgccgccttctgcCGTCCCCTGCACTTCAat                                                                                                                              <  2:326053‑M2/18‑1 (MQ=255)                                  ggcgctgatatggtttcatgccttcgctcctcatcttacttttctacagacaaaaaaaaggcgactcatcagtcgccttaaaaatcagtttgccagcgccgccttctgcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCg                                                                                                              <  1:101001‑M2/34‑1 (MQ=255)                                  gccgctgatatggtttcatgccttcgctcctcatcttacttttctacagacaaaaaaaaggcgactcatcagtcgccttaaaaatcagtttgccagcgccgccttctgcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCg                                                                                                              <  1:202626‑M2/34‑1 (MQ=255)                                     gctgatatggtttcatgccttcgctcctcatcttacttttctacagacaaaaaaaaggcgactcatcagtcgccttaaaaatcagtttgccagcgccgccttctgcCGTCCCCTGCACTTCAATCATGCGCCCGTCTtcggcc                                                                                                           >  1:457353‑M2/107‑141 (MQ=255)                                            tggtatcatgccttcgctccgcatcttacgtgtctacagacaaaaaaaaggcgactcagcagtcgccttaaaaatcagtctgccagcgctgcgttcggcTGGCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTa                                                                                                    <  2:300737‑M2/44‑1 (MQ=255)                                             ggtttcatgcctgcgctcctcatcttacttttctacagacaaaaaaaaggcgactcatcagtcgccttaaaaatcagtttgccagcgccgccttctgcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTAc                                                                                                   <  1:377094‑M2/45‑1 (MQ=255)                                               tttcatgccttcgctcctcatcttacttttctacagacaaaaacaaggcgcctcatcagtcgccttaaaaatcagtttgccagcgccgccttctgcCGTCCCCTGCACTTCAATGATGCGCCCGTCCTCCGGCACCACCACGt                                                                                                 >  2:158451‑M2/97‑143 (MQ=255)                                                       cttcgcacctcatcttacttctctacagacaaaaaaaaggcgactcatcagtcgccttaaaaatcagtttgccagcgccgccttctgcCGTCCCCTGCACTTCAAAGAAGCGCCCGTCTTCGGTCATCACTACGTTCATGTCg                                                                                         <  2:609204‑M2/55‑1 (MQ=255)                                                             cctcatcttacttttctacagacaaaaaaaaaggcgactcatcagtcgccttaaaaatcagtttgccagcgccgccttctgcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGtcggtacct                                                                                   >  2:173007‑M2/83‑139 (MQ=255)                                                                taatcttactgttctacagacaaaaaaaaggcgactcatcagtcgccttaaaaatcagtttgccagcgccgccttctgcCGTCCCCTGCACTTCCATGATGCGCCCGACTTCGGTCCTCACTACGCTCATGTCGGTCTCTGCg                                                                                >  2:174073‑M2/80‑143 (MQ=255)                                                                   tcttacttttctacagacaaaaaaaaggcgactcatcagtcgccttaaaaatcagtttgccagcgccgccttctgcCGTCCACTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTGTCTGCGGCa                                                                             >  1:362791‑M2/77‑143 (MQ=255)                                                                    cttacttttctacagacaaaaaaaaggcgactcatcagtcgccttaaaaatcagtttgccagcgccgccttctgcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCag                                                                            >  1:424829‑M2/76‑143 (MQ=255)                                                                     ttccttttctacaggcaaaaaaaaggcgactcatcagtcgccttaaaaatcagtctgccagcgccgccttctgcCGGCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCaga                                                                           <  1:13247‑M2/69‑1 (MQ=255)                                                                         ttttctacagacaaaaaaaaggcgactcatcagtcgccttaaaaatcagtttgccagcgccgccttctgcCGTCTCCTGCACTTCAATGATCCGCCCGTCTTCCGTCACCGCTACGTTCCTTGCGGTCcccgcggccgccact                                                                       >  1:204227‑M2/71‑136 (MQ=255)                                                                         ttttctacagacaaaaaaaaggcgactcatcagtcgccttaaaaatcagtttgccagcgccgccttctgcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCCGCGGCAGCGtcc                                                                       >  1:272409‑M2/71‑142 (MQ=255)                                                                         ttttctacagacaaaaaaaaggcgactcatcagtcgccttaaaaatcagtttgccagcgccgccttctgcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGGTCATGTCGCTCTCTGCcgcagaggcg                                                                       >  2:345649‑M2/71‑140 (MQ=255)                                                                         ttttctacagacaaaaaaaaggcgactcatcagtcgccttaaaaatcagtttgccagcgccgccttctgcCGTCCCCTGCACTTCAATGATGCGCCCGTCTCCGGTCCTCACTACGTTCATGTCGCGCTCTGCGGCAGAGTCt                                                                       >  2:377732‑M2/71‑143 (MQ=255)                                                                           ttctacagacaaaaaaaaggcgactcatcagtcgccttaaaaatcagtttgccagcgccgccttctgcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGCCTTc                                                                     >  1:529669‑M2/69‑143 (MQ=255)                                                                             ctacctacaaaaaaaaggcgactcatcagtcgccttaaaaagcagtttgccagcgccgccttctgcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCaa                                                                   <  1:450975‑M2/77‑1 (MQ=255)                                                                             ctacagacaaaaaaaaggcgactcatcagtcgccttaaaaatcagtttgccagcgccgccttctgcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTCCATGTCGGTCTCTGCGGCAGGGTCTTCaa                                                                   >  1:155194‑M2/67‑143 (MQ=255)                                                                              tacagacaaaaaaaaggcgactcatcagtcgccttaaaaatcagtttgccagcgccgccttctgcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTtcacc                                                                  >  2:304666‑M2/66‑141 (MQ=255)                                                                                 agacaaaaaaagggcgacgcatcgggcgccttaaaaatcagtttgccagcgccgccttctgcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTa                                                               <  2:230119‑M2/81‑1 (MQ=255)                                                                                      tttttaaggcgactgatgagtcgccttaaaaatcagtttgccagcgccgccttctgcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGCCATCACTACGCTCATGTCGGTCTCTGCGGCGGTGTCTTCAACGTACTCCa                                                          >  1:97942‑M2/58‑143 (MQ=255)                                                                                      ataaaaaggcgcctcatcagtcgccttaatgatccgtttgccagcgccgccttctgcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCa                                                          <  2:72215‑M2/86‑1 (MQ=255)                                                                                         aaacggcgcctcatctctcgccttaaacatcagtttgccagcgccgccttctgcCGTCCCCTGCACTTCAATGCTGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCAGAt                                                       <  2:318677‑M2/89‑1 (MQ=255)                                                                                          aaaggcgattcatcagtcgccttaaaaatcagtttgccagcgccgccttctgcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCAGATc                                                      <  2:537119‑M2/90‑1 (MQ=255)                                                                                          aaaggcgactcatcagtcgtgttagagttcagtttgccagcgccgccttctgcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCAGATc                                                      <  1:519086‑M2/90‑1 (MQ=255)                                                                                                gactcatcagtcgccttaaaaatcagtttgccagcgccgccttctgcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGCATTCAACGCATTCCAGATCGCAAAc                                                >  2:53127‑M2/48‑143 (MQ=255)                                                                                                         gtcgccttaaaaatcagtttgccagcgccgccttctgcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCCGAGTCTTCAACGTATTCCAGATCGCAAACCGCTTcccc                                       >  1:166901‑M2/39‑143 (MQ=255)                                                                                                               ttaaaaatcagtttgccagcgccgccttctgcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCAGATCGCAAACCGCTTCGCCGTTcac                                 >  1:403868‑M2/33‑143 (MQ=255)                                                                                                                               cagcgccgccttctgcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCAGATCGCAAACCGCTTCGCCGTTCACCATTCCGCCAGAAACt                 >  1:657375‑M2/17‑143 (MQ=255)                                                                                                                                     cgctttctgcTGTCCCCTGCACTTCAATGATGCGCACGTATTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCAGATCGCAAACCGCTTCGCCGTTCACAATTCCGACAGAAACTGCGGCt           <  1:287303‑M2/133‑1 (MQ=255)                                                                                                                                     cgccttctgcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGGCCTCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCAGAGCGCAAACCGCTGCGCCGTCTCCACTTCCGCCAGAACCTGGGgcg           >  2:290416‑M2/11‑142 (MQ=255)                                                                                                                                       ccttctgcCGTCCCATGCACTTCAGTGATGGTCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCAGATCGCAAACCGCTTCGCCGTTCACAATTCCGACAGAAACTGCGGCTAc         <  2:217929‑M2/135‑1 (MQ=255)                                                                                                                                              cCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCAGATCGCACACCGCTTCGCCGTTCACAATTCCGACAGAAACTGCGGCTACCATCCCt  >  2:646324‑M2/2‑143 (MQ=255)                                                                                                                                               |                                                                                                                                              ACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATTCCCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCAGATCGCAAACCGCTTCGCCGTTCACAATTCCGACAGAAACTGCGGCTACCATCCCT  >  NC_000913/3815799‑3816082

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: ‑