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

AACACCTGTACCGCAGGCGCTGGCTGGCCGCTTTGCTGGGCCTGGAAGTGGGCGGCCATGAAGCCATCGCGCAGGGCGCGGGCTACACCGGCCAGCGGGCCTTCCTGAGGCAGCAGCAGGGCGATCTTGGTCAGCGGCTGGCTGGCCAGCTCCTTGAGCTTGGTCAGCGGCTGGCTGGCCAGCTCCTTGAGCTTGGTCAGGGCCAGCGGCAGTTGCTGGGCGGCAGGGTGGTCCGGGTGCTGTTTGACCCAGG  >  NC_002947/1515644‑1515896                                                                                                                             |                                                                                                                                 aaCACCTGTACCGCAGGCGCTGGCTGGCCGCTTTGCTGGGCCTGGAAGTGGGCGGCCATGAAGCCATCGCGCAGGGCGCGGGCTACACCGGCCAGCGGGCCTTCCTGAGGCAGCAGCAGGGCGATcttggtcagcg                                                                                                                       <  1:113201‑M2/136‑12 (MQ=255)   cacCTGTACCGCAGGCGCTGGCTGGCCGCTTTGCTGGGCCTGGAAGTGGGCGGCCATGAAGCCATCGCGCAGGGCGCGGGCTACACCGGCCAGCGGGCCTTCCTGAGGCAGCAGCAGGGCGATcttggtcagcggc                                                                                                                     >  1:376169‑M2/1‑123 (MQ=255)                 gcgcTGGCTGGCCGCTTTGCTGGGCCTGGAAGTGGGCGGCCATGAAGCCATCGCGCAGGGCGCGGGCTACACCGGCCAGCGGGCCTTCCTGAGGCAGCAGCAGGGCGATcttggtcagcggctggctggccagctc                                                                                                       <  2:76676‑M2/136‑28 (MQ=255)                 gcgcTGGCTGGCCGCTTTGCTGGGCCTGGAAGTGGGCGGCCATGAAGCCATCGCGCAGGGCGCGGGCTACACCGGCCAGCGGGCCTTCCTGAGGCAGCAGCAGGGCGATcttggtcagcggctggctggccagctc                                                                                                       <  1:192692‑M2/136‑28 (MQ=255)                 gcgcTGGCTGGCCGCTTTGCTGGGCCTGGAAGTGGGCGGCCATGAAGCCATCGCGCAGGGCGAGGGCTACACCGGCCAGCGGGCCTTCCTGAGGCAGCAGCAGGGCGATcttggtcagcggctggctggccagctc                                                                                                       <  1:142309‑M2/136‑28 (MQ=255)                      ggctgggCGCTTTGCTGGGCCTGGAAGTGGGCGGCCATGAAGCCATCGCGCAGGGCGCGGGCTACACCGGCCAGCGGGCCTTCCTGAGGCTGCAGCAGGGCGATcttggtcagcggctggctggccagctccttga                                                                                                  >  1:478323‑M2/1‑104 (MQ=255)                      ggctggcCGCTTTGCTGGGCCTGGAAGTGGGCGGCCATGAAGCCATCGCGCAGGGCGCGGGCTACACCGGCCAGCGGGCCTTCCTGAGGCAGCAGCAGGGCGATcttggtcagcggctggctggccagctccttga                                                                                                  >  2:387058‑M2/1‑104 (MQ=255)                                  tGCTGGGCCTGGAAGTGGGCGGCCATGAAGCCATCGCGCAGGGCGCGGGCTACACCGGCCAGCGGGCCTTCCTGAGGCAGCAGCAGGGCGATctgggtcagcggctggctggccagctccttgagcttggtcaggg                                                                                      >  1:8751‑M2/1‑92 (MQ=255)                                            ggAAGTGGGCGGCCATGAAGCCATCGCGCAGGGCGCGGGCTACACCGGCCAGCGGGCCTTCCTGAGGCAGCAGCAGGGCGATcttggtcagcggctggctggccagctccttgagcttggtcagggccagcggcag                                                                            <  1:192873‑M2/136‑55 (MQ=255)                                            ggAAGTGGGCGGCCATGAAGCCATCGCGCAGGGCGCGGGCTACACCGGCCAGCGGGCCTTCCTGAGGCAGCAGCAGGGCGATcttggtcagcggctggctggccagctccttgagcttggtcagggccagcggcag                                                                            <  1:69740‑M2/136‑55 (MQ=255)                                                      ggcCATGAAGCCATCGCGCAGGGCGCGGGCTACACCGGCCAGCGGGCCTTCCTGAGGCAGCAGCAGGGCGATcttggtcagcggctggctggccagctccttgagcttggtcagggccagcggcagttgctgggcg                                                                  <  1:82747‑M2/136‑65 (MQ=255)                                                      ggcCATGAAGCCATCGCGCAGGGCGCGGGCTACACCGGCCAGCGGGCCTTCCTGAGGCAGCAGCAGGGCGATcttggtcagcggctggctggccagctccttgagcttggtcagggccagcggcagttgctgggcg                                                                  <  1:1820‑M2/136‑65 (MQ=255)                                                          aTGAAGCCATCGCGCAGGGCGCGGGCTACACCGGCCAGCGGGCCTTCCTGAGGCAGCAGCAGGGCGATcttggtcagcggctggctggccagctccttgagcttggtcagggccagcggcagttgctgggcggcag                                                              >  2:232202‑M2/1‑68 (MQ=255)                                                                  agcgcgcaGGGCGCGGGCTACACCGGCCAGCGGGCCTTCCTGCGGCAGCAGCAGGGCGATcttggtcagcggctggctggccagctccttgagcttggtcagggccagcggcagttgctgggcggcagggtggtcc                                                      >  1:334099‑M2/3‑60 (MQ=255)                                                                                              aGCGGGCCTTCCTGAGGCAGCAGCAGGGCGATcttggtcagcggctggctggccagctccttgagcttggtcagggccagcggcagttgctgggcggcagggtggtccgggtgctgtttgacccaggcttcgatgg                          <  2:383892‑M2/136‑105 (MQ=255)                                                                                              aGCGGGCCTTCCTGAGGCAGCAGCAGGGCGATcttggtcagcggctggctggccagctccttgagcttggtcagggccagcggcagttgctgggcggcagggtggtccgggtgctgtttgacccaggcttcgatgg                          <  2:392126‑M2/136‑105 (MQ=255)                                                                                               gCGGGCCTTCCTGAGGCAGCAGCAGGGCGATcttggtcagcggctggctggccagctccttgagcttggtcagggccagcggcagttgcggggcggcagggtggtccgggtgctgtttgacccaggcttcgatggc                         <  1:254630‑M2/136‑106 (MQ=255)                                                                                                                     cagGGCGATcttggtcagcggctggctggccagctccttgagcttggtcagggccagcggcagttgctgggcggcagggtggtccgggtgctgtttgacccaggcttcgatggccgcttg                   <  1:22808‑M2/120‑112 (MQ=255)                                                                                                                     cagGGCGATcttggtcagcggctggctggccagctccttgagcttggtcagggccagcggcagttgctgggcggcagggtggtccgggtgctgtttgacccaggcttcgatggccgcttg                   >  2:22808‑M2/1‑9 (MQ=255)                                                                                                                      agGGCGATcttggtcagcggctggctggccagctccttgagcttggtcagggccagcggcagttgctgggcggcagggtggtccgggtgctgtttgacccaggcttcgatggccgcttgctgctgctcaagggtgc  <  1:417555‑M2/136‑129 (MQ=255)                                                                                                                             |                                                                                                                                 AACACCTGTACCGCAGGCGCTGGCTGGCCGCTTTGCTGGGCCTGGAAGTGGGCGGCCATGAAGCCATCGCGCAGGGCGCGGGCTACACCGGCCAGCGGGCCTTCCTGAGGCAGCAGCAGGGCGATCTTGGTCAGCGGCTGGCTGGCCAGCTCCTTGAGCTTGGTCAGCGGCTGGCTGGCCAGCTCCTTGAGCTTGGTCAGGGCCAGCGGCAGTTGCTGGGCGGCAGGGTGGTCCGGGTGCTGTTTGACCCAGG  >  NC_002947/1515644‑1515896

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