Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
JC JC	NC_000913	3,486,010	IS2 (–) +5 bp :: Δ32 bp +GTCTTAATACTAGTTTTTAGACTATCACTTATTTAAGTGATATTGGTTGTCTGGAGATTCAGGGGGCCAGTCTA	intergenic (‑192/‑106)	yhfA ← / → crp	OsmC family protein YhfA/DNA‑binding transcriptional dual regulator CRP

New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913	3186128 =	NA (NA)	11 (0.660) +GTCTTAATACTAGTTTTTAGACTATCACTTATTTAAGTGATATTGGTTGTCTGGAGATTCAGGGGGCCAGTCTA	10/348	0.7	100%	noncoding (33/1331 nt)	IS2	repeat region
	?	NC_000913	3486010 =	0 (0.000)					intergenic (‑192/‑110)	yhfA/crp	OsmC family protein YhfA/DNA‑binding transcriptional dual regulator CRP
*	?	NC_000913	= 3187426	NA (NA)	27 (1.140)	24/496	0.1	100%	noncoding (1331/1331 nt)	IS2	repeat region
	?	NC_000913	= 3486014	0 (0.000)					intergenic (‑196/‑106)	yhfA/crp	OsmC family protein YhfA/DNA‑binding transcriptional dual regulator CRP

TTTCAAAGCTCGTATGGAATCCCCAAGCAAAGATGCCACCCCGGGCAACGTTTACGCTCAGGCTGACTACAAGATTGCTTACGAGTAATCGCATTTGGCCGGGTTCGCCCGGCCACATCATTAACGGATTAATGATAAGTGGATCAGATGTATAAAAAATTAAAATTAACCACAATAAGCGAATTGATTAAAAATATTTATTGTTCATTATCCGTTATTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGTTTTTAGACTAGTCATTGGAGAACAGATGATTGATGTCTTAGGGCCGGAGAAACGCAGACGGCGTACCACACAGGAAAAGATCGCAATTGTTCAGCAGAGCTTTGAACCGGGGATGACGGTCTCCCTCGTTGCCCGGCAACATGGTGTAGCAGCCAGCCAGTTATTTCTCTGGCGTAAGC  >  NC_000913/3185878‑3186340                                                                                                                                                                                                                                                           |                                                                                                                                                                                                                     tttCAAAGCTCGTATGGAATCCCCAAGCACTGATGCCACCCCGGGCAACGTTTACGCTCAGGCTGACTCCAAGATTGCTTACGAGTAATCGCATTTGGCCGGGTTCGCCCGGCCACATCATTAACGGATTAATGATAAGTGGATCAGATGTATAAAAAATTAAAATTAACCACAATAAGCGAATTGATTAAAAATATTTATTGTTCATTATCCGTTATTAGACTGGCCCCCTGAATCTCCAGACAACCAat                                                                                                                                                                                                                      >  2:112146/1‑251 (MQ=255)     aaatcTCGTATGGAATCCACAAGCAGAGTTGCCCCCAAGGGCAACGTTTACGCTCAGGCTGACTACAAGATTGCTTACGAGAAATCGCCTTTGGCCGGGTTCGCCCGGCCACATCATTAACGGATTAATGATAAGTGGATCAGATGTATAAACAATTAAAATTAACCACAATAAGCGAATTGATTAAACATATTTATTGTTCATTATCCGTTATTAGACTGGCCCCCTGACTCTCCAGACAACCAATATc                                                                                                                                                                                                                   >  2:110676/5‑250 (MQ=255)                  aaTCCCCAAGCAAAGATGCCACCCCGGGCAACGTTTACGCTCAGGCTGACTACAAGATTGCTTACGAGTAATCGCATTTGGCCGGGTTCGCCCGGCCACATCATTAACGGATTAATGATAAGTGGATCAGATGTATAAAAAATTAAAATTAACCACAATAAGCGAATTGATTAAAAATATTTATTGTTCATTATCCGTTATTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTg                                                                                                                                                                                                      >  1:103766/1‑250 (MQ=255)                             aaaGATGCCACCCCGGGCAACGTTTACGCTCAGGCTGACTACAAGATTGCTTACGAGTAATCGCATTTGGCCGGGTTCGCCCGGCCACATCATTAACGGATTAATGATAAGTGGATCAGATGTATAAAAAATTAAAATTAACCACAATAAGCGAATTGATTAAAAATATTTATTGTTCATTATCCGTTATTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAAt                                                                                                                                                                                           >  2:127115/1‑250 (MQ=255)                                            ggCAACGTTTACGCTCAGTCTGACTACAAGATTGCTTACGAGTAATCGCATTTGGCCGGGTTCGGACGGCCACATCATGAACGGATTAATGATAAGTGGATCAGATGTATAAAAAAATAAAATTAACCACAATAAGCGAATTGATTAAAAATATTTATTGTTCATTATCCGTTATTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGTTTTTAGTCTa                                                                                                                                                                           <  2:123956/251‑1 (MQ=255)                                                                  aCTACAAGATTGCTTACGAGTAATCGCATTTGGCCGGGTTCGCCCGGCCACATCATTAACGGATTAATGATAAGTGGATCAGATGTATAAAAAATTAAAATTAACCACAATAAGCGAATTGATTAAAAATATTTATTGTTCATTATCCGTTATTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGTTTTTAGACTAGTCATTGGAGAACAGAtgattg                                                                                                                                                     >  1:84698/1‑251 (MQ=255)                                                                           ttGCTTACGAGTAATCGCATTTGGCCGGGTTCGCCCGGCCACATCATTAACGGATTAATGATAAGTGGATCAGATGTATAAAAAATTAAAATTAACCACAATAAGCGAATTGATTAAAAATATTTATTGTTCATTATCCGTTATTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGTTTTTAGACTAGTCATTGGAGAACAGATGATTGATGTCtt                                                                                                                                              >  2:53351/1‑249 (MQ=255)                                                                                 aCGAGTAATCGCATTTGGCCGGGTTCGCCCGGCCACATCATTAACGGATTAATGATAAGTGGATCAGATGTATAAAAAATTAAAATTAACCACAATAAGCGAATTGATTAAAAATATTTATTGTTCATTATCCGTTATTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGTTTTTAGACTAGTCATTGGAGAACAGATGATTGATGTCTTAGGGCCgg                                                                                                                                      <  1:179340/251‑1 (MQ=255)                                                                                                              cGGCCACATCATTAACGGCTTAATGATACGTGGATCAGATGTATAAAAACTTAAAATTAACCACAATAAGCGAATTGATTAAAAATATTTATTGTTCATTATCCGTTATTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGTTTTTAGACTAGTCATTGGAGAACAGATGATTGATGTCTTAGGGCCGGAGAAACGCAGACGGCGTACCACACAGGa                                                                                                          >  2:109858/1‑250 (MQ=255)                                                                                                                  cacctaattcaCGGATTACTTATAAGAGGGTCAGCTGTATAAAGAATTAAACTTAACAACAATAAGCGAATTGTTTGGCAATATTTATTGTTCATTATCCGCTATTAAACTGGCACACTGACTCTAAAGACAACAAATATCACTTAAATAAGTGATAGTCTTAATGCTAGTTTTTAGACTAGTAATTGGAGCACAGATGATTGATGTCTTAGGGCCGAAGAAACGCAGAGGGCGTAACACAAAGGaaa                                                                                                        >  2:217723/11‑248 (MQ=12)                                                                                                                     atcatTAACGGATTAATGATAAGTGGATCAGATGTATAAAAAATTAAAATTAACCACAATAAGCGAATTGATTAAAAATATTTATTGTTCATTATCCGTTATTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGTTTTTAGACTAGTCATTGGAGAACAGATGATTGATGTCTTAGGGCCGGAGAAACGCAGACGGCGTACCACACAGGaaa                                                                                                        >  2:7756/1‑245 (MQ=255)                                                                                                                           aaCGGATTAATGATAAGTGGATCAGATGTATAAAAAATTAAAATTAACCACAATAAGCGAATTGATTAAAAATATTTATTGTTCATTATCCGTTATTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGTTTTTAGACTAGTCATTGGAGAACAGATGATTGATGTCTTAGGGCCGGAGAAACGCAGACGGCGTACCACACAGGAAAAGATCGCAATTg                                                                                            >  2:40405/1‑251 (MQ=255)                                                                                                                                          aGTGGATCAGATGTATAAAAAATTAAAATTAACCACAATAAGCGAATTGATTAAAAATATTTATTGTTCATTATCCGTTATTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGTTTTTAGACTAGTCATTGGAGAACAGATGATTGATGTCTTAGGGCCGGAGAAACGCAGACGGCGGGCCACACAGGAAAAGATCGCAATTGTTCAGCAGAGCTTTg                                                                             <  1:22462/251‑1 (MQ=255)                                                                                                                                                           aaaaaTTAAAATTAACCACAATAAGCGAATTGATTAAAAATATTTATTGTTCATTATCCGTTATTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGTTTTTAGACTAGTCATTGGAGAACAGATGATTGATGTCTTAGGGCCGGAGAAACGCAGACGGCGTACCACACAGGAAAAGATCGCAATTGTTCAGCAGAGCTTTGAACCGGGGATGACGGtc                                                            >  1:106134/1‑251 (MQ=255)                                                                                                                                                               attaaaattaaCCACAATAAGCGAATTGATTAAAAATATTTATTGTTCATTATCCGTTATTAGACTGGCCCCCTGAATCTCCATACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGTTTTTAGACTAGTCATTGGAGAACAGATGATTGATGTCTTAGGGCTGGAGAAACGCAGACGGCGTACCACACAGGAAAAGATCGCAATTGTTCAGCAGAGCTTTGAACCGGGGATGACGGTCTccc                                                        >  1:207717/1‑251 (MQ=255)                                                                                                                                                                  aaaattaaCCACAATAAGCGAATTGATTAAAAATATTTATTGTTCATTATCCGTTATTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGTTTTTAGACTAGTCATTGGAGAACAGATGATTGATGTCTTAGGGCCGGAGAAACGCAGACGGCGTACCACACAGGAAAAGATCGCAATTGTTCAGCAGAGCTTTGAACCGGGGATGACGGTCTCCCTCg                                                     >  2:115943/1‑251 (MQ=255)                                                                                                                                                                          ccACAATAAGCGAATTGATTAAAAATATTTATTGTTCATTATCCGTTATTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGTTTTTAGACTAGTCATTGGAGAACAGATGATTGATGTCTTAGGGCCGGAGAAACGCAGACGGCGTACCACACAGGAAAAGATCGCAATTGTTCAGCAGAGCTTTGAACCGGGGATGACGGTCTCCCTCGTTGCCCgg                                             <  1:167195/251‑1 (MQ=255)                                                                                                                                                                           cacaATAAGCGAATTGATTAAAAATATTTATTGTTCATTATCCGTTATTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGTTTTTAGACTAGTCATTGGAGAACAGATGATTGATGTCTTAGGGCCGGAGAAACGCAGACGGCGTACCACACAGGAAAAGATCGCAATTGTTCAGCAGAGCTTTGAACCGGGGATGACGGTCTCCCTCGTTGCCCGGc                                            >  2:228649/1‑251 (MQ=255)                                                                                                                                                                                             tAAAAATATTTATTGTTCATTACCCGTTATTAGACTGGCACCCTAAATCTACAGACAACTAATATCACTTAAATAAGTGAAAGTCTTAATACTAGTTTTTAGACTAGACATTGGAGAACAGATGATTGATGCATTAGGGCCGGAGCAAAGCAGACGGCGTACCACACAGGAAAAGATCGAAATAGTTCAGCAGAGCTTTGAACCGGGGATGACGGACTCCTTAGTTGCCCGGCAACATGGTGTatcggcc                           >  2:57003/1‑244 (MQ=14)                                                                                                                                                                                               aaaaTATTTATTGTTCATTATCCGTTATTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGTTTTTAGACTAGTCATTGGAGAACAGATGATTGATGTCTTAGGGCCGGAGAAACGCAGACGGCGTACCACACAGGAAAAGATCGCAATTGTTCAGCAGAGCTTTGAACCGGGGATGACGGTCTCCCTCGTTGCCCGGCAACATGGTGTAGcagccagc                        >  1:199210/1‑251 (MQ=255)                                                                                                                                                                                                          tGTTCATTATCCGTTATTAGACTGGCCCCCTGCATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGTTTTTAGACTAGTCATTGGAGAACAGATGATTGATGTCTTAGGGCCGGAGAAACGCAGACGGCGTACCACACAGGAAAAGATCGCAATTGTTCAGCAGAGCTTTGAACCGGGGATGACGGTCTCCCTCGTTGCCCGGCAACATGGTGTAGCAGCCAGCCAGTTATTtct             >  2:135973/1‑251 (MQ=37)                                                                                                                                                                                                          tGTTCATTATCCGTTATTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGTTTTTAGACTAGTCATTGGAGAACAGATGATTGATGTCTTAGGGCCGGAGAAACGCAGACGGCGTACCACACAGGAAAAGATCGCAATTGTTCAGCAGAGCTTTGAACCGGGGATGACGGTCTCCCTCGTTGCCCGGCAACATGGTGTAGCAGCCAGCCAGTTATTtc              <  1:125809/250‑1 (MQ=37)                                                                                                                                                                                                                ttATCCGTTATTAGACTGGCCCACTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGTTTTTAGACTAGTCATTGGAGAACAGATGATTGATGTCTTAGGGCCGGAGAAACGCAGACGGCGTACCACACAGGAAAAGATCGCAATTGTTCAGCAGAGCTTTGAACCGGGGATGACGGTCTCCCTCGTTGCCCGGCAACATGGTGTAGCAGCCAGCCAGTTATTTCTCTGGCg       >  2:31597/1‑251 (MQ=18)                                                                                                                                                                                                                     cGTTATTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGTTTTTAGACTAGTCATTGGAGAACAGATGATTGATGTCTTAGGGCCGGAGAAACGCAGACGGCGTACCACACAGGAAAAGATCGCAATTGTTCAGCAGAGCTTTGAACCGGGGATGACGGTCTCCCTCGTTGCCCGGCAACATGGTGTAGCAGCCAGCCAGTTATTTCTCTGGCGTAAGc  <  1:199500/251‑1 (MQ=31)                                                                                                                                                                                                                                                           |                                                                                                                                                                                                                     TTTCAAAGCTCGTATGGAATCCCCAAGCAAAGATGCCACCCCGGGCAACGTTTACGCTCAGGCTGACTACAAGATTGCTTACGAGTAATCGCATTTGGCCGGGTTCGCCCGGCCACATCATTAACGGATTAATGATAAGTGGATCAGATGTATAAAAAATTAAAATTAACCACAATAAGCGAATTGATTAAAAATATTTATTGTTCATTATCCGTTATTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGTTTTTAGACTAGTCATTGGAGAACAGATGATTGATGTCTTAGGGCCGGAGAAACGCAGACGGCGTACCACACAGGAAAAGATCGCAATTGTTCAGCAGAGCTTTGAACCGGGGATGACGGTCTCCCTCGTTGCCCGGCAACATGGTGTAGCAGCCAGCCAGTTATTTCTCTGGCGTAAGC  >  NC_000913/3185878‑3186340

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