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)	28 (1.010)	25/488	0.2	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)	32 (1.160)	28/486	0.2	100%	noncoding (1258/1258 nt)	IS3	repeat region
	?	NC_000913	2468298 =	0 (0.000)					coding (85/921 nt)	yfdH	bactoprenol glucosyl transferase

GGTTGCTGCTACGATAATGCCTGCGTGGAAAGCTTCTTTCATTCGCTGAAAGTGGAATGTATCCATGGAGAACACTTTATCAGCCGGGAAATAATGCGGGCAACGGTGTTTAATTATATCGAATGTGATTACAATCGGTGGCGGCGGCACAGTTGGTGTGGCGGCCTCAGTCCGGAACAATTTGAAAACAAGAACCTCGCTTAGGCCTGTGTCCATATTACGTGGGTAGGATCAGTCATCCCTTCGTGGATCATTTTGACTGCCGCCTGCGCTATTCTGTTTTTTGGCGCGCTGGCGAGCTGAAATCGCTCTTTAAACCCGACTTCCTGGGTCTCGGTTTCACTATTACCAGACATTATTTTCGCCGCACCGCCATGAAAGCGCGTAACCACGCCTTTTTGTTCGAGAAAGCGCAAATCGGCACGGATTGTCGCTTCCGAGGCAGCAAATACTCCCGCCAGATCCTGAAC  >  NC_000913/2171195‑2171664                                                                                                                                                                                                                                          |                                                                                                                                                                                                                                             ggTTGCTGCTACGATAATGCCTGCGTGGAAAGCTTCTTTCATTCGCTGAAAGTGGAATGTATCCATGGAGAACACTTTATCAGCCGGGAAATAATGCGGGCAACGGTGTTTAATTATATCGAATGTGATTACAATCGGTGGCGGCGGCACAGTTGGTGTGGCGGCCTCAGTCCGGAACAATTTGAAAACAAGAACCTCGCTTAGGCCTGTGTCCATATTACGTGGGTAGGATCAGTCATCCCTTCGTGGAt                                                                                                                                                                                                                             <  2:237716/251‑1 (MQ=37)                                                    gTGGAATGTATCCATGGAGAACACTTTATCAGCCGGGAAATAATGCGGGCAACGGTGTTTAATTATATCGAATGTGATTACAATCGGTGGCGGCGGCACAGTTGGTGTGGCGGCCTCAGTCCGGAACAATTTGAAAACAAGAACCTCGCTTAGGCCTGTGTCCATATTACGTGGGTAGGATCAGTCATCCCTTCGTGGATCATTTTGACTGCCGCCTGCGCTATTCTGTTTTTTGGCGCGCTGGCGAGCTg                                                                                                                                                                          <  2:253025/251‑1 (MQ=255)                                                      ggAATGTATCCATGGAGAACAATTTATCAGCCGGGAAATAATGCGGGCAACGGTGTTTAATTATATCGAATGTGATTACAATCGGTGGCGGCGGCACAGTTGGTGTGGCGGCCTCAGTCCGGAACAATTTGAAAACAAGAACCTCGCTTAGGCCTGTGTCCATATTACGTGGGTAGGATCAGTCATCCCTTCGTGGATCATTTTGACTGCCGCCTGCGCTATTCTGTTTTTTGGCGCGCTGGCGAGCTGa                                                                                                                                                                         <  1:193663/250‑1 (MQ=255)                                                                cATGGAGAACACTTTATCAGCCGGGAAATAATGCGGGCAACGGTGTTTAATTATATCGAATGTGATTACAATCGGTGGCGGCGGCACAGTTGGTGTGGCGGCCTCAGTCCGGAACAATTTGAAAACAAGAACCTCGCTTAGGCCTGTGTCCATATTACGTGGGTAGGATCAGTCATCCCTTCGTGGATCATTTTGACTGCCGCCTGCGCTATTCTGTTTTTTGGCGCGCTGGCGCGCTGAAATCGCTCttt                                                                                                                                                              >  2:17818/1‑251 (MQ=255)                                                                     agaACACTTTATCAGCCGGGAAATAATGCGGGCAACGGTGTTTAATTATATCGAATGTGATTACAATCGGTGGCGGCGGCACAGTTGGTGTGGCGGCCTCAGTCCGGAACAATTTGAAAACAAGAACCTCGCTTAGGCCTGTGTCCATATTACGTGGGTAGGATCAGTCATCCCTTCGTGGATCATTTTGACTGCCGCCTGCGCTATTCTGTTTTTTGGCGCGCTGGCGAGCTGAAATCGCTCTTTAAAcc                                                                                                                                                         >  2:49966/1‑251 (MQ=255)                                                                                      gggAAATAATGCGGGCAACGGTGTTTAATTATATCGAATGTGATTACAATCGGTGGCGGCGGCACAGTTGGTGTGGCGGCCTCAGTCCGGAACAATTTGAAAACAAGAACCTCTCTTAGGCCTGTGTCCATATTACGTGGGTAGGATCAGTCATTCCTTCGGGTATCATTGGGATTGCCGCCTGCGCTATTCTGTTTTTTGGCGCGCTGGCTAGCTGAAATCGCTCTCTAAACCCTACGTCCGGGGTCTCg                                                                                                                                        <  2:138747/251‑1 (MQ=17)                                                                                         aaaTAATGCGGGCCAAGGTGTTTAATTATATCGAATGTGATTACAATCGGTGGCGGCGGCACAGTTGGTGTGGCGGCCTCAGTCCGGAACAATTTGAAAACACGAACCTCGCTTAGGCCTGTGTCCATATTACGTGGGTAGGATCAGTCCTCCCTTCGTGGATCCTTTTGAATGCCGCCTGCGCCATTCTGTTTTTTGGCGCGCTGGCGAGCTGAAATCGCTCTTTTAACCCCACTTCCTGGGACTCGGtt                                                                                                                                     >  2:248640/1‑251 (MQ=18)                                                                                                               tAATTATATCGAATGCGATTACAAGCGGAGGCGGCGGCACAGTTGGTGTGGCGGCCTCAGTCCGGAACAATTTGAAAACAAGAACCTCGCGTAGGCCTGTGTCCATATGACGTGGGTATGATCAGTCATCCCTTCGTGGATCATTTTGACGGCCGCCTGCGCTATTCTGTTTTCTGGCGCGCTGGCGAGCTGAAATCGCTCGTTAAACCCGACTTCCTGGGGCTCGGTTTCAATATTCCCAGACattatt                                                                                                                <  2:139583/250‑1 (MQ=25)                                                                                                                aaTTATATCGAATGTGATTACAATCGGTGGCGGCGGCAAAGTTGGTGTGGCGGCCTCAGTCCGGAACAATTTGAAAACAAGAACCTCGCTTAGGCCTGTGTCCATATTACGTGGGTAGGATCAGTCATACCTTCGTGGATCATTTTGACAGCCGCCTGCGCTATTCTGTTTTTTGGCGCGCTGGCGAGCTGTAATCGCTCTTTCAACCCGACTTCCTGGGTCTCGGTTTCACTATTACAAGACATTAtttt                                                                                                              >  1:122230/1‑251 (MQ=255)                                                                                                                  ttATATCGAATGTGATTACAATCGGTGGCGGCGGCACAGTTGGTGTGGCGGCCTCAGTCCTGAACAATTTGAAAACAAGATCCTCGCTTAGGCCTGTGTCCATATTACGTGGTTAGGATCAGTCATCCCTTCGTGGATCATTTTGACTTCCGCCTTCGCTATTCTGGTTTGGGGCGCGCTTGCAAGCTGAAATCGCTCtt                                                                                                                                                               <  1:234726/200‑1 (MQ=18)                                                                                                                  ttATATCGAATGTGATTACAATCGGTGGCGGCGGCACAGTTGGTGTGGCGGCCTCAGTCCGGAACAATTTGAAAACCAGGACCTCGCTTAGGCCTGTGTCCATATTACGTGGGTAGGATCAGTCATCACTTTGTGGATCTTTTTGACTGCCGCCTGCGCTTTTCTGTTTTTTGGCGCGCTGGCGAGCTGAAATCGCTCtt                                                                                                                                                               >  2:234726/1‑200 (MQ=21)                                                                                                                         gAATGTGATTACAATCGGTGGCGGCGGCACAGTTGGTGTGGCGGCCTCAGTCCGGAACAATTTGAAAACAAGAACCTCGCTTAGGCCTGTGTCCATATTACGTGGGTAGGATCAGTCATCCCTTCGAGGATCATTTTGACTGCCGCCTGCGCTATTCTGTTTTTTGGCGCGCTGGCGAGCTGAAATCGCTCTTTAAACCCGACTTCCTGGGTCTCGGTTTCACTATTACCAGACATTATTTTCGCCGCAcc                                                                                                     <  2:195482/251‑1 (MQ=255)                                                                                                                                                                                  cAATTTGAAAACAAGAACCTCGCTTAGGCCTGTGTCCATATTACGGGGGTAGGATCAGTCATCCCTTCGTGGATCATTTTGACTGCCGCCTGCGCTATTCTGTTTTTTGGCGCGCTGGCGAGCTGAAATCGCTCTTTAAACCCGACTTCCTGGGTCTCGGTTTCACTATTACCAGACATTATTTTCGCCGCACCGCCATGAAAGCGCGTAACCACGCCTTTTTGTTCGAGAAAGCGCAAATCGGCACGGAt                                            <  1:17818/251‑1 (MQ=255)                                                                                                                                                                                    aTTTGAAAACAAGAACCTCGCTTAGGCCTGTGTCCATATTACGTGGGTAGGATCAGTCATCCCTTCGTGGATCATTTTGACTGCCGCCTGCGCTATTCTGTTTTTTGGCGCGCTGGCGAGCTGAAATCGCTCTTTAACCCCGACTTCCTGGGTCTCGGTTTCACTATTACCAGACATTATTTTCGCCGCACCGCCATGAAAGCGCGTAACCACGCCTTTTTGTTCGAGAAAGCGCAAATCGGCACGGATTg                                          <  1:220801/251‑1 (MQ=255)                                                                                                                                                                                                      cGCTTAGGCCTGTGTCCATATTACGTGGGTAGGATCAGTCATCCCTTCGTGGATCATTTTGACTGCCGCCTGCGCTATTCTGTTTTTTGGCGCGCTGGCGAGCTGAAATCGCTCTTTAAACCCGACTTCCTGGGTCTCGGTTTCACTATTACCAGACATTATTTTCGCCGCACCGCCATGAAAGCGCGTAACCACGCCTTTTTGTTCGAGAAAGCGCAAATCGGCACGGATTGTCGCTTCCGAGGCAGCaa                        <  1:49966/251‑1 (MQ=255)                                                                                                                                                                                                                            aCGTGGGTAGGATCAGTCATCCCTTCGTGGATCATTTTGACTGCCGCCTGCGCTATTCTGTTTTTTGGCGCGCTGGCGAGCTGAAATCGCTCTTTAAACCCGACTTCCTGGGTCTCGGTTTCACTATTCCCAGACATTATTTTCGCCGCACCGCCATGAAAGCGCGTAACCACGCCTTTTTGTTCGAGAAAGCGCAAATCGGCACGGATTGTCGCTTCCGAGGCAGCAAATACTCCCGCCAGATCCTGAAc  <  1:167297/251‑1 (MQ=255)                                                                                                                                                                                                                                          |                                                                                                                                                                                                                                             GGTTGCTGCTACGATAATGCCTGCGTGGAAAGCTTCTTTCATTCGCTGAAAGTGGAATGTATCCATGGAGAACACTTTATCAGCCGGGAAATAATGCGGGCAACGGTGTTTAATTATATCGAATGTGATTACAATCGGTGGCGGCGGCACAGTTGGTGTGGCGGCCTCAGTCCGGAACAATTTGAAAACAAGAACCTCGCTTAGGCCTGTGTCCATATTACGTGGGTAGGATCAGTCATCCCTTCGTGGATCATTTTGACTGCCGCCTGCGCTATTCTGTTTTTTGGCGCGCTGGCGAGCTGAAATCGCTCTTTAAACCCGACTTCCTGGGTCTCGGTTTCACTATTACCAGACATTATTTTCGCCGCACCGCCATGAAAGCGCGTAACCACGCCTTTTTGTTCGAGAAAGCGCAAATCGGCACGGATTGTCGCTTCCGAGGCAGCAAATACTCCCGCCAGATCCTGAAC  >  NC_000913/2171195‑2171664

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