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	34 [0]	[0] 34	[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)	33 (1.330)	31/264	0.0	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

ACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATTCCCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCG  >  NC_000913/3815717‑3815995                                                                                                                                              |                                                                                                                                          acacCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCc                                                                                                                                          <  2:387610‑M1/143‑2 (MQ=255)    ccTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGCTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGTCTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtc                                                                                                                                       <  1:359320‑M1/143‑5 (MQ=255)     cTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcc                                                                                                                                      >  1:388122‑M1/1‑138 (MQ=255)      tgcttgcaAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtccc                                                                                                                                     <  1:476166‑M1/143‑7 (MQ=255)            ctAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcac                                                                                                                               >  1:319114‑M1/1‑131 (MQ=255)                  gcAAATTCAGTAAGCTGGCGCTGATATGGTTTCATGACTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaat                                                                                                                         <  2:327402‑M1/143‑19 (MQ=255)                  gcAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCCGTTTGCCAGCTCCGCCTTCTGCcgtcccctccacttcaat                                                                                                                         >  1:128896‑M1/1‑125 (MQ=255)                     aaTTCAATAACCTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgat                                                                                                                      <  2:188483‑M1/143‑22 (MQ=255)                     aaTTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgat                                                                                                                      <  2:131463‑M1/143‑22 (MQ=255)                       ttCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTACAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttccaatgatg                                                                                                                    >  2:18407‑M1/1‑120 (MQ=255)                       ttCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcggcccctgcacttcaacaatgc                                                                                                                    >  2:139093‑M1/1‑120 (MQ=255)                              aaCTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtc                                                                                                             <  1:30554‑M1/143‑31 (MQ=255)                                cTGGGGCGGATATTGTGTCGTGCCTTCCCTCCTCATGTTACTTTTCTACAGGAATAAAAAAGGCGGCTCATCAGTCGCCTTAAAAATCTGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcctgtctt                                                                                                           <  2:40451‑M1/143‑33 (MQ=255)                                    ccctgCTATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggt                                                                                                       <  2:201136‑M1/141‑37 (MQ=255)                                        gATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatc                                                                                                   >  1:100949‑M1/1‑103 (MQ=255)                                                      ccTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactacgttcatgtc                                                                                     <  1:42227‑M1/143‑55 (MQ=255)                                                            ctccgcATCTTACTTTTCGACAGACAAAAACAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactacgttcatgtcggtctc                                                                               <  2:63386‑M1/143‑61 (MQ=255)                                                               ctcATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactacgttcatgtcgggctctgc                                                                            >  1:111974‑M1/1‑80 (MQ=255)                                                                tcATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcgcgtcttcggtcagcactacgttcatgtaggtctctgca                                                                           >  1:15572‑M1/1‑79 (MQ=255)                                                                         ttttCTACAGAC‑AAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactacgttcatgtcggtctctgcggcagagcctt                                                                 >  1:319153‑M1/1‑69 (MQ=255)                                                                          cgtatacagACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactacgttcatgtcggtctctgcggcagagtctt                                                                 <  1:84593‑M1/139‑75 (MQ=255)                                                                                   aCAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatccctacgcccatgtcggtctctgcggcagaggcttcaacgtgtt                                                        >  2:468507‑M1/1‑60 (MQ=255)                                                                                   aCAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactacgttcatgtcggtctctgcggcagagtcttcaacgtatt                                                        <  1:36170‑M1/143‑84 (MQ=255)                                                                                    cAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactacgttcatgtcggtctctgcggcagagtcttcaacgtattc                                                       <  1:397275‑M1/143‑85 (MQ=255)                                                                                     aaaaaaaaGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcagtcatcactacgttcatgtcggtctctgcggcagagtcatcaacgtattcc                                                      >  1:370304‑M1/1‑58 (MQ=255)                                                                                        aaaaaGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactacgtccatgtcggtctctgcgccccagtctccatcgtattccaga                                                   >  2:362129‑M1/1‑55 (MQ=255)                                                                                                     atcaGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccggctaaggtcaaccaacagtgtatgtcggtctctgcggcagagtcttcaacgtattctagatcgcaagcggctt                                      >  1:111209‑M1/1‑42 (MQ=255)                                                                                                      tcaGTCGCCTTAAAAATCAGTGTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtctgcgggcagcactacgttcatgtcggtctctgcggcagagtcttcaacgtattccagatcgcaaaccgcttc                                     <  1:278190‑M1/143‑103 (MQ=255)                                                                                                       caGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactaccttcatgtcggtctctgcggcagagtcttcaacgtattccagatcgcaaaccgcttcg                                    <  1:237689‑M1/143‑104 (MQ=255)                                                                                                                  ataatcAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactacgttcatgtcggtctctgcagcagagtcttcaacgtattccagatcgcaaaccgcttcgccgttcacaat                         <  1:70853‑M1/141‑115 (MQ=255)                                                                                                                   aaaTCAGTTTGCCAGCGCCGGCTTATGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcacgacggtcatgccggtctctgcggcagagtcttcaacgtattcccgatcgcaaaccgcttcgccgttcacaatt                        <  1:313482‑M1/143‑116 (MQ=255)                                                                                                                       cAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcataactacaatcatgtaggtctctgcgacagagccttcaacgtattccagatagcaaaccgcttcgccgttcacaattccaa                    >  1:476560‑M1/1‑24 (MQ=255)                                                                                                                                  cgccgcCTTCTCCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactacgttcatgtcggtctctgcggcagagtcttcaacgtattccagatcgcaacccgcttcgccattcacaattccgacagaaactgcg         <  2:355770‑M1/143‑131 (MQ=255)                                                                                                                                         ttCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactacgctcatgtcggtctctgcggcagagtcttcaacgtattccagatcgcaaaccgcttcgccgttcacaattccgacagaaactgcggctacca  >  1:260947‑M1/1‑6 (MQ=255)                                                                                                                                              |                                                                                                                                          ACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATTCCCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCG  >  NC_000913/3815717‑3815995

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