Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
JC JC	NC_000913	1,093,481	IS2 (+) +5 bp	coding (606‑610/1359 nt)	dgcT →	putative diguanylate cyclase DgcT

New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913	381260 =	NA (NA)	37 (0.860)	32/482	0.4	97.3%	noncoding (1/1331 nt)	IS2	repeat region
	?	NC_000913	= 1093485	1 (0.020)					coding (610/1359 nt)	dgcT	putative diguanylate cyclase DgcT
*	?	NC_000913	1093481 =	1 (0.020)	34 (0.790)	33/486	0.3	97.1%	coding (606/1359 nt)	dgcT	putative diguanylate cyclase DgcT
	?	NC_000913	= 1651548	NA (NA)					noncoding (706/706 nt)	IS2	repeat region

CGGAGTCCGGAGAGTAACGGAATAGCAGAGAGCTTCGTGAAAACGATAAAGCGTGACTACATCAGTATCATGCCCAAACCAGACGGGTTAACGGCAGCAAAGAACCTTGCAGAGGCGTTCGAGCATTATAACGAATGGCATCCGCATAGTGCGCTGGGTTATCGCTCGCCACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAGGGGCAAATCCACAGGTGAACTACGCTCCTCTGTTTGTTACGCAATAAAAACTGGTGTTTTCGACTATGCAAAACAGTTTCCCTCCTCACGCAATCTGGAAAAATTTGGTGAGGCCCGACAAGATTTAACCATAAAAGAACTGGCTGAAAAATTTCTGGCACTGAAAGAAACTGAAGTCGCCAAAACATCACTCAACACATACCGTGCCGTCATCAAAAATATCCTGAGCATAATCGGTGAAAAAAATCTTGCCTCATCG  >  NC_000913/1651304‑1651796                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                         cGGAGTCCGGAGAGTAACGGAATAGCAGAGAGCTTCGTGAAAACGATAAAGCGTGACTACATCAGTATCATGCCCAAACCAGACGGGTTAACGGCAGCAAAGAACCTTTCAGAGGCGTTCGAGCATTATAACGAATGGCATCCGCATAGTGCGCTGGGTTATCGCTCGCCACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAGGGGCAAATCCACAGGTg                                                                                                                                                                                                                                                    <  1:293873/251‑1 (MQ=14)  ggAGTCCGGAGAGTAACGGAATAGCAGAGAGCTTCGTGAAAACGATAAAGCGTGACTACATCAGTATCATGCCCAAACCAGACGGGTTAACGGCAGCAAAGAACCTTGCAGAGGCGTTCGAGCATTATAACGAATGGCATCCGCATAGTGCGCTGGGTTATCGCTCGCCACGGGAATATCTGCGGCAGCGGGATTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAGGGGCAAATCCACAGGTg                                                                                                                                                                                                                                                    >  2:363194/1‑250 (MQ=17)              gTAAAGGAATATCATAGGGCTTCGTGAAAACGAGAAAGCGTGTCTACAGCAGTATCATGCCCAATCCCGGCGGGTTAACGGCAGCAAAGACCCTTGCAGAGGCGTTCGAGCATTATAGCGAATGGCATCCGCATAGTGCGCTGGGTTATCTCTCGCCACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAGCAGATGTCGGGAAATATAGGTGCAAATCCACAGGTGAACTACGctcctc                                                                                                                                                                                                                                       <  2:266919/251‑1 (MQ=11)                aaCGGAATAGCAGATAGCTTCGTGAAAACGATAAAGCTGGACTACATCAGTATCATGCCCAAACCAGACGGGTTAACGGCAGCAAAGAACCTTGCAGAGGCGTTCGAGCATTATAACGAATGGCATCCGCATAGTGCGCTGGGTTATCGCTCGCCACGGTAATATCTGCGGCAGCGGGCGTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAGGGGCAAATCCACAGGTGAACTACGCTCCTCtg                                                                                                                                                                                                                                     <  2:362973/251‑1 (MQ=255)                        aGCAGAGAGCTTCGTGAAAACGATAAAGCGTGACTACATCAGTATCATGCCCAAACCAGACGGGTTAACGGCAGCAAATAACCTTGCAGAGGCGTTCGTGCATTATAACGAATGGCATCCGCATAGTGCGCTGGGTGAACGCTCTCCACGGGAATAGCTGCGGCAGCGTGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATATGGGCAAATCCACAGGTGAACTACGCTCCTCTGTTTGTTAc                                                                                                                                                                                                                             <  1:219980/251‑1 (MQ=14)                         gCAGAGAGCTTCGTGAAAACGATAAAGCGTGACTACATCAGTATCATGCCCAAACCAGACGGGTTAACGGCAGCAAAGAACCTTGCAGAGGCGTTCGAGCATTATAACGAATGGCATCCGCATAGTGCGCTGGGTTATCGCTCGCCACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAGGGGCAAATCCACAGGTGAACTACGCTCCTCTTTTTGTTACg                                                                                                                                                                                                                            <  2:427222/251‑1 (MQ=255)                             agagCTTCGTGAAAACGAGAAAGCGGGACTAAATCAGTATCATGGCCAAACCATACGGGCTATCGGCAGCAAAGAACCTTGCCGAGGCGTTCGAGCATTATAACGAATGGCATGCGCATAGTGCGCTGGGTTATCGCTCGCCACGGGAATATCTGCGGCAGCGGGCTTGTAATTTGTTAAGTGATAACGGATGTCTGGAAATATAGGGGCAAATCCACAGGTGAACTACGCTCCTCTGTTTGATTCGCAAt                                                                                                                                                                                                                        <  2:220998/251‑1 (MQ=14)                                gctgtgTGAAAGCGATAATGCGTGACTACATCAGTATCATGCCCAAACCCGACGGGTTAACGGCAGCAAAGTACCTTGCAGAGGCGTTCGCGCATTATAACGAATGGCATCCGCAAAGTGCGCTAGGTTGTCGCTCGCCACGGGAATATCTGCGGAAGCGGGCTGTTCATGTGTTAAGTGATAACAGATGTCTGGAAATATAGGAGTAAAACCACAGGTGAACTACGCTCCTCTGTTTTTTCCGCAATaaa                                                                                                                                                                                                                     <  2:25433/246‑1 (MQ=14)                                  ttCGTGAAAACGATAAAGCGTGACTACATCAGTATCATGCCCAAACCAGACGGGTTAACGGCAGCAAAGAACCTTGCAGAGGCGTTCGAGCATTATAACGAATGGCATCCGCATAGTGCGCTGGGTTATCGCTCGCCACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAGGGGCAAATCCACAGGTGAACTACGCTCCTCTGTTTGTTACGCAATaaaa                                                                                                                                                                                                                    >  1:228576/1‑250 (MQ=255)                                          aaCGATAAAGCGTGACTACATCAGTATCATGCCCAAACCAGACGGGTTAACGGCAGCAAAGAACCTTGCAGAGGCGTTCGAGCATTATAACGAATGGCATCCGCATAGTGCGCTGGGTTATCGCTCGCCACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAGGGGCAAATCCACAGGTGAACTACGCTCCTCTGTTTGTTACGCAATAAAAACTGGTGtt                                                                                                                                                                                                           <  1:343924/251‑1 (MQ=255)                                                   gCGTGACTACATCAGTATCATGCCCAAACCAGACGGGTTAACGGCAGCAAAGAACCTTGCAGAGGCGTTCGAGCATTATAACGAATGGCATCCGCATAGTGCGCTGGGTTATCGCTCGCCACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAGGGGCAAATCCACAGGTGAACTACGCTCCTCTGTTTGTTACGCAATAAAAACTGGTGTTTTCGACTAt                                                                                                                                                                                                  <  2:55386/251‑1 (MQ=255)                                                             aTCAGTATCATGCCCAAACCAGACGGGTTAACGGCAGCAAAGAACCTTGCAGAGGCGTTCGAGCATTATAACGAATGGCATCCGCATAGTGCGCTGGTTTATCGCTCGCCACGGGAATATCTTCGGCAGCGTGCTTGTAATGGGTTAAGTGATAACAGATGTCTTGAAATATAGGGGCAAATCCACAGGTGAACTACGCTCctct                                                                                                                                                                                                                                      <  2:213411/205‑1 (MQ=255)                                                             aTCAGTATCATGCCCAAACCAGACGGGTTAACGGCAGCAAAGAACCTTGCAGAGGCGTTCGAGCATTATAACGAATGGCATCCGCATAGTGCGCTGGGTTATCGCTCGCCACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAGGGGCAAATCCACAGGTGAACTACGCTCctct                                                                                                                                                                                                                                      >  1:213411/1‑205 (MQ=40)                                                                             aaCCAGACGGGTTAACGGCAGCAAAGAACCTTGCAGAGGCGTTCGAGCATTATAACGAATGGCATCCGCATAGTGCGCTGGGTTATCGCTCGCCACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAGGGGCAAATCCACAGGTGAACTACGCTCCTCTGTTTGTTACGCAATAAAAACTGGTGTTTTCGACTATGCAAAACAGTTTCCCTCCTCACGCaa                                                                                                                                                                        <  1:363194/251‑1 (MQ=255)                                                                                                    aaGAACCTTGCAGAGGCGTTCGAGCATTATAACGAATGGCATCCGCATAGTGCGCTGGGTTATCGCTCGCCACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAGGGGCAAATCCACAGGTGAACTACGCTCCTCTGTTTGTTACGCAATAAAAACTGGTGTTTTCGACTATGCAAAACAGTTTCCCTCCTCACGAAATCTGGAAAAATTTGGTGAGGccc                                                                                                                                                 <  1:47344/251‑1 (MQ=255)                                                                                                            tGCAGAGGCGTTCGAGCATTATAACGAATGGCATCCGCATAGTGCGCTGGGTTATCGCTCGCCACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAGGGGCAAATCCACAGGTGAACTACGCTCCTCTGTTTGTTACGCAATAAAAACTGGTGTTTTCGACTATGCAAAACAGTTTCCCTCCTCACGCAATCTGGAAAAATTTGGTGAGGCCCGACAAg                                                                                                                                           <  1:41669/249‑1 (MQ=255)                                                                                                             gCAGAGGCGTTCGAGCATTATAACGAATGGCATCCGCATAGTGCGCTGGGTTATCGCTCGCCACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAGGGGCAAATCCACAGGTGAACTACGCTCCTCTGTTTGTTACGCAATAAAAACTGGTGTTTTCGACTATGCAAAACAGTTTCCCTCCTCACGCAATCTGGAAAAATTTGGTGAGGCCCGACAAGAtt                                                                                                                                        <  1:417559/251‑1 (MQ=255)                                                                                                             gCAGAGGCGTTCGAGCATTATAACGAATGGCATCCGCATAGTGCGCTGGGTTATCGCTCGCCACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAGGGGCAAATCCACAGGTGAACTACGCTCCTCTGTTTGTTACGCAATAAAAACTGGTGTTTTCGACTATGCAAAACAGTTTCCCTCCTCACGCAATCTGGAAAAATTTGGTGAGGCCCGACAAGAt                                                                                                                                         <  2:228576/250‑1 (MQ=255)                                                                                                              cAGAGGCGTTCGAGCATTATAACGAATGGCATCCGCATAGTGCGCTGGGTTATCGCTCGCCACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAGGGGCAAATCCACAGGTGAACTACGCTCCTCTGTTTGTTACGCAATAAAAACTGGTGTTTTCGACTATGCAAAACAGTTTCCCTACTCACGCAATCTGGAAAAATTTGGTGAGGCCCGACAAGAttt                                                                                                                                       <  1:381347/251‑1 (MQ=255)                                                                                                                                              ccGCATAGTGCGCTGGGTTATCGCTCGCCACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAGGGGCAAATCCACAGGTGAACTACGCTCCTCTGTTTGTTACGCAATAAAAACTGGTGTTTTCGACTATGCAAAACAGTTTCCCTCCTCACGCAATCTGGAAAAATTTGGTGAGGCCCGACAAGATTTAACCATAAAAGAACTGGCTGAAAAATTTCTgg                                                                                                       >  2:394524/1‑251 (MQ=255)                                                                                                                                                    aGTGCGCTGGGTTATCGCTCGCCACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAGGGGCAAATCCACAGGTGAACTACGCTCCTCTGTTTGTTACGCAATAAAAACTGGTGTTTTCGACTATGCAAAACAGTTTCCCTCCTCACGCAATCTGGAAAAATTTGGTGAGGCCCGACAAGATTTAACCATAAAAGAACTGGCTGAAAAATTTCTGGCACTg                                                                                                  <  1:132435/250‑1 (MQ=255)                                                                                                                                                        cgcTGGGGTATCGCTCGACCCGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAGGGGCAAATCCACAGGTGAACTACGCTCCTCTGTTTGTTACGCAATAAAAACTGGTGGTTTCGACTATGCAAAACAGTTTCCCTCCTCACGCAATCTGGAAAAATTTGGTGAGGCCCGACAAGATTTAACCATAAAAGAACTGGCTGAAAAATTTCTGGCACTgaaaga                                                                                             <  2:455420/251‑1 (MQ=255)                                                                                                                                                                  gcgctcgcCACGGGAAGATCTTCGGCTGCTGGCTTGTAATGGGGTAAGTGATAACAGAGGGCGTGAAATATAGGGGCAAATCCACAGGGGATCTACGCTCCTCTGTTTTTTACGCAATAAAAACTGGTGTTTTCGATTATGCAAAACTGTTTCCCTCCTCACGCAATCTGGAAAAATTTGGTGAGGCCCGACAAGATATAACCATAAAAGAACAGGCTGAAAAATTTCTGGCACTGAAAGAAACTGAAgct                                                                                   <  2:387855/250‑3 (MQ=255)                                                                                                                                                                       cgcCACGGGAATATCTTCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAGGGGCAAATCCACAGGTGAACTACGCTCCTCTGTTTGTTACGCAATAAAAACTGGTGTTTTCGACTATGCAAAACAGTTTCCCTCCTCACGCAATCTGGAAAAATTTGGTGAGGCCCGACAAGATTTAACCATAAAAGAACTGGCTGAAAAATTTCTGGCACTGAAAGAAACTGAAGTCGcc                                                                                <  1:195763/249‑1 (MQ=255)                                                                                                                                                                       cgcCACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAGGGGCAAATCCACAGGTGAACTACGCTCCTCTGTTTGTTACGCAATAAAAACTGGTGTTTTCGACTATGCAAAACAGTTTCCCTCCTCACGCAATCTGGAAAAATTTGGTGAGGCCCGACAAGATTTAACCATAAAAGAACTGGCTGAAAAATTTCTGGCACTGAAAGAAACTGAAGTCGCCaa                                                                              <  1:356393/251‑1 (MQ=255)                                                                                                                                                                                                 gCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAGGGGCAAATCCACAGGTGAACTACGCTCCTCTGTTTGTTACGCAATAAAAACTGGTGTTTTCGACTATGCAAAACAGTTTCCCTCCTCACGCAATCTGGAAAAATTTGGTGAGGCCCGAAAAGATTTAACCATAAAAGAACTGGCTGAAAAATTTCTGGCACTGAAAGAAACTGAAGTCGCCAAAACATCACTCAACACAAACCGTGCCg                                                    >  2:11792/1‑251 (MQ=255)                                                                                                                                                                                                    tGTAAGGGGTTAAGTGATAAAAGATGTTTGGAAATATAGGGGCAAATCCACAGGTGAACTACGCTCCTCTGTTTGTTACGCAATAAAAACTGGTGTTTACGACTATGCAAAACAGTTTCCCTCCTCACGCAATCTGGAAAAATTTGGTGAGGCCCGACAAGATTTAACCATAAAAGAACTGGCTGAAAAATTTCTGGCACTGAAAGAAACTGAAGTCGCCAAAACATCACTCAACACATACCGTGCCGtc                                                  <  2:157339/250‑1 (MQ=255)                                                                                                                                                                                                          gggTTAAGTGATAACAGATGTCTGGAAATATAGGGGCAAATCCACAGGTGAACTACGCTCCTCTGTTTGTTACGCAATAAAAACTGGTGTTTTCGACTATGCAAAACAGTTTCCCTCCTCACGCAATCTGGAAAAATTTGGTGAGGCCCGACAAGATTTAACCATAAAAGAACTGGCTGAAAAATTTCTGGCACTGAAAGAAACTGAAGTCGCCAAAACATCACTCAACACATACCGTGCCGTCATCaaaa                                           <  1:304460/251‑1 (MQ=255)                                                                                                                                                                                                          gggTTAAGTGATAACAGATGTCTGGAAATATAGGGGCAAATCCACAGGTGAACTACGCTCCTCTGTATGTTACGCAATAAAAACTGGTGTTTTCGACTATGCAAAACAGTTTCCCTCTTCACGCAATCTGGAACAATTTGGTGAGGCACGACAAGATTTAACCATAAAAGGACTGGCTGAAAAATTTCTGGCACTGAACAAAACTGAAATCGCCAAAGCAAAAATCAACACATACAGTGCCGACCTCaaa                                            >  2:92316/1‑250 (MQ=255)                                                                                                                                                                                                                                  gAAATATAGGGGCAAATCCACAGGTGAACTACGCTCCTCTGTTTGTTACGCAATAAAAACTGGTGTTTTCGACTATGCAAAACAGTTTCCCTCCTCACGCAATCTGGAAAAATTTGGTGAGGCCCGACAAGATTTAACCATAAAAGAACTg                                                                                                                       >  1:110838/1‑151 (MQ=255)                                                                                                                                                                                                                                  gAAATATAGGGGCAAATCCACAGGTGAACTACGCTCCTCTGTTTGTTACGCAATAAAAACTGGTGTTTTCGACTATGCAAAACAGTTTCCCTCCTCACGCAATCTGGAAAAATTTGGTGAGGCCCGACAAGATTTAACCATAAAAGAACTg                                                                                                                       <  2:110838/151‑1 (MQ=255)                                                                                                                                                                                                                                                   ccACAGGTTATCTACGCTCCTCTGTTTGTTACGCAATAAAAACTGGTGTTTTCGACTATGCAAAACAGTTTCCCTCCTCACGCAATCTGGAAAAATTTGGTGAGGCCCGACAAGATTTAACCATAAAAGAACTGGCTGAAAAATTTCTGGCACTGAAAGAAACTGAAGTCGCCAAAACATCACTCAACACATACCGTGCCGTCATCAAAAATATCCTGAGCATAATCGGTGAAAAAAATCTTGCCTCATCg  <  2:470250/251‑1 (MQ=255)                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                         CGGAGTCCGGAGAGTAACGGAATAGCAGAGAGCTTCGTGAAAACGATAAAGCGTGACTACATCAGTATCATGCCCAAACCAGACGGGTTAACGGCAGCAAAGAACCTTGCAGAGGCGTTCGAGCATTATAACGAATGGCATCCGCATAGTGCGCTGGGTTATCGCTCGCCACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAGGGGCAAATCCACAGGTGAACTACGCTCCTCTGTTTGTTACGCAATAAAAACTGGTGTTTTCGACTATGCAAAACAGTTTCCCTCCTCACGCAATCTGGAAAAATTTGGTGAGGCCCGACAAGATTTAACCATAAAAGAACTGGCTGAAAAATTTCTGGCACTGAAAGAAACTGAAGTCGCCAAAACATCACTCAACACATACCGTGCCGTCATCAAAAATATCCTGAGCATAATCGGTGAAAAAAATCTTGCCTCATCG  >  NC_000913/1651304‑1651796

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.