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	26 [0]	[0] 25	[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)	24 (0.640)	22/264	0.8	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

AAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATTCCCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCAGATCGCAAACCGCTTCGCCGTTCACAATTCCGACAGAAACTGCGGCTAC  >  NC_000913/3815807‑3816075                                                                                                                                       |                                                                                                                                       ttgctaagcgcaaattcaataaactggcgctgatatggtttcatgccttcgctcctcatcttacgtttctacagacaaaaaaaaggcgactcatcagtcgccttaaaaatcagtttgccagcgccgccttctgcCGTCCCCTg                                                                                                                                <  2:10439‑M2/9‑1 (MQ=255)           caaattcaatacaccggcgctgatatggtttcatgccttcgcccctcatcttacttttctacagacaaaaaaaaggcgactcatcagtcgccttaaaaatcagtttgccagcgccgccttctgcCGTCCCCTGCACTTCAatg                                                                                                                      <  1:522987‑M2/19‑1 (MQ=255)             aattcaataaactggcgctgatatggtttcatgccttcgctcctcatcttacttttctacagacaaaaaaaaggcgactcatcagtcgccttaaaaatcagtttgccagcgccgccttctgcCGTCCCCTGCACTTCAatgat                                                                                                                    <  1:651028‑M2/21‑1 (MQ=255)                tcaacacactggtgctgatatggtgtcatgccttcgctcctcatcttacttctctacaggcaaaaaacaggcgactcatcagtcgccttaaaaatcagttggccagcgccgccttctgcCGTCCCCTGCACTTCAATGATgcg                                                                                                                 <  2:547561‑M2/24‑1 (MQ=255)                  aataaactggcgctgatatggtttcatgccttcgctcctcatcttacttttcgacagacaaaaaaaaggcgactcatcagtcgccttaaaaatcagtttgccagcgccgccttctgcCGTCCCCTGCACTTCAATGATGCGcc                                                                                                               <  1:649586‑M2/26‑1 (MQ=255)                   ataaactggcgctgatatggtttcatgccttcgctcctcatcttacttttctacagacaaaaaaaaggcgactcatcagtcgccttaaaaatcagtttgccagcgccgccttctgcCTTCCCCTGCACTTCCATGATGCAccc                                                                                                              >  1:46021‑M2/117‑143 (MQ=255)                   ataaactggcgctgatatggtttcatgccttcgctcctcatcttacttttctacagacaaaaaaaaggcgactcatcagtcgccttaaaaatcagtttgccagcgccgccttcagcAGTCCCCAGAACTTCAATGATGCGccc                                                                                                              >  2:616091‑M2/117‑143 (MQ=255)                          ggcgctgatatggtttcaggcctgcgctcctcatcttactttgctacagacaaaaaaaaggcgactcatcagtcgccttaaaaatcagtttgccagcgccgccttctgcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCg                                                                                                       <  1:138890‑M2/34‑1 (MQ=255)                                        ttcatgccttcgcccctcatcttacttttctacagacaaaaaaaaggcgactcatcagtcgccttaaaaatcagtttgccagcgccgccttctgcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGtt                                                                                         <  2:262617‑M2/48‑1 (MQ=255)                                             gccttcgctcctcatcttacttttctacagacaaaaaaaaggcgactcatcagtcgccttaaaaatcagtttgccagcgccgccttctgcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGt                                                                                    >  1:464436‑M2/91‑143 (MQ=255)                                                                 ttttctacagacaaaaaaaaggcgactcatcagtcgccttaaaaatcagtttgccagcgccgccttctgcCCTCCCCTGCACTTCAATGCTGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCCCTGCGGCCGAGTCt                                                                >  2:444869‑M2/71‑143 (MQ=255)                                                                      tacagacaaaaaaaaggcgactcatcagtcgccttaaaaatcagtttgccagcgccgccttctgcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAAc                                                           >  1:536633‑M2/66‑143 (MQ=255)                                                                       acagacagaaaacaggcgactcatcagtcgccctaaaaaccagtttgccagcgccgccttctgcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACg                                                          <  1:11719‑M2/79‑1 (MQ=255)                                                                         agacaaaaaaaaggcgactcatcagtcgccttaaaaatcagtttgccagcgccgccttctgcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTa                                                        <  1:149731‑M2/81‑1 (MQ=255)                                                                          gacaaaaaaaaggcgactcatcagtcgccttaaaaatcagtttgccagcgccgccttctgcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTACTCAACGTAt                                                       >  2:92765‑M2/62‑143 (MQ=255)                                                                             aaaaaaaaggcgactcatcagtcgccttaaaaatcagtttgccagcgccgccttctgcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGGCGCcaacgtcgccg                                                    >  2:488758‑M2/59‑138 (MQ=255)                                                                                 aaagggcgactcatcagtcgccttacaaatcagtttgccagcgccgccttctgcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCAGAt                                                <  2:462722‑M2/89‑1 (MQ=255)                                                                                 aaaaggcgactcatcagtcgccttaaaaatcagtttgccagcgccgccttctgcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCAGAt                                                <  1:391556‑M2/89‑1 (MQ=255)                                                                                        gactcatcagtcgccttaaaaatcagtttgccagcgccgccttctgcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTCTTCCAGATCgcaacc                                         >  1:513323‑M2/48‑141 (MQ=255)                                                                                         actcatcagtcgccttaaaaatcagtttgccagcgccgccttctgcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCAGATCGCAAAcc                                        >  1:250703‑M2/47‑143 (MQ=255)                                                                                           tcatcagtcgccttaaaaatcagtttgccagcgccgccttctccCCGCCGCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTGCGTTCAGGTCGGTCTCTGCGGCAGAGCCTGCAACGTATTCCGGATCGCGAACCGc                                      >  2:398463‑M2/45‑143 (MQ=255)                                                                                                      cttaaaaatcagtttgccagcgccgccttctgcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCAGATCGCACACCGCTTCGCCGTTca                           >  1:97973‑M2/34‑143 (MQ=255)                                                                                                         aaaaatcaggttgccagcgccgccttctgcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCAGATCGCAAACCGCTTCGCCGTTCACaa                        <  1:492306‑M2/113‑1 (MQ=255)                                                                                                            agtccgtttgccagcgccgcctgctgcCGTCCCCTGCACTTCAATGACGCGCCCGTCTTCGGTCATCACTCCGCTCATGTCGGTCTCTGCGGCAGCGTCTTCAACGTATTCCAGATCGCAAACCGCTTCGCCGTTCACAATTc                     <  2:645305‑M2/116‑1 (MQ=255)                                                                                                                               ccttctgcCGTCCCCTGCACTGCATTGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCAGATCGCAAACCGCTTCGCCGTTCACAATTCCGACAGAAACTGCGGCTAc  <  1:450859‑M2/135‑1 (MQ=255)                                                                                                                                       |                                                                                                                                       AAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATTCCCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCAGATCGCAAACCGCTTCGCCGTTCACAATTCCGACAGAAACTGCGGCTAC  >  NC_000913/3815807‑3816075

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