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	41 [0]	[0] 40	[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)	35 (0.980)	32/264	0.1	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

AAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATTCCCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCAGATCGCAAACCGCTTCGCCGTTCACAATTCCGACAGAAACTGCGGCTACCATCCCTT  >  NC_000913/3815804‑3816083                                                                                                                                          |                                                                                                                                               ggcttgctagcctctagtttaatacactggccctggtgtggtttcatgccttcgctcctgatcttacttttctacagacaaaaaaaaggcgactcatcagtcgccttaaaaatctgtttgccagcgccgccttctgcCGTccc                                                                                                                                           <  2:554402‑M2/6‑1 (MQ=255)     tgctaagcgcaaatccaataaactggcgctgatatggtttcatgccgtcgctcctcatcttacttttctacagacaaaaaaaaggcgactcatcagtcgccttaaaaatcagtttgccagcgccgccttctgcCGTCCCCTGc                                                                                                                                       <  1:611076‑M2/10‑1 (MQ=255)        taagcgcaaagtcaataaactggcgctgatatggtttcatgccttcgctcctcatcttacttttctacagacaaaaaaaaggcgactcatcagtcgccttaaaaatcagtttgccagcgccgccttccgcCGTCCCCTgcaat                                                                                                                                    <  2:398563‑M2/13‑3 (MQ=255)         aagcgcaaattcaataaactggcgctgatatggtttcatgccttcgctcctcatcttacttttctacagacaaaaaaaagacgactcatcagtcgccttaaaaatcagtttgccagcgccgccttctgcCgtcccctcgtctt                                                                                                                                   >  1:187805‑M2/130‑137 (MQ=255)         aagcgcaaaattacgaaactggcgagaatatagatacaagccgtcgctcctcatcttacttatctacagccaaaaaaaaggcgactcatcagtcgccttaaaaatcagtttgccagcgccgccttctgcCGTCCCCTGCACtt                                                                                                                                   <  1:146149‑M2/14‑1 (MQ=255)          agcgcaaattcaataaactggcgctgatatggtttcatgccttcgctcctcatcttacttttctacagacagaaaaaaggcgactcatcagtcgccttacaaacacgtttcccaacgccgccttctgcCGGCCCCTGCcctcg                                                                                                                                  >  1:491898‑M2/129‑141 (MQ=255)            ctcaaattcaataaactggcgctgatatggtttcatgccttcgctcctcatcttacttttctacagacaaaaaaaaggcgactcatcagtcgccttaaaaatcagtttgccagcgccgccttctgcCGTCCCCTGCACTTCaa                                                                                                                                <  2:719991‑M2/17‑1 (MQ=255)            cgtaaattcaataaactgacgctgatatggttgcaggccttcgctcctcatcttacttttctacagacaaaaaaaaggcgactcatcagtcgccttaaaaatcagtttgccagcgccgccttctgcCGTCCCCTGCACTTCaa                                                                                                                                <  1:88101‑M2/17‑1 (MQ=255)                aattgaataaactggcgctgatatggtttcatgccttcgctcctcatcttacttttctacagacaaaaaaaaggcgactcatcagtcgccttaaaaatcagtttgccagcgccgccttctgcCGTCCCCTGCACTTCAatgat                                                                                                                            <  2:330128‑M2/21‑1 (MQ=255)                 attcaataacctggcgctgatatggtttcatgccttcgcacctcatcttacttttctacagacaaaaaaaaggcgactcatcagtcgccttaaaaatcagtttgccagcgccgccttctgcCGTCCCCTGCACTTCAatgatg                                                                                                                           <  1:476394‑M2/22‑1 (MQ=255)                               agcagagatggtgacatgccttcgctccacctcttacttttcgacagacaaaaaaaaggcgactcatcagtcgtcttaaaactcagtttgccagcgccgccttctgcCGACCCCTGCACTTCAATGATGCGCCCGTCTTCGGt                                                                                                             <  1:287345‑M2/36‑1 (MQ=255)                                         ggttgatgtcttcggtcttcatcttatgtttctacagacaaaaaaaaggcgactcatcagtcgccttaaaattcagtttgccagcgccgccttctgcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACg                                                                                                   <  2:204870‑M2/46‑1 (MQ=255)                                                      gctcctcatcttatttttctacatactaaaaataggcgtctcatcagtcgccttaaaaatcagtttgccagcgccgcctgctgcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGtct                                                                                      <  1:341959‑M2/59‑1 (MQ=255)                                                      gctcctcatcttacttttctacagacaaaaaaaaggcgactcatcagtcgccttaaaaatcagtttgccagcgccgccttctgcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTAAGATCATGACGGtct                                                                                      >  2:417161‑M2/85‑143 (MQ=255)                                                       gccctcagattacttttctacagacagaaaaaatgcgtctcttcagtcgccttaaaaatcattttgctagccccgccttctgcCGTCCCCTGCACTTCAATGATGCGCTCGTCTTCGGTCATCACTACGTTCATGTCGGtctc                                                                                     <  2:319882‑M2/60‑1 (MQ=255)                                                       ctcctcatcttacttttctacagacaaaaaaaaggcgactcatcagtcgccttaaaaatcagtttgccagcgacgccttctgcCGTCCCCTGCACTTCAATGCTGCGCCCGTCTTCGGTCATCACTACGATCATGTCGGtctc                                                                                     >  2:234070‑M2/84‑143 (MQ=255)                                                        tcctcatcttacttttctacagacaaaaaacacgcgactcatcagtcgccttaaaaatcagtttgccagcgccgccttctgcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTctct                                                                                    <  2:363856‑M2/61‑1 (MQ=255)                                                                     ttgctacagacaaaaaaaaggcgactcatcagtcgccttaaaaatcagtttgcaagcgccgccttctgcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCtt                                                                       <  2:683693‑M2/74‑1 (MQ=255)                                                                        ttacatacaaaaaaaaggcgattcttcagtcgctttaataatcattttgccagcgccgcattctgcCGGCCCCTGCACTGCAATGATGCGCCCGGCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCaa                                                                    <  1:205416‑M2/77‑1 (MQ=255)                                                                        ctacagacaaaaaaaaggctactcatcagtcgccttaaaaatcagtttgccagcgccgccttctgcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCaa                                                                    <  1:718811‑M2/77‑1 (MQ=255)                                                                         tacagacaaaaaaaaggcgactcatcagtcgccttaaaaatcagtttgccagcgccgccttctgcCGTTCGCTGCACTTCAATGTTGCGCCCGTTTTCGGTCATCACTACGTTCATGTAGTGCTCCACGGCAGAGACTtcacc                                                                   >  1:146784‑M2/66‑141 (MQ=255)                                                                              acaaaaaaaaggcgactcttcagtcgccttaaaaatcagtttgccagcgccgccttctgcCGTCCCCTGCACTTCAATGAAGCGCCCCTCTTCGGTCATCACTAAGTTCATGTCGGTCTCTGCGGCAAAgactccaaacaact                                                              >  2:231906‑M2/61‑137 (MQ=255)                                                                                aaaaaaaaggcgactcatcagtcgccttaaaaatcagtttgccagcgccgccttctgcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACGACGTTCATGTCGGTCTCTGCAGCAGAGTCGTCAACGAATTcc                                                            >  2:237899‑M2/59‑143 (MQ=255)                                                                                 aaaaaaaggcgactcatcagtcgccttaaaaatcagtttgccagcgccgccttctgcCGTCCCCTGCCCTTCAATGATGCGCCCGTCTTCGGTAAACACTACAATCATACCGGACTCTAGGGCAAAGTCTTCAAAGTATTCCa                                                           >  1:421611‑M2/58‑143 (MQ=255)                                                                                    ataaggcgactcatcagtcgggttaaagagcagatcggcggcgcctcgtcctgcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCAGAt                                                        <  1:328270‑M2/89‑1 (MQ=255)                                                                                            atccatcagtcgccttaaaaatcagtttgcctgcgccaccttctgcCGTCCCCTGCACTTCATTGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTTCAGATCGCAAAcc                                                <  2:93318‑M2/97‑1 (MQ=255)                                                                                                    gtcgccttaaaaatcagtttgccagcgccgacttctgcCGTCCCCTGCACTTCAATGATGCGCCCGTCATCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCAGATCGCAAAACGCAACGcc                                        >  2:28179‑M2/39‑143 (MQ=255)                                                                                                     tcgccataaaaagcagtgtgccagcgccgccttctgcCGGCGCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCAGATCGCAAACCGCTTCGCCg                                       <  1:690338‑M2/106‑1 (MQ=255)                                                                                                            aaaaatcagtttgccagcgccgccttctgcCGTCCCCTGCATTTCTATGATTCGCCCGTCTTCGTTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCAGATCGCAAACCGCTTCGCCGTTCACaa                                <  1:265995‑M2/113‑1 (MQ=255)                                                                                                             aaaatcagtttgccagcgccgccttctgcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCCCTACGTTCATGTCGGTCTCTGCGGCAGAGCCTCCAACGTACTCCAGCCCGCAAACCCCTTCGCCGTCCACAAt                               >  2:490615‑M2/30‑143 (MQ=255)                                                                                                                  cagtttgccagcgcccctttctgcCGTCTCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCAGATCGCAAACCGCTTCGCCGTTCACAATTCCGa                          <  2:189418‑M2/119‑1 (MQ=255)                                                                                                                     tttgccagcgccgccttctgcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCAGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCAGATCGCAAACCGCTTCGCCGTTCACAATTCCGACAg                       >  1:376062‑M2/22‑143 (MQ=255)                                                                                                                       tgccagcgccgccttctgcCGTCCCCTGCACTTCAATGATGCGCCCGTCTCCGGTCATCACTACGTTCATGACGGTCTCTGCGGCAGAGTCTTCAACGTATTCCAGATCGCAAACCGCTTCGCCGCTCACAATTACGCCAGaa                     >  2:296328‑M2/20‑143 (MQ=255)                                                                                                                         ccagcgccgccttctgcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCAGATCGCAAACCGCTTCGCCGTTCACACTTCCGACAGAAAc                   >  2:456515‑M2/18‑143 (MQ=255)                                                                                                                             cgccgccttctgcCGTCGCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCAGATCGCAAACCGCTTCGCCGTTCACAATTCCGACAGAAACTGCg               <  2:612193‑M2/130‑1 (MQ=255)                                                                                                                                cgccttctgcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCAGATCGCAAACCGCTTCGCCGTTCACAATTCCGACAGAAACTGCGGCt            <  1:591892‑M2/133‑1 (MQ=255)                                                                                                                                cgccttctgcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGCTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTACAGATAGCAAACCACTTAGCCGATCACAATTCCGACAGAAACTgcggat            >  1:677540‑M2/11‑141 (MQ=255)                                                                                                                                    ttctgcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCAGAGCGCAAACCGCTTCGCCGTTCACAATTCCGACAGAAACTGCGGCTACCa        <  1:71595‑M2/137‑1 (MQ=255)                                                                                                                                         cCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGCGTCTTCAACGTATTCCAGATCGCAAACCGCTTCGCCGTTCACAATTCCGACAGAAACTGCGGCTACCATCCCt   >  1:601738‑M2/2‑143 (MQ=255)                                                                                                                                          cGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCCGCGTCTTCCACGTATTCCAGATCGCAAACCACTTAGCCGTTCACAATTACGCCAGAAACTGCGGCTACCCTCCCtt  >  1:7614/1‑143 (MQ=255)                                                                                                                                          |                                                                                                                                               AAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATTCCCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCAGATCGCAAACCGCTTCGCCGTTCACAATTCCGACAGAAACTGCGGCTACCATCCCTT  >  NC_000913/3815804‑3816083

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