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	29 [0]	[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	0 (0.000)	28 (0.740)	24/264	0.6	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

AAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATTCCCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCAGATCGCAAACCGCTTCGCCGTTCACAATTCCGACAGAAACTGCGGCTACCATCCC  >  NC_000913/3815804‑3816081                                                                                                                                          |                                                                                                                                             tgcttgctaagcgcaaattcaataaactggcgctgatatggtttcatgccgtcgctcctcatcttacttttcgacagacaaaaaaaaggcgactcatcagtcgccttaaaaatcagtttgccagcgccgccttctgcCGTccc                                                                                                                                         <  2:549038‑M2/6‑1 (MQ=255)   cttgcgaagcgcaaatgcaataagctggcgctgatatggtttcatgccttcgctcctcatcttacttttctacagacaaaaaaaaggcgactcatcagtcgccttaaaaatcagtttgccagcgccgccttctgcCGTCCCCt                                                                                                                                       <  1:213683‑M2/8‑1 (MQ=255)           gcgcaaattcaataaactggcgctgatatggtttcatgccttcgctcctcatcttacttttctacagacaaaaaaaaggcgactcatcagtcgccttaaaaatcagtttgccagagccgccttctggCGGCCACTGCCCTTCa                                                                                                                               >  1:115288‑M2/128‑143 (MQ=255)                aattcaacaatccggcgctgatatggtttcatgccttcgctcctcatcttacttttctacagacaaaaaaaaggcgactcatcagtcgccttaaaaatcagtttgccagcgccgccttctgcCGTCCCCTGCACTTCAatgat                                                                                                                          <  2:53450‑M2/21‑1 (MQ=255)                aatgcaatcagctggcgctgatatggtttcatgccttcgctcctcatcttacttttctacagacaaaaaaaaggcgactcatcagtcgccttaaaaatcagtttgccagcgccgccttctgcCGTCCCCTGCACTTCAatgat                                                                                                                          <  2:78621‑M2/21‑1 (MQ=255)                  ttcaataaactggcgctgatatggtttcatgccttcgctcctcatcttacttttctacagacaaaaaaaaggcgactcatcagtcgccttaaaaatcagtttgccagcgccgccttctgcCGTCCCCTGCACTTCACTGATgc                                                                                                                        >  1:390710‑M2/121‑143 (MQ=255)                    caataaactggcgctgatatggtttcatgccttcgctcctcatcttacttttctacagacaaaaaaaaggcgactcatcagtcgccttaaaaatcagtttgccagcgccgccttctgcCGTCCCCTGCACTTCAATGATgcgc                                                                                                                      <  1:375182‑M2/25‑1 (MQ=255)                    caagaaactggcgagaatatggtttcatgccttcggtcctcatcatactttactaaagacaaaaaaaaggcgactcataagtcgccttaaaaatcagtttgccagcgccgccttctgcCGTCCCCTGCACTTCAATGATgcgc                                                                                                                      <  2:137997‑M2/25‑1 (MQ=255)                     aataaactggcgctgatatggtttcatgccttcgctcctggtgttacttttctacagacaaaaaaaaggcgactcatcagtcgccttaaaaatcagtttgccagcgccgccttctgcCGTCCCCTGCACTTCAATGATGCGcc                                                                                                                     >  1:172457‑M2/118‑143 (MQ=255)                      ataaactggcgctgatatggtttcatgccttcgctcctcatcttacttttctacagacaaaaaaaagacgactcatcagtcgccttaaaaatcagtttgccagcgccgccttctgcCGCCCCCTGCACCTCAATGatgcgcac                                                                                                                    >  1:138187‑M2/117‑141 (MQ=255)                                   gatatggtttcatgccttcgctcctcatcttacttttctacagacaaaaaaaaggcgactcatcagtcgccttaaaaatcagtttgccagcgccgccttctgcCGTCCCCTGCACTTCACTGATGCGCCCGTCTTCGGtcatc                                                                                                       >  1:67321‑M2/104‑143 (MQ=255)                                              atgccttcgctcctcatcttacttttctacagacaaaaaaaaggcgactcatcagtcgccttaaaaatcagtttgccagcgccgccttctgcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCAt                                                                                            <  1:628626‑M2/51‑1 (MQ=255)                                                gccttcgctcctcatcttacttttctacagacaaaaaaaaggcgactcatcagtcgccttcaaaatcagtttgccagagccgcctcccgcCGTTCCCGGCCCTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCTTGt                                                                                          >  1:270888‑M2/91‑143 (MQ=255)                                                  ctgcgcttctcatctgacttgacgacagacaaaaaaaaggcgactcagcagtcgccttaaagatcagtttgccagcgccgccttctgcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCg                                                                                        <  2:448391‑M2/55‑1 (MQ=255)                                                      gctcctcatctgacttttctacagacaaaaaaaaggcgactcatcagtcgccttaaaaatcagtttgccagcgccgccttctgcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGtct                                                                                    <  2:15845‑M2/59‑1 (MQ=255)                                                          gtcaggttactgttctacagacaaaaaaaaggcggctcatcagtcgccttaaaaatcagtttgccagcgccgccttctgcCGTCCCCTTCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGc                                                                                <  1:402735‑M2/63‑1 (MQ=255)                                                                ttacttttctacagacaaaaaaaaggcgactcatcagtcgccttaaaaatcagtttgccagcgccgccttctgcCGGCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCaga                                                                          <  1:10474‑M2/69‑1 (MQ=255)                                                                    ttttctacagacaaaaaaaaggcgactcatcagtcgccttaaaaatcagtttgccagcgccgccttctgcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCGCTGCGGCAGAGTCt                                                                      >  1:97806‑M2/71‑143 (MQ=255)                                                                      ttctacagacaaaaaaaaggcgactcatcagtcgccttaaaaatcagtttgccagcgccgccttctgcCGTCCCCTGCACTGCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCACAGTCTTc                                                                    >  2:482518‑M2/69‑143 (MQ=255)                                                                         tacagacaaaaaaaaggcgactcatcagtcgccttaaaaatcagtttgccagcgccgccttctgcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAAc                                                                 >  1:267956‑M2/66‑143 (MQ=255)                                                                              acaaaaaaaaggcgactcatcagtcgccttaaaaatcagtttgccagcgccgccttctgcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTAtt                                                            >  1:363226‑M2/61‑143 (MQ=255)                                                                                atttttaaggcgactgatgagtcgccttaaaaatcagtttgccagcgccgccttctgcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTcc                                                          >  2:594293‑M2/59‑143 (MQ=255)                                                                                    aaaaggcgactcatcagtcgccttaaaaatcagtttgccagcgccgccttctgcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCAGAt                                                      <  1:53096‑M2/89‑1 (MQ=255)                                                                                              tcatcagtcgccttaaaaatcagtttgccagcgccgccttctgcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCAAGATCGCAAACCGc                                            >  1:555423‑M2/45‑143 (MQ=255)                                                                                              tcatcagtcgccttaaaaatcagtttgccagcgccgccttctgcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAAAGTCTTCAACGTATTCCAGATCGAAAACGGc                                            >  1:426489‑M2/45‑143 (MQ=255)                                                                                                          ttaaaaatcagtttgccagcgccgccttctgcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCAGATCGCAAACCGCTTCGCCGttccc                                >  1:89470‑M2/33‑141 (MQ=255)                                                                                                            aaaaatcagtttgccagcgccgccttctgcCGTCCCCTGCACTTCAATGATGCGCCCGTCCTCGGCCATCACTACCTCCACGTCGGTCTCTTCGGCAGAGCCTTCACCGTATTCCAGCTCTCAAACCGCTTCGCCCTTCCCaa                              >  2:712222‑M2/31‑143 (MQ=255)                                                                                                              aaatcagtttgccagcgccgccttctgcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCCTGCCGCGCACGGCGGCAGAGTCGTCAACGTATGCCACATCGCAAACCGCTCCGCCGTCCACAAtt                            >  1:336266‑M2/29‑143 (MQ=255)                                                                                                                                        gcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGATCTCTGCGGCAGAGTCTTCAACGTATTCCCAATCGCACACCGCTTCGCCGTTCCCAATTCCGACAGAAACTGCGGCTAGCAGccc  >  2:406248‑M2/3‑143 (MQ=255)                                                                                                                                          |                                                                                                                                             AAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATTCCCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCAGATCGCAAACCGCTTCGCCGTTCACAATTCCGACAGAAACTGCGGCTACCATCCC  >  NC_000913/3815804‑3816081

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