Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
JC JC	NC_000913	2,468,298	IS3 (+) +3 bp	coding (85‑87/921 nt)	yfdH →	bactoprenol glucosyl transferase

New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913	2170171 =	NA (NA)	23 (0.730)	20/486	0.5	100%	noncoding (1/1258 nt)	IS3	repeat region
	?	NC_000913	= 2468300	0 (0.000)					coding (87/921 nt)	yfdH	bactoprenol glucosyl transferase
*	?	NC_000913	= 2171428	NA (NA)	26 (0.830)	24/484	0.3	100%	noncoding (1258/1258 nt)	IS3	repeat region
	?	NC_000913	2468298 =	0 (0.000)					coding (85/921 nt)	yfdH	bactoprenol glucosyl transferase

CAATACCATCAGCACCGACAAACATTAAATCGGCATTAATATCCTGTAATGAACGCTCGGCAATGGAACCATGCATCGAGCGCGTTTTATGGCGGACCGTGCCACCACAGACGACCAGAGTAATGTCTTTATTTTCGGAAAGGGCAAACGCTGCCGGGAGACTGTTGGTAATCACCGTGATATTTTTGGCGGTCATTAATCCTTCGGCGATAAGCATTGTCGTACTTCCGCTGTCGAGAATAACGGTCTGATCTTACCCAGCAATAGTGGACACGCGGCTAAGTGAGTAAACTCTCAGTCAGAGGTGACTCACATGACAAAAACAGTATCAACCAGTAAAAAACCCCGTAAACAGCATTCGCCTGAATTTCGCAGTGAAGCCCTGAAGCTTGCTGAACGCATCGGTGTTACTGCCGCAGCCCGTGAACTCAGCCTGTATGAATCACAACTCTACAACTGGCGCAGTAAACAGCAAAATCAGCAGACGTCTTC  >  NC_000913/2169923‑2170414                                                                                                                                                                                                                                                         |                                                                                                                                                                                                                                                    cAATACCATCAGCACCGACAAACATTAAATCGGCATTAATATCCTGTAATGAACGCTCGGCAATGGAACCATGCATCGAGCGCGTTTTATGGCGGACCGTGCCACCACAGACGACCAGAGTAATGTCTTTATTTTCGGAAAGGGCAAACGCTGCCGGGAGACTGTTGGTAATCACCGTGATATTTTTGGCGGTCATTAATCCTTCGGCGCTAAGCATTGTCGTACTTCCGCTGTCGAGAATAACGGTCTg                                                                                                                                                                                                                                                    >  1:13618/1‑250 (MQ=255)           aGCACCGACAAACATTAAATCGGCATTAATATCCTGTAATGAACGCTCGGCACTGGAACCATGCATCGAGCGCGTTTTATGGCGGACCGTGCCACCACAGACGACCAGAGTAATGTCTTTATTTTCGGAAAGGGCAAACGCTGCCGGGAGACTGTTGGTAATCACCGTGATATTTTTGGCGGTCATTAATCCTTCGGCGATAAGCATTGTCGTACTTCCGCTGTCGAGAATAACGGTCTGATCTTAccc                                                                                                                                                                                                                                           >  1:157680/1‑249 (MQ=255)                   cAAACATTAAATCGGCATTAATATCCTGTAATGAACGCTCGGCAATGGAACCATGCATCGAGCGCGTTTTATGGCGGACCGTGCCACCACAGACGACCAGAGTAATGTCTTTATTTTCGGAAAGGGCAAACGCTGCCGGGAGACTGTTGGTAATCACCGTGATATTTTTGGCGGTCATTAATCCTTCGGCGATAAGCATTGTCGTACTGCCGCTGTCGAGAATAACGGTCTGATCTTACCCAGCAATAGTg                                                                                                                                                                                                                                 <  2:156087/251‑1 (MQ=255)                    aaaCATTAAAACGGCATTAATATCCTGTAATGAACGCTCGGCAATGGAACCATGCATCGAGCGCGTTTTATGGCGGACCGTGCCACCACAGACGACCAGAGTAATGTCTTTATTTTCGGAAAGGGCAAACGCTGCCGGGAGACTGTGGGGAATCACCGTGAGAGTTGTGGCGGTCATTAATCCTTCGGCGATAAGCATTGTCGGACTGCCGCTGTCGGGAATAGCGGTCGGATCTTAGCCAGGAAGAGTgg                                                                                                                                                                                                                                <  2:199604/251‑1 (MQ=255)                                                aaTGAACGCTCGGCAATGGAACCATGCATCGAGCGCGTTTTATGGCGGCCCGTGCCACCACAGACGACCAGAGTAATGTCTTTATTTTCGGAAAGGGCAAACGCTGCCGGGAGACTGTTGGTAATCACCGTGATATTTTTGGCGGTCATTAATCCTTCGGCGATAAGCATTGTCGTACTTCCGCTGTCGAGAATAACGGTCTGATCTTACCCAGCAATAGTGGCCACGCGGCTAAGTGAGTAAACtctc                                                                                                                                                                                                      >  1:121861/1‑249 (MQ=255)                                                    aaCGCTCGGCAATGGAACCATGCATCGAGCGCGTTTTATGGCGGACCGTGCCACCACAGACGACCAGAGTAATGTCTTTATTTTCGGAAAGGGCAAACGCTGCCGGGAGACTGTTGGTAATCACCGTGATATTTTTGGCGGTCATTAATCCTTCGGCGATACGCATTGTCGTACTTCCGCTGTCGAAAATAACGGTCTGATCTTACCCAGCAATAGTGGACACGCGGCTAAGTGAGTAAACTCtcagtca                                                                                                                                                                                                 >  1:42575/1‑250 (MQ=255)                                                                   aaCCATGCATCGAGCGCCTTTTATGGCGGACCGTGCCACCACAGACGACCAGAGTAATGTCTTTATTTTCGGAAAGGGCAAACGCTGCCGGGAGACTGTTGGTAATCACCGTGATATTTTTGGCGGTCATTAATCCTTCGGCGATAAGCATTGTCGTACTTCCGCTGTCGAGAATAACGGTCTGATCTTACCCAGCAATAGTGGACACGCGGCTAAGTGAGTAAACTCTCAGTCAGAGGTTACTCACATGa                                                                                                                                                                                 <  1:264071/251‑1 (MQ=255)                                                                                                aCCGTGCCACCACAGACGACCAGAGTAATGTCTTTATTTTCGGAAAGGGCAAACGCTGCCGGGAGACTGTTGGTAATCACCGTGATATTTTTGGCGGTCATTAATCCTTCGGCGATAAGCATTGTCGTACTTCCGCTGTCGAGAATAACGGTCTGATCTTACCCAGCAATAGTGGACACGCGGCTAAGTGAGTAAACTCTCAGTCAGAGGTGACTCACATGACAAAAACAGTATCAACCAGTAAAAAAccc                                                                                                                                                    <  2:163482/251‑1 (MQ=255)                                                                                                 ccGTGCCACCACAGACGACCAGAGTAATGTCTTTATTTTCGGAAAGGGCAAACGCTGCCGGGAGACTGTTGGTAATCACCGTGATATTTTTGGCGGTCATTAATCCTTCGGCGATAAGCATTGTCGTACTTCCGCTGTCGAGAATAACGGTCTGATCTTACCCAGCAATAGTGGACACGCGGCTAAGTGAGTAAACTCTCAGTCAGAGGTGACTCACATGACACAAACAGTATCAACCCGTAAAAAAcccc                                                                                                                                                   >  2:295265/1‑251 (MQ=255)                                                                                                                       agTAATGTCTTTATTTTCGGAAAGGGCAAACGCTGCCGGGAGACGGTTGGTAATCACCGTGATATTTTTGGCGGTCATTAATCCTTCGGCGATAAGCATTGTCGTACTTCCGCTGTCGAGAATAACGGTCTGATCTTACCCAGCAATAGTGGACACGCGGCTAAGTGAGTAAACTCTCAGTCAGAGGTGACTCACATGACAAAAACAGTATCAACCAGTAAAAAACCCCGTAAACAGCATTCGCCTGAAtt                                                                                                                             <  2:121861/251‑1 (MQ=255)                                                                                                                                     tttCGGAAAGGGCAAACGCTGATGGGAGAGTTTTGGTAATCACCGTGATCTTTATGGCGGTCATTAATCCTTCGGAGATAAACATTGTCGTACTTCCGCTGTCGAGAATAAAGGTCTGATATTACCCAGCAATAGTGGACACGCGGCTAAGTGAGTCAACTCTCAGTCAGAGGTGACTCACATGACAAAAACAGTATCAACCAGTCACAAACCCCGTAAACATCATTCGCCTGAATTTCGAAGTGAAGccc                                                                                                               >  2:200533/1‑251 (MQ=18)                                                                                                                                             aGGGCAAACGCTGCCGGGAGACTGGTGGTAATCAACGTGATATTTTTGGCGGTCATTAATCCTTCGGCGATAAGCATTGTCGTACTTCCGCTGTCGAGAATAACGGTCTGATCTTACCCAGCAATAGTGGACACGCGGCTAAGTGAGTAAACTCTCAGTCAGAGGTGACTCACATGACAAAAACAGTATCAACCAGTAAAAAACCCCGTAAACAGCATTCGCCTGAATTTCGCAGTGAAGACCTGAAGCtt                                                                                                       <  2:42575/251‑1 (MQ=255)                                                                                                                                                  aaaCGCTGCCGGGAGACTGTTGGTAATCACCGTGATATTTTTGGCGGTCATTAATCCTTCGGCGATAAGCATTGTCGTACTTCCGCTGTCGAGAATAACGGTCTGATCTTACCCCGCAATAGTGGAAACGCGGCTAAGTGAGTAAACTCTCAGTCAGAGGTGACTAACATGACAAAAACAGTATCAACCAGTAAAAAACCCCTTAAACAGCATTAGCCTGAATTTCGCAGTGAAGCCCTGAAGCTTGCTg                                                                                                   >  2:135704/1‑250 (MQ=255)                                                                                                                                                  aaaCGCTGCCGGGAGACTGTTGGTAATCACCGAGATATTTTTGGCGGTCATTAATCCTTCGGCGATAAGCATTGTCGTACTTCCGCTGTCGAGAATAACGGTCTGATCTTACCCAGCAATAGTGGACACGCGGCTAAGTGAGTAAACTCTCAGTCAGAGGTGACTCACATGACAAAAACAGAATCAACCAGTAAAAAACCCAGTAAACAGCATTCGCCTGAATTTCGCAGTGAAACCCTGAAGCTTGCTg                                                                                                   >  2:198478/1‑250 (MQ=255)                                                                                                                                                            gggAGACTGTTGGTAATCACCGTGATATTTTTGGCGGTCATTAATCCTTCGGCGATAAGCATTGTCGTACTTCCGCTGTCGAGAATCACGGTCTGATCTTACCCAGCAATAGTGGACACGCGGCTAAGTGAGTAAACTCTCAGTCAGAGGTGACTCACATGACAAAAACAGTATCAACCAGTAAAAAACCCAGTAAACAGCATTCGCCTGAATTTCGCAGTGAAGCCCTGAAGCTTGCTGAACGCATCggg                                                                                        >  1:300508/1‑250 (MQ=255)                                                                                                                                                              gagaCTGTTGGTAATCACCGTGATATTTTTGGCGGTCATTAATCCTTCGGCGATAAGCATTGTCGTACTTCCGCTGTCGAGAATAACGGTCTGATCTTACCCAGCAATAGTGGACACGCGGCTAAGTGAGTAAACTCTCAGTCAGAGGTGACTCACATGACAAAAACAGTATCAACCAGTAAAAAACCCCGTAAACAGCATTCGCCTGAATTTCGCAGTGAAGCCCTGAAGCTTGCTGAACGCATCGgtgt                                                                                      >  1:266117/1‑251 (MQ=255)                                                                                                                                                                    gTTGGTAATCACCGTGATATTTTTGGCGGTCATTAATCCTTCGGCGATAAGCATTGTCGTACTTCCGCTGTCGAGAATAACGGTCTGATCTTACCCAGCAATAGTGGACACGCGGCTAAGTGAGTAAACTCTCAGTCAGAGGTGACTCACATGACAAAAACAGTATCAACCAGTAAAAAACCCCGTAAACAGCATTCGCCTGAATTTCGCAGTGAAGCCCTGAAGCTTGCTGAACGCATCGGTGTTACTGc                                                                                <  2:13618/251‑1 (MQ=255)                                                                                                                                                                          aaTCACCGTGATATTTATGGCGGTCATTAATCCTTCGGCGATAAGCATTGTCGTACTTCCGCTGTCGAGAATAACGGTCTGATCTTACCCAGCAATAGTGGACACGCGGCTAAGTGAGTAAACTCTCAGTCAGAGGTGACTCACATGACAAAAACAGTATCAACCAGTAAAAAACCCCGTAAACAGCATTCGCCTGAATTTCGCAGTGAAGCACTGAAGCTTGCTGAACGCATCGGTGTTACTGACGCAGc                                                                          >  2:93595/1‑251 (MQ=255)                                                                                                                                                                             cACCGTGATATTTTTGGCGGTCATTAATCCTTCGGCGATAAGCATTGTCGTACTTCCGCTGTCGAGAATAACGGTCTGATCTTACCCAGCAATAGTGGACACGCGGCTAAGTGAGTAAACTCTCAGTCAGAGGTGCCTCACATGACAAAAACAGTATCAACCAGTAAAAAACCCCGTAAACAGCATTCGCCTGAATTTCGCAGTGAAGCCCTGAAGCTTGCTGAACGCATCGGTGTTCCTGCCGCAGCCCg                                                                       >  1:197566/1‑251 (MQ=255)                                                                                                                                                                             cACCGTGATATTTTTGGCGGTCATTAATCCTTCGGCGATAAGCATTGTCGTACTTCCGCTGTCGAGAATAACGGTCTGATCTTACCCAGCAATAGTGGACACGCGGCTAAGTGAGTAAACTCTCAGTCAGAGGTGACTCACATGACAAAAACAGTATCAACCAGTAAAAAACCCCGTAAACAGCATTCGCCTGAATTTCGCAGTGAAGCCCTGAAGCTTGCTGAACGCATCGGTGTTACTGCCGCAGCCCg                                                                       <  2:6963/251‑1 (MQ=255)                                                                                                                                                                                            ggCGGTCATTAATCCTTCGGCGATAAGCATTGTCGTACTTCCGCTGTCGAGAATAACGGTCTGATCTTACCCAGCAATAGTGGACACGCGGCTAAGTGAGTAAACTCTCAGTCAGAGGTGACTCACATGACAAAAACAGTATCAACCAGTAAAAAACCCCGTAAACAGCATTCGCCTGAATTTCGCAGTGAAGCCCTGAAGCTTGCTGAACGCATCGGTGTTACTGCCGCAGCCCGTGAACTCAGCCTGt                                                         >  1:3092/1‑250 (MQ=255)                                                                                                                                                                                                                                    ccGCTGTCGAGAATAACGGTCTGATCTTACCCAGCAATAGTGGACACGCGGCTAAGTGAGTAAACTCTCAGTCAGAGGTGACTCACATGACAAAAACAGTATCAACCAGTAAAAAACCCCGTAAACAGCATTCGCCTGAATTTCGCAGTGAAGCCCTGAAGCTTGCTGAACGCATCGGTTTTACTGCCGCAGCCCGTGAACTCAGCCTGTATGAATCACAACTCTACAACTGGCGCAGTAAACAGCaaaa                 <  1:52292/250‑1 (MQ=38)                                                                                                                                                                                                                                      gCTGTCGAGAATAACGGTCTGATCTTACCCAGCAATAGTGGACCCGCGGCTAAGTGAGTAAACTCTCAGTCAGAGGCGACTCACATGACAAAAACAGTATCAACCAGTAAAAAACCCCGTAAACAGCATTCGCCTGAATTTCGCAGTGAAGCCCTGAAGATTGCTGAACGCATCGGTGTTACTGCCGCAGCCCGTGAACTCAGCCTGTATGAATCACAACTATACAACTGGCGCAGTAATCAGCAAAATc               >  2:241605/1‑250 (MQ=255)                                                                                                                                                                                                                                                  aaCGGTCTGATCTTACCCAGCAATAGTGGACACGCGGCTAAGTGAGTAAACTCTCAGTCAGAGGTGACTCACATGACAAAAACAGTATCAACCAGTAAAAAACCCCGTAAACAGCATTCGCCTGAATTTCGCAGTGAAGCCCTGAAGCTTGCTGAACGCATCGGAGTTACTGCCGCAGCACGTGAACTCAGCCGGTATGAGTCACAACTCTACAACTGGCGCAGTAACCAGCAAAATCAGCAGACGtcttc  >  1:6561/1‑251 (MQ=9)                                                                                                                                                                                                                                                         |                                                                                                                                                                                                                                                    CAATACCATCAGCACCGACAAACATTAAATCGGCATTAATATCCTGTAATGAACGCTCGGCAATGGAACCATGCATCGAGCGCGTTTTATGGCGGACCGTGCCACCACAGACGACCAGAGTAATGTCTTTATTTTCGGAAAGGGCAAACGCTGCCGGGAGACTGTTGGTAATCACCGTGATATTTTTGGCGGTCATTAATCCTTCGGCGATAAGCATTGTCGTACTTCCGCTGTCGAGAATAACGGTCTGATCTTACCCAGCAATAGTGGACACGCGGCTAAGTGAGTAAACTCTCAGTCAGAGGTGACTCACATGACAAAAACAGTATCAACCAGTAAAAAACCCCGTAAACAGCATTCGCCTGAATTTCGCAGTGAAGCCCTGAAGCTTGCTGAACGCATCGGTGTTACTGCCGCAGCCCGTGAACTCAGCCTGTATGAATCACAACTCTACAACTGGCGCAGTAAACAGCAAAATCAGCAGACGTCTTC  >  NC_000913/2169923‑2170414

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 13 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑
Reads not counted as support for junction
read_name Not counted due to insufficient overlap past the breakpoint.
read_name Not counted due to not crossing MOB target site duplication.