Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
MC JC	NC_000913	1,095,503	Δ10,045 bp	IS3‑mediated	ycdU–[ymdA]	ycdU, serX, ghrA, ycdX, ycdY, ycdZ, csgG, csgF, csgE, csgD, csgB, csgA, csgC, [ymdA]

Missing coverage evidence...
			seq id	start	end	size	←reads	reads→	gene	description
*	*	÷	NC_000913	1094415–1095503	1105547	10045–11133	27 [25]	[1] 84	[insF4]–[ymdA]	[insF4],insE4,ycdU,serX,ghrA,ycdX,ycdY,ycdZ,csgG,csgF,csgE,csgD,csgB,csgA,csgC,[ymdA]

New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913	315229 =	NA (NA)	83 (0.870) +T	61/596	0.3	98.8%	noncoding (1/1255 nt)	IS3	repeat region
	?	NC_000913	1105548 =	1 (0.010)					coding (135/312 nt)	ymdA	uncharacterized protein YmdA

AATGCGTTAAACCCAAATGAGATCTATCGACTGGCTGATAACGCACTTTACGAGGCGAAAGAGACCGGGCGTAATAAGGTGGTTGTGAGGGATGTGGTGAATTTTTGTGAGTCACCATAAAGCGGCATTTTGATCCTACCCACGTAATATGGACACAGGCCTAAGCGAGGTTCTTGTTTTCAAATTGTTCCGGACTGAGGCCGCCACACCAACTGTGCCGCCGCCACCGATTGTAATCACATTCGATATAATTAAACACCGTTGCCCGCATTATTTCCCGGCTGATAAAGTGTTCTCCATGGATACATTCCACTTTCAGCGAATGAAAGAAGCTTTCCACGCAGGCATTATCGTAGCAGCAACCTTTTGCGCTCATACTTCCACGCAGATTATGCCGCTTCAGTTGCGCCTGATAATCTGCTGAACAG  >  NC_000913/1094115‑1094542                                                                                                                                                                                                                                                                                                            |                                                                                                                                 aaTGCGTTAAACCCAAATGAGATCTATCGACTGGCTGATAACGCACTTTACGAGGCGAAAGAGACCGGGCGTAATAAGGTGGTTGTGAGGGATGTGGTGAATTTTTGTGAGTCNCCATAAAGCGGCATTTTGATCCTACCCACGTAATATGGACACAGGCCTAAGCGAGGTTCTTNTTTTCAAATTGTTCCGGACTGAGGCCGCCACACCAACTGTGCCGCCGCCACCGATTGTAATCACATTCGATATAATTAAACACCGTTGCCCGCATTATTTCCCGGCTGATAAAGTGTTCTCCAt                                                                                                                                  <  2:615413/300‑1 (MQ=255) aaTGCGTTAAACCCAAATGAGATCTATCGACTGGCTGATAACGCACTTTACGAGGCGAAAGAGACCGGGCGTAATAAGGTGGTTGTGAGGGATGTGGTGAATTTTTGTGAGTCACCATAAAGCGGCATTTTGATCCTACCCACGTAATATGGACACAGGCCTAAGCGAGGTTCTTGTTTTCAAATTGTTCCGGACTGAGGCCGCCACACCAACTGTGCCGCCGCCACCGATTGTAATCACATTCGATATAATTAAACACCGTTGCCCGCATTATTTCCCGGCTGATAAAGTGTTCTCCAt                                                                                                                                  <  1:486845/300‑1 (MQ=255)      gTTAAACCCAAATGAGATCTATCGACTGGCTGATAACGCACTTTACGAGGCGAAAGAGACCGGGCGTAATAAGGTGGTTGTGAGGGATGTGGTGAATTTTTGTGAGTCACCATAAAGCGGCATTTTGATCCTACCCACGTAATATGGACACAGGCCTAAGCGAGGTTCTTGTTTTCAAATTGTTCCGGACTGAGGCCGCCACACCAACTGTGCCGCCGCCACCGATTGTAATCACATTCGATATAATTAAACACCGTTGCCCGCATTATTTCCCGGCTGATAAAGTGTTCTCCATGGATa                                                                                                                             <  1:358013/300‑1 (MQ=255)          aaCCCAAATGAGATCTATCGACTGGCTGATAACGCACTTGACGAGGCGAAAGAGACCGGGCGTAATAAGGTGGTTGTGAGGGATGTGGTGAATTTTTGTGAGTCANCATAAAGCGGCATTTTGATCCTACCCACGTAATATGGACACACGCCTAAGCGAGGTTCTTGTTTTCAAATTGTTCCGGACTGAGGCCGCCACACCAACTGTGCCGCCGCCACCGATTGTAATCACATTCGATATAATTAAACACCGTTGCCCGCATTATTTCCCGGCTGATAAAGTGTTCTCCATGGATACAtt                                                                                                                         <  2:661130/300‑1 (MQ=255)                                  gCTGATAACGCACTTTACGAGGCGAAAGAGACCGGGCGTAATAAGGTGGTTGTGAGGGATGTGGTGAATTTTTGTGAGTCACCATAAAGCGGCATTTTGATCCTACCCACGTAATATGGACACAGGCCTAAGCGAGGTTCTTGTTTTCAAATTGTTCCGGACTGAGGCCGCCACACCAACTGTGCCGCCGCCACCGATTGTAATCACATTCGATATAATTAAACACCGTTGCCCGCATTATTTCCCGGCTGATAAAGTGTTCTCCATGGATACATTCCACTTTCAGCGAATGAAAGAAGc                                                                                                 >  2:587607/1‑300 (MQ=255)                                   cTGATAACGCACTTTACTAGGCGAAAGAGACCGGGCGTAATAAGGTGGTTGTGAGGGATGTGGTGAATTTTTGTGAGTCACCATAAAGCGGCATTTTGATCCTACCCACGTAATATGGACACAGGCCTAAGCGAGGTTCTTGTTTTCAAATTGTTCCGGACTGAGGCCGCCACACCAACTGTGCCGCCGCCACCGATTGTAATCACATTCGATATAATTAAACACCGTTGCCCGCATTATTTCCCGGCTGATAAAGTGTTCTCCATGGATACATTCCACTTTCAGCGAATGAAAGAAGCt                                                                                                <  2:343318/300‑1 (MQ=255)                                      aTAACGCACTTTACGAGGCGAAAGAGACCGGGCGTAATAAGGTGGTTGTGAGGGATGTGGTGAATTTTTGTGAGTCACCATAAAGCGGCATTTTGATCCTACCCACGTAATATGGACACAGGCCTAAGCGAGGTTCTTGTTTTCAAATTGTTCCGGACTGAGGCCGCCACACCAACTGTGCCGCCGCCACCGATTGTAATCACATTCGATATAATTAAACACCGTTGCCCGCATTATTTCCCGGCTGATAAAGTGTTCTCCATGGATACATTCCACTTTCAGCGAATGAAAGAAGCTTTc                                                                                             >  1:573986/1‑300 (MQ=255)                                                     aGGCGAAAGAGACCGTGCGTAATAAGGTGGTTGTGAGGGATGTGGTGAATTTTTGTGAGTCACCATAAAGCGGCATTTTGATCCTACCCACTTAATATGGACACAGGCCTAAGCGAGGTTCTTGTTTTCAAATTGTTCCGGACTGAGGCCGCCACACCAACTGTGCCGCCGCCACCGATTGTAATCACATTCGATATAATTAAACACCGTTGCCCGCATTATTTCCCGGCTGATAAAGTGTTCTCCATGGATACATTCCACTTTCAGCGAATGAAAGAAGCTTTCCACGCAGGCATTATc                                                                              <  2:410156/300‑1 (MQ=255)                                                     aGGCGAAAGAGACCGGGCGTAATAAGGTGGTTGTGAGGGATGTGGTGAATTTTTGTGAGTCACCATAAAGCGGCATTTTGATCCTACCCACGTAATATGGACACAGGCCTAAGCGAGGTTCTTGTTTTCAAATTGTTCCGGACTGAGGCCGCCACACCAACTGTGCCGCCGCCACCGATTGTAATCACATTCGATATAATTAAACACCGTTGCCCGCATTATTTCCCGGCTGATAAAGTGTTCTCCATGGATACATTCCACTTTCAGCGAATGAAAGAAGCTTTCCACGCAGGCATTATc                                                                              <  1:43597/300‑1 (MQ=255)                                                     aGGCGAAAGAGACCGGGCGTAATAAGGTGGTTGTGAGGGATGTGGTGAATTTTTGTGAGTCACCATAAAGCGGCATTTTGATCCTACCCACGTAATATGGACACAGGCCNAAGCGAGGTTCTTGTTTTCAAATTGTTCCGGACTGAGGCCGCCACACCAACTGTGCCGCCGCCACCGATTGTAATCACATTCGATATAATTAAACACCGTTGCCCGCATTATTTCCCGGCTGATAAAGTGTTCTCCATGGATACATTCCACTTTCAGCGAATGAAAGAAGCTTTCCACGCAGGCATTATc                                                                              <  1:532509/300‑1 (MQ=255)                                                       gCGAAAGAGACCGGGCGTAATAAGGTGGTTGTGAGGGATGTGGTGAATTTTTGTGAGTCACCATAAAGCGGCATTTTGATCCTACCCACGTAATATGGACACAGGCCTAAGCGAGGTTCTTGTTTTCAAATTGTTCCGGACTGAGGCCGCCACACCAACTGTGCCGCCGCCACCGATTGTAATCACATTCGATATAATTAAACACCGTTGCCCGCATTATTTCCCGGCTGATAAAGTGTTCTCCATGGATACATTCCACTTTCAGCGAATGAAAGAAGCTTTCCACGCAGGCATTATCGt                                                                            <  1:363729/300‑1 (MQ=255)                                                         gAAAGAGACCGGGCGTAATAAGGTGGTTGTGAGGGATGTGGTGAATTTTTGTGAGTCACCATAAAGCGGCATTTTGATCCTACCCACGTAATATGGACACAGGCCTAAGCGAGGTTCTTGTTTTCAAATTGTTCCGGACTGAGGCCGCCACACCAACTGTGCCGCCGCCACCGATTGTAATCACATTCGATATAATTAAACACCGTTGCCCGCATTATTTCCCGGCTGATAAAGTGTTCTCCATGGATACATTCCACTTTCAGCGAATGAAAGAAGCTTTCCACGCAGGCATTATCGTag                                                                          >  1:87112/1‑300 (MQ=255)                                                                aCCGGGCGTAATAAGGTGGTTGTGAGGGATGTGGTGAATTTTTGTGAGTCACCATAAAGCGGCATTTTGATCCTACCCACGTAATATGGACACAGGCCTAAGCGAGGTTCTTGTTTTCAAATTGTTCCGGACTGAGGCCGCCACACCAACTGTGCCGCCGCCACCGATTGTAATCACATTCGATATAATTAAACACCGTTGCCCGCATTATTTCCCGGCTGATAAAGTGTTCTCCATGGATACATTCCACTTTCAGCGAATGAAAGAAGCTTTCCACGCAGGCATTATCGTAGCAGCAAc                                                                   <  1:725271/300‑1 (MQ=255)                                                                    ggCGTAATAAGGTGGTTGTGAGGGATGTGGTGAATTTTTGTGAGTCACCATAAAGCGGCATTTTGATCCTACCCACGTAATATGGACACAGGCCTAAGCGAGGTTCTTGTTTTCAAATTGTTCCGGACTGAGGCCGCCACACCAACTGTGCCGCCGCCACCGATTGTAATCACATTCGATATAATTAAACACCGTTGCCCGCATTATTTCCCGGCTGATAAAGTGTTCTCCATGGATACATTCCACTTTCAGCGAATGAAAGAAGCTTTCCACGCAGGCATTATCGTAGCAGCAACCttt                                                               >  1:475349/1‑300 (MQ=255)                                                                         aagaaGGTGGTTGTGAGGGATGTGGTGAATTTTTGTGAGTCACCATAAAGCGGCATTTTGATCCTACCCACGTAATATGGACACAGGCCTAAGCGAGGTTCTTGTTTTCAAATTGTTCCGGACTGAGGCCGCCACACCAACTGTGCCGCCGCCACCGATTGTAATCACATTCGATATAATTAAACACCGTTGCCCGCATTATTTCCCGGCTGATAAAGTGTTCTCCATGGATACATTCCACTTTCAGCGAATGAAAGAAGCTTTCCACGCAGGCATTATCGTAGCAGCAACCTTTTgcgc                                                          <  2:590117/297‑1 (MQ=255)                                                                                  gtTGTGAGGGATGTGGTGAATTTTTGTGAGTCACCATAAAGCGGCATTTTGATCCTACCCACGTAATATGGACACAGGCCTAAGCGAGGTTCTTGTTTTCAAATTGTTCCGGACTGAGGCCGCCACACCAACTGTGCCGCCGCCACCGATTGTAATCACATTCGATATAATTAAACACCGTTGCCCGCATTATTTCCCGGCTGATAAAGTGTTCTCCATGGATACATTCCACTTTCAGCGAATGAAAGAAGCTTTCCACGCAGGCATTATCGTAGCAGCAACCTTTTGCGCTCATACTTc                                                 >  1:480223/1‑300 (MQ=255)                                                                                                 gtgAATTTTTGTGAGTCACCATAAAGCGGCATTTTGATCCTACCCACGTAATATGGACACAGGCCTAAGCGAGGTTCTTGTTTTCAAATTGTTCCGGACTGAGGCCGCCNCACCAACTGTGCCGCCGCCACCGATTGTAATCACATTCGATATAATTAAACACCGTTGCCCGCATTATTTCCCGGCTGATAAAGTGTTCTCCATGGATACATTCCACTTTCAGCGAATGAAAGAAGCTTTCCACGCAGGCATTATCGTAGCAGCAACCTTTTGCGCTCATACTTCCACGCAGATTATGcc                                  <  1:545196/300‑1 (MQ=255)                                                                                                 gtgAATTTTTGTGAGTCACCATAAAGCGGCATTTTGATCCTACCCACGTAATATGGACACAGGCCTAAGCGAGGTTCTTGTTTTCAAATTGTTCCGGACTGAGGCCGCCACACCAACTGTGCCGCCGCCACCGATTGTAATCACATTCGATATAATTAAACACCGTTGCCCGCATTATTTCCCGGCTGATAAAGTGTTCTCCATGGATACATTCCACTTTCAGCGAATGAAAGAAGCTTTCCACGCAGGCATTATCGTAGCAGCAACCTTTTGCGCTCATACTTCCACGCAGATTATGcc                                  <  2:687144/300‑1 (MQ=255)                                                                                                 gtgAATTTTTGTGAGTCACCATAAAGCGGCATTTTGATCCTACCCACGTAATATGGACACAGGCCTAAGCGAGGTTCTTGTTTTCAAATTGTTCCGGACTGAGGCCGCCACACCAACTGTGCCGCCGCCACCGATTGTAATCACATTCGATATAATTAAACACCGTTGCCCGCATNATTTCCCGGCTGATAAAGTGTTCTCCATGGATACATTCCACTTTCAGCGAATGAAAGAAGCTTTCCACGCAGGCATTATCGTAGCAGCAACCNTTTGCGCTCATACTTCCACGCAGATTATGcc                                  <  2:548295/300‑1 (MQ=255)                                                                                                     aTTTTTGTGAGTCACCATAAAGCGGCATTTTGATCCTACCCACGTAATATGGACACAGGCCTAAGCGAGGTTCTTGTTTTCAAATTGTTCCGGACTGAGGCCGCCACACCAACTGTGCCGCCGCCACCGATTGTAATCACATTCGATATAATTAAACACCGTTGCCCGCATTATTTCCCGGCTGATAAAGTGTTCTCCATGGATACATTCCACTTTCAGCGAATGAAAGAAGCTTTCCACGCAGGCATTATCGTAGCAGCAACCTTTTGCGCTCATACTTCCACGCAGATTATGCCGCtt                              >  1:378890/1‑300 (MQ=40)                                                                                                           gtgAGTCACCATAAAGCGGCATTTTGATCCTACCCACTTAATATGGACACAGGCCTAATCGAGGTTCTTGTTTTCAAATTGTTCCGGACTGAGGCCGCCACACCAACTGTGCCGCCGCCACCGATTGTAATCACATTCGATATAATTAGACACCGTTGCCCGCATTATTTCCCGGCTGATAAAGTGTTCTCCATGGATACATTCCACTTTCAGCGAATGAAAGAAGCTTTCCACGCAGGCATTATCGTAGCAGCAACCTTTTGCGCTCATACTTCCACGCAGATTATGCCGCTTCAGTTg                        <  2:308130/300‑1 (MQ=255)                                                                                                               gTCACCATAAAGCGGCATTTTGATCCTACCCACGTAATATGGACACAGGCCTAAGCGAGGTTCTTGTTTTCAAATTGTTCCGGACTTAGGCCGCCACACCAACTGTGCCGCCGCCACCGATTGTAATCACATTCGATATAATGAAACACCGTTGCCCGCATTATTTCCCGGCTGATAAAGTGTTCTCCATGGATACATTCCACTTTCAGCGAATGAAAGAAGCTTTCCACGCAGGCATTATCGTAGCAGCAACCTTTTGCGCTCATACTTCCACGCAGATTATGCCGCTTCAGTTGCGcc                    <  2:7533/300‑1 (MQ=37)                                                                                                                   ccATAAAGCGGCATTTTGATCCTACCCACGTAATATGGACACAGGCCTAAGCGAGGTTCTTGTTTTCAAATTGTTCCGGACTGAGGCCGCCACACCAACTGTGCCGCCGCCACCGATTGTAATCACATTCGATATAATTAAACACCGTTGCCCGCATTATTTCCCGGCTGATAAAGTGTTCTCCATGGATNCATTCCACTTTCAGCGAATGAAAGAAGCTTTCCACGCAGGCATTATCGTAGCAGCAACCTTTTGCGCTCATACTTCCACGCAGATTATGCCGCTTCAGTTGCGCCTGAt                >  1:560687/1‑300 (MQ=25)                                                                                                                   ccATAAAGCGGCATTTTGATCCTACCCACGTAATATGGACACAGGCCTAAGCGAGGTTCTTGTTTTCAAATTGTTCCGGACTGAGGCCGCCACACCAACTGTGCCGCCGCCACCGATTGTAATCACATTCGATATAATTAAACACCGTTGCCCGCATTATTTCCCGGCTGATAAAGTGTTCTCCATGGATACATTCCACTTTCAGCGAATGAAAGAAGCTTTCCACGCAGGCATTATCGTAGCAGCAACCTTTTGCGCTCATACTTCCACGCAGATTATGCCGCTTCAGTTGCGCCTGAt                >  1:88723/1‑300 (MQ=35)                                                                                                                               tgtttGATCCTACCCACGTAATATGGACACATGCCTAAGCTAGGTTCTTGTTTTCAAATTGTTCCGGACTGAGGCCGCCACACCAACTGTGCCGCCGCCACCGATTGTAATCACATTCGATATAATTAAACACCGTTGCCCGCATTATTTCCCGGCTGATAAAGTGTTCTCCATGGATACATTCCACTTTCAGCGAATGAAAGAAGCTTTCCACGCAGGCATTATCGTAGCAGCAACCTTTTGCGCTCATACTTCCACGCAGATTATGCCGCTTCAGTTGCGCCTGATAATCTGCTGAAc    <  2:584539/298‑1 (MQ=11)                                                                                                                               aTTTTGATCCTACCCACGTAATATGGACACAGGCCTAAGCGAGGTTCTTGTTTTCAAATTGTTCCGGACTGAGGCCGCCACACCAACTGTGCCGCCGCCACCGATTGTAATCACATTCGATATAATTAAACACCGTTGCCCGCATTATTTCCCGGCTGATAAAGTGTTCTCCATGGATACATTCCACTTTCAGCGAATGAAAGAAGCTTTCCACGCAGGCATTATCGTAGCAGCAACCTTTTGCGCTCATACTTCCACGCAGATTATGCCGCTTCAGTTGCGCCTGATAATCTGCTGAAc    <  1:422084/300‑1 (MQ=31)                                                                                                                                 tttGATCCTACCCACGTAATATGGACACAGGCCTAAGCGAGGTTCTTGTTTTCAAATTGTTCCGGACTGAGGCCGCCACACCAACTGTGCCGCCGCCACCGATTGTAATCACATTCGATATAATTAAACACCGTTGCCCGCATTATTTCCCGGCTGATAAAGTGTTCTCCATGGAGACATTCCACTTTCAGCGAATGAAAGAAGCTTTCCACGCATGCATTATCGTAGCAGCAACCTTTTGCGCTCATACTTCCACGCAGATTATGCCGCTTCAGTTGCGCCTGATAATCTGCTGAACAg  <  2:334970/300‑1 (MQ=11)                                                                                                                                                                                                                                                                                                            |                                                                                                                                 AATGCGTTAAACCCAAATGAGATCTATCGACTGGCTGATAACGCACTTTACGAGGCGAAAGAGACCGGGCGTAATAAGGTGGTTGTGAGGGATGTGGTGAATTTTTGTGAGTCACCATAAAGCGGCATTTTGATCCTACCCACGTAATATGGACACAGGCCTAAGCGAGGTTCTTGTTTTCAAATTGTTCCGGACTGAGGCCGCCACACCAACTGTGCCGCCGCCACCGATTGTAATCACATTCGATATAATTAAACACCGTTGCCCGCATTATTTCCCGGCTGATAAAGTGTTCTCCATGGATACATTCCACTTTCAGCGAATGAAAGAAGCTTTCCACGCAGGCATTATCGTAGCAGCAACCTTTTGCGCTCATACTTCCACGCAGATTATGCCGCTTCAGTTGCGCCTGATAATCTGCTGAACAG  >  NC_000913/1094115‑1094542

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