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	30 [0]	[0] 30	[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)	27 (1.010)	24/264	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

AAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATTCCCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCAGATCGCAAACCGCTTCGCCGTTCACAATTCCGACAGAAACTGCGGCTACCATCCCT  >  NC_000913/3815801‑3816082                                                                                                                                             |                                                                                                                                              acctacttgctaagcgcaaattcaataaactggcgctgatatggtttcatgccttcgctcctcatcttacttttctacagacaaaaaaaaggcgactcatcagtcgccttaaaaatcagtttgccagcgccgccttctgcCGt                                                                                                                                             >  1:79741‑M2/141‑143 (MQ=255)               cgcaaattcactaacctggcccggatatggtttcatgccttcgctcctcatcttacttttctacagccaaaaaaaatgcgactcatcagtcgccttaaaaatcagtttgccagcgccgccttctgcCGTCCCCTGCACTTCaa                                                                                                                               <  1:324936‑M2/17‑1 (MQ=255)                gcaaattcaataaactggcgctgatatggtttcatgccttcgctcctcatcttacttttctacagacaaaaaaaaggcgactcatcagtcgccttaaaaatcagtttgccagcgccgccttctgcCGTCCCCTGCACTTCAat                                                                                                                              <  1:481604‑M2/18‑1 (MQ=255)                gcaaattcaataaactggcgctgatatggtttcatgccttcgctcctcatcttacttttctacagacaaaaaaaaggcgactcatcagtcgccttaaaaatcagtttgccagcgccgccttctgcCGTCCCCTGCACTTCAat                                                                                                                              >  1:259857‑M2/126‑143 (MQ=255)                 caaattcaataaactggcgctgatatggtttcatgccttcgctcctcatcttacttttctacagacaaaaaaaaggcgactcatcagtcgccttaaaaatcagtttgccagcgccgccttctgcCGTCCCCTGCACTTCAatg                                                                                                                             <  1:325450‑M2/19‑1 (MQ=255)                     ttcaataaactggcgctgatatggtttcatgccttcgctcctcatcttacttttctacagacaaaaaaaaggcgactcatcagtcgccttaaaaatcagtttgccagcgccgccttctgcCGTCCCCTGCCCTTCAATGATgc                                                                                                                         >  1:479785‑M2/121‑143 (MQ=255)                          taaactggcgctgatatggtttcatgccttcgctcctcatcttacttttctacagacaaaaaaaaggcgactcatcagtcgccttaaaaatcagtttgccagcgccgccttctgcCGTCCCCTGCCCTTAGATGATGCGACCg                                                                                                                    >  1:156468‑M2/116‑143 (MQ=255)                           aaactggcgctgatatggtttcatgccttcgctcctcatcttacttttctacagacaaaaaaaaggcgactcatcagtcgccttaaaaatcagtttgccagcgccgccttctgcCGTCCCCTGCACTTCAATGATGCGCCCGt                                                                                                                   <  2:428179‑M2/29‑1 (MQ=255)                                gccgctgatatggtttcatgccttcgctcctcatcttacgtttctacagacaaaaaaaatgcgactcatcagtcgccttaaaaatcagtttgccagcgccgccttctgcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCg                                                                                                              <  1:207567‑M2/34‑1 (MQ=255)                                  cgctgatatgggttcatgccttcgctcctcttcttacttttctacagacaaaaagaaggcgactcatcagtcgccttaaaaatcagtttgccagcgccgccttctgcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGt                                                                                                            <  2:120922‑M2/36‑1 (MQ=255)                                     tgatatggtttcatgccttcgctcctcatcttacctttctacagacaaaaaaagggcgactcatcagtcgccttaaaaatcagtttgccagcgccgccttctgcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGtcat                                                                                                         <  2:245511‑M2/39‑1 (MQ=255)                                     agatatggtttcatgcctccgctcctcatcttacttttctacagacaaaaaaaaggcgactcatcagtcgccttaaaaatcagtttgccagcgccgccttctgcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGtcat                                                                                                         <  1:436797‑M2/39‑1 (MQ=255)                                                          ctcctcatcttacttttctacagacaaaaaaaaggcgactcatcagtcgccttaaaaatcagtttgccagcgccgccttctgcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGtctc                                                                                    >  1:146188‑M2/84‑143 (MQ=255)                                                                  cttacttttctacagacaaaaaaaaggcgactcatcagtcgccttaaaaatcagtttgccagcgccgccttctgcCGTCCCCTGCACTTCAATGATGCGCCCGTCTGCGGGCACCACTACGTTCATGCAGGTCTCTGCCGCag                                                                            >  2:69048‑M2/76‑143 (MQ=255)                                                                     acttttctacagacaaaaaaaaggcgactcatcagtcgccttaaaaatcagtttgccagcgccgccttctgcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGt                                                                         >  1:366899‑M2/73‑143 (MQ=255)                                                                        tttctacagacaaaaaaaaggcgactcatcagtcgccttaaaaatcagtttgccagcgccgccttctgcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCtt                                                                      <  1:66223‑M2/74‑1 (MQ=255)                                                                           ctacagccgaaaaaccggcgactcatcagtcgccttaaaaatcagcttgccagcgccgccttctgcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCaa                                                                   <  1:313113‑M2/77‑1 (MQ=255)                                                                           ctacagacaaaaaaaaggcgactcagcagtcgccttaaaaatcagtttgccagcgccgccttctgcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCaa                                                                   <  1:368026‑M2/77‑1 (MQ=255)                                                                                gacaaaaaaaaggcgactcatcagtcgccttaaaaatcagtttgccagcgccgccttctgcCGTCCCCTGCACTTCAATGATGCGCCCGTCTCCGGTCATCACTACGGTCATTTCGGTCTCTGCGGCAGAGTCTTCAACGTAt                                                              <  1:501161‑M2/82‑1 (MQ=255)                                                                                                   atcagtcgccttaaaaatcagtttgccagcgcagccttctgcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTCCGTTCATGTCGGTCTCTGACGCCGAGTCTCCAACGTATCCCCGATCGCAAACcgcct                                           >  2:273768‑M2/43‑141 (MQ=255)                                                                                                               aaaaatcagtttgccagcgccgccttctgcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCAGATCGCAAACCGCTTCGCCGTTCACaa                               <  1:172126‑M2/113‑1 (MQ=255)                                                                                                                 aaatcagtttgccagcgccgccttctgcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCAGATCGCAAACCGCTTCGCCGTTCacaagt                             >  2:118982‑M2/29‑141 (MQ=255)                                                                                                                 aaatcagtttgccagcgccgccttctgcCGTCCCCTGCACCTCAATGATGCGCCCGCCTCCGGTCATCACTAGGTTCATGTCGGTCTCTGCGGCCGAGTCTTCAACGTATTCCAGATCGCGGACCGCTTCGCCGTTCCCAAtt                             >  1:9978‑M2/29‑143 (MQ=255)                                                                                                                     cagtgtgccagcgccgccttctgcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGGCTCTGCGGCAGAGTCTTCAACGTATTCCAGATCGCAAACCGCTTCGCCGTTCACAATTCCGa                         <  2:87897‑M2/119‑1 (MQ=255)                                                                                                                      agtttgccagcgccgccttctgcCGTCCCCTGCACTTCAATGAGGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCAGATCGCAAACCGCTTCGCCGTTCACAATTCCGAc                        <  2:286492‑M2/120‑1 (MQ=255)                                                                                                                             cagcgccgccttctgcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCCTCACTACGTTCATGTCGGTCTCTGCGGCCGAGTCTTCCACGTATCCCACATCGCAAACCGCTTCGCCGATCACAATTCCGACAGAAACt                 >  1:314278‑M2/17‑143 (MQ=255)                                                                                                                                 gccgccttctgcCGTCCCCTGCACTCCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCGGCTGCAGAGTCTTCAACGTATTCCAGTTCGCAAACCGCTTCGCCGTTCACAATTCCGACAGAAACTGCgg             <  1:442476‑M2/131‑1 (MQ=255)                                                                                                                                       ttctgcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGCGTCTTCAACGTATTCCAGATCGCAAACCGCTTCGCCGTTCACAATTCCGACAGAAACTGCGGCTACCa       >  1:477111‑M2/7‑143 (MQ=255)                                                                                                                                        tctgcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCAGATCGCAAACCGCTTCGCCGTTCACAATTCCGACAGAAACTGCGGCTACAAt      >  1:345615‑M2/6‑143 (MQ=255)                                                                                                                                            cCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCAGATCGCAAACCGCTTCGCCGTTCACAATTCGGACAGAACCTGCGGCTACCAtcccc  >  2:47570‑M2/2‑142 (MQ=255)                                                                                                                                             |                                                                                                                                              AAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATTCCCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCAGATCGCAAACCGCTTCGCCGTTCACAATTCCGACAGAAACTGCGGCTACCATCCCT  >  NC_000913/3815801‑3816082

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