Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
MC JC	NC_002947	1,515,769	Δ33 bp	100%	pseudogene (757‑789/1851 nt)	PP_1325 ←	putative lipoprotein

Missing coverage evidence...
			seq id	start	end	size	←reads	reads→	gene	description
*	*	÷	NC_002947	1515769	1515801	33	20 [0]	[0] 20	PP_1325	putative lipoprotein

New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_002947	1515802 =	0 (0.000)	12 (1.040)	9/160	NT	100%	pseudogene (756/1851 nt)	PP_1325	putative lipoprotein
	?	NC_002947	= 1515768	0 (0.000)					pseudogene (790/1851 nt)	PP_1325	putative lipoprotein

TGCTGGGCCTGGAAGTGGGCGGCCATGAAGCCATCGCGCAGGGCGCGGGCTACACCGGCCAGCGGGCCTTCCTGAGGCAGCAGCAGGGCGATCTTGGTCAGCGGCTGGCTGGCCAGCTCCTTGAGCTTGGTCAGCGGCTGGCTGGCCAGCTCCTTGAGCTTGGTCAGGGCCAGCGGCAGTTGCTGGGCGGCAGGGTGGTCCGGGTGCTGTTTGACCCAGGCTTCGATGGCCGCTTGCTGCTGCTCCAGGGTGC  >  NC_002947/1515677‑1515929                                                                                                                              |                                                                                                                                aacacctgtaccgcaggcgctggctggccgctttgctgggcctggaagtgggcggccatgaagccatcgcgcagggcgcgggctacaccggccagcgggccttcctgaggcagcagcagggcgatCTTGGTCAGCg                                                                                                                       <  1:113201‑M1/11‑1 (MQ=255)   cacctgtaccgcaggcgctggctggccgctttgctgggcctggaagtgggcggccatgaagccatcgcgcagggcgcgggctacaccggccagcgggccttcctgaggcagcagcagggcgatCTTGGTCAGCggc                                                                                                                     >  1:376169‑M1/124‑136 (MQ=255)                 gcgctggctggccgctttgctgggcctggaagtgggcggccatgaagccatcgcgcagggcgcgggctacaccggccagcgggccttcctgaggcagcagcagggcgatCTTGGTCAGCGGCTGGCTGGCCAGCTc                                                                                                       <  2:76676‑M1/27‑1 (MQ=255)                 gcgctggctggccgctttgctgggcctggaagtgggcggccatgaagccatcgcgcagggcgcgggctacaccggccagcgggccttcctgaggcagcagcagggcgatCTTGGTCAGCGGCTGGCTGGCCAGCTc                                                                                                       <  1:192692‑M1/27‑1 (MQ=255)                 gcgctggctggccgctttgctgggcctggaagtgggcggccatgaagccatcgcgcagggcgagggctacaccggccagcgggccttcctgaggcagcagcagggcgatCTTGGTCAGCGGCTGGCTGGCCAGCTc                                                                                                       <  1:142309‑M1/27‑1 (MQ=255)                      ggctgggcgctttgctgggcctggaagtgggcggccatgaagccatcgcgcagggcgcgggctacaccggccagcgggccttcctgaggctgcagcagggcgatCTTGGTCAGCGGCTGGCTGGCCAGCTCCTTGa                                                                                                  >  1:478323‑M1/105‑136 (MQ=255)                      ggctggccgctttgctgggcctggaagtgggcggccatgaagccatcgcgcagggcgcgggctacaccggccagcgggccttcctgaggcagcagcagggcgatCTTGGTCAGCGGCTGGCTGGCCAGCTCCTTGa                                                                                                  >  2:387058‑M1/105‑136 (MQ=255)                                  tgctgggcctggaagtgggcggccatgaagccatcgcgcagggcgcgggctacaccggccagcgggccttcctgaggcagcagcagggcgatCTGGGTCAGCGGCTGGCTGGCCAGCTCCTTGAGCTTGGTCAggg                                                                                      >  1:8751‑M1/93‑136 (MQ=255)                                            ggaagtgggcggccatgaagccatcgcgcagggcgcgggctacaccggccagcgggccttcctgaggcagcagcagggcgatCTTGGTCAGCGGCTGGCTGGCCAGCTCCTTGAGCTTGGTCAGGGCCAGCGGCAg                                                                            <  1:192873‑M1/54‑1 (MQ=255)                                            ggaagtgggcggccatgaagccatcgcgcagggcgcgggctacaccggccagcgggccttcctgaggcagcagcagggcgatCTTGGTCAGCGGCTGGCTGGCCAGCTCCTTGAGCTTGGTCAGGGCCAGCGGCAg                                                                            <  1:69740‑M1/54‑1 (MQ=255)                                                      ggccatgaagccatcgcgcagggcgcgggctacaccggccagcgggccttcctgaggcagcagcagggcgatCTTGGTCAGCGGCTGGCTGGCCAGCTCCTTGAGCTTGGTCAGGGCCAGCGGCAGTTGCTGggcg                                                                  <  1:82747‑M1/64‑1 (MQ=255)                                                      ggccatgaagccatcgcgcagggcgcgggctacaccggccagcgggccttcctgaggcagcagcagggcgatCTTGGTCAGCGGCTGGCTGGCCAGCTCCTTGAGCTTGGTCAGGGCCAGCGGCAGTTGCTGggcg                                                                  <  1:1820‑M1/64‑1 (MQ=255)                                                          atgaagccatcgcgcagggcgcgggctacaccggccagcgggccttcctgaggcagcagcagggcgatCTTGGTCAGCGGCTGGCTGGCCAGCTCCTTGAGCTTGGTCAGGGCCAGCGGCAGTTGCTGGGCGGCAg                                                              >  2:232202‑M1/69‑136 (MQ=255)                                                                  agcgcgcagggcgcgggctacaccggccagcgggccttcctgcggcagcagcagggcgatCTTGGTCAGCGGCTGGCTGGCCAGCTCCTTGAGCTTGGTCAGGGCCAGCGGCAGTTGCTGGGCGGCAGGGTGGTcc                                                      >  1:334099‑M1/61‑136 (MQ=255)                                                                                              agcgggccttcctgaggcagcagcagggcgatCTTGGTCAGCGGCTGGCTGGCCAGCTCCTTGAGCTTGGTCAGGGCCAGCGGCAGTTGCTGGGCGGCAGGGTGGTCCGGGTGCTGTTTGACCCAGGCTTCGATgg                          <  2:383892‑M1/104‑1 (MQ=255)                                                                                              agcgggccttcctgaggcagcagcagggcgatCTTGGTCAGCGGCTGGCTGGCCAGCTCCTTGAGCTTGGTCAGGGCCAGCGGCAGTTGCTGGGCGGCAGGGTGGTCCGGGTGCTGTTTGACCCAGGCTTCGATgg                          <  2:392126‑M1/104‑1 (MQ=255)                                                                                               gcgggccttcctgaggcagcagcagggcgatCTTGGTCAGCGGCTGGCTGGCCAGCTCCTTGAGCTTGGTCAGGGCCAGCGGCAGTTGCGGGGCGGCAGGGTGGTCCGGGTGCTGTTTGACCCAGGCTTCGATGGc                         <  1:254630‑M1/105‑1 (MQ=255)                                                                                                                     cagggcgatCTTGGTCAGCGGCTGGCTGGCCAGCTCCTTGAGCTTGGTCAGGGCCAGCGGCAGTTGCTGGGCGGCAGGGTGGTCCGGGTGCTGTTTGACCCAGGCTTCGATGGCCGCTtg                   <  1:22808‑M1/111‑1 (MQ=255)                                                                                                                     cagggcgatCTTGGTCAGCGGCTGGCTGGCCAGCTCCTTGAGCTTGGTCAGGGCCAGCGGCAGTTGCTGGGCGGCAGGGTGGTCCGGGTGCTGTTTGACCCAGGCTTCGATGGCCGCTtg                   >  2:22808‑M1/10‑120 (MQ=255)                                                                                                                      agggcgatCTTGGTCAGCGGCTGGCTGGCCAGCTCCTTGAGCTTGGTCAGGGCCAGCGGCAGTTGCTGGGCGGCAGGGTGGTCCGGGTGCTGTTTGACCCAGGCTTCGATGGCCGCTTGCTGCTGCTCAAGGGTGc  <  1:417555‑M1/128‑1 (MQ=255)                                                                                                                              |                                                                                                                                TGCTGGGCCTGGAAGTGGGCGGCCATGAAGCCATCGCGCAGGGCGCGGGCTACACCGGCCAGCGGGCCTTCCTGAGGCAGCAGCAGGGCGATCTTGGTCAGCGGCTGGCTGGCCAGCTCCTTGAGCTTGGTCAGCGGCTGGCTGGCCAGCTCCTTGAGCTTGGTCAGGGCCAGCGGCAGTTGCTGGGCGGCAGGGTGGTCCGGGTGCTGTTTGACCCAGGCTTCGATGGCCGCTTGCTGCTGCTCCAGGGTGC  >  NC_002947/1515677‑1515929

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