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)	16 (0.740)	16/490	0.4	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)	25 (1.150)	23/488	0.1	100%	noncoding (1258/1258 nt)	IS3	repeat region
	?	NC_000913	2468298 =	0 (0.000)					coding (85/921 nt)	yfdH	bactoprenol glucosyl transferase

ACCATCAGCACCGACAAACATTAAATCGGCATTAATATCCTGTAATGAACGCTCGGCAATGGAACCATGCATCGAGCGCGTTTTATGGCGGACCGTGCCACCACAGACGACCAGAGTAATGTCTTTATTTTCGGAAAGGGCAAACGCTGCCGGGAGACTGTTGGTAATCACCGTGATATTTTTGGCGGTCATTAATCCTTCGGCGATAAGCATTGTCGTACTTCCGCTGTCGAGAATAACGGTCTGATCTTACCCAGCAATAGTGGACACGCGGCTAAGTGAGTAAACTCTCAGTCAGAGGTGACTCACATGACAAAAACAGTATCAACCAGTAAAAAACCCCGTAAACAGCATTCGCCTGAATTTCGCAGTGAAGCCCTGAAGCTTGCTGAACGCATCGGTGTTACTGCCGCAGCCCGTGAACTCAGCCTGTATGAATCACAACTCTACAACTGGCGCAGTAAACAGCAAAATCAGCAGACGTCTT  >  NC_000913/2169927‑2170413                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                   aCCATCAGCACCGACAAACATTAAATCGGCATTAATATCCTGTAATGAACGCTCGGCAATGGAACCATGCATCGAGCGCGTTTTATGGCGGACCGTGCCACCACAGACGACCAGAGTAATGTCTTTATTTTCGGAAAGGGCAAACGCTGCCGGGAGACTGTTGGTAATCACCGTGATATTTTTGGCGGTCATTAATCCTTCGGCGATAAGCATTGTCGTACTTCCGCTGTCGAGAATAACGGTCTGATCtt                                                                                                                                                                                                                                              <  1:119763/251‑1 (MQ=255)             gACAAACATTAAATCGGCATTAATATCCTGTAATGAACGCTCGGCAATGGAACCATGCATCGAGCGCGTTTTATGGCGGACCGTGCCACCACAGACGCCCAGAGCAATAACTTTATTTTCGGAAAGGGCAACCGCTGCCGGGAGACTGTTGGTAATCACCGTGATATTTTTAGCGGTCATTAATCCTTCGGCGATAAGCATTGTCGTACTTCCGCTGTCGAGAATAACGGTCTGATCTTACCCAGCAAt                                                                                                                                                                                                                                    >  1:19525/1‑249 (MQ=255)             gACAAACATTAAATCGGCATTAATATCCTGTAATGAACGCTCGGCAATGGAACCATGCATCGAGCGCGTTTTATGGCGGACCGTGCCACCACAGACGACCAGAGTAATGTCTTTATTTTCGGAAAGGGCAAACGCTGCCGGGAGACTGTTGGTAATCACCGTGATATTTTTGGCGGTCATTAATCCTTCGGCGATAAGCATTGTCGTAATTCCGCTGTCGCGCATAACGGTCTGATCTTACCCAGCAAt                                                                                                                                                                                                                                    >  2:57720/1‑249 (MQ=255)                    aTTAAATCGGCATTAATATCCTGTAATGAACGCTCGGCAATGGAACCATGCATCGAGCGCGTTTTATGGCGGACCGTGCCACCACAGACGACCAGAGTAATGTCTTTATTTTCGGAAAGGGCAAACGCTGCCGGGAGACTGTTGGTAATCACCGTGATATTTTTGGCGGTCATTAATCCTTCGGCGATAAGCATTGTCGTACTTCCGCTGTCGAGAATAACGGTCTGATCTTACCCAGCAATAGTGGacac                                                                                                                                                                                                                           <  2:134865/251‑1 (MQ=255)                     ttAAATCGGCATTAATATCCTGTAATGAACGCTCGGCAATGGAACCATGCATCGAGCGCGTTTTATGGCGGACCGTGCCACCACAGACGACCAGAGTAATGTCTTTATTTTCGGAAAGGGCAAACGCTGCCGGGAGACTGTTGGTAATCACCGTGATATTTTTGGCGGTCATTAATCCTTCGGCGATAAGCATTGTCGTACTTCCGCTGTCGAGAATAACGGTCTGATCTTACCCAGCAATAGTGGacacc                                                                                                                                                                                                                          >  2:83995/1‑250 (MQ=255)                               aTTAATATCCTGTAATGAACGCTCGGCAATGGAACCATGCATCGAGCGCGTTTTATGGCGGACCGTGCCACCACAGACGACCAGAGTAATGTCTTTATTTTCGGAAAGGGCAAACGCTGCCGGGAGACTGTTGGTAATCACCGTGATATTTTTGGCGGTCATTAATCCTTCGGCGATAAGCATTGTCGTACTTCCGCTGTCGAGAATAACGGTCTGATCTTACCCAGCAATAGTGGACACGCGGCTAAGTg                                                                                                                                                                                                                <  1:225399/251‑1 (MQ=255)                                         tGTAATGAACGCTCGGCAATGGAACCATGCATCGAGCGCGTTTTATGGCGGACCGTGCCACCACAGACGACCAGAGTAATGTCTTTATTTTCGGAAAGGGCAAACGCTGCCGGGAGACTGTTGGTAATCACCGTGATATTTTTGGCGGTCATTAATCCTTCGGCGATAAGCATTGTCGTACTTCCGCTGTCGAGAATAACGGTCTGATCTTACCCAGCAATAGTGGACACGCGGCTAAGTGAGTAAActct                                                                                                                                                                                                      >  1:222271/1‑251 (MQ=255)                                                  cGCTCGGCAATGGAACCATGCATCGAGCGCGTTTTATGGCGGACCGTGCCACCACAGACGACCAGAGTAATGTCTTTATTTTCGGAAAGGGCAAACGCTGCCGGGAGACTGTTGGTAATCACCGTGATATTTTTGGCGGTCATTAATCCTTCGGCGATAAGCATTGTCGTACTTCCGCTGTCGAGAATAACGGTCTGATCTTACCCAGCAATAGTGGACACGCGGCTAAGTGAGTAAACTCTCAGTCagag                                                                                                                                                                                             <  2:96145/251‑1 (MQ=255)                                                                     gCATCGAGCGCGTTTTATGGCGGACCGTGCCACCACAGACGACCAGAGTAATGTCTTTATTTTCGGAAAGGGCAAACGCTGCCGGGAGACTGTTGGTAATCACCGTGATATTTTTGGCGGTCATTAATCCTTCGGCGATAAGCATTGTCGTACTTCCGCTGTCGATAATAACGTTCTGATCTTACCCAGCAATAGTGGACACGCGGCTAAGTGAGTAAACTCTCAGTCAGAGGTGACTCACATGACaaaa                                                                                                                                                                           <  1:211650/250‑1 (MQ=255)                                                                      cATCGAGCGCGTTTTATGGCGGACCGTGCCACCACAGACGACCAGAGTAATGTCTTTATTTTCGGAAAGGGCAAACGCTGCCGGGAGACTGTTGGTAATCACCGTGATATTTTTGGCGGTCATTAATCCTTCGGCGATAAGCATTGTCGTACTTCCGCTGTCGAGAATAACGGTCTGAGCTTACCCAGCAATAGTGGACACGCGGCTAAGTGAGTAAACTCTCAGTCAGAGGTTACTCACATGACAAAAAc                                                                                                                                                                         <  2:137251/251‑1 (MQ=255)                                                                      cATCGAGCGCGTGTTATGGCGGACCGTGCCACCACAGACGACCAGAGTAATGTCTTTATTTTCGGAAAGGGCAAACGCTGCCGGGAGACTGTTGGTAATCACCGTGATATTTTTGGCGGTCATTAATCCTTCGGCGATAAGCATTGTCGTACTTCCGCTGTCGAGAATAACGGTCTGATCTTACCCAGCAATAGTGGACACGCGGCTAAGTGAGTAAACTCTCAGTCAGAGGTGACTCACATGACAAAAAc                                                                                                                                                                         <  1:207775/251‑1 (MQ=255)                                                                              gcgTTTTATGGCGGACCGTGCCACCACAGACGACCAGAGTAATGTCTTTATTTTCGGAAAGGGCAAACGCTGCCGGGAGACTGTTGGTAATCACCGTGATATTTTTGGCGGTCATTAATCCTTCGGCGATAAGCATTGTCGTACTTCCGCTGTCGAGAATAACGGTCTGATCTTACCCAGCAATAGTGGACACGCGGCTAAGTGAGTAAACTCTCAGTCAGAGGTGACTCACATGACAAAAACAGTATCaa                                                                                                                                                                 <  1:164325/251‑1 (MQ=255)                                                                              gcgTTTTATGGCGGACCGTGAAACAACAGACGACCAGAGTAATGTCTTTATTTTCGGAAAGGGCAAACGCCGCCGGGAGACTGTTGGTAATCACCGTGATATTTTTGGCGGTCATTAATCCTTCGGCGATAAGCATTGTCGTACTTCCGCTGTCGAGAATAACGGTCTGATCTTACCCAGCAATAGTGGACACGCGGCTAAGTGAGTAAACTCTCAGTCAGAGGTGACTCACATGACAAAAACAGTaaca                                                                                                                                                                  >  2:211573/1‑247 (MQ=255)                                                                                  tttATGGCGGACCGTGCCACCACAGACGACCAGAGTAATGTCTTTATTTTCGGAAAGGGCAAACGCTGCCGGGAGACTGTTGGTAATCACCGTGATATTTTTGGCGGTCATTAATCCTTCGGCGATAAGCATTGTCGTACTTCCGCTGTCGAGAATAACGGTCTGATCTTACCCAGCAATAGGGGACACGCGGCTAAGTGAGTAAACTCTCAGTCAGAGGTGACTCACATGACAAAAACAGTATCAACCa                                                                                                                                                              <  1:190115/250‑1 (MQ=255)                                                                                       ggCGGACCGTGCCACCACAGACGACCAGAGTAATGTCTTTATTTTCGGAAAGGGCAAACGCTGCCGGGAGACTGTTGGTAATCACCGTGATATTTTTGGCGGTCATTAATCCTTCGGCGATAAGCATTGTCGTACTTCCGCTGTCGAGAATAACGGTCTGATCTTACCCAGCAATAGTGGACACGCGGCTAAGTGAGTAAACTCTCAGTCAGAGGTGACTCACATGACAAAAACAGTATCAACCAGTaaa                                                                                                                                                         >  2:104444/1‑250 (MQ=255)                                                                                                 gCCACCACAGACGACCAGAGTAATGTCTTTATTTTCGGAAAGGGCAAACGCTGCCGGGAGACTGTTGGTAATCACCGTGATATTTTTGGCGGTCATTAATCCTTCGGCGATAAGCATTGTCGTACTTCCGCTGTCGAGAATAACGGTCTGATCTTACCCAGCAATAGTGGACACGCGGCTAAGTGAGTAAACTCTCAGTCAGAGGTGACTCACATGACAAAAACAGAATCAACCAGTAAAAAACCCCGTaa                                                                                                                                              <  1:71179/251‑1 (MQ=255)                                                                                                                    gTAATGTCTTTATTTTCGGAAAGGGCAAACGCTGCCGGGAGACTGTTGGTAATCACCGTGATATTTTTGGCGGTCATTAATCCTTCGGCGATAAGCATTGTCGTACTTCCGCTGTCGAGAATAACGGTCTGATCTTACCCAGCAATAGTGGACACGCGGCTAAGTGAGTAAACTCTCAGTCAGAGGTGACTCACATGACAAAAACAGTATCAACCAGTAAAAAACCCCGTAAACAGCATTCGCCTGAAttt                                                                                                                           >  1:5228/1‑251 (MQ=255)                                                                                                                                                aCGCTGCCGGGAGACTGTTGGTAATCACCGTGATATTTTTGGCGGTCATTAATCCTTCGGCGATAAGCATTGTCGTACTTCCGCTGTCGAGAATAACGGTCTGATCTTACCCAGCAATAGTGGACACGCGGCTAAGTGAGTAAACTCTCAGTCAGAGGTGACTCACATGACAAAAACAGTATCAACCAGTAAAAAACCCCGTAAACAGCATTCGCCTGAATTTCGCAGTGAAGCCCTGAAGCTTGCTGAAc                                                                                               <  1:73901/251‑1 (MQ=255)                                                                                                                                                   ctgacGGGAGACTGTTGGTAATCACCGTGATATTTTTGGCGGTCATTAATCCTTCGGCGATAAGCATTGTCGTACTTCCGCTGTCGAGAATAACGGTCTGATCTTACCCAGCAATAGTGGACACGCGGCTAAGTGAGTAAACTCTCAGTCAGAGGTGACTCACATGACAAAAACAGTATCAACCAGTCAAAAACCCCGTAAACAGCATTCGCCTGAATTTCGCAGTGAAGCCCTGAAGCTTGCTGAACGc                                                                                             >  2:90565/5‑250 (MQ=255)                                                                                                                                                         ggAGACTGTTGGTAATCACCGTGATATTTTTGGCGGTCATTAATCCTTCGGCGATAAGCATTGTCGTACTTCCGCTGTCGAGAATAACGGTCTGATCTTACCCAGCAATAGTGGACACGCGGCTAAGTGAGTAAACTCTCAGTCAGAGGTGACTCACATGACAAAAACAGTATCAACCAGTAAAAAACCCCGTAAACAGCATTCGCCTGAATTTCGCAGTGAAGCCCTGAAGCTTGCTGAACGCATCGgtg                                                                                      >  1:50604/1‑251 (MQ=255)                                                                                                                                                                    gTAATCACCGTGATATTTTTGGCGGTCATTAATCCTTCGGCGATAAGCATTGTCGTACTTCCGCTGTCGAGAATAACGGTCTGATCTTACCCAGCAATAGTGGACACGCGGCTAAGTGAGTAAACTCTCAGTCAGAGGTGACTCACATGACAAAAACAGGATCAACCAGTAAAAAACCCCGTAAACAGCATTCGCCTGAATTTCGCAGTGAAGCCTTGAAGCTTGCTGAACGCATCGGTGTTACTGCCGCa                                                                           <  1:210431/251‑1 (MQ=255)                                                                                                                                                                        tCACCGTGATATTTTTGGCGGTAATTAATCCTTCGGCGATAAGCATTGTCGTACTTCCGCTGTCGAGAATCACGGTCTGATCTTACACAGCAATAGTGGACACGCGGCTAAGTGAGTAAACTCTCAGTCAGAGGTGACTCACATGACAAAAACAGTATCAACCAGTAAAAAACCCCGTAAACAGCATTCGCCTGAATTTCGCAGTGAAGCCCTGAAGCTTGCTGAACGCATCGGTGTTACTGCCGCAGccc                                                                       >  2:205114/1‑251 (MQ=255)                                                                                                                                                                                                    tCCTTCGGCGATAAGCATTGTCGTACTTCCGCTGTCGAGAATAACGGTCTGATCTTACCCAGCAATAGTGGACACGCGGCTAAGTGAGTAAACTCTCAGTCAGAGGTGACTCACATGACAAAAACAGTATCAACCAGTAAAAAACCCCGTAAACAGCATTCGCCTGAATTTCGCAGTGAAGCCCTGAAGCTTGCTGAACGCATCGGTGTTACTGCCGCAGCCCGTGAACTCAGCCTGTATGAATCACAAct                                           >  1:118660/1‑251 (MQ=255)                                                                                                                                                                                                          ggCGATAAGCATTGTCGTACTTCCGCTGTCGAGAATAACGGTCTGATCTTACCCAGCAATAGTGGACACGCGGCTAAGTGAGTAAACTCTCAGTCAGAGGTGACTCACATGACAAAAACAGTATCAACCAGTAAAAAACCCCGTAAACAGCATTCGCCTGAATTTCGCAGTGAAGCCCTGAAGCTTGCTGAACGCATCGGTGTTACTGCCGCAGCCCGTGAACTCAGCCTGTATGAATCACAACTCTACa                                      >  1:113755/1‑250 (MQ=255)                                                                                                                                                                                                                            aCTTCCGCTGTCGAGAATAACGGTCTGATCTTACCCAGCAATAGTGGACACGCGGCTAAGTGAGTAAACTCTCAGTCAGAGGTGACTCACATGACAAAAACAGTATCAACCAGTAAAAAACCCCGTAAACAGCATTCGCCTGAATTTCGCAGTGAAGCCCTGAAGCTTGCTGAACGCATCGGTGTTACTGCCGCAGCCCGTGAACTCAGCCTGTATGAATCACAACTCTACAACTGGCGCAGTAAACAGcc                   >  2:200428/1‑250 (MQ=40)                                                                                                                                                                                                                                        gagaATAACGGTCTGATCTTACCCAGCAATAGTGGACACGCGGCTAAGTGAGTAAACTCTCAGTCAGAGGTGACTCACATGACAAAAACAGTATCAACCAGTAAAAAACCCCGTAAACAGCATTCGCCTGAATTTCGCAGTGAAGCCCTGAAGCTTGCTGAACGCATCGGTGTTACTGCCGCAGCCCGTGAACTCAGCCTGTATGAATCACAACTCTACAACTGGCGCAGTAAACAGCAAAATCAGCAGAc       >  2:122264/1‑251 (MQ=34)                                                                                                                                                                                                                                        gagaATAACGGTCTGATCTTACCCAGCAATAGTGGACACGCGGCTAAGTGAGTAAACTCTCAGTCAGAGGTGACTCACATGACAAAAACAGTATCAACCAGTAAAAAACCCCGTAAACAGCATTCGCCTGAATTTCGCAGTGAAGCCCTGAAGCTTGCTGAACGCATCGGTGTTACTGCCGCAGCCCGTGAACTCAGCCTGTATGAATCACAACTCTACAACTGGCGCAGTAAACAGCAAAATCAGCAGAc       >  1:28154/1‑251 (MQ=34)                                                                                                                                                                                                                                         agaATAACGGTCTGATCTTACCCAGCAATAGTGGACACGCGGCTAAGTGAGTAAACTCTCAGTCAGAGGTGACTCACATGACAAAAACAGTATCAACCAGTAAAAAACCCCGTAAACAGCATTCGCCTGAATTTCGCAGTGAAGCCCTGAAGCTTGCTGAACGCATCGGTGTTACTGCCGCAGCCCGTGAACTCAGCCTGTATGAATCACAACTCTACAACTGGCGCAGTAAACAGCAAAATCAGCAGACg      <  2:74599/251‑1 (MQ=33)                                                                                                                                                                                                                                             tAACGGTCTGATCTTACCCAGCAATAGTGGACACGCGGCTAAGTGAGTAAACTCTCAGTCAGAGGTGACTCACATGACAAAAACAGTATCAACCAGTAAAAAACCCCGTAAACAGCATTCGCCTGAATTTCGCAGTGAAGCCCTGAAGCTTGCTGAACGCATCGGTGTTACTGCCGCAGCCCGTGAACTCAGCCTGTATGAATCACAACTCTACAACTGGCGCAGTAAACAGCAAAATCAGCAGACGtctt  >  2:86974/1‑251 (MQ=32)                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                   ACCATCAGCACCGACAAACATTAAATCGGCATTAATATCCTGTAATGAACGCTCGGCAATGGAACCATGCATCGAGCGCGTTTTATGGCGGACCGTGCCACCACAGACGACCAGAGTAATGTCTTTATTTTCGGAAAGGGCAAACGCTGCCGGGAGACTGTTGGTAATCACCGTGATATTTTTGGCGGTCATTAATCCTTCGGCGATAAGCATTGTCGTACTTCCGCTGTCGAGAATAACGGTCTGATCTTACCCAGCAATAGTGGACACGCGGCTAAGTGAGTAAACTCTCAGTCAGAGGTGACTCACATGACAAAAACAGTATCAACCAGTAAAAAACCCCGTAAACAGCATTCGCCTGAATTTCGCAGTGAAGCCCTGAAGCTTGCTGAACGCATCGGTGTTACTGCCGCAGCCCGTGAACTCAGCCTGTATGAATCACAACTCTACAACTGGCGCAGTAAACAGCAAAATCAGCAGACGTCTT  >  NC_000913/2169927‑2170413

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 24 ≤ ATCG/ATCG < 35 ≤ 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.