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	31 [0]	[0] 31	[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)	30 (0.780)	26/262	0.5	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

ACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATTCCCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCAGATCGCAAACCGCTTCGCCGTTCACAATTCCGACAGAAACTGCGGCTACCATCC  >  NC_000913/3815799‑3816080                                                                                                                                               |                                                                                                                                            acccctgcttgctaagcgcaaattcaataacctggcgctgatatggtttcatgccttcgctcctcatcttacttttctacagaccaaaaaaaggcgactcatcagtcgccttaaaaatcagtttgccagcgccgccttctgcc                                                                                                                                             <  1:569650‑M2/1‑1 (MQ=255) acacctgcttgctaagcgcaaattcaataaactggcgctgatatggtttcatgccttcgctcctcatcttacttttctacagacaaaaaaaaggcgactcatcagtcgccttaaaaatcagtttgccagcgccgccttctgcc                                                                                                                                             >  1:598720‑M2/143‑143 (MQ=255)         ttgctaagcgcaaagtcaataaactggcgctgatatggtttcatgccttcgctcctcatcttacttttctacagacaaaaaaaaggcgactcatcagtcgccttaaaaatcagtttgccagcgccgccttctgcCGTCCCCTg                                                                                                                                     <  2:553095‑M2/9‑1 (MQ=255)         ttgcaaagcgcaaattcaagaaactggcgctgatatggtgtcatgccttcgctcctcatcttacttttcgacagacaaaaaaaaggcgactcatcagtcgccttaaaaatcagtttgccagcgccgccttctgcCGTCCCCTg                                                                                                                                     <  2:553754‑M2/9‑1 (MQ=255)                  gcaaattcaataaactggcgctgatatggtttcatgcattcgctcctcatcttacttttctacagacaaaaaaaaggcgactcatcagtcgccttaaacatcagtttgccagcgccgccgtctgcTGTCCCATGCACTccact                                                                                                                            >  1:121492‑M2/126‑141 (MQ=255)                       ttcaataaactggcgctgatatggtttcatgccttcgctcctcatcttacttttctacagacaaaaaaaaggcgactcatcagtcgccttaaaaatcagtttgccagcgccgccttctgcCGTCCCCTGCACTTCAATGATgc                                                                                                                       <  2:194708‑M2/23‑1 (MQ=255)                        tgaatagactggcgcttatgtgtgttcaggccttggctcctcgtcttatttttctacagacaaaaaaagggcgactcatcagtcgccttaaaaatcagtttgccagcgccgccttcttcCGTCCCCTGCACTTCAATGATgcg                                                                                                                      <  2:644860‑M2/24‑1 (MQ=255)                        tcaataaactggcgctgatatggtttcatgccttcgctcctcatcttacttttctacagacaaaaaaaaggcgactcatcagtcgccttaaaaatcagtttgccagcgccgccttctgcCGTCCCCTGCACTTCCATGCTgcg                                                                                                                      >  1:437658‑M2/120‑143 (MQ=255)                                            tggttccatgccttcgctcctcatcttacttttctacagacaaaaaaaaggcgactcatcagtcgccttaaaaatcagtttgccagcgccgccttctgcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTa                                                                                                  <  2:389071‑M2/44‑1 (MQ=255)                                                  catgccttcgctcctcatcttacttttctacagacaaaaaaaaggcgactcatcagtcgccttaaaactcagtttgccagcgccgccttctgcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCa                                                                                            <  2:70343‑M2/50‑1 (MQ=255)                                                   atgccttcgctcctcatcttacttttctacagacaaaaaaaaggcgactcatcagtcgccttaaaaatcagtttgccagcgccgccttctgcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGCCATCACTACGATCAt                                                                                           >  2:605846‑M2/93‑143 (MQ=255)                                                     gccttcgctcctcatctgacttttctacagacaaaaaaaaggcgactcatcagtcgccttaaaaatcagtttgccagcgccgccttctgcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGt                                                                                         <  1:513691‑M2/53‑1 (MQ=255)                                                             tcctcatcttacttttctacagacaaaaaaaaggcgactcatcagtcgccttaaaaatcagtttgccagcgccgccttctgcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGggctcg                                                                                 >  1:21916‑M2/83‑142 (MQ=255)                                                             tcctcatcttacttttctacagacaaaaaaaaggcgactcatcagtcgccttaaaaatcagtttgccagagccgccttctgcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTctct                                                                                 >  2:413343‑M2/83‑143 (MQ=255)                                                                  atcttacttttctacagacaaaaaaaaggcgactcatcagtcgccttaaaaatcagtttgccagcgccgccttctgcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGc                                                                            <  2:513898‑M2/66‑1 (MQ=255)                                                                   tcttacttttctacagacaaaaaaaaggcgactcatcagtcgccttaaaaatcagtttgccagcgccgccttctgcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCa                                                                           >  1:298160‑M2/77‑143 (MQ=255)                                                                         ttttctacagacaaaaaaaaggcgactcatcagtcgccttaaaaatcagtttgccagcgccgccttctgcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCt                                                                     >  1:743390‑M2/71‑143 (MQ=255)                                                                             cgacagacaaaagaaaggcgacgcatcagtcgccttaaaaatcagtttgccagcgccgccttctgcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCaa                                                                 <  1:119552‑M2/77‑1 (MQ=255)                                                                              tacagacaaaaaaaaggcgactcatcagtcgccttaaaaatcagtttgccagcgccgccttctgcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAAc                                                                >  1:393227‑M2/66‑143 (MQ=255)                                                                              tacagacaaaaaaaaggcgactcatcagtcgccttaaaaatcagtttgccagcgacgccttctgcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAAc                                                                >  1:396771‑M2/66‑143 (MQ=255)                                                                              tacagacaaaaaaaaggagactcatcattcgccttcacgtttgtcgcgccagcgccgccttctgcCGTCCCCTGTACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAAc                                                                <  1:551047‑M2/78‑1 (MQ=255)                                                                                   acaaaaaaaaggcgactcatcagtcgccttaaaaatcagtttgccagcgccgccttctgcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTAtt                                                           >  1:594991‑M2/61‑143 (MQ=255)                                                                                                   tcatcagtcgccctaaaaatcagtttgccagcgccgccttctgcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCACCGTGATCCATGTGGCAaaccac                                           >  2:306942‑M2/45‑141 (MQ=255)                                                                                                     atcagtcgccttaaaactcagtctgccagcgccgccttctgcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCAGATCGCAAACCGCtt                                         <  2:584115‑M2/101‑1 (MQ=255)                                                                                                       cggtcgcgttaaaagtcagttggccagcgccgccttctgcCGTCCCCTGCACTGCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCAGATCGCAAACCGCTTCg                                       <  2:301278‑M2/103‑1 (MQ=255)                                                                                                            gccttaaaaatcagtttgccagcgccgccttctgcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCAGATCGCAAACCGCTTCGCCGtt                                  <  1:532837‑M2/108‑1 (MQ=255)                                                                                                              cttaaaaatcagtttgccagcgccgccttctgcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGTGTATTCAACGTATTCCAGATCGCACACCGCTTCGCCGTAca                                >  1:64048‑M2/34‑143 (MQ=255)                                                                                                                            tgccagcgccgccttctgcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCAGATCGTAAACCGCTTCGTCGTTCACAATTCCGACAGaa                  <  2:396771‑M2/124‑1 (MQ=255)                                                                                                                            tgccagcgccgccttctgcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCAGATCGCAAACCGCTTCGCCGTTCACAATTCCGACAGaa                  <  2:393227‑M2/124‑1 (MQ=255)                                                                                                                              ccagcgccgccttctgcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTAGTCCAGATCGCAAACCGCTTCGCCGTTCACAATTCCgacagcacc                >  1:181130‑M2/18‑139 (MQ=255)                                                                                                                                            tgcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTGCATGTCGGTCGCTGCGGCAGCGTCTTCAACGTATTCCAGATCGCAAACCGCTTCGCCGTTCACACTTCCGACAGACACTGCGcctcccctcc  >  1:94909‑M2/4‑143 (MQ=255)                                                                                                                                               |                                                                                                                                            ACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATTCCCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCAGATCGCAAACCGCTTCGCCGTTCACAATTCCGACAGAAACTGCGGCTACCATCC  >  NC_000913/3815799‑3816080

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