Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
MC JC	NC_000913	574,454	Δ137 bp	IS5‑mediated	nmpC	nmpC

Missing coverage evidence...
			seq id	start	end	size	←reads	reads→	gene	description
*	*	÷	NC_000913	574454	575597–574586	133–1144	116 [0]	[25] 26	nmpC–[insH2]	nmpC,[insH2]

New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913	= 574453	0 (0.000)	117 (1.100)	77/592	0.1	100%	intergenic (+114/+76)	quuD/nmpC	prophage antitermination protein Q homolog QuuD/putative outer membrane porin NmpC
	?	NC_000913	2101749 =	NA (NA)					noncoding (1195/1195 nt)	IS5	repeat region

GGTTGTACCAATGCTGCATGCAGTGAATGCGTAGCATGGTTTCCAGCGGATAAGGTCGCCGGCCATTACCAGCCTTGGGGTAAAACGGCTCGATGACTTCCACCATGTTTTGCCATGGCAGAATCTGCTCCATGCGGGACAAGAAAATCTCTTTTCTGGTCTGACGGCGCTTACTGCTGAATTCACTGTCGGCGAAGGTAAGTTGATGACTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCTTAGTGAAGTCATTTTTGTCAAGCAGGTTGATTTTGTAATCAACGAAAGTAGACATATTTTTGTTGAAGTAATAGGTTGCACCTACATCAACATATTTGACTAAGTCCTGAT  >  NC_000913/575504‑575897                                                                                               |                                                                                                                                                                                                                                                                                                            ggTTGTACCAATGCTGCATGCAGTGAATGCGTAGCATGGTTTCCAGCGGATAAGGTNGCCGGCCATTACCAGCCTTGGGGTAAAACGGCTCGATGACTTCCACCATGTTTTGCCATGGCAGAATCTGCTCCATGCGGGACAAGAAAATCTCTTTTCTGGTCTGACGGCGCTTACTGCTGAATTCACTGTCGGCGAAGGTAAGTTGATGACTCATGATGAACCCTGTTCTATGGCTCCAGATGACANACATGATCTAATATCAGGGACTTGTTCGCACCTTCCTTAGTGAAGTCATTTTTg                                                                                                >  2:793209/1‑300 (MQ=21) ggTTGTACCAATGCTGCATGCAGTGAATGCGTAGCATGGTTTCCAGCGGATAAGGTCGCCGGCCATTACCAGCCTTGGGGTAAAACGGCTCGATGACTTCCACCATGTTTTGCCATGGCAGAATCTGCTCCATGCGGGACAAGAAAATCTCTTTTCTGGTCTGACGGGGCTTACTGCTGAATTCACTGTCGGCGAAGGTAAGTTGATGACTCATGATGAACCCTGTTCTATGGCTCCAGATGACANACATGATCTCATATCAGGGACTTGTTCGCACCTTCCTTAGTGAAGTCATTTTTg                                                                                                >  2:625538/1‑300 (MQ=21) ggTTGTACCAATGCTGCATGCAGTGAATGCGTAGCATGGTTTCCAGCGGATAAGGTCGCCGGCCATTACCAGCCTTGGGGTAAAACGGCTCGATGACTTCCACCATGTTTTGCCATGGCAGAATCTGCTCCATGCGGGACAAGAAAATCTCTTTTCTGGTCTGACGGCGCTTACTGCTGAATTCACTGTCGGCGAAGGTAAGTTGATGACTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCTTAGTGAAGTCATTTTTg                                                                                                >  1:263824/1‑300 (MQ=34)           aTGCTGCATGCAGTGAATGCGTAGCATGGTTTCCAGCGGATAAGGTCGCCGGCCATTACCAGCCTTGGGGTAAAACGGCTCGATGACTTCCACCATGTTTTGCCATGGCAGAATCTGCTCCATGNGGGACAAGAAAATCTCTTTTCTGGTCTGACGGCGCTTACTGCTGAATTCACTGTCGGCGAANGTAAGTTGATGACTCATGATGAACCCTGTTCTCTGGCTCCAGATGACAAACATGATCTCATATCCGGGACTTGTTCGCACCTTCCTTAGTGAAGTAATTTTTGTCAAGCAGGt                                                                                      >  2:671846/1‑300 (MQ=255)                     caGTGAATGCGTAGCATGGTTTCCAGCGGATAAGGTCGCCGGCCATTACCAGCCTTGGGGTAAAACGGCTCGATGACTTCCACCATGTTTTGCCATGGCAGAATCTGCTCCATGCGGGACACGAAAATCTCTTTTCTGGTCTGACGGCGCTTACTGCTGAATTCACTGTCGGCGAAGGTAAGTTGATGACTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCTTAGTGAAGTCATTTTTGTCAAGCAGGTTGATTTTGTa                                                                            >  1:196603/1‑300 (MQ=255)                     caGTGAATGCGTAGCATGGTTTCCAGCGGATAAGGTCGCCGGCCATTACCAGCCTTGGGGTAAAACGGCTCGATGACTTCCACCATGTTTTGCCATGGCAGAATCTGCTCCATGCGGGACAAGAAAATCTCTTTTCTGGTCTGACGGCGCTTACTGCTGAATTCACTGTCGGCGAAGGTAAGTTGATGACTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCTTAGTGAAGTCATTTTTGTCAAGCAGGTTGATTTTGTa                                                                            >  1:183802/1‑300 (MQ=255)                       gTGAATGCGTAGCATGGTTTCCAGCGGATAAGGTCGCCGGCCATTACCAGCCTTGGGGTAAAACGGCTCGATGACTTCCACCATGTTTTGCCATGGCAGAATCTGCTCCATGCGGGACAAGAAAATCTCTTTTCTGGTCTGACGGCGCTTACTGCTGAATTCACTGTCGGCGAAGGTAAGTTGATGACTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCTTAGTGAAGTCATTTTTGTCAAGCAGGTTGATTTTGTAAt                                                                          >  1:437626/1‑300 (MQ=255)                                      ggTTTCCAGCGGATAAGGTCGCCGGCCATTACCAGCCTTGGGGTAAAACGGCTCGATGACTTCCACCATGTTTTGCCATGGCAGAATCTGCTCCATGCGGGACAAGAAAATCTCTTTTCTGGTCTGACGGCGCTTACTGCTGAATTCACTGTCGGCGAAGGTAAGTTGATGACTCATGATGAACCCTGTTCTATNGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCTTAGTGAAGTCATTTTTGTCAAGCAGGTTGATTTTGTAATCAACGAAAGTGGACa                                                           >  2:755522/1‑300 (MQ=255)                                       gTTTCCAGCGGATAAGGTCGCCGGCCATTACCAGCCTTGGGGTAAAACGGCTCGATGACTTCCACCATGTTTTGCCATGGCAAAATCTGCTCCATGCGGGACAAGAAAATCTCTTTTCTGGTCTGACGGCGCTTACTGCTGAATTCACTGTCGGCGAAGGTAAGTTGATGACTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCTTAGTGAAGTCATTTTTGTCAAGCAGGTTGATTTTGTAATCAACGAAAGTAGACat                                                          >  1:222603/1‑300 (MQ=255)                                                gggtAAGGTCGCCGGCCATTACCAGCATTGGGGTAAAACGGCTCGATGACTTCCACCATGTTTTGCCATGGCAGAATCTGCTCCATGCGGGACAAGAAAATCTCTTTTCTGGTCTGACGGCGCTTACTGCTGAATTCACTGTCGGCGAAGGTAAGTTGATGACTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCTTAGTGAAGTCATTTTTGTCAAGCAGGTTGATTTTGTAATCAACGAAAGTAGATATATTTttgtt                                                 >  2:264256/4‑300 (MQ=255)                                                ggATAAGGTCGCCGGCCATTACCAGCCTTGGGGTAAAACGGCTCGATGACTTCCACCATGTTTTGCCATGGCAGAATCTGCTCCATGCGGGACAAGAAAATCTCTTTTCTGGTCTGACGGCGCTTACTGCTGAATTCACTGTCGGCGAAGGTAAGTTGATGACTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCTTAGTGAAGTCATTTTTGTCAATCAGGTTGATTTTGTAATCAACGAAAGTAGACATATTTttgtt                                                 >  1:742863/1‑300 (MQ=255)                                                          gCCGGCCATTACCAGCCTTTGGGTAAAACGGCTCGATGACTTCCACCATTTTTTGCCATGGCAGAATCTGCTCCATGCGGGACAAGAAAATCTCTTTTCTGGTCTGACGGCGCTTACTGCTGAATTCACTGTCGGCGAAGGTAAGTTGATGACTCATGATGAACCCTGTTCTATGNCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCTTAGTGAAGTCATTTTTGTCAAGCAGGTTGATTTTGTAATCAANGAAAGTAGACATATTTTTGTTGAAGTAATAg                                       <  2:613788/300‑1 (MQ=255)                                                            cGGCCATTACCAGCCTGGGGGTAAAACGGCTCGATGACTTCCACCATGTTTTGCCATGGCAGAATCTGCTCCATGCGGGACAAGAAAATCTCTTTTCTGGTCTGACGGCGCTTACTGCTGAATTCACTGTCGGCGAAGGTAAGTTGATGACTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCTTAGTGAAGTCATTTTTGTCAAGCAGGTTGATTTTGTAATCAACGAAAGTAGACATATTTTTGTTGAAGTAATAGGt                                     <  1:231395/300‑1 (MQ=255)                                                              gCCATTACCAGCCTTGGGGTAAAACGGCTCGATGACTTCCACCATGTTTTGCCATGGCAGAATCTGCTCCATGCGGGACAAGAAAATCTCTTTTCTGGTCTGACGGCGCTTACTGCTGAATTCACTGTCGGCGAAGGTAAGTTGATGACTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCTTAGTGAAGTCATTTTTGTCAAGCAGGTTGATTTTGTAATCAACGAAAGTAGACATATTTTTGTTGAAGTAATAGGTTg                                   >  1:756283/1‑300 (MQ=255)                                                              gCCATTACCAGCCTTGGGGTAAAACGGCTCGATGACTTCCACCATGTTTTGCCATGGCAGAATCTGCTCCATGCGGGACAAGAAAATCTCTTTTCTGGTCTGACGGCGCTTACTGCTGAATTCACTGTCGGCGAAGGTAAGTTGATGACTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCTTAGTGAAGTCATTTTTGTCAAGCAGGTTGATTTTGTAATCAACGAAAGTAGACATATTTTTGTTGAAGTAATAGGTTg                                   >  2:294792/1‑300 (MQ=255)                                                                 aTTACCAGCCTTGGGGTAAAACGGCTCGATGACTTCCACCATGTTTTGCCATGGCAGAATCTGCTCCATGCGGGACAAGAAAATCTCTTTTCTGGTCTGACGGCGCTTACTGCTGAATTCACTGTCGGCGAAGGTAAGTTGATGACTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCTTAGTGAAGTCATTTTTGTCAAGCAGGTTGATTTTGTAATCAACGAAAGTAGACATATTTTTGTTGAAGTAATAGGTTGCAc                                <  1:443281/300‑1 (MQ=255)                                                                  ttACCAGCCTTGGGGTAAAACGGCTCGATGANTTCCACCATGTTTTGCCATGGCAGAATCTGCTCCATGCGGGACAAGAAAATCTCTTTTCTGGTCTGACGGCGCTTACTGCTGAATTCACTGTNGGCGAAGGTAAGTTGATGACTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCTTAGTGAAGTCATTTTTGTCAAGCAGGTTGATTTTGTAATCAACGAAAGTAGACATATTTTTGTTGAAGTAATAGGTTGCAcc                               >  2:606951/1‑300 (MQ=255)                                                                  ttACCAGCCTTGGGGTAAAACGGCTCGATGACTTCCACCATGTTTTGCCATGGCAGAATCTGCTCCATGCGGGACAAGAAAATCTCTTTTCTGGTCTGACGGCGCTTACTGCTGAATTCACTGTCGGCGAAGGTAAGTTGATGACTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTAATATCAGGGACTTGTTCGCACCTTCATTAGTGAAGTCATTTTTGTCAAGCAGGTTGATTTTGTAATCAACGAAAGTAGACATATTTTTGTTGAAGTAACAGGTTGCAcc                               >  2:143415/1‑300 (MQ=255)                                                                       aGCCTTGGGGTAAAACGACTCGATGACTTCCACCATGTTTTGCCATGGCAGAATCTGCTCCATGCGGGACAAGAAAATCTCCTTTCTGGTCTGACGGCGCTTACTGCTGAATTCACTGTCTGCGAAGGTAAGTTGATTACCCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCTTAGTGAAGTCATTTTTGTCAAGCAGGTTGATTTTGTAATCAACGAAAGTAGACATATTTTTGTTGAAGTAATAGGTTGCACCTACAt                          <  2:36973/300‑1 (MQ=255)                                                                        gCCTTGGGGTAAAACGGCTCGATGACTTCCACCATGTTTTGCCATGGCAGAATCTGCTCCATGCGGGACAAGAAAATCTCTTTTCTGGTCTGACGGCGCTTACTGCTGAATTCACTGTCGGCGAAGGTAAGTTGATGACTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCTTAGTGAAGTCATTTTTGTCAAGCAGGTTGATTTTGTAATCAACGAAAGTAGACATATTTTTGTTGAAGTAATAGGTTGCACCTACATc                         <  1:458031/300‑1 (MQ=255)                                                                        gCCTGGGGGTAAAACGGCTCGATGACTTCCACCATGTTTTGCCATGGCAAAATCTGCTCCATGCGGGACAAGAAAATCTCTTTTCTGGTCTGACGGCGCTTACTGCTGAATTCACTGTCGTCGAAGGTAAGTTGATGACTCACGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATAGCAGGGACTTGTTCGCACCTTCCTTAGTGAAGTCATTTTTGTCAAGCAGGTTGATTTTGTAATCAACGAAAGTAGACATATTTTTGTTGAAGTAATAGGTTGCACCTACATc                         <  2:21476/300‑1 (MQ=255)                                                                         ccTTGGGGTAAAACGGCTCGATGACTTCCAACATGTTTTGCCATGGCAGAATCTGCTCAATGCGGGACAAGAAAATCTATTTTCTGGTCTGACGGCGCTTACTGCTGAATTCACTGTCGGCGAAGGTAAGTTGATGACTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTAATATCAGGGACTTGTTCGCACCTTCCTTAGTGAAGTCATTTTTGTCAAGCAGGTTGATTTTGTAATCAACGAAAGTAGACATATTTTTGTTGAAGTAATAGGTTGCACCTACATCa                        >  2:96002/1‑300 (MQ=255)                                                                              gggTAAAACGGCTCGATGACTTCCACCATGTTTTGCCATGGCAGAATCTGCTCCATGCGGGACAAGAAAATCTCTTTTCTGGTCTGACGGCGCTTACTGCTGAATTCACTGTCGGCGAAGGTAAGTTGATGACTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCTTAGTGAAGTCATTTTTGTCAAGCAGGTTGATTTTGTAATCAACGAAAGTAGACATATTTTTGTTGAAGTAATAGGTTGCACCTACATCAACata                   >  1:490248/1‑300 (MQ=255)                                                                                  aaaaCGGCTCGAGGGTTTCCACCATGTTTTGCCATGTCAGAACCTGCGCCATGCGGGACAAGAAAATCTCTTTTCTGGTCTGACGGCGCTTACTGCTGAATTCACTGTCTGCGAAGGTAAGGTGGTTACTCATGATGAACCTTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCATGGACTTGTTCGCACCTTTTTTAGTGAAGTCATTGGTGTCAAGCAGGTTGATTTTGTAATCAACGAAAGGAGACATATTTTTGTTGAAGTAATAGGTTGCACCTACATCGACATATTTg               <  2:139445/300‑1 (MQ=18)                                                                                    aaCGGCTCGATGACTTCCACCATGTTTTGCCATGGCAGAATCTGCTCCATGCGGGACAAGAAAATCTCTTTTCTGGTCTGACGGCGCTTACTGCTGAATTCACTGTCGGCGAAGGTAAGTTGATGACTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCTTAGTGAAGTCATTTTTGTCAAGCAGGTTGATTTTGTAATCAACGAAAGTAGACATATTTTTGTTGAAGTAATAGGTTGCACCTACATCAACATATTTGAc             <  1:368732/300‑1 (MQ=255)                                                                                               gACTTCCACCATGTTTTGCCATGGCAGAATCTGCTCCATGCGGGACAAGAAAATCTCTTTTCTGGTCTGACGGCGCTTACTGCTGAATTCACTGTCGGCGAAGGTAAGTTGATGACTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCTTAGTGAAGTCATTTTTGTCAAGCAGGTTGATTTTGTAATCAACGAAAGTAGACATATTTTTGTTGAAGTAATAGGTTGCACCTACATCAACATATTTGACTAAGTCCTGAt  <  1:125394/300‑1 (MQ=255)                                                                                               |                                                                                                                                                                                                                                                                                                            GGTTGTACCAATGCTGCATGCAGTGAATGCGTAGCATGGTTTCCAGCGGATAAGGTCGCCGGCCATTACCAGCCTTGGGGTAAAACGGCTCGATGACTTCCACCATGTTTTGCCATGGCAGAATCTGCTCCATGCGGGACAAGAAAATCTCTTTTCTGGTCTGACGGCGCTTACTGCTGAATTCACTGTCGGCGAAGGTAAGTTGATGACTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCTTAGTGAAGTCATTTTTGTCAAGCAGGTTGATTTTGTAATCAACGAAAGTAGACATATTTTTGTTGAAGTAATAGGTTGCACCTACATCAACATATTTGACTAAGTCCTGAT  >  NC_000913/575504‑575897

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