Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
JC JC	NC_000913	332,367	IS2 (–) +5 bp		[pdeL]	[pdeL]

New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913	332367 =	0 (0.000)	12 (0.620)	11/278	0.4	100%	intergenic (+491/‑4)	yahV/pdeL	protein YahV/DNA‑binding transcriptional regulator/c‑di‑GMP phosphodiesterase PdeL
	?	NC_000913	= 1469240	NA (NA)					noncoding (1/1331 nt)	IS2	repeat region
*	?	NC_000913	= 332371	0 (0.000)	21 (1.100)	18/274	0.1	100%	coding (1/1089 nt)	pdeL	DNA‑binding transcriptional regulator/c‑di‑GMP phosphodiesterase PdeL
	?	NC_000913	= 382590	NA (NA)					noncoding (1331/1331 nt)	IS2	repeat region

CCGGAGAGTAACGGAATAGCAGAGAGCTTCGTGAAAACGATAAAGCGTGACTACATCAGTATCATGCCCAAACCAGACGGGTTAACGGCAGCAAAGAACCTTGCAGAGGCGTTCGAGCATTATAACGAATGGCATCCGCATAGTGCGCTGGGTTATCGCTCGCCACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAGGGGCAAATCCAATTACCTATCAGGCAGTTTGCGCCAATCACTACCGGACCCTGAATAACGGCACCATGACAAATACGGGTATTCGCGCCGATAACAACTTCACCTGCAGTTTCATCAATAATTACCTGGTCTGCAATCTGAATATTTTTCTCTTTAACATTAAGCTTCAGATAGTTTTTCAGTTCCTGTGACAGAAGGTTGGGAATAAGCGGTTGACGCTTCCACCAGTGCGGAATAG  >  NC_000913/382352‑382817                                                                                                                                                                                                                                               |                                                                                                                                                                                                                                    ccGGAGAGTAACGGAATAGCAGAGAGCTTCGTGAAAACGATAAAGCGTGACTACATCAGTATCATGCCCAAACCAGACGGGTTAACGGCAGCAAAGAACCTTGCAGAGGCGTTCGAGCATTATAACGAATGGCATCCGCATAGTGCGCCGGGTTATCGCTCGCCACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAGGGGCAAATCCAAttgtttattgg                                                                                                                                                                                                                         >  1:489137/1‑242 (MQ=2)  cGGAGAGTAACGGAATAGCAGAGAGCTTCGTGAAAACGATAAAGCGTGACTACATCAGTATCATGCCCAAACCAGACGGGTTAACGGCAGCAAAGAACCTTGCAGAGGCGTTCGAGCATTATAACGAATGGCAGCCGCATAGTGCGCTGGGTTATCGCTCGCCACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAAGATATGGGCAAATCCAATTACCTATCAgg                                                                                                                                                                                                                        <  1:163234/251‑1 (MQ=21)  cGGAGAGTAACGGAATAGCAGAGAGCTTCGTGAAAACGATAAAGCGTGACTACATCAGTATAATGCCCAAACCAGACGGGTTAACGGCAGCAAAGACCCTTGCAGAGGCGTTCGAGCATTATAACGAATGGCATCCGCATAGTGCGCTGGGTTATCGCTCGCCACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAGGGGCAAATCCAATTACCTATCAgg                                                                                                                                                                                                                        >  1:276470/1‑251 (MQ=21)                               gTGAAAACGATAAAGAGTGACTACATCAGTATCATGCCCAAACCAGACGGGTTAACGGCAGCAAAGAACCTTGCAGAGGCGTTCGAGCATTATAACGAATGGCATCCGCATAGTGCGCTGGGTTATCGCTCGCCACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAGGGGCAAATCCAATTACCTATCAGGCAGTTTGCGCCAATCACTACCGGAcc                                                                                                                                                                                              >  1:154123/1‑248 (MQ=255)                                                                                           gcaAAGAACCTTGCAGAGGCGTTCGAGCATTATAACGAATGGCATCCGCATAGTGCGCTGGGTTACCGCTCGCCACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAGGGGCAAATCCAATTACCTATCAGGCAGTTTGCGCCAATCACTACCTGACCCTGAATAACGGCACCATGACAAATACGGGTATTCGCGCCGATAACAACTTCAc                                                                                                                                         <  1:189293/241‑1 (MQ=255)                                                                                                  aCCTTGCAGAGGCGTTCGAGCATTATAACGAATGGCATCCGCATAGTGCGCTGGGTTATCGCTCGCCACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCGGGAAATATAGGGGCAAATCCAATTACCTATCAGGCAGTTTGCGCCAATCACTACCGGACCCTGAATAACg                                                                                                                                                                                    <  1:191878/191‑1 (MQ=255)                                                                                                              ngTTCGAGCATTATAACGAATGGCATCCGCACAGTGCGCTGGGTTATCGCTCGCCACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAGGGGCAAATCCAATTACCTATCAGGCAGTTTGCGCCAATCACTACCGGACCCTGAATAACGG‑ACCATGACAAATACGGGTATTCGCGCCGATAACAACTTCACCTGCAGTTTCATCAATAATTACCTGGTCTg                                                                                                           >  1:130670/2‑251 (MQ=255)                                                                                                               gTTCGAGCATTATAACGAATGGCAGCCGCATAGTGCGCTGGGTTATCGCTCGCCACGGGAATATCGGCGGCAGCGGGCTTGTAAGGGGTAAAGTGATAACAGATGTCGGGAAATATAGGGGCAAATCCAATTACCTATCAGGCAGTTTGCGCCAATCACTACCGGACCCTg                                                                                                                                                                                           <  1:440322/171‑1 (MQ=18)                                                                                                                                   ggCATCCGCATAGTGCGCTGGGTTATCGCTCGCCACGGGAAGATCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAGGGGCAAATCCAATTACCTATCAGGCAGTTTGCGCCAATCACTACCGGACCCTGAATAACGGCACCATGACAAATACTGTTATTCGCTCCGATAACAACTTCACCTGCAGTTTCATCAATAATTACCTGGTCTGCAATCTGAATATTTTtctc                                                                                        <  1:329318/250‑1 (MQ=255)                                                                                                                                                    cTGGGTTATCGCTCGCCACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAGGGGCAAATCCAATTACCTATCAGGCAGTTTGCGCCAATCACTACCGGACCCTGAATAACGGCACCAGTACAAATACGGGTAGTCGCGcc                                                                                                                                                       <  1:140980/170‑1 (MQ=255)                                                                                                                                                      gggTTATCGCTCGCCACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAGGGGCAAATCCAATTACCTATCAGGCAGTTTGCGCCAATCACAACCGGACCCTGAATAACGGCCCCATGACAAATACGGGTATTCGCGCCGATAACAACTTCACCTGCAGTTTCATCAATAATTACCTGGTCTGCAATCTGAATATTTTTCTCTTTAACATTAAGCTTc                                                                        >  1:189777/1‑247 (MQ=255)                                                                                                                                                        gTTATCGCTCGCCACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAGGGGCAAATCCAATTACCTATCAGGCAGTTTGCGCCAATCACTACCg                                                                                                                                                                                                  >  1:234532/1‑123 (MQ=255)                                                                                                                                                                            tagcggCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATTTCTGGAAAGATAGGGGCAAATCCAATTACCGATCAGTCAGTTTGCGCCAATCACTACCGGACCCTGAATAACGGCACCATGACAAATACggg                                                                                                                                                                 <  1:488805/131‑1 (MQ=25)                                                                                                                                                                                               gTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAGGGGCAAATCCAATTACCTATCAGGCAGTTTGCGCCAATCACTACCGGACCCTGAATAACGGCACCATGACAAATACGGTTTTTCGCTCCGATAACAACGTCACCTGCAGTTTCATCAATAGTTACCTGGTCTGCAATCTTAATAGTTTTCTCTTGAACATTAAGCTTCAGATAGTTTTTCAGTTCCTGTGACAGAAGGTTGGGAATAAGCgg                            <  1:273346/250‑1 (MQ=255)                                                                                                                                                                                                             gATAACAGATGTCTGGAAATATAGGGGCAAATCCAATTACCTATCAGGCAGTTTGCGCCAATCACTACCGGACCCTGAATAACGGCACCATGACAAATACGGGTATTCGCGCCGATAACAACTTCACCTGCAGTTTCATCAATAATTACCTGGTCTGCAATCTGAATATTTTTCTCTTTAACATTAAGCTTCAGATAGTTTTTCAGTTCCTGTGACAGAAGGTTGGGAATAAGCGGTTGACGCTTccacc              >  1:390331/1‑250 (MQ=255)                                                                                                                                                                                                                        tCTGGAAATATAGGGGCAAATCCAATTACCTATCAGGCAGTTTGCGCCAATCACTACCGGACCCTGAATAACGGCACCATGACAAATACGGGTATTCGCGCCGATAACAACTTCan                                                                                                                                         <  1:310805/116‑2 (MQ=255)                                                                                                                                                                                                                        tCTGGAAATATAGGGGCAAATCCAATTACCTATCAGGCAGTTTGCGCCAATCACTACCGGACCCTGAATAACGGCACCATGACAAATACGGGTATTCGCGCCGATAACAACTTCAc                                                                                                                                         >  2:310805/1‑116 (MQ=255)                                                                                                                                                                                                                         atGGAAATAAAGGGGCAAAACCCAATACCTATAAGGCAGTTTGCGCCAATCACTACCGGACCCTGAATAACGGCACCATGACAAATACGGGTAGTCGCGCCGATAACAACTTCACCTGCAGTTTCATCAATAATTACCTGGTCTGCAATCTGAATATTTTTCTCTTTAAAATTAAGCTTCAGATATGTTTTCAGTTCCTGGGACAGAAGGTTGGGAATCAGCGGTTGACGCTTCCACCAGTGCGTAATan  <  1:320596/249‑2 (MQ=255)                                                                                                                                                                                                                                               |                                                                                                                                                                                                                                    CCGGAGAGTAACGGAATAGCAGAGAGCTTCGTGAAAACGATAAAGCGTGACTACATCAGTATCATGCCCAAACCAGACGGGTTAACGGCAGCAAAGAACCTTGCAGAGGCGTTCGAGCATTATAACGAATGGCATCCGCATAGTGCGCTGGGTTATCGCTCGCCACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAGGGGCAAATCCAATTACCTATCAGGCAGTTTGCGCCAATCACTACCGGACCCTGAATAACGGCACCATGACAAATACGGGTATTCGCGCCGATAACAACTTCACCTGCAGTTTCATCAATAATTACCTGGTCTGCAATCTGAATATTTTTCTCTTTAACATTAAGCTTCAGATAGTTTTTCAGTTCCTGTGACAGAAGGTTGGGAATAAGCGGTTGACGCTTCCACCAGTGCGGAATAG  >  NC_000913/382352‑382817

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