Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
MC JC	NC_000913	1,394,106	Δ1,938 bp	IS5‑mediated	[mppA]–ynaI	[mppA], ynaI

Missing coverage evidence...
			seq id	start	end	size	←reads	reads→	gene	description
*	*	÷	NC_000913	1394106	1397224–1396040	1935–3119	68 [3]	[24] 27	[mppA]–insH1	[mppA], ynaI, insH1

New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913	= 1394105	3 (0.050)	65 (1.010)	49/198	0.3	95.6%	coding (879/1614 nt)	mppA	murein tripeptide (L‑ala‑gamma‑D‑glutamyl‑meso‑DAP) transporter subunit
	?	NC_000913	1396044 =	NA (NA)					noncoding (1/1195 nt)	IS5	repeat region

ACGGGAGAAAAAATCGGCTCAAACATGAAGAAATGAAATGACTGAGTCAGCCGAGAAGAATTTCCCCGCTTATTCGCACCTTCCCTAACTAATCAATGCGTTGATTGTAAATCCAGCTAAGAGGTGAGGTTTTCAGAGCAGACAACGGTGAAATGTCATGGTATTGTTAC  >  NC_000913/1397155‑1397324                                                                       |                                                                                                    aagagagaAAAAATCGGCTCAAACATGAAGAAATGAAATGACTGAGTCAGCCGAGAAGAATTTCCCCGCTTATTCGCACCTTCCCTAACTAATCAATg                                                                          <  2:437312/94‑1 (MQ=255) aCGGGAGAAAAAATCGGCTCAAACATGAAGAAATGAAATGACTGAGTCAGCCGAGAAGAATTTCCCCGCTTATTCGCACCTTCCCTAACTAATCAATGCg                                                                        >  1:446820/1‑100 (MQ=255)                ggCTCAAACATGAAGACATCAATTTACTGAGTGAGCCGAGAAGAATGTCCCCGCTTATTCGCGCCTTCCCTATCTCATCAATGCGTGGATTGTAAATCCa                                                         <  2:563534/100‑1 (MQ=17)                 gCTCAAACATGAAGAAATGAAATGACTGAGTCAGCCGAGAAGAATTTCCCCGCTTATTCGCACCTTCCCTAACTAATCAATGCGTTGATTGTAAATCCAg                                                        <  1:710331/100‑1 (MQ=255)                  cTCAAACATGAAGAAATGAAATGACTGAGTCAGCCGAGAAGAATTTCCCCGCTTATTCGCACCTTCCCTAACTAATCAATGCGTTGATTGTAAATCCAGc                                                       >  1:381725/1‑100 (MQ=255)                   tCAAACATGAAGAAATGAAATGACTGAGTCAGCCGAGAAGAATTTCCCCGCTTATTCGCACCTTCCCTAACTAATCAATGCGTTGATTGTAAATCCAGCt                                                      <  2:723026/100‑1 (MQ=255)                        cATGAAGAAATGAAATGACTGAGTCAGCCGAGAAGAATTTCCCCGCTTATTCGCACCTTCCCTAACTAATCAATGCGTTGATTGTAAATCCAGCTAagag                                                 >  2:1562029/1‑100 (MQ=255)                         aTGAAGAAATGAAATGACTGAGTCAGCCGAGAAGAATTTCCCCGCTTATTCGCACCTTCCCTAACTAATCAATGCGTTGATTGTAAATCCAGCTAAgagg                                                <  2:464427/100‑1 (MQ=255)                          tGAAGAAATGAAATGACTGAGTCAGCCGAGAAGAATTTCCCCGCTTATTCGCACCTTCCCTAACTAATCAATGCGTTGATTGTAAATCCAGCTAAgaggt                                               >  2:525935/1‑100 (MQ=255)                           gaagaaATGAAATGACTGAGTCAGCCGAGAAGAATTTCCCCGCTTATTCGCACCTTCCCTAACTAATCAATGCGTTGATTGTAAATCCAGCTAAgaggtg                                              <  1:691845/100‑1 (MQ=255)                             agaaATGAAATGACTGAGTCAGCCGAGAAGAATTTCCCCGCTTATTCGCACCTTCCCTAACTAATCAATGCGTTGATTGTAAATCCAGCTAAgaggtgag                                            >  2:275677/1‑100 (MQ=255)                             agaaATGAAATGACTGAGTCAGCCGAGAAGAATTTCCCCGCTTATTCGCACCTTCCCTAACTAATCAATGCGTTGATTGTAAATCCAGCTAAgaggtgag                                            >  2:491204/1‑100 (MQ=255)                              gaaatgaaatGACTGAGTCAGCCGAGAAGAATTTCCCCGCTTATTCGCACCTTCCCTAACTAATCAATGCGTTGATTGTAAATCCAGCTAAgaggtgagg                                           <  1:312682/100‑1 (MQ=255)                                   gaaatgaCTGAGTCAGCCGAGAAGAATTTCCCCGCTTATTCGCACCTTCCCTAACTAATCAATGCGTTGATTGTAAATCCAGCTAAGAGGTGAGGTTTTc                                      >  2:164541/1‑100 (MQ=255)                                     aatgaCTGAGTCAGCCGAGAAGAATTTCCCCGCTTATTCGCACCTTCCCTAACTAATCAATGCGTTGATTGTAAATCCAGCTAAGAGGTGAGGTTTTCag                                    <  1:264910/100‑1 (MQ=255)                                        gaCTGAGTCAGCCGAGAAGAATTTCCCCGCTTATTCGCACCTTCCCTAACTAATCAATGCGTTGATTGTAAATCCAGCTAAGAGGTGAGGTTTTCAGAGc                                 <  1:809484/100‑1 (MQ=255)                                        gaCTGAGTCAGCCGAGAAGAATTTCCCCGCTTATTCGCACCTTCCCTAACTAATCAATGCGTTGATTGTAAATCCAGCTAAGAGGTGAGGTTTTCAGAGc                                 <  1:1487908/100‑1 (MQ=255)                                         aCTGAGTCAGCCGAGAAGAATTTCCCCGCTTATTCGCACCTTCCCTAACTAATCAATGCGTTGATTGTAAATCCAGCTAAGAGGTGAGGTTTTCAGAGCa                                <  1:952640/100‑1 (MQ=255)                                                cAGCCGAGAAGAATTTCCCCGCTTATTCGCACCTTCCCTAACTAATCAATGCGTTGATTGTAAATCCAGCTAAGAGGTGAGGTTTTCAGAGCAGACAACg                         >  1:112233/1‑100 (MQ=255)                                                 aGCCGAGAAGAATTTCCCCGCTTATTCGCACCTTCCCTAACTAATCAATGCGTTGATTGTAAATCCAGCTAAGAGGTGAGGTTTTCAGAGCAGACAACgg                        <  2:245087/100‑1 (MQ=255)                                                    cGAGAAGAATTTCCCCGCTTATTCGCACCTTCCCTAACTAATCAATGCGTTGATTGTAAATCCAGCTAAGAGGTGAGGTTTTCAGAGCAGACAACGGTGa                     <  2:808470/100‑1 (MQ=255)                                                             tttCCCCGCTTATTCGCACCTTCCCTAACTAATCAATGCGTTGATTGTAAATCCAGCTAAGAGGTGAGGTTTTCAGAGCAGACAACGGTGAAATGTCATg            >  1:89265/1‑100 (MQ=255)                                                                 cccGCTTATTCGCACCTTCCCTAACTAATCAATGCGTTGATTGTAAATCCAGCTAAGAGGTGAGGTTTTCAGAGCAGACAACGGTGAAATGTCATGGTAt        >  2:175134/1‑100 (MQ=255)                                                                 cccGCTTATTCGCACCTTCCCTAACTAATCAATGCGTTGATTGTAAATCCAGCTAAGAGGTGAGGTTTTCAGAGCAGACAACGGTGAAATGTCATGGTAt        >  1:241099/1‑100 (MQ=255)                                                                       tATTCGCACCTTCCCTAACTAATCAATGCGTTGATTGTAAATCCAGCTAAGATGTGAGGTTTTCAGAGCAGACAACAGTGAAATGTCATGGTATTGTTAc  <  2:1253469/100‑1 (MQ=255)                                                                       tATTCGCACCTTCCCTAACTAATCAATGCGTTGATTGTAAATCCAGCTAAGAGGTGAGGTTTTCAGAGCAGACAACGGTGAAATGTCATGGTATTGTTAc  >  1:39645/1‑100 (MQ=255)                                                                       tATTCGCACCTTCCCTAACTAATCAATGCGTTGATTGTAAATCCAGCTAAGAGGTGAGGTTTTCAGAGCAGACAACGGTGAAATGTCATGGTATTGTTAc  >  1:1520772/1‑100 (MQ=255)                                                                       |                                                                                                    ACGGGAGAAAAAATCGGCTCAAACATGAAGAAATGAAATGACTGAGTCAGCCGAGAAGAATTTCCCCGCTTATTCGCACCTTCCCTAACTAATCAATGCGTTGATTGTAAATCCAGCTAAGAGGTGAGGTTTTCAGAGCAGACAACGGTGAAATGTCATGGTATTGTTAC  >  NC_000913/1397155‑1397324

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