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	45 [0]	[0] 45	[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)	40 (0.950)	36/262	0.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

AAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATTCCCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCAGATCGCAAACCGCTTCGCCGTTCACAATTCCGACAGAAACTGCGGCTACCAT  >  NC_000913/3815802‑3816078                                                                                                                                            |                                                                                                                                          cctgcttgctaagcgcaaattcaataaactggcgctgatatggtttcatgccttcgctcctcatcttacttttctacagacaaaaaaaaggcgactcatcagtcgccttaaaaatcagtttgccagcgccgccttctgcCGTc                                                                                                                                        <  1:550678‑M2/4‑1 (MQ=255) cctgctcgctaagcgcaaattcaataaactggcgctgatatggtttcatgccttcgctcctcatcttacttttctacagacaaaaaaaaggcgactcatcaggcgccttaaaaatcagtttgccagcgccgccttctgcCGTc                                                                                                                                        <  1:543168‑M2/4‑1 (MQ=255)  ctgcttgctaagcgcaaattcaataaactggcgctgatatggtttcatgccttcgctcctcatcttacttttctcccgacaaaaaacaggcgactcatcagtcgccttaaaaatcactttgccagcgccgccctctgcCGTcc                                                                                                                                       >  1:647539‑M2/139‑143 (MQ=255)   tgcttgctaagcgcaaattcaataaactggcgctgatatggtttcatgccttcgctcctcatcttacttttctacagacaaaaaaaaggcgactcatcagtcgccttaaaaatcagtttgccagcgccgccttctgcCGTccc                                                                                                                                      <  1:382153‑M2/6‑1 (MQ=255)          taagcgcaaattcaataaactggcgctgatatggtttcatgccttcgctcctcatcttacttttctacagacaaaaaaaaggcgcctcatcagtcgccttaacaatcagtttgccagcgccgccttctgcCGTCCCCTGCACt                                                                                                                               >  1:16213‑M2/131‑143 (MQ=255)                         taaactggcgctgatatggtttcatgccttcgctcctcatcttacttttctacagacaaaaacaaggcgcctcatcagtcgccttaaaaataagtttgccagcgccgccttctgcCGTCCCCCGCACTTCAATGATCCGCCCg                                                                                                                >  2:815967‑M2/116‑143 (MQ=255)                           aactggcgctgatatggtttcatgccttcgctcctcatcttactcttctaccgtcagaaaaaaggcgactcatcagtcgccttaaaaatcagtttgccagcgccgccttctgcCGTCCCCTGCACTTCAATGATGCGCCCGTc                                                                                                              <  1:692155‑M2/30‑1 (MQ=255)                            actggcgctgatatggtttcatgccttcgctcctcatcttacttttctacagacaaacaaaaggcgactcatccgtcgccttaaaactcagtttgccagcgccgccttctgcCGTCCCCTGCACTGCAATGATGCCcccggcc                                                                                                             >  1:147472‑M2/113‑140 (MQ=255)                              tggcgctgatatggattcatgccttcgctcctcatcttacttttctacagacaaaaaaaaggcgactcatcagtcgccttaaaaatcagtatgcaagcgccgcctgctgcCGACACCTGCACTTCAACGAACAGCCTGTCTTc                                                                                                           >  2:549‑M2/111‑143 (MQ=255)                                              catgccttcgctcctcctcttacttttctacagacgaaaaaaaaggcggctcatcagtcgccttaaaaatcagtttgccagcgccgccttctgcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTc                                                                                           <  2:620299‑M2/49‑1 (MQ=255)                                                  gccttcgctcctcatcttacttttctacagacaaaaaaaaggcgactcatcagtcgccttaaaaatcagtatgccagcgccgccttctgcCGTCCCCTGCACTTCAATGATGCGCCCGTCTCCGGTGCTCActacattactgt                                                                                       >  2:827950‑M2/91‑138 (MQ=255)                                                     ttcgctcctcatcttacttttctacagacaaaaaaaaggcgactcatcagtcgccttaaaaatcagtttgccagcgccgccttctgcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCgg                                                                                    >  1:458378‑M2/88‑143 (MQ=255)                                                      tcgctcctcatcttacttttctacagacaaaaaaaaggcgactcatcagtcgccttaaaaatcagtttgccagcgccgccttctgcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGt                                                                                   <  1:493315‑M2/57‑1 (MQ=255)                                                           cctcatcttacttttctacagccaaaaaaaaggcgactcatcagtcgccttaaaaatcagtttgccagcgccgccttctgcCGTACCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTg                                                                              <  1:276512‑M2/62‑1 (MQ=255)                                                             tcatcttacttttctacagacaaaaaaaaggcgactcatcagtcgccttaaaaatcagtttgccagcgccgccttctgcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCg                                                                            <  1:798741‑M2/64‑1 (MQ=255)                                                              catcttacttttctacagacaaaaaaaaggcgactcatcagtcgccttaaaaatcagtttgccagcgccgccttctgcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCgg                                                                           >  1:504407‑M2/79‑143 (MQ=255)                                                                tcgcacttttccacagccaaaacaaaggcgactcatcagtcgccttaaaaatcagtttgccagcgccgcctgctgcCGTCCCCCGCACTTCAATGATGCGCCCGTCTTCGGGCATCACTACGTTCATGTCGGTCTCTGCGGCa                                                                         <  1:676034‑M2/67‑1 (MQ=255)                                                                  ttacctttctacagacaaaaaaaaggcgactcatcagtcgccttaaaaagcagtttgccagcgccgccttctgcCGTCCCCTGCACTTCAATGATGCGCCCGTCGTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCaga                                                                       <  1:9382‑M2/69‑1 (MQ=255)                                                                     cgttgctaccgacaaaaaaaaggcgactcatcagtcgccttaaaaatcagtttgccagcgccgccttctgcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTc                                                                    <  1:791839‑M2/72‑1 (MQ=255)                                                                       tttctacagacaaaaaaaaggcgactcatcagtcgccttaaaaatcagtttgccagcgccgccttctgcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCtt                                                                  <  1:53832‑M2/74‑1 (MQ=255)                                                                           tacagacaaaaaaaaggcgactcatcagtcgccttaaaaatcagtttgccagcgccgccttctgcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCCCTCCCCCCACGTCCGTCCCCCCGCCAGAGTCTCCAAc                                                              >  2:105665‑M2/66‑143 (MQ=255)                                                                                acaaaaaaaaggcgactcatcagtcgccttaaaaatcagtttgccagcgccgccttctgcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTAtt                                                         <  1:580722‑M2/83‑1 (MQ=255)                                                                                acaaaaaaaaggcgactcatcagtcgccttaaaaatcagtttgccagcgccgccttctgcCGTCCCCTGCACTTCAATGATCCGCCCGTCTTCGGTCATCACTACGTTCATGTCAGCCACTACGCCAAAGAATTAAACGTAtt                                                         >  1:808786‑M2/61‑143 (MQ=255)                                                                                  aaaaaaaaggcgactcatcagtcgccttaaaaatcagtttgccagcgccgccttctgcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCAc                                                                                                  >  2:739271‑M2/59‑100 (MQ=255)                                                                                  aaaaaaaaggcgactcatcagtcgccttaaaaatcagtttgccagcgccgccttctgcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCAc                                                                                                  <  1:739271‑M2/42‑1 (MQ=255)                                                                                  aaaaaaaaggcgactcatcagtcgccttaaaaatcagtttgccagcgccgccttctgcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGtattac                                                       >  1:685463‑M2/59‑141 (MQ=255)                                                                                    aaaaaaggcgactcatcagtcgccttaaaaatcagtttgccagcgccgccttctgcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCAg                                                     >  1:313662‑M2/57‑143 (MQ=255)                                                                                    aaaaaaggcgactcatcagtcgccttaaaaatcagtttgccagcgccgccttctgcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTCCAACGTATTCCAg                                                     >  1:473197‑M2/57‑143 (MQ=255)                                                                                            cgactcatcagtcgccttaaaaatcagtttgccagcgccgccttctgcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCAGATCGCaaa                                             <  2:823077‑M2/95‑1 (MQ=255)                                                                                                tcatcagtcgccttaaaaatcagtttgccagcgccgccttctgcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTGAACGGATTCCAGACCGCAAACCGc                                         >  2:398254‑M2/45‑143 (MQ=255)                                                                                                  atcagtcgccttaaaaatcagtttgccagcgccgccttctgcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCAGATCGCAAACCGctg                                       >  2:67352‑M2/43‑142 (MQ=255)                                                                                                   tcaggggcctgaaacatcagtttgccagcgccgccttcggcCGTCCGCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCAGATCGCAAACCGCTTc                                      <  1:204898‑M2/102‑1 (MQ=255)                                                                                                    cagtcgccttaaaaatcagtttgccagcgccgccttctgcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCAGATCGCAAACCGCTTCg                                     >  1:558724‑M2/41‑143 (MQ=255)                                                                                                          ccttaaaaatcagtttgccagcaccgccttctgcGGTCCCCAGCACGGCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCAGATCGCAAACCGCTTCGCCGTTc                               <  1:305526‑M2/109‑1 (MQ=255)                                                                                                                     cgtttgccagcgccgccttctgcCGTCCCGTGCACTTCAATGACGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCAGATCGCAAACCGCTTCGCCGTTCACAATTCCGAc                    <  2:647539‑M2/120‑1 (MQ=255)                                                                                                                     cgtttcccatcgccgcctcctgcCGTCCCCTGCACCTCAAGTATGCGCCCGTCTTCGGTCATCACTACGCTCATGTCGGTCTCTGCCGCAGAGTCTTCAACGTACTCCAGATCGCAAACCGCTTCGCCGTTCACAATTCCGAc                    <  2:76311‑M2/120‑1 (MQ=255)                                                                                                                      gtttgccagcgccgccttctgcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCAGATCGCAAACCGCTTCGCCGTTCACAATTCCGACa                   <  2:132006‑M2/121‑1 (MQ=255)                                                                                                                         tgccagcgccgccttctgcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCAGATCGCAAACCGCTTCGCCGTTCACACTTCCGACAGaa                >  1:510283‑M2/20‑143 (MQ=255)                                                                                                                            cagcgccgccttctgcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCAGATCGCAAACCGCTTCGCCGTTCACAATTCCGACAGAAACt             >  1:601796‑M2/17‑143 (MQ=255)                                                                                                                                gccgccttctgcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCAGATCGCAAACCGCTTCGCCGTTCACAATTCCGACAGAAACTGCgg         <  1:199474‑M2/131‑1 (MQ=255)                                                                                                                                  cgccttctgcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTGGGTCTCTCCGGCAGAGTCTTCAACGTATTCCAGATCGCAAACCGCTTCGCCGGTCACAATTCCGACACAAACCGCGgcc       >  2:426809‑M2/11‑142 (MQ=255)                                                                                                                                  cgccttctgcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCAGATCGCAAACCGCTTCGCCGTTCACAATTCCGACACAAACTGCGGCt       >  2:498240‑M2/11‑143 (MQ=255)                                                                                                                                   gccttctgcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGCCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCAGATCGCAAACAGCTTCGCCGTTCACCATCCCGACAGCAACGGCGGCTa      >  2:570865‑M2/10‑143 (MQ=255)                                                                                                                                    ccttctgcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCAGATCGCAAACCGCTTCGCCGTTCACAATTCCGCCAGAAACTGCGGCTAc     >  1:826776‑M2/9‑143 (MQ=255)                                                                                                                                       tctgcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCAGATCGCAAACCGCTTCGCCGTTCACAATTCCGACAGAAACTGCGGCTACCAt  <  1:674245‑M2/138‑1 (MQ=255)                                                                                                                                            |                                                                                                                                          AAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATTCCCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCAGATCGCAAACCGCTTCGCCGTTCACAATTCCGACAGAAACTGCGGCTACCAT  >  NC_000913/3815802‑3816078

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