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	23 [0]	[0] 23	[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)	21 (0.880)	18/492	0.3	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

CGTAAACTCGCCAAACTTTAACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATTCCCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCAGATCGCAAACCGCTTCGCCGTTCACAATTCCGACAGAAACTGCGGCTACCATCCCTTTCATCGGATTGGTTTTCAGCTTGCCGTTTTCCACCAGCTTCTGTAGCGCATCTACCAGCGCCACGCAGGCACCCGTAATCGACGCGGTACGCGTGCCACCATCAGCCTGAAGCACGTCGCAGTCCAGCG  >  NC_000913/3815697‑3816212                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                                                aacgtccttacattgccagatcgcgcccggtattaaacagcccggcgttgaataaataggggcttttgcgcccggatttcagcggaaactcgccaaactttaacacctgcttgctaagcgcaaattcaataaactggcgctgatatggtttcatgccttcgctcctcatcttacttttctacagacaaaaaaaaggcgactcatcagtcgccttaaaaatcagtttgccagcgccgccttctgcCGTCCCCTGCACTTCAATGATGCGCCCGTCt                                                                                                                                                                                                                                                   <  1:107011‑M2/31‑1 (MQ=255)                                     caagcccggcgttgaagaaataggggcttttgcgcccggatttcagcttaaactcgccaaactttaacacctgcttgctaagcgcaaattcaataaactggcgctgatatggtttcatgccttcgctcctcatcttacttttctacagacaaaaaaaaggcgactcatcagtcgccttaaaaatcagtttgccagcgccgccttctgcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCa                                                                                                                                                                                                               <  1:71325‑M2/67‑1 (MQ=255)                                     acagcccggcgttgaagaaataggggcttttgcgcccggatttcagcgtaaactcgccaaactttaacacctgcttgctaagcgcaaattcaataaactggcgctgatatggtttcatgccttcgctcctcatcttacttttctacagacaaaaaaaaggcgactcatcagtcgccttaaaaatcagtttgccagcgccgccttctgcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCa                                                                                                                                                                                                               <  1:164046‑M2/67‑1 (MQ=255)                                       agcccggcgttgaagaaataggggcttttgcgcccggatttcagcgtaaactcgccaaactttaacacctgcttgctaagcgcaaattcaataaactggcgctgatatggtttcatgccttcgctcctcatcttacttttctacagacaaaaaaaaggcgactcatcagtcgccttaaaaatcagtttgccagcgccgccttctgcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCaga                                                                                                                                                                                                             <  1:106420‑M2/69‑1 (MQ=255)                                       agcccggcgttgaagaaataggggcttttgcgcccggatttcagcgtaaactcgccaaactttaacacctgcttgctaagcgcaaattcaataaactggcgctgatatggtttcatgccttcgctcctcatcttacttttctacagacaaaaaaaaggcgactcatcagtcgccttaaaaatcagtttgccagcgccgccttctgcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCag                                                                                                                                                                                                              >  2:74361‑M2/207‑274 (MQ=255)                                             gcgttgaagaaataggggcttttgcgcccggatttcagcgtaaactcgccaaactttaacacctgcttgctaagcgcaaattcaataaactggcgctgatatggtttcatgccttcgctcctcatcttacttttctacagacaaaaaaaaggcgactcatcagtcgccttaaaaatcagtttgccagcgccgccttctgcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCtt                                                                                                                                                                                                        <  1:1375‑M2/74‑1 (MQ=255)                                              ctttgaagaaatgggggcttttgctaccggatttcaacgaaaacacgccaaactttaacacctgcttgctaagcgcaaattcaataaactggcgctgatatggtttcatgccttcgctcctcatcttacttttctaccgacaaaaacaaggcgactcatcagtcgccttaaaaatcagtttgccagcgccgccttctgcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCa                                                                                                                                                                                                      <  2:143598‑M2/76‑1 (MQ=255)                                                                 tttgcgcccggatttcagcgtaaactcgccaaactttaacacctgcttgctaagcgcaaattcaatacactggcgctgatatggtttcatgccttcgctcctcatcttacttttctacagacaaaaaaaaggcgactcatcagtcgccttaaaaatcaatttgccagcgccgccttctgcAGTCCCCTGCACTTCAATGCTGCGCCCGTCTTCGGTCATCACTACGGTCATGTAGGTCTCTGCGGTAGAGTCTTCAacgtattccatcacgaaaa                                                                                                                                                                                   >  1:135590‑M2/181‑266 (MQ=255)                                                                      gcccggatttcagcgtaaactcgccaaactttaacaccctcttgctaagcgcaaatgcaataaactggcgctgatatggtttcatgccttcgctcctcatcttacttttctacagacaaaaaaaaggcgactcagcagtcgccttaaaaatcagtttgccagcgccgccttctgcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCGTCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCAGATCGCAAACCGCt                                                                                                                                                                              <  2:155122‑M2/100‑1 (MQ=255)                                                                                                             gcttgctaagcgcaaattcaataaactggcgctgatatggtttcatgccttcgctcctcatcttacttttctacagacaaaaaaaaggcgactcatcagtcgccttaaaaatcagtttgccagcgccgccttctgcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCAGATCGCAAACCGCTTCGCCGTTCACAATTCCGACAGAAACTGCGGCTACCATc                                                                                                                                       >  2:27865‑M2/137‑275 (MQ=255)                                                                                                                                          cgctgatatggtttcatgccttcgctcctcatcttacttttctacagacaaaaaaaaggcgactcatcagtcgccttaaaaatcagtttgccagcgccgccttctgcCGTCCCCTGCACTTCAATCATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACTTATTCCAGATCGCAACCCGCTTCGCAGTTCACAATTCCGACAGAAACGGCGGCTAACAACCCTTTCAACGACTTGGTTTTCACCTTGcc                                                                                                          >  2:185727‑M2/108‑275 (MQ=255)                                                                                                                                                  tgttttcatgccttcgctcctcatcttacttttctacagacaaaaaaaaggcgactcatcagtcgccttaaaaatcagtttgccagcgccgccttctgcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCAGATCGCAAACCGCTTCGCCGTTCACAATTCCGACAGAAACTGCGGCTACCATCCCTTTCATCGGATTGGTTTTCAGCTTGCCGTTTTcca                                                                                                  <  1:194196‑M2/176‑1 (MQ=255)                                                                                                                                                      ttcatgccttcgctcctcatcttacttttctaccgacaaaaaaaaggcgactcatccgtcgccttaaaaatcagtttgccagcgccgccttctgcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCAGATCGCAAACCGCTTCGCCGTTCACAATTCCGACAGAAACTGCGGCTACCATCCCTTTCATCGGATTGGTTTTCAGCTTGCCGTTTTCCACCAg                                                                                              <  2:33736‑M2/180‑1 (MQ=255)                                                                                                                                                                               ttttctacagacaaaaaaaaggcgactcatcagtcgccttaaaaatcagtttgccagcgccgccttctgcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCAGATCGCAAACCGCTTCGCCGTTCACAATTCCGACAGAAACTGCGGCTACCATCCCTTTCATCGGATTGGTTTTCAGCTTCCCGTTTTCCACCAGCTTCTGTAGCGCATCtcccacccc                                                                      >  2:101057‑M2/71‑270 (MQ=255)                                                                                                                                                                                  tctacagacaaaaaaaaggcgactcatcagtcgccttaaaaatcagtttgccagcgccgccttctgcCGTCCCTTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCAGATCGCAAACCGCTTCGCCGTTCACAATTCCGACAGAAACTGCGGCTACCATCCCTTTCATCGGATTGGTTTTCAGCTTGCCGTTTTCCACCAGCTTCTGTCGCCCATCTCCCAGcgcccc                                                                   >  2:96874‑M2/68‑272 (MQ=255)                                                                                                                                                                                         acaaaaaaaaggcgactcatcagtcgccttaaaaatcagtttgccagcgccgccttctgcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCAGATCGCAAACCGCTTCGCCGTTCACAATTCCGACAGAAACTGCGGCTACCATCCCTTTCATCGGCTTGGTTTTCAGCTTGCCGTTTTCCACCAGCTTCTGTAGCGCATTTACCAGCGCCACCCAggccc                                                           >  2:224059‑M2/61‑273 (MQ=255)                                                                                                                                                                                                        ctcatcaggcgccgtaaaaagccgtttgccaccgcccccgtctgcCGTACCCTGCACTTCAATGTTGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCAGATCGCAAACCGCTTCGCCGTTCACAATTCCGACAGAAACTGCGGCTACCATCCCTTTCATCGGATTGGTTTTCAGCTTGCCGTTTTCCACCAGCTTCTGTAGCGCATCTACCAGCGCCACGCAGGCACCCGTAATCGAcg                                               <  1:95146‑M2/227‑1 (MQ=255)                                                                                                                                                                                                                  gccttaaaaatcagtttgccagcgccgccttctgcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCAGATCGCAAACCGCTTCGCCGTTCACAATTCCGACAGAAACTGCGGCTACCATCCCTTTCATCGGATTGGTTTTCAGCTTGCCGTTTTCCACCAGCTTCTGTAGCGCATCTACCAGCGCCACGCAGGCACCCGTAATCGACGCGGTACGCGTGcc                                  >  2:170315‑M2/36‑275 (MQ=255)                                                                                                                                                                                                                  gccttaaaaatcagtttgccagcgccgccttctgcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCAGATCGCAAACCGCTTCGCCGTTCACAATTCCGACAGAAACTGCGGCTACCATCCCTTTCATCGGATTGGTTTTCAGCTTGCCGTTTTCCACCAGCTTCTGTAGCGCATCTACCAGCGCCACGCAGGCACCCGTAATCGACGCGGTACGCGTGcc                                  >  2:230522‑M2/36‑275 (MQ=255)                                                                                                                                                                                                                  gccttaaaaatcagtttgccagcgccgccttctgcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCAGATCGCAAACCGCTTCGCCGTTCACAATTCCGACAGAAACTGCGGCTACCATCCCTTTCATCGGATTGGTTTTCAGCTTGCCGTTTTCCACCAACTTCTGTAGCGCATCTACCAGCGCCACCCCGGCACCCGTAATCGACGCGGTACGCGGGcc                                  >  2:106046‑M2/36‑275 (MQ=255)                                                                                                                                                                                                                                   gccagcgccgccttctgcCGTCCCCTGCTCTTCAATGCTGCGCCCGTCTTCGGTCATCACTCCGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCAGATCGCAAACCGCTTCGCCGTTCACAATTCCGACAGAAACTGCGGCTACCAGCCCTTTCATCGGATTGGTTTTCAGCTTGCCGTTTTCCACCAGCTTCTGTAGCGCATCTACCAGCGCCACGCAGGCACCCGTAATCGACGCGGTACGCGTGCCACCATCAGCCTGAAGCa                 <  1:74362‑M2/257‑1 (MQ=255)                                                                                                                                                                                                                                          cccccttcggcGGTCCCCTGCACGGGCGCGCTGCGCCCGTCTTGGGTCATCACTACGTTCATGTCGGTCTCTGCGGCGGAGTCTTCAACGTATTCCAGATCGCAAACCGGTTCGCCGTTCACAATTCCGACAGAGACTGCGGCTACCATCCCTTTCATCGGATTGGTTTTCAGCTTGCCGTTTTCCACCAGCTTCTGTAGCGCATCTACCAGCGCCACGCAGGCACCCGTAATCGACGCGGTACGCGTGCCACCATCAGCCTGAAGCACGTCGCa          <  2:126552‑M2/264‑1 (MQ=255)                                                                                                                                                                                                                                                  tgcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCAGATCGCAAACCGCTTCGCCGTTCACAATTCCGACAGAAACTGCGGCTACCATCCCTTTCATCGGATTGGTTTTCAGCTTGCCGTTTTCCACCAGCTTCTGTAGCGCATCTACCAGCGCCACGCAGGCACCCGTAATCGACGCGGTACGCGTGCCACCATCAGCCTGAAGCACGTCGCAGTCCAGCg  <  1:230524‑M2/272‑1 (MQ=255)                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                                                CGTAAACTCGCCAAACTTTAACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATTCCCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCAGATCGCAAACCGCTTCGCCGTTCACAATTCCGACAGAAACTGCGGCTACCATCCCTTTCATCGGATTGGTTTTCAGCTTGCCGTTTTCCACCAGCTTCTGTAGCGCATCTACCAGCGCCACGCAGGCACCCGTAATCGACGCGGTACGCGTGCCACCATCAGCCTGAAGCACGTCGCAGTCCAGCG  >  NC_000913/3815697‑3816212

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