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

ACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATTCCCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTC  >  NC_000913/3815719‑3816000                                                                                                                                            |                                                                                                                                               acCTACTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgt                                                                                                                                             >  1:79741‑M1/1‑140 (MQ=255)               cgcAAATTCACTAACCTGGCCCGGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGCCAAAAAAAATGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaa                                                                                                                               <  1:324936‑M1/143‑18 (MQ=255)                gcAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaat                                                                                                                              <  1:481604‑M1/143‑19 (MQ=255)                gcAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaat                                                                                                                              >  1:259857‑M1/1‑125 (MQ=255)                 cAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatg                                                                                                                             <  1:325450‑M1/143‑20 (MQ=255)                     ttCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcccttcaatgatgc                                                                                                                         >  1:479785‑M1/1‑120 (MQ=255)                          tAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcccttagatgatgcgaccg                                                                                                                    >  1:156468‑M1/1‑115 (MQ=255)                           aaaCTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgt                                                                                                                   <  2:428179‑M1/143‑30 (MQ=255)                                gccgctgATATGGTTTCATGCCTTCGCTCCTCATCTTACGTTTCTACAGACAAAAAAAATGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcg                                                                                                              <  1:207567‑M1/141‑35 (MQ=255)                                  cgcTGATATGGGTTCATGCCTTCGCTCCTCTTCTTACTTTTCTACAGACAAAAAGAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggt                                                                                                            <  2:120922‑M1/143‑37 (MQ=255)                                     tGATATGGTTTCATGCCTTCGCTCCTCATCTTACCTTTCTACAGACAAAAAAAGGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcat                                                                                                         <  2:245511‑M1/143‑40 (MQ=255)                                     agaTATGGTTTCATGCCTCCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcat                                                                                                         <  1:436797‑M1/142‑40 (MQ=255)                                                          ctcctcATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactacgttcatgtcggtctc                                                                                    >  1:146188‑M1/1‑83 (MQ=255)                                                                  cTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtctgcgggcaccactacgttcatgcaggtctctgccgcag                                                                            >  2:69048‑M1/1‑75 (MQ=255)                                                                     aCTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactacgttcatgtcggtctctgcggcagagt                                                                         >  1:366899‑M1/1‑72 (MQ=255)                                                                        tttCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactacgttcatgtcggtctctgcggcagagtctt                                                                      <  1:66223‑M1/143‑75 (MQ=255)                                                                           cTACAGCCGAAAAACCGGCGACTCATCAGTCGCCTTAAAAATCAGCTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactacgttcatgtcggtctctgcggcagagtcttcaa                                                                   <  1:313113‑M1/143‑78 (MQ=255)                                                                           cTACAGACAAAAAAAAGGCGACTCAGCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactacgttcatgtcggtctctgcggcagagtcttcaa                                                                   <  1:368026‑M1/143‑78 (MQ=255)                                                                                gACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtctccggtcatcactacggtcatttcggtctctgcggcagagtcttcaacgtat                                                              <  1:501161‑M1/143‑83 (MQ=255)                                                                                                   atcaGTCGCCTTAAAAATCAGTTTGCCAGCGCAGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactccgttcatgtcggtctctgacgccgagtctccaacgtatccccgatcgcaaaccgcct                                           >  2:273768‑M1/1‑42 (MQ=255)                                                                                                               aaaaaTCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactacgttcatgtcggtctctgcggcagagtcttcaacgtattccagatcgcaaaccgcttcgccgttcacaa                               <  1:172126‑M1/143‑114 (MQ=255)                                                                                                                 aaaTCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactacgttcatgtcggtctctgcggcagagtcttcaacgtattccagatcgcaaaccgcttcgccgttcacaagt                             >  2:118982‑M1/1‑28 (MQ=255)                                                                                                                 aaaTCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacctcaatgatgcgcccgcctccggtcatcactaggttcatgtcggtctctgcggccgagtcttcaacgtattccagatcgcggaccgcttcgccgttcccaatt                             >  1:9978‑M1/1‑28 (MQ=255)                                                                                                                     cAGTGTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactacgttcatgtcgggctctgcggcagagtcttcaacgtattccagatcgcaaaccgcttcgccgttcacaattccga                         <  2:87897‑M1/143‑120 (MQ=255)                                                                                                                      aGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgaggcgcccgtcttcggtcatcactacgttcatgtcggtctctgcggcagagtcttcaacgtattccagatcgcaaaccgcttcgccgttcacaattccgac                        <  2:286492‑M1/143‑121 (MQ=255)                                                                                                                             cAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcctcactacgttcatgtcggtctctgcggccgagtcttccacgtatcccacatcgcaaaccgcttcgccgatcacaattccgacagaaact                 >  1:314278‑M1/1‑16 (MQ=255)                                                                                                                                 gccgccTTCTGCcgtcccctgcactccaatgatgcgcccgtcttcggtcatcactacgttcatgtcggtctcggctgcagagtcttcaacgtattccagttcgcaaaccgcttcgccgttcacaattccgacagaaactgcgg             <  1:442476‑M1/143‑132 (MQ=255)                                                                                                                                       ttCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactacgttcatgtcggtctctgcggcagcgtcttcaacgtattccagatcgcaaaccgcttcgccgttcacaattccgacagaaactgcggctacca       >  1:477111‑M1/1‑6 (MQ=255)                                                                                                                                        tCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactacgttcatgtcggtctctgcggcagagtcttcaacgtattccagatcgcaaaccgcttcgccgttcacaattccgacagaaactgcggctacaat      >  1:345615‑M1/1‑5 (MQ=255)                                                                                                                                            ccgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactacgttcatgtcggtctctgcggcagagtcttcaacgtattccagatcgcaaaccgcttcgccgttcacaattcggacagaacctgcggctaccatcccc  >  2:47570‑M1/1‑1 (MQ=255)                                                                                                                                            |                                                                                                                                               ACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATTCCCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTC  >  NC_000913/3815719‑3816000

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