Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
JC JC	NC_000913	3,960,601	IS2 (+) +5 bp	intergenic (‑788/‑72)	ppiC ← / → rep	peptidyl‑prolyl cis‑trans isomerase C/ATP‑dependent DNA helicase Rep

New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913	3960601 =	0 (0.000)	46 (1.130)	44/488	0.1	100%	intergenic (‑788/‑76)	ppiC/rep	peptidyl‑prolyl cis‑trans isomerase C/ATP‑dependent DNA helicase Rep
	?	NC_000913	= 4499511	NA (NA)					noncoding (1331/1331 nt)	IS2	repeat region
*	?	NC_000913	= 3960605	0 (0.000)	52 (1.280)	48/488	0.1	100%	intergenic (‑792/‑72)	ppiC/rep	peptidyl‑prolyl cis‑trans isomerase C/ATP‑dependent DNA helicase Rep
	?	NC_000913	4498181 =	NA (NA)					noncoding (1/1331 nt)	IS2	repeat region

AACAGTAATCATCCCGGGCAAATGCCCGGGAATTATTCTAGGATTATTTTCTTTGTTAAAAAAGACAAACGGTATTAACTGATGTATTTACTATTTACCGCTCCCTGCAGGGGGGATTATCAGAATCCGCGGTTATTCATTATTTATATACATTCAATAAAAAAGTAACCTCATAAAAAATCCTGTACTCACAGGGTGATGCGGGATAACTGACGTAAGGAGTGCTGCAGGGATTTGATGGTCCCTTGTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGTTTTTAGACTAGTCATTGGAGAACAGATGATTGATGTCTTAGGGCCGGAGAAACGCAGACGGCGTACCACACAGGAAAAGATCGCAATTGTTCAGCAGAGCTTTGAACCGGGGATGACGGTCTCCCTCGTTGCCCGGCAACATGGTGTAGCAGCCAGCCAGTTATTTCTCTGGCGTAAGCAAT  >  NC_000913/4497933‑4498428                                                                                                                                                                                                                                                         |                                                                                                                                                                                                                                                        aaCAGTAATCATCCCGGGCAAATGCCCGGGAATTATTCTAGGATTATTTTCTTTGTTAAAAAAGACAAACGGTATTAACTGATGTATTTACTATTTACCGCTCCCTGCAGGGGGGATTATCAGAATCCGCGGTTATTCATTATTTATATACATTCAATAAAAAAGTAACCTCATAAAAAATCCTGTACTCACAGGGTGATGCGGGATAACTGACGTAAGGAGTGCTGCAGGGATTTGATGGTCCCTTGTAg                                                                                                                                                                                                                                                       <  2:54124/251‑1 (MQ=255) aaCAGTAATCATCACGGGCAAATGCCCGGGAATAATTCTAGGATTATTTTCTTTGTTAAAAAAGACAAACGGTATTAACTGATGTATTTACTATTTACCGCTCCCTGCAGGGGGGATTATCAGAATCCGCGGTTATTCATTATTTATATACATTCAATACAAACGTAACCTCATAAAAAATCCTGTACTCACAGGGTGATGCGGGATAACTGACGTAAGGAGTGCTGCAGGGATTTGATGGTCCCTTGt                                                                                                                                                                                                                                                         >  2:366916/1‑249 (MQ=255)   cAGTAATCATCCCGGGCAAATGCCCGGGAATTATTCTCGGATTATTTTCTTTGTTAAAAAAGACAAACGGTATTAACTGATGTATTTACTATTTACCGCTCCCTGCAGGGGGGATTATCAGAATCCGCGGTTATTCATTATTTATATACATTCAATAAAAAAGTAACCTCATAAAAAATCCTGTACTAACCGGGTGATGCGGGATAACTGCAGTAAGGAGTGCTGCAGGGATTTGATGGTCCCTTGTAGAc                                                                                                                                                                                                                                                     >  2:184445/1‑251 (MQ=255)      tAATCATCCCGGGCAAATGCCCGGGAATTATTCTAGGATTATTTTCTTTGTTAAAAAAGACAAACGGTATTAACTGATGTATTTACTATTTACCGCTCCCTGCAGGGGGGATTATCAGAATCCGCGGTTATTCATTATTTATATACATTCAATAAAAAAGTAACCTCATAAAAAATCCTGTACTCACAGGGTGATGCGGGATAACTGACGTAAGGAGTGCTGCAGGGATTTGATGGTCCCTTGTAGACTgg                                                                                                                                                                                                                                                  >  2:6603/1‑251 (MQ=255)                                gttaaacTAGGATGATTTTATTTGTAACAGATGATAAACGGTATTAACTGATGTATTTAATATTTACGGCTACCGGCAGAGGGGCTTATCAGAATTAGCTGTTACTCATAATTTATATACATTCAATAAAAAAGTAACCTCATAAAAAATACTGTAATCACAGGGTGATGCGGGATAAATACCGTAAGGAGAGCTGACGGGATTTGATAGTCTCTTGTAGACTGACCTCCTGAATCTCCAGACAACCAat                                                                                                                                                                                                                         >  2:312194/7‑250 (MQ=255)                                      cTAGGATTATTTTCTTTGTTAAAAAAGACAAACGGCATTAACTGATGTATTTACTATTTACCGCTCCCTGCAGGGGGGATTATCAGAATCCGCGGTTATTCATTATTTATATACATTCAATAAAAAAGTAACCTCATAAAAAATCCTGTACTCACAGGGTGATGCGGGATAACTGACGTAAGGAGTGCTGCAGGGATTTGATGGTCCCTTGTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACtt                                                                                                                                                                                                                  >  1:43124/1‑251 (MQ=255)                                                                          aTTAACTGATGTATTTACTATTTACCGCTCCCTGCAGGGGGGATTATCAGAATCCGCGGTTATTCATTATTTATATACATTCAATAAAAAAGTAACCTCATAAAAAATCCTGTACTCACAGGGTGATGCGGGATAACTGACGTAAGGAGTGCTGCAGGGATTTGATGGTCCCTTGTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGTTTTTAGACt                                                                                                                                                                               >  2:354680/1‑250 (MQ=255)                                                                                                              ggggggATTATCAGAATCCGCGGTTATTCATTATTTATATACATTCAATAAAAAAGTAACCTCATAAAAAATCCTGTACTCACAGGGTGATGCGGGATAACTGACGTAAGGAGTGCTGCAGGGATTTGATGGTCCCTTGTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACGTAAATAAGTGATAGTCTTAATGCTAGTTTTTAGACTAGTCATTGGAGAACAGATGATTGATGTCTTATGGCCg                                                                                                                                          <  2:27695/251‑1 (MQ=255)                                                                                                                          aGAATCCGCGGTTATTCATTATTTATATACATTCAATAAAAAAGTAACCTCATAAAAAATCCTGTACTCACAGGGTGATGCGGGATAACTGACGTAAGGAGTGCTGCAGGGATTTGAGGGTCCCTTGTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGTTTTTAGACTAGTCATTGGAGAACAGATGATTGATGTCTTAGGGCCGGAGAAACGCAGa                                                                                                                              <  1:366916/251‑1 (MQ=255)                                                                                                                                         tCATTCTTTATATACATTCAATAAAAAAGTAACCTCATAAAAAATCCTGTACTCACAGGGTGATGCGGGATAACTGACGTAAGGAGTGCTGCAGGGATTTGATGGTCCCTTGTAGACTGGCCCCATGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGTTTTTAGACTAGTCATTGGAGAACAGATGATTGATGTCTTAGGACCGGAGAAACGCAGACGGCGTACcacac                                                                                                                 >  2:365747/1‑249 (MQ=255)                                                                                                                                                         ttCAATAAAAAAGTAACCTCATAAAAAATCCTGTACTCACAGGGTGATGCGGGATAACTGACGTAAGGAGTGCTGCAGGGATTTGATGGTCCCTTGTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGTTTTTAGACTAGTCATTGGAGAACAGATGATTGATGTCTTAGGGCCGGAGAAACGCAGACGGCGTACCACACAGGAAAAGATCGCAATTg                                                                                               >  2:223689/1‑251 (MQ=255)                                                                                                                                                                          cTCATAAAAAATCCTGTACTCACAGGGTGATGCGGGATAACTGACGTAAGGAGTGCTGCAGGGATTTGATGGTCCCTTGTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGTTTTTAGACTAGTCATTGGAGAACAGATGATTGATGTCTTAGGGCCGGAGAAACGCAGACGGCGTACCACACAGGAAAAGATCGCAATTGTTCAGCAGAGCTTTGa                                                                               >  1:71504/1‑250 (MQ=255)                                                                                                                                                                               aaaaaaTCCTGTACTCACAGGGTGATGCGGGATAACTGACGTAAGGAGTGATGCAGGGATTTGATGGTCCCTTGTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGTTTTTAGACTAGTCATTGGAGAACAGATGATTGATGTCTTAGGGCCGGAGAAACGCAGACGGCGTACCACACAGGAAAAGATCGCAATTGTTCAGCAGAGCTTTGAACCggg                                                                         >  1:86865/1‑251 (MQ=255)                                                                                                                                                                                 aaaaTCCTGTACTCACAGGGTGATGCGGGATAACTGACGTAAGGAGTGCTGCAGGGATTTGATGGTCCCTTGTAGACTGGCCCCCTGCATCTCCAGACAACCAATATCACTTAAATAAGAGATAGTCTTAATACTAGTTTTTAGACTAGTCATTGGAGAACAGATGATTGATGTCTTAGGGCCGGAGAAACGCAGACGGCGTACCACACAGGAAAAGATCGCAATTGTTCAGCAGAGCTTTGAACCgggg                                                                        >  1:73256/1‑250 (MQ=255)                                                                                                                                                                                     tCCTGTACTCACAGGGTGATGCGGGATAACTGACGTAAGGAGTGCTGCAGGGATTTGATGGTCCCTTGTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGTTTTTAGACTAGTCATTGGAGAACAGATGATTGATGTCTTAGGGCCGGAGAAACGCAGACGGCGTACCACACAGGAAAAg                                                                                                         >  2:357212/1‑213 (MQ=255)                                                                                                                                                                                     tCCTGTACTCACAGGGTGATGCGGGATAACTGACGTAAGGAGTGCTGCAGGGATTTGATGGTCCCTTGTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGTTTTTAGACTAGTCATTGGAGAACAGATGATTGATGTCTTAGGGCCGGAGAAACGCAGACGGCGTACCACACAGGAAAAg                                                                                                         <  1:357212/213‑1 (MQ=255)                                                                                                                                                                                      ccTGTACTCACAGGGTGATGCGGGATAACTGACGTAAGGAGTGCTGCAGGGATTTGATGGTCCCTTGTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATaa                                                                                                                                                                                                            <  1:53771/113‑1 (MQ=255)                                                                                                                                                                                      ccTGTACTCACAGGGTGATGCGGGATAACTGACGTAAGGAGTGCTGCAGGGATTTGATGGTCCCTTGTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGTTTTTAGACTAGTCATTGGAGAACAGATGATTGATGTCTTAGGGCCGGAGAAACGCAGACGGCGTACCACACAGGAAAAGATCGCAATTGTTCAGCAGAGCTTTGAACCGGGGATGACg                                                                  >  1:319325/1‑251 (MQ=255)                                                                                                                                                                                      ccTGTACTCACAGGGTGATGCGGGATAACTGACGTAAGGAGTGCTGCAGGGATTTGATGGTACCTTGTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATaa                                                                                                                                                                                                            >  2:53771/1‑113 (MQ=255)                                                                                                                                                                                            cTCACAGGGTGATGCGGGATAACTGACGTAAGGAGTGCTGCAGGGATTTGATGGTCCCTTGTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGTTTTTAGACTAGTCATTGGAGAACAGATGATTGATGTCTTAGGGCCGGAGAAACGCAGACGGCGTACCACACAGGAAAAGATCGCAATTGTTCAGCAGAGCTTTGAACCGGGGATGACGGTCTcc                                                            >  1:101256/1‑251 (MQ=255)                                                                                                                                                                                                caGGGTGGTGCGGGATAACTGACGTAAGTAGTGCTTCAGATATTTGATGGTCCCTTGTAGACTGGCCCCCTTAATCGCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGTTTTGAGACTAGTCATTGGAGAACAGATGAGGGATGTCTTAGTGCCGGAGAAACGCAGACGGCGTACCGCACAGGAAAAGATCGCAATTGTTCAGCAGATCTTTGAACCGTGGATGACGGTCTCCCTCg                                                        <  2:265820/251‑1 (MQ=14)                                                                                                                                                                                                   ggTGATGCGGGATAACTGACGTAAGGAGTGCTGCAGGGATTTGATGGTCCCTTGTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATTCTAGTTTTTAGACTAGTCATTGGAGAACAGATGATTGATGTCTTAGGGCCGGTGAAACGCAGACGGCGTACCACACAGGAAAAGATCGCAATTGTTCAGCAGTGCTTTGAACCGGGGATTACGGTCTCCCTCGTTg                                                     <  1:340586/251‑1 (MQ=255)                                                                                                                                                                                                    gTGATGCGGGATAACTGACGTAAGGAGTGCTGCAGGGATTTGATGGTCCCTTGTAGACTGGCCCCCTGAATCTCCAGACAACAAATATCACTTAAATAAGAGATAGTCTTAATACTAGTTTTTAGACTAGTCATTGGAGAACAGATGATTGATGTCTTAGGGCCGGAGAAAAGCAGACGGCGTACCACACAGGAAAAGATCGCAATTGTTCAGCAGAGCTTTGAACCGGGGATGACGGTCTCCCTCGTTGc                                                    >  1:262172/1‑251 (MQ=255)                                                                                                                                                                                                                aaCTGACGTAAGGAGTGCTGCAGGGATTTGATGGTCCCTTGTAGACTGGCCCCCTGAATCTCCAGACAACCAATAGCACTTAAATAAGTGATAGTCTTAATACTAGTTTTTAGACTAGTCATTGGAGAACAGATGATTGATGTCTTAGGGCCGGAGAAACGCAGACGGCGTACCACACAGGAAAAGATCGCAATTGTTCAGCAGAGCTTTGAACCGGGGATGACGGTCTCCCTCGTTGCCCGGCAACATgg                                        <  2:101256/251‑1 (MQ=255)                                                                                                                                                                                                                     aCGTAAGGAGTGCTGCAGGGAGTTGATGGTCCCTTGGAGACTGGCCCTCTGAATCGCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACGAGTTTTTAGACTAGTCATTGGAGAACAGATGATTGATGTCTTAGGGCCGGAGAAACGCAGACGGCGTGCCAGACAGGAAAAGTTCGCAATTGTTTGGCATAGGTTTGAACCGTGGATGTCGGTCTCCCTCGTTGCCCGGCAACATGGTGTag                                   <  1:248313/251‑1 (MQ=14)                                                                                                                                                                                                                       gTAAGGAGTGCTGCAGGGATTTGATGGTCCCTTGTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGTTTTTAGACTAGTCATTGGAGAACAGATGATTGATGTCTTAGGGCCGGAGAAAAGCAGACGGCGTACCACACAGGAAAAGATCGCAATTGTTCAGCAGAGCTTTGAACAGGGGATGACGGTCTCCCTCGTTGCCCGGCAACATGGTGTagc                                  >  1:360388/1‑250 (MQ=255)                                                                                                                                                                                                                        tAAGGAGTGCTGCAGGGATTTGATGGTCCCTTGTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGTTTTTAGACTAGTCATTGGAGAACAGATGATTGATGTCTTAggg                                                                                                                                             >  1:320171/1‑142 (MQ=255)                                                                                                                                                                                                                        tAAGGAGTGCTGCAGGGATTTGATGGTCCCTTGTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGTTTTTAGACTAGTCATTGGAGAACAGATGATTGATGTCTTAggg                                                                                                                                             <  2:320171/142‑1 (MQ=255)                                                                                                                                                                                                                                             gATGGTCCCTTGTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGTTTTTAGACTAGTCATTGGAGAACAGATGATTGATGTCTTAGGGCCGGAGAAACGCAGACGGCGTACCACACAGGAAAAGATCGCAATTGTTCAGCAGAGCTTTGAACCGGGGATGACGGTCTCCCTCGTTGCCCGGCAACATGGTGTAGCAGCCAGCCAGTTATTTCTCTgg            <  1:254599/250‑1 (MQ=35)                                                                                                                                                                                                                                                       tGTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGTTTTTAGACTAGTCATTGGAGAACAGATGATTGATGTCTTAGGGCCGGAGAAACGCAGACGGCGTACCACACAGGAAAAGATCGCAATTGTTCAGCAGAGCTTTGAACCGGGGATGACGGTCTCCCTCGTTGCCCGGCAACATGGTGTAGCAGCCAGCCAGTTATTTCTCTGGCGTAAGCAAt  >  1:202517/1‑250 (MQ=11)                                                                                                                                                                                                                                                       tGTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGTTTTTAGACTAGTCATTGGAGAACAGATGATTGATGTCTTAGGGCCGGAGAAACGCAGACGGCGTACCACACAGGAAAAGATCGCAATTGTTCAGCAGAGCTTTGAACCGGGGATGACGGTCTCCCTCGTTGCCCGGCAACATGGTGTAGCAGCCAGCCAGTTATTTCTCTGGCGTAAGCAAt  >  1:330321/1‑250 (MQ=11)                                                                                                                                                                                                                                                       tGTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACGTAAATAAGTGATAGTCTTAATACTAGTTTTTAGACTAGTCATTGGAGAACAGATGATTGATGTCTTAGGGCCGGAGAAACGCAGACGGCGTACCACACAGGAAAAGATCGCAATTGTTCAGCAGAGCTTTGAACCGGGGATGACGGTCTCCCTCGTTGCCCGGCAACATGGTGTAGCAGCCAGCCAGTTATTTCTCTGGCGTAAGCaa   <  2:33959/249‑1 (MQ=11)                                                                                                                                                                                                                                                        gTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCTGTTAAATAAGTGATAGTCTTAATACTAGGTGTTAGACTAGTCATTGGAGAACAGATGATTGATGTCTTAGGGCCGGAGAAACGCAGACGGCGTACCCCACAGGAAAAGATCGCAACTGGTCAGCAGAGCTTTGAACGGTGGATGACGTTCTCCCTCGTTTCCCGGCAGCATGGTGGAGCAGCCAGCCAGTTATTTCTCGTGCGTAAg      <  1:35577/245‑1 (MQ=11)                                                                                                                                                                                                                                                        gTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCGCTTAAATAAGTGATAGTCTTAATACTAGTGTTTAGACTAGTCATTGGAGAACAGATGATTGATGTCTTAGGGCCGGAGAAACGCAGACGGCGTACCACACAGGAAAAGATCGCAATTGTTCAGCAGAGCTTTGAACCGGGGATGACGGTCTCCCTCGTTGCCCGGCAACATGGTGTAGCAGCCAGCCAGTTATTTCTCTGGCGTAAGCa    <  2:143372/247‑1 (MQ=11)                                                                                                                                                                                                                                                        gTAGACTGGCCCACTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGTTTTTAGACTAGTCATTGGAGAACAGATGATGGATGTCTTAGGGTCGGAGAAACGCAGACGGCGTACCACACAGGAAAAGATCGCAATTGTTCAGCAGATCTTTGAACCGGGGATGACGGTCTCCCTCGTTGCCCGGCAACATGGTGTAGCAGCCAGCCAGTTATTTCTCTGGCGTAAGc     <  2:209713/246‑1 (MQ=2)                                                                                                                                                                                                                                                         |                                                                                                                                                                                                                                                        AACAGTAATCATCCCGGGCAAATGCCCGGGAATTATTCTAGGATTATTTTCTTTGTTAAAAAAGACAAACGGTATTAACTGATGTATTTACTATTTACCGCTCCCTGCAGGGGGGATTATCAGAATCCGCGGTTATTCATTATTTATATACATTCAATAAAAAAGTAACCTCATAAAAAATCCTGTACTCACAGGGTGATGCGGGATAACTGACGTAAGGAGTGCTGCAGGGATTTGATGGTCCCTTGTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGTTTTTAGACTAGTCATTGGAGAACAGATGATTGATGTCTTAGGGCCGGAGAAACGCAGACGGCGTACCACACAGGAAAAGATCGCAATTGTTCAGCAGAGCTTTGAACCGGGGATGACGGTCTCCCTCGTTGCCCGGCAACATGGTGTAGCAGCCAGCCAGTTATTTCTCTGGCGTAAGCAAT  >  NC_000913/4497933‑4498428

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.