Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
MC JC	NC_000913	2,948,255	Δ13 bp	coding (744‑756/1254 nt)	amiC ←	N‑acetylmuramoyl‑L‑alanine amidase C

Missing coverage evidence...
			seq id	start	end	size	←reads	reads→	gene	description
*	*	÷	NC_000913	2948255	2948267	13	33 [0]	[0] 33	amiC	N‑acetylmuramoyl‑L‑alanine amidase C

New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913	= 2948254	0 (0.000)	32 (1.180)	27/272	0.1	100%	coding (757/1254 nt)	amiC	N‑acetylmuramoyl‑L‑alanine amidase C
	?	NC_000913	2948268 =	0 (0.000)					coding (743/1254 nt)	amiC	N‑acetylmuramoyl‑L‑alanine amidase C

TTTGCCGCAGTACTGGTTGCACCTTTGGTTGAGAGCGCAAACACAGAGGAACCGCTCGGCTGACGACTGGTAAAGGCGTCGGCATGGATAGAGACAAACAGGTCAGCACGCTGTTTCTGGGCTTTTGCTACGCGCACTTGCAATGGAATGAAGATGTCTTCATTGCGCGTCATGTACACCTTCATATTGCCCTCTTTCTCGATCAGAGAGCGCAAACGGCGAGCTATTTGCAATACCACATCTTTTTCGCGCGTTTTGTATTTCCCCACCGCAC  >  NC_000913/2948133‑2948406                                                                                                                                        |                                                                                                                                           ttgtgccagatattttgccgcagtactggttgcacctttggttgagagcgcaaacacagaggaaccgctcggctgacgactggtaaaggcgtcggcatggatagagacaaacaggtcagcacgctgtttctgggcACTTg                                                                                                                                        >  1:199779‑M2/136‑140 (MQ=255) ttgtgccagatattttgccgcagtactggttgcacctttggttgagagcgcaaacacagaggaaccgctcggctgacgactggtaaaggcgtcggcatggatagagacaaacaggtcagcacgctgtttctgggcACTTg                                                                                                                                        >  1:31426‑M2/136‑140 (MQ=255)   gtgccagatattttgccgcagtactggttgcacctttggttgagagcgcaaacacagaggaaccgctcggctgacgactggtaaaggcgtcggcatggatagagacaaacaggtcagcacgctgtttctgggcACTTGCa                                                                                                                                      <  2:66972‑M2/7‑1 (MQ=255)    tgccagatattttgccgcagtactggttgcacctttggttgagagcgcaaacacagaggaaccgctcggctgacgactggtaaaggcgtcggcatggatagagacaaacaggtcagcacgctgtttctgggcACTTGCaa                                                                                                                                     <  1:132366‑M2/8‑1 (MQ=255)    tgccagatattttgccgcagtactggttgcacctttggttgagagcgcaaacacagaggaaccgctcggctgacgactggtaaaggcgtcggcatggatagagacaaacaggtcagcacgctgtttctgggcACTTGCaa                                                                                                                                     <  2:404240‑M2/8‑1 (MQ=255)    tgccagatattttgccgcagtactggttgcacctttggttgagagcgcaaacacagaggaaccgctcggctgacgactggtaaaggcgtcggcatggatagagacaaacaggtcagcacgctgtttctgggcACGTGCaa                                                                                                                                     >  2:450967‑M2/133‑140 (MQ=255)              tttgccgcagtactggttgcacctttggttgagagcgcaaacacagaggaaccgctcggctgacgactggtaaaggcgtcggcatggatagagacaaacaggtcagcacgctgtttctgggcACTTGCAATGGAATGAAg                                                                                                                           >  2:64088‑M2/123‑140 (MQ=255)                           tggttgcacctttggttgagagcgcaaacacagaggaaccgctcggctgacgactggtaaaggcgtcggcatggatagagacaaacaggtcagcacgctgtttctgggcACTTGCAATGGAATGAAGATGTCTTCATTgc                                                                                                              <  2:199779‑M2/31‑1 (MQ=255)                              ttgcacctttggttgagagcgcaaacacagaggaaccgctcggctgacgactggtaaaggcgtcggcatggatagagacaaacaggtcagcacgctgtttctgggcACTTGCAATGGAATGAAGATGTCTTCATTGcgcg                                                                                                           <  1:64088‑M2/34‑1 (MQ=255)                                  acctttggttgagagcgcaaacacagaggaaccgctcggctgacgactggtaaaggcgtcggcatggatagagacaaacaggtcagcacgctgtttctgggcACTTGCAATGGAATGAAGATGTCTTCATTGCGCGTCAt                                                                                                       >  1:231435‑M2/103‑140 (MQ=255)                                  acctttggttgagagcgcaaacacagaggaaccgctcggctgacgactggtaaaggcgtcggcatggatagagacaaacaggtcagcacgctgtttctgggcACTTGCAATGGAATGAAGATGTCTTCATTGCGCGTCAt                                                                                                       >  1:453305‑M2/103‑140 (MQ=255)                                  acctttggttgagagcgcaaacacagaggaaccgctcggctgacgactggtaaaggcgtcggcatggatagagacaaacaggtcagcacgctgtttctgggcACTTGCAATGGAATGAAGATGTCTTCATTGCGCGTCAt                                                                                                       >  2:459268‑M2/103‑140 (MQ=255)                                       tggttgagagcgcaaacacagaggaaccgctcggctgacgactggtaaaggcgtcggcatggatagagacaaacaggtcagcacgctgtttctgggcACTTGCAATGGAATGAAGATGTCTTCATTGCGCGTCATGTaca                                                                                                  <  1:288232‑M2/43‑1 (MQ=255)                                            gagagcgcaaacacagaggaaccgctcggctgacgactggtaaaggcgtcggcatggatagagacaaacaggtcagcacgctgtttctgggcACTTGCAATGGAATGAAGATGTCTTCATTGCGCGTCATGTACACCTTc                                                                                             <  2:453305‑M2/48‑1 (MQ=255)                                              gagcgcaaacacagaggaaccgctcggctgacgactggtaaaggcgtcggcatggatagagacaaacaggtcagcacgctgtttctgggcACTTGCAATGGAATGAAGATGTCTTCATTGcgcg                                                                                                           >  1:390837‑M2/91‑124 (MQ=255)                                              gagcgcaaacacagaggaaccgctcggctgacgactggtaaaggcgtcggcatggatagagacaaacaggtcagcacgctgtttctgggcACTTGCAATGGAATGAAGATGTCTTCATTGcgcg                                                                                                           <  2:390837‑M2/34‑1 (MQ=255)                                                    aaacacagaggaaccgctcggctgacgactggtaaaggcgtcggcatggatagagacaaacaggtcagcacgctgtttctgggcACTTGCAATGGAATGAAGATGTCTTCATTGCGCGTCATGTACACCTTCATATTGcc                                                                                     <  2:168307‑M2/56‑1 (MQ=255)                                                          agaggaaccgctcggctgacgactggtaaaggcgtcggcatggatagagacaaacaggtcagcacgctgtttctgggcACTTGCAATGGAATGAAGATGTCTTCATTGCGCGTCATGTACACCTTCATATTGCCCTCttt                                                                               <  2:112825‑M2/62‑1 (MQ=255)                                                                      cggctgacgactggtaaaggcgtcggcatggatagagacaaacaggtcagcacgctgtttctgggcACTTGCAATGGAATGAAGATGTCTTCATTGCGCGTCATGTACACCTTCATATTGCCCTCTTTCTCGATCagaga                                                                   >  1:441537‑M2/67‑140 (MQ=255)                                                                       ggctgacgactggtaaaggcgtcggcatggatagagacaaacaggtcagcacgctgtttctgggcACTTGCAATGGAATGAAGATGTCTTCATTGCGCGTCATGTACACCTTCATATTGCCCTCTTTCTCGATCagagag                                                                  <  2:31426‑M2/75‑1 (MQ=255)                                                                                 tggtaaaggcgtcggcatggatagagacaaacaggtcagcacgctgtttctgggcACTTGCAATGGAATGAAGATGTCTTCATTGCGCGTCATGTACACCTTCATATTGCCCTCTTTCTCGATCAGAGAGCGCAAACGGc                                                        <  2:441537‑M2/85‑1 (MQ=255)                                                                                   gtaaaggcgtcggcatggatagagacaaacaggtcagcacgctgtttctgggcCCTTGCAATGGAATGAAGATGTCTTCATTGCGCGTCATGTACACCTTCATATTGCCCTCTTTCTCGATCAGAGAGCGCAAACGGCGa                                                      >  1:322702‑M2/54‑140 (MQ=255)                                                                                    taaaggcgtcggcatggatagagacaaacaggtcagcacgctgtttctgggcACTTGCAATGGAATGAAGATGTCTTCATTGCGCGTCATGTACACCTTCATATTGCCCTCTTTCTCGATCAGAGAGCGCAAACGGCGAg                                                     <  1:459268‑M2/88‑1 (MQ=255)                                                                                            tcggcatggatagagacaaacaggtcagcacgctgtttctgggcACTTGCAATGGAATGAAGATGTCTTCATTGCGCGTCATGTACACCTTCATATTGCCCTCTTTCTCGATCAGAGAGCGCAAACGGCGAGCTATTTGc                                             >  1:264609‑M2/45‑140 (MQ=255)                                                                                             tggcatggatagagtgaaacaggttagcacgttgtttctgggcACTTGCAATGGAATGAAGATGTCTTCATTGCGCGTCATGTACACCTTCATATTGCCCTCTTTCTCTATCAGAGAGCGCAAACGGTGAGCTATTTGCa                                            <  2:322702‑M2/97‑1 (MQ=255)                                                                                             cggcatggatagagacaaacaggtcagcacgctgtttctgggcACTTGCAATGGAATGAAGATGTCTTCATTGCGCGTCATGTACACCTTCATATTGCCCTCTTTCTCGATCAGAGAGCGCAAACGGCGAGCt                                                   <  1:256772‑M2/90‑1 (MQ=255)                                                                                             cggcatggatagagacaaacaggtcagcacgctgtttctgggcACTTGCAATGGAATGAAGATGTCTTCATTGCGCGTCATGTACACCTTCATATTGCCCTCTTTCTCGATCAGAGAGCGCAAACGGCGAGCt                                                   >  2:256772‑M2/44‑133 (MQ=255)                                                                                                   ggatagagacaaacaggtcagcacgctgtttctgggcACTTGCAATGGAATGAAGATGTCTTCATTGCGCGTCATGTACACCTTCATATTGCCCTCTTTCTCGATCAGAGAGCGCAAACGGCGAGCTATTTGCAATACca                                      >  2:243887‑M2/38‑140 (MQ=255)                                                                                                              aacaggtcagcacgctgtttctgggcACTTGCAATGGAATGAAGATGTCTTCATTGCGCGTCATGTACACCTTCATATTGCCCTCTTTCTCGATCAGAGAGCGCAAACGGCGAGCTATTTGCAATACCACATCTTTTTcg                           >  2:165332‑M2/27‑140 (MQ=255)                                                                                                                  ggtcagcacgctgtttctgggcACTTGCAATGGAATGAAGATGTCTTCATTGCGCGTCATGTACACCTTCATATTGCCCTCTTTCTCGATCAGAGAGCGCAAACGGCGAGCTATTTGCAATACCACATCTTTTTcgcgcg                       <  1:243887‑M2/118‑1 (MQ=255)                                                                                                                   gtcagcacgctgtttctgggcACTTGCAATGGAATGAAGATGTCTTCATTGCGCGTCATGTACACCTTCATATTGCCCTCTTTCTCGATCAGAGAGCGCAAACGGCGAGCTATTTGCAATACCACATCTTTTTCGCGCGt                      >  2:426815‑M2/22‑140 (MQ=255)                                                                                                                             tgtttctgggcACTTGCAATGGAATGAAGATGTCTTCATTGCGCGTCATGTACACCTTCATATTGCCCTCTTTCTCGATCAGAGAGCGCAAACGGCGAGCTATTTGCAATACCACATCTTTTTCGCGCGTTTTGTATTTc            >  1:134916‑M2/12‑140 (MQ=255)                                                                                                                                       cACTTGCAATGGAATGAAGATGTCTTCATTGCGCGTCATGTACACCTTCATATTGCCCTCTTTCTCGATCAGAGAGCGCAAACGGCGAGCTATTTGCAATACCACATCTTTTTCGCGCGTTTTGTATTTCCCCACCGCAc  >  1:131489‑M2/2‑140 (MQ=255)                                                                                                                                        |                                                                                                                                           TTTGCCGCAGTACTGGTTGCACCTTTGGTTGAGAGCGCAAACACAGAGGAACCGCTCGGCTGACGACTGGTAAAGGCGTCGGCATGGATAGAGACAAACAGGTCAGCACGCTGTTTCTGGGCTTTTGCTACGCGCACTTGCAATGGAATGAAGATGTCTTCATTGCGCGTCATGTACACCTTCATATTGCCCTCTTTCTCGATCAGAGAGCGCAAACGGCGAGCTATTTGCAATACCACATCTTTTTCGCGCGTTTTGTATTTCCCCACCGCAC  >  NC_000913/2948133‑2948406

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 22 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑