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 [1]	[0] 29	[rph]–[rph]	[rph],[rph]

New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913	= 3815858	1 (0.020)	27 (0.670)	21/264	1.0	98.3%	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

AAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATTCCCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCAGATCGCAAACCGCTTCGCCGTTCACAATTCCGACAGAAACTGCGGCTACCATCC  >  NC_000913/3815803‑3816080                                                                                                                                           |                                                                                                                                            ctgcgtgctaagcgcaaattcaatatactcgcgctgatatggtttcatgccttcgctcctcatcttacttttctacagacaaaaaaaaggcgactcatcagtcgccttaaaaatcagtttgccagcgccgccttctgcCGTcc                                                                                                                                         <  1:66355‑M2/5‑1 (MQ=255)       gctaagcgcacattcaataacctggcgctgatatggtttcatgccttcgctcctcatcttacttttctacagacaaaaaaaaggcgactcatcagtcgccttaaatatcagtttgccagcgccgccttctgcCGTCCCCTGCa                                                                                                                                   <  1:687540‑M2/11‑1 (MQ=255)       gctaagcgcaaattcaataaactggcgctgatatggtttcatgccttcgctcctcatcttacttttctacagacaaaaaaaaggcgactcatcagtcgccttaaaaatcagtttgccagcgccgccttctgcCGTCCCCTGCa                                                                                                                                   <  1:677044‑M2/11‑1 (MQ=255)       gctaagcgcaaattcaataaactggcgctgatatggtttcatgccttcgctcctcatcttacttttctacagacaaaaaaaaggcgactcatcagtcgccttaaaaatcagtttgccagcgccgccttctgcCGTCCCCTGCa                                                                                                                                   <  1:676336‑M2/11‑1 (MQ=255)       gctaagcgcaaattcaataaactggcgctgatatgctttcatgccttcgctcctcatcttacttttctacagacaaaaaaaaggcgactcatcagtcgccttaaaaatcagtttgccagcgccgccttctgcCGTCCCCTGCa                                                                                                                                   <  2:679652‑M2/11‑1 (MQ=255)       ctaaacgcaaattgaataaactggcgttgatatggtgtcatgcattcactcctcatcttacatatcaacagacaaaaaaaaaggcgactcatcagtcgccttaaaaatcagtttgccagcgccgccttctgcCGTCCCCTGCa                                                                                                                                   <  2:467845‑M2/11‑1 (MQ=255)         taagcgcaaattcaataaactggcgctgatatggtttcatgccttcgctcctcatcttacttttctacagacaaaaaaaaggcgactcatcagtcgccttaaaactcagtttgccagcgccgccttctgcCGTCCCCTGCACt                                                                                                                                 <  2:687071‑M2/13‑1 (MQ=255)             cgcaaattcaatacactggcgctgatatggtttcatgccttcgctcctcatcttacttttctacagacaaaaaaaaggcgactcatcagtcgccttaaaaatcagtttgccagcgccgccttctgcCGTCCCCTGCACTTCaa                                                                                                                             <  1:594596‑M2/17‑1 (MQ=255)              gcaaattcaataaactggcgctgatatggtttcatgccttcgctcctcatcttacttttctacagacaaaaaaacggcgactcatcagtcgccttaaaaatcagtttgccagcgccgccttctgcCGTCCCCTGCCCTTCAat                                                                                                                            >  1:194797‑M2/126‑143 (MQ=255)                aaattcaataaactggcgctgatatggtttcatgccttcgctcctcatcttacttttctacagacaaaaaaaaggcgactcatcagtcgccttaaaaatcagattgccagcgccgccttctgaCGTCCCCTGCActtcatgca                                                                                                                          >  2:416040‑M2/124‑139 (MQ=255)                    tcaataaactggcgctgatatggtttcatgccttcgctcctcatcttacttttctacagacaaaaaaaaggcgactcatcagtcgccttaaaaatcagtttgccagcgccgccttctgcCGTCCCCTGCACTTCAATGATgcg                                                                                                                      <  2:645478‑M2/24‑1 (MQ=255)                                 gctgatatggtttcatgccctcgctcctcaacttgcttttctacagacaaaaaaaaggcgactcatcagtcgccttaaaaatcagtttgccagcgccgccttctgcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGtc                                                                                                         <  2:451331‑M2/37‑1 (MQ=255)                                      tattgtttcatgcctttgctcctcatctttcttttctacagacaaaaaaatggcgactcatcagtcgccttaaaaatcagtttgccagcgccgccttctgcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCAc                                                                                                    <  2:487034‑M2/42‑1 (MQ=255)                                           tttcatgccttcactcctcatcttacttttctacagacaaaaaaaaggcgactcatcagtcgccttagaaaccagtttgccagcgccgccttctgcCGTCCCCTGCACTTCACTGACGCTCCCGTTTTCGAACATCACTACGt                                                                                               >  1:314826‑M2/97‑143 (MQ=255)                                             tcatgccttcgctcctcatcttacttctctacagacaaaaaaaaggcgactcatcagtcgccttaaaaatcagtttgccagcgccgccttctgcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTc                                                                                             <  1:299286‑M2/49‑1 (MQ=255)                                             gcatgccttcgctcgtcatcttacttttctacacacaaaaaaaaggcgactcatcagtcgccttaaaaatcagtttgccagcgccgccttctgcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTc                                                                                             <  2:194797‑M2/49‑1 (MQ=255)                                                               tcttacttttctacagacaaaaaaaaggcgactcatcagtcgccttaaaaatcagtttgccagcgccgccttctgcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCa                                                                           >  1:573549‑M2/77‑143 (MQ=255)                                                                          tacagacaaaaaaaaggcgactcatcagtcgccttaaaaatcagtttgccagcgccgccttctgcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAgagtctcccac                                                                >  2:87700‑M2/66‑138 (MQ=255)                                                                                            gactcatcagtcgccttaaaaatcagtttgccagcgccgccttctgcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCAGATCGCAAAc                                              <  1:91452‑M2/96‑1 (MQ=255)                                                                                                 atcagtcgccttaaaaatcagtttgccagcgccgccttctgcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCATCGTATTCCCGCTTGCAAACCGCtt                                         >  1:677687‑M2/43‑143 (MQ=255)                                                                                                  tcagtcgccttaaaaatcagtttgccagcgccgccttctgcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCAGATCGCAAACCGCTTc                                        <  2:406828‑M2/102‑1 (MQ=255)                                                                                                  tcagtcgccttaaaaatcagtgtgccagcgccgccttctgcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCAGATCGCAAACCGCTTc                                        <  1:432886‑M2/102‑1 (MQ=255)                                                                                                        gccttaaaaatcagtttgccagcgccgccttctgcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCAGATCGCAAACCGCTTCGCCGtt                                  <  2:573549‑M2/108‑1 (MQ=255)                                                                                                         ccttaaaagtcaggttgccagcgccgccttctgcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCAGATCGCAAACCGCTTCGCCGTTc                                 <  2:472608‑M2/109‑1 (MQ=255)                                                                                                                    agtttgccagcgccgccttctgcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCAGATCGCAAACCGCTTCGCCGTTCACAATTccgcc                      >  1:419439‑M2/24‑141 (MQ=255)                                                                                                                            agcgccgccttctgcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCAGATCGCAAACCGCTTCGCCGTTCACAATTCCGACAGAAACTg              >  1:656408‑M2/16‑143 (MQ=255)                                                                                                                              cgccgcctgctgcCGTCCCCTGCACTGCAGTGGTGCGCCCGTCTTCGGTCATCACTACGTGCATGGCGGTCTCTGCGGCTGTGTCTTCAACGTATGCCAGATCGCAAACCGCTTCGCCGTTCACAATTCCGACAGAAACTGCg            <  2:130188‑M2/130‑1 (MQ=255)                                                                                                                                   ccttctgcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCAGATCGCAAACCGCTTCGCCGTTCCCAATTCCGACAGAAAATACGGGTAc       >  1:337477‑M2/9‑143 (MQ=255)                                                                                                                                        tgcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCAGATCGCAAACCGCTTCGCCGGTCACAATTCCGACAGAACCTGCGGCTACCATcc  >  1:168285‑M2/4‑143 (MQ=255)                                                                                                                                           |                                                                                                                                            AAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATTCCCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCAGATCGCAAACCGCTTCGCCGTTCACAATTCCGACAGAAACTGCGGCTACCATCC  >  NC_000913/3815803‑3816080

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