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	32 [0]	[0] 33	[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)	29 (0.600)	25/262	1.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

ACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATTCCCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGT  >  NC_000913/3815717‑3815997                                                                                                                                              |                                                                                                                                            aaacctGCGTGCTAAGCGCAAATTCAATAATCTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCc                                                                                                                                            <  2:167364‑M1/141‑2 (MQ=255)     cTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcc                                                                                                                                        >  2:928633‑M1/1‑138 (MQ=255)             tAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTCCTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAATTTACCAGCGCCGCCTTCTGCcgtccactgcact                                                                                                                                >  1:480827‑M1/1‑130 (MQ=255)                 cgaaaattcAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaa                                                                                                                            <  2:291457‑M1/140‑18 (MQ=255)                  gcAAATTCAATAAGCTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaat                                                                                                                           <  1:560531‑M1/143‑19 (MQ=255)                     aaTTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTAGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgat                                                                                                                        >  2:911256‑M1/1‑122 (MQ=255)                          aaTAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcc                                                                                                                   >  1:438777‑M1/1‑117 (MQ=255)                          aaTAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCAACTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcc                                                                                                                   <  1:202467‑M1/143‑27 (MQ=255)                                cTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtctt                                                                                                             <  2:530269‑M1/143‑33 (MQ=255)                                 tGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttc                                                                                                            <  1:334998‑M1/143‑34 (MQ=255)                                  ggCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAACAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcg                                                                                                           <  2:215078‑M1/143‑35 (MQ=255)                                       tGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcat                                                                                                      <  2:239581‑M1/143‑40 (MQ=255)                                       tGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcat                                                                                                      <  1:48481‑M1/143‑40 (MQ=255)                                                      ccTTCGCTCCTCAGCTTACTTTTCTGCCGACAAAAAAAAGCCGACGCGTCTGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactacgttcatgtc                                                                                       <  1:353887‑M1/143‑55 (MQ=255)                                                        ttCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactacgctcatgtcgg                                                                                     >  2:545328‑M1/1‑87 (MQ=255)                                                         tCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcctcactacgttcatgtcggt                                                                                    >  2:889611‑M1/1‑86 (MQ=255)                                                                     tttctttTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactacgttcatgtcggtctctgcggcaga                                                                        <  2:361134‑M1/140‑70 (MQ=255)                                                                     ttACGTTTCTACAGACCATAAAAAGGCGACTCATCAGTCGCCTTAAAACTCAGTTTGCCAGCGCCGCCTGCTGCcgtcccctgcacttcaatgatgcgcccgccttcggtcatcactacgttcatgtcggtctctgcggcaga                                                                        <  2:253979‑M1/143‑70 (MQ=255)                                                                                    cAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactacgttcatgtcggtctctgcggcagagtcttcaacgtattc                                                         <  1:748837‑M1/143‑85 (MQ=255)                                                                                                 aCTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactacgttcatgtcggtctctgcggcagagtctcaaacgtattccagatctcacacg                                            >  1:668173‑M1/1‑46 (MQ=255)                                                                                                     atcaGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactacgttcatgtcggtctctgcggcagcgcgttcaacgtattccagatcgcaaaccgctt                                        >  2:462391‑M1/1‑42 (MQ=255)                                                                                                     atcaGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactacgttcatgtcggtctctgcggcagagtcttcaacgtattccagatcgcaaaccgctt                                        >  2:796998‑M1/1‑42 (MQ=255)                                                                                                               ttAAAAGTCAGTTTGCCAGCGCCGCCTTCTGCcgccccctgcacttcaatgatgcgcccgtcttcggtcatcactacgttcatgtcggtctctgcggcagagtcttcaacgtattccagatcgcaaaccgcttcgccgttcac                              <  1:754112‑M1/143‑112 (MQ=255)                                                                                                                    aaTCAGTTTGCCAGCGCCGCCTTCTGCcgccccctgcacttcaatgatgcgcccgtcttcggtcatcacgacgttcatgtcggtctctgcggcagagtcttcaacgtattccagatcgcaaaccgcttcgccgttcacaattc                         <  2:438777‑M1/143‑117 (MQ=255)                                                                                                                     aTCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactacgttcatgtcggtctctgcggcagcgtcttcaacgtattccagctcgcacaccgcttcgccgttcacaattcc                        >  1:198203‑M1/1‑26 (MQ=255)                                                                                                                          tttGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactacgttcatgtcggtctctgcggcagagtcttcaacgtattccagatcgcaaaccgcttcgccattcacaattccgacag                   >  1:119676‑M1/1‑21 (MQ=255)                                                                                                                                aGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactacgttcatgtcggtctctgcggcagagtcttcaacgtattccagatcgcaaaccgcttcgccgttcacaattccgacagaaactg             >  2:345301‑M1/1‑15 (MQ=255)                                                                                                                                aGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactacgttcatgtcggtctctgcggcagagtcttcaacgtattccagatcgcaaaccgcttcgccgttcacaattccgacagaaactg             >  1:566461‑M1/1‑15 (MQ=255)                                                                                                                                 gcgcCGCCTTCTGCcgtcccctgcactgcaatgatgcgcccgtcttcggtcatcactacgttcatgtcggtctctgcggcagagtcttcaacgtattccagatcgcaaaccgcttcgccgttcacaattccgacagaaactgc            <  2:480827‑M1/143‑130 (MQ=255)                                                                                                                                       ccTTCTGCcgtcccctgcacttcaaggatgcgcccgtctgcggtcatcactacgttcatgtcggtctctgcggcagagtcttcaacgtattccagatcgcaaaccgcttcgccgttcacaattccgacagaaactgcggctac      <  1:789711‑M1/143‑136 (MQ=255)                                                                                                                                         ttCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactacgttcatgtcggtctctgcggcagagtcttcaacgtattccagatcgcaaaccgcttcgccgttcacaattccgacagaaactgcggctacca    <  1:796998‑M1/143‑138 (MQ=255)                                                                                                                                           cTGCcgtcccctgctcctccccgatgcgcccgtcttcggtcagcactacgttcatgtcggtctctgcggcagagtcttcaacgtattccagatcgcaaaccgcttcgccgttcacaattccgacagaaactgcggctaccatc  <  2:452243‑M1/143‑140 (MQ=255)                                                                                                                                              |                                                                                                                                            ACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATTCCCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGT  >  NC_000913/3815717‑3815997

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