Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
MC JC	CP000730	2,144,161	4 bp→TGGTTTATCT	coding (480‑483/1254 nt)	sdrH ←	Ser‑Asp rich fibrinogen/bone sialoprotein bindin protein SdrH

Missing coverage evidence...
			seq id	start	end	size	←reads	reads→	gene	description
*	*	÷	CP000730	2144161	2144164	4	22 [0]	[0] 17	sdrH	Ser‑Asp rich fibrinogen/bone sialoprotein bindin protein SdrH

New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	CP000730	= 2144160	0 (0.000)	17 (0.540) +TGGTTTATCT	13/250	1.8	100%	coding (484/1254 nt)	sdrH	Ser‑Asp rich fibrinogen/bone sialoprotein bindin protein SdrH
	?	CP000730	2144165 =	0 (0.000)					coding (479/1254 nt)	sdrH	Ser‑Asp rich fibrinogen/bone sialoprotein bindin protein SdrH

TCGAGCCACCAGAATGATTGGAATCCCCAGGTTGATCGGGATCTGGTGACGGGTTAGGATTTGGCTTATTAGGGTCTGGTTTTGGATTCGGATTTGGCTTTGGTTTATCTGGATCTGGTTTTGGATCCGGATTTGGCTTATTAGGGTCTGGTTTTGGATCTGGATTTGGTTTCGGGTTATCTGGATCTGGTTTTGGATCCGGATTTGGCT  >  CP000730/2144022‑2144231                                                                                                                                           |                                                                        tCGAGCCACCAGAATGATTGGAATCCCCAGGTTGATCGGGATCTGGTGACGGGTTAGGATTTGGCTTATTAGGGTCTGGTTTTGGATTCGGATTTGGCTTTGGTTTATCTGGATCTGGTTTTGGATCCGGATTTGGCTTtg                                                                       >  1:5248‑M1/1‑139 (MQ=255)  cGAGCCACCAGAATGATTGGAATCCCCAGGTTGATCGGGATCTGGTGACGGGTTAGGATTTGGCTTATTAGGGTCTGGTTTTGGATTCGGATTTGGCTTTGGTTTATCTGGATCTGGTTTTGGATCCGGATTTGGCTTtgg                                                                      <  1:208104‑M1/141‑4 (MQ=255)        accaGAATGATTGGAATCCCCAGGTTGATCGGGATCTGGTGACGGGTTAGGATTTGGCTTATTAGGGTCTGGTTTTGGATTCGGATTTGGCTTTGGTTTATCTGGATCTGGTTTTGGATCCGGATTTGGCTTtggtttatc                                                                <  1:257039‑M1/141‑10 (MQ=255)             aaTGATTGGAATCCCCAGGTTGATCGGGATCTGGTGACGGGTTAGGATTTGGCTTATTAGGGTCTGGTTTTGGATTCGGATTTGGCTTTGGTTTATCTGGATCTGGTTTTGGATCCGGATTTGGCTTtggtttatctgggt                                                           >  1:219668‑M1/1‑127 (MQ=255)                       aTCCCCAGGTTGATCGGGATCTGGTGACGGGTTAGGATTTGGCTTATTAGGGTCTGGTTTTGGATTCGGATTTGGCTTTGGTTTATCTGGATCTGGTTTTGGATCCGGATTTGGCTTtggtttatctgggtctggttttgg                                                 >  1:289374‑M1/1‑117 (MQ=255)                         ccccAGGTTGATCGGGATCTGGTGACGGGTTAGGATTTGGCTTATTAGGGTCTGGTTTTGGATTCGGATTTGGCTTTGGTTTATCTGGATCTGGTTTTGGATCCGGATTTGGCTTtggtttatctgggtctggttttggat                                               <  1:191847‑M1/141‑27 (MQ=255)                           ccAGGTTGATCGGGATCTGGTGACGGGTTAGGATTTGGCTTATTAGGGTCTGGTTTTGGCTTCGGATTTGGTTTTGGTTTATCTGGATCTGGTTTTGGATCCGGATTTGGCTTtggtttatctgggtctgggtttggatcc                                             >  2:271765‑M1/1‑113 (MQ=255)                              ggTTGATCGGGATCTGGTGACGGGTTAGGATTTGGCTTATTAGGGTCTGGTTTTGGATTCGGATTTGGCTTTGGTTTATCTGGATCTGGTTTTGGATCCGGATTTGGCTTtggtttatctgggtctggttttggatccggt                                          >  2:161366‑M1/1‑110 (MQ=255)                                ttGATCGGGATCTGGTGACGGGTTAGGATTTGGCTTATTAGGGTCTGGTTTTGGATTCGGATTTGGCTTTGGTTTATCTGGATCTGGTTTTGGATCCGGATTTGGCTTtggtttatctgggtctggttttggatccggttt                                        <  1:33102‑M1/141‑34 (MQ=255)                                ttGATCGGGATCTGGTGACGGGTTAGGATTTGGCTTATTAGGGTCTGGTTTTGGATTCGGATTTGGCTTTGGTTTATCTGGATCTGGTTTTGGATCCGGATTTGGCTTtggtttatctgggtctggttttggatccggttt                                        <  2:202415‑M1/141‑34 (MQ=255)                                 tGATCGGGATCTGGTGACGGGTTAGGATTTGGCTTATTAGGGTCTGGTTTTGGATTCGGATTTGGCTTTGGTTTATCTGGATCTGGTTTTGGATCCGGATTTGGCTTtggtttatctgggtctggttttggatccggtttt                                       <  1:327699‑M1/141‑35 (MQ=255)                                     cGGGATCTGGTGACGGGTTAGGATTTGGCTTATTAGGGTCTGGTTTTGGATTCGGATTTGGCTTTGGTTTATCTGGATCTGGTTTTGGATCCGGATTTGGCTTt                                                                        >  1:24852‑M1/1‑103 (MQ=255)                                     cGGGATCTGGTGACGGGTTAGGATTTGGCTTATTAGGGTCTGGTTTTGGATTCGGATTTGGCTTTGGTTTATCTGGATCTGGTTTTGGATCCGGATTTGGCTTt                                                                        <  2:24852‑M1/104‑2 (MQ=255)                                           cTGGTGACGGGTTAGGATTTGGCTTATTAGGGTCTGGTTTTGGATTCGGATTTGGCTTTGGTTTATCTGGATCTGGTTTTGGATCCGGATTTGGCTTtggtttatctgggtctggttttggatccggttttggtttcgggt                             <  2:92742‑M1/141‑45 (MQ=255)                                             ggTGACGGGTTAGGATTTGGCTTATTAGGGTCTGGTTTTGGATTCGGATTTGGCTTTGGTTTATCTGGATCTGGTTTTGGATCCGGATTTGGCTTtggtttatctgggtctggttttggatccggttttggtttcgggttt                           >  2:228004‑M1/1‑95 (MQ=255)                                                 aCGGGTTAGGATTTGGCTTATTAGGGTCTGGTTTTGGATTCGGATTTGGCTTTGGTTTATCTGGATCTGGTTTTGGATCCGGATTTGGCTTtggtttatctgggtctggttttggatccggttttggtttcgggttatctg                       >  1:394833‑M1/1‑91 (MQ=255)                                                  cGGGTTAGGATTTGGCTTATTAGGGTCTGGTTTTGGATTCGGATTTGGCTTTGGTTTATCTGGATCTGGTTTTGGATCCGGATTTGGCTTtggtttatctgggtctggtttagatccggttttggtttcgggttatctgga                      <  2:161156‑M1/141‑52 (MQ=255)                                                   gggTTAGGATTTGGCTTATTAGGGTCTGGTTTTGGATTCGGATTTGGCTTTGGTTTATCTGGATCTGGTTTTGGATCCGGATTTGGCTTtggtttatctgggtctggttttggatccggttttggtttcgggttatctgga                     >  1:86978‑M1/1‑89 (MQ=255)                                                      ttAGGATTTGGCTTATTAGGGTCTGGTTTTGGATTCGGATTTGGCTTTGGTTTATCTGGATCTGGTTTTGGATCCGGATTTGGCTGtgggttatcggggtctggttttggatccggttttggtttcgggttatctggatct                  >  1:236664‑M1/1‑86 (MQ=255)                                                            tttGGCTTATTAGGGTCTGGTTTTGGATTCGGATTTGGCTTTGGTTTATCTGGATCTGGTTTTGGATCCGGATTTGGCTTtggtttatctgggtctggttttggatccggttttggtttcgggttatctggatctggttta            >  1:354086‑M1/1‑80 (MQ=255)                                                                  ttattaGGGTCTGGTTTTGGATTCGGATTTGGCTTTGGTTTATCTGGATCTGGTTTTGGATCCGGATTTGGCTTgggtttatctgggtctggttttggatccggttttggtttcgggttatctggatctggtttaggattc      >  1:337314‑M1/1‑74 (MQ=255)                                                                      taGGGTCTGGTTTTGGATTCGGATTTGGCTTTGGTTTATCTGGATCTGGTTTTGGATCCGGATTTGGCTTtggtttatctgggtctggttttggatccggttttggtttcgggttatctgtatctggtttaggattcggac  <  1:42489‑M1/141‑72 (MQ=255)                                                                                                                                           |                                                                        TCGAGCCACCAGAATGATTGGAATCCCCAGGTTGATCGGGATCTGGTGACGGGTTAGGATTTGGCTTATTAGGGTCTGGTTTTGGATTCGGATTTGGCTTTGGTTTATCTGGATCTGGTTTTGGATCCGGATTTGGCTTATTAGGGTCTGGTTTTGGATCTGGATTTGGTTTCGGGTTATCTGGATCTGGTTTTGGATCCGGATTTGGCT  >  CP000730/2144022‑2144231

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