Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
JC JC	NC_000913	3,560,081	IS1 (+) +9 bp	intergenic (+215/+533)	rtcR → / ← glpG	DNA‑binding transcriptional activator RtcR/rhomboid protease GlpG

New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913	3560081 =	1 (0.040)	53 (1.830)	46/484	0.0	98.1%	intergenic (+215/+541)	rtcR/glpG	DNA‑binding transcriptional activator RtcR/rhomboid protease GlpG
	?	NC_000913	= 3584195	NA (NA)					noncoding (768/768 nt)	IS1	repeat region
*	?	NC_000913	= 3560089	1 (0.040)	30 (1.040)	27/484	0.2	96.7%	intergenic (+223/+533)	rtcR/glpG	DNA‑binding transcriptional activator RtcR/rhomboid protease GlpG
	?	NC_000913	3583428 =	NA (NA)					noncoding (1/768 nt)	IS1	repeat region

CCCTTTGACGTGGTGATATGGATGACGGATGGCTGGCCGCTGTATGAATCCCGCCTGAAGGGAAAGCTGCACGTAATCAGCAAGCGATATACGCAGCGAATTGAGCGGCATAACCTGAA‑TCTGAGGCAGCACCTGGCACGGCTGGGACGGAAGTCGCTGTCGTTCTCAAAATCGGTGGAGCTGCATGACAAAGTCATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCGACATTAAAAGTAACTTTTTTGAAATTAGAGATAATATTGGCACAGAAAATATATTCAGTTCGTTTAGTGATTTCTCTTCATTCTTAAACGAAATCATGGACAGTTGTTCATAATAATCAAAAGGAAACAAAAATCATTGTAAACCATGCCATCTTTTTAGAGATTACCAGCAGAAAATTCCTTGTTCGAAATTTAGCAGTGCAAGCTATGGAAAAACTTACAACGGAACTACACTCCTTA  >  NC_000913/3583950‑3584436                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                  cccTTAGACGTGGTGATATGGATGACGGATGGCTGGCCGCTGTATGAATCCCGCCTGAAGGGAAAGCGGCACGTAATCAGCAAGCGATATACGCAGCGAATTGAGCGGCATAACCTGAA‑TCTGAGGCAGCACCTGGCACGGCTGGGACGGAAGTCGCTGTCGTTCTCAAAATCGGTGGAGCTGCATGACAAAGTCATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCGACa                                                                                                                                                                                                                                               <  1:331159/250‑1 (MQ=14)        aCGTGGTGATATGGATGGCGGATGGCTGGCCGCTGTATGAATCCCGCCTGAAGGGAAAGCTGCACGTAATCAGCAAGCGATATACGCAGCGAATTGAGCGGCATAACCTGAA‑TCTGAGGCAGCACCTGGCACGGCTGGGACGGAAGTCGCTGTCGTTCTCAAAATCGGTGGAGCTGCATGACAAAGTCATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCGACATTAAAAGt                                                                                                                                                                                                                                       <  1:159445/251‑1 (MQ=21)            ggtgATATGGATGACGGATGGCTGGCCGCTGTATGAATCCCGCCTGAAGGGAAAGCTGCACGTAATCAGCAAGCGATATACGCAGCGAATTGAGCGGCATAACCTGAA‑TCTGAGGCAGCACCTGGCACGGCTGGGACGGAAGTCGCTGTCGTTCTCAAAATCGGTGGAGCTGCATGACAAAGTCATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATGACCGACATTAAAAGTAAc                                                                                                                                                                                                                                    <  1:98077/250‑1 (MQ=37)               gATATGGATGACGGATGGCTGGCCGCTGTATGAATCCCGCCTGAAGGGAAAGCTGCACGTAATCAGCAAGCGATATACGCAGCGAATTGAGCGGCATAACCTGAA‑TCTGAGGCAGCACCTGGCACGGCTGGGACGGAAGTCGCTGTTGTTCTCAAAATCGGTGGAGCTGCATGACAAAGTCATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCGAc                                                                                                                                                                                                                                                >  1:81136/1‑235 (MQ=14)               gATATGGATGACGGATGGCTGGCCGCTGTATGAATCCCGCCTGAAGGGAAAGCTGCACGTAATCAGCAAGCGATATACGCAGCGAATTGAGCGGCATAACCTGAA‑TCTGAGGCAGCACCTGGCACGGCTGGGACGGAAGTCGCTGTCGTTCTCAAAATCGGTGGAGCTGCATGACAAAGGCATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCGACATTAAAAGTACCtttt                                                                                                                                                                                                                                >  2:273067/1‑251 (MQ=25)               gATAGGGATGACGGAAGGCTGGCCGCTGTATGAATCCCGCCTGAAGGGAAAGCTGCACGTAATCAGCAAGCGATATACGCAGCGAATTGAGCGGCATAACCTGAA‑TCTGAGGCAGCACCTGGCACGGCTGGGACGGAAGTCGCTGTTGTTCTCAAAATCGGTGGAGCTGCATGACAAAGTCATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCGAc                                                                                                                                                                                                                                                <  2:81136/235‑1 (MQ=9)                  agggagtaCGGAGGGCTGGCCGCTGTATGAATCCCGCCTGAAGGGAAAGCTGCACGTAATCAGCAAGCGATATACGCAGCGAATTGAGCGGCATAACCTGAA‑TCTGAGGCAGCACCTGGCACGGCTGGGACGGAAGTCGCTGTCGTTCTCAAAATCGGTGGAGCTGCATGACAAAGTCATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCGACATTAAAAGTAACTTTTTTg                                                                                                                                                                                                                             <  2:121137/244‑1 (MQ=255)                            gagggctggcAGCTGTAGGAATCCCGCCTGAAGGGAAAGCTGCACGTAATCAGCAAGCGATATACGCAGCGAATTGAGAGGCATAACCTGAA‑TCTGAGGCAGCACCTGGCACGGCTGGGACGGAAGTCGCTGTCGTTCTCAAAATCGGTGGAGCTGCAGGACAAAGTCATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCGACATTAAAAGTAACTTTTTTGAAATTAgaga                                                                                                                                                                                                                   <  2:123145/248‑1 (MQ=255)                                   ggccgctttagaaaTCCCGCCTGAAGGGAAAGCTGTACGTAATCAGCAAGCGATATACGCAGCGAATTGAGCGGCATAACCTGAA‑TCTGAGGCAGCACCTGGCACGGCTGGGACGGAAGTCGCTGTCGTTCTCAAAATCGGTGGAGCTGCATGACAAAGTCATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCGACATTAAAAGTAACTTTTTTGAAATTAGAGATAATAtt                                                                                                                                                                                                            <  2:315155/239‑1 (MQ=255)                                                       cTGAAGGGAAAGCTGCACGTAATCAGCAAGCGATATACGCAGCGAATTGAGCGGCATAACCTGAA‑TCTGAGGCAGCACCTGGCACGGCTGGGACGGAAGTCGCTGTCGTTCTCAAAATCGGTGGAGCTGCATGACAAAGTCATCGGGCATTATATGAACATAAAACACTATCAATAAGTTGGAGTCATTACCGACATTAAAAGTAACTTTTTTGAAATTAGAGATAATATTGGCACAGAAAATATATTc                                                                                                                                                                                          >  1:165402/1‑249 (MQ=255)                                                                  ccTGCACGTAATCAGCAAGCGATATACGCCGCGAATTGAGCGGCATAACCAGAA‑TCTGAGGCAGCACCTGGCACGGCTGGGACGGAAGTCGCTGTCGTTCTCAAAATCGGTGGAGCTGCATGCCAAAGTCATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCGACATTAAAAGTACCTTTTTTGAAATTAGAGATAATATTGGCACAGAAAATATATTAAGTTCGTTTAGTg                                                                                                                                                                             >  2:42365/2‑251 (MQ=255)                                                                         gTAATCAGCAAGCGATATACGCAGCGAATTGAGCGGCATAACCTGAA‑TCTGAGGCAGCACCTGGCACGGCTGGGACGGAAGTCGCTGTCGTTCTCAAAATCGGTGGAGCTGCATGACAAAGTCATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCGACATTAAAAGTAACTTTTTTGAAATTAGAGATAATATTGGCACAGAAAATATATTCAGTTCGTTTAGGGATTtctc                                                                                                                                                                      >  2:5983/1‑251 (MQ=255)                                                                              cAGCAAGCGATATACGCAGCGAATTGAGCGGCATAACCTGAA‑TCTGAGGCAGCACCTGGCACGGCTGGGACGGAAGTCGCTGTCGTTCTCAAAATCGGTGGAGCTGCATGACAAAGACATCGGGCATTATCTGAACATAAAACACTATCAATAAGGTGGAGTCATTACCGACATTAAAAGTAACTTTTTTGTAATTAGAGATACTATTGGCACAGAAAATATATTCAGTTCGTTTAGGGATTTCTCTTCAt                                                                                                                                                                 >  2:272996/1‑251 (MQ=255)                                                                                         ataCGCAGCGAATTGAGCGGCATAACCTGAA‑TCTGAGGCAGCACCTGGCACGGCTGGGACGGAAGTCGCTGTCGTTCTCAAAATCGGTGGAGCTGCATGACAAAGTCATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCGACATTAAAAGTAACTTTTTTGAAATTAGAGATAATATTGGCACAGAAAATATATTCAGTTCGTTTAGTGATTTCTCTTCATTCTTAAACGa                                                                                                                                                       >  1:10710/1‑250 (MQ=255)                                                                                                 cGAATTGAACGGCAAAACATGTA‑TCTGAGGCAGCTCCTGGCACGGTTGGGACGGAATACGCGGTCGTTCTAAAAATCGGTTGATCTGCATGACAAAGTCATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCGACATTAAAAGTAACTGTTTTGAAATTAGAGATAATATTGGCACAGAAAATATATTCAGTTCGTTTAGTGATTTCTCTTCATTCTTAAACGAAATCATGGa                                                                                                                                              <  2:181885/251‑1 (MQ=21)                                                                                                   aaTTGAGCGGCATAACCTGAA‑TCTGAGGCAGCACCTGGCACGGCTGGGACGGAAGTCGCTGTCGTTCTCAAAATCGGTGGAGCTGCATGACAAAGTTATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCGACATTAAAAGTAACTTTTTTGAAATTAGAGATAATATTGGCACAGAAAATATATTCAGTTCGTTTAGTGATTTCTCTTCATTCTTAAACGAAATCATGGACa                                                                                                                                            <  1:241699/251‑1 (MQ=255)                                                                                                         gcgggctAACCTGAAGTC‑GAGGCAGCACCCGGCACGTCTGGGACGGAAGTCGTTGTCGATCACAAATTCGGTGGAGCAGCATGACAAAGTCATCTGGAATTATCTGAAAAAAAAACACTATCAATAAGTTGGAGTCATTACCGACATTAAAAGTAACTTTTTTGAAATTAGAGATAATATTGGCACAGCAAATATATTCAGTTCGTTTAGTGATTTCTCTTCATTCTTAAACGAAATCATGGCCAgttgtt                                                                                                                                      <  2:110424/245‑1 (MQ=21)                                                                                                                                        ggCACGGCTGGGACGGAAGTCGCTGTCGTTCTCAAAATCGGTGGAGCTGCATGACAAAGTCATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCGACATTAAAAGTAACTTTTTTGAAATTAGAGATAATATTGGCACAGAAAATATATTCAGTTCGTTTAGTGATTTCTCTTCATTCTTAAACGAAATCATGGACAGTTGTTCATAATAATCAAAAGGAAACAAAAATCAt                                                                                                         <  1:5983/250‑1 (MQ=255)                                                                                                                                           ccGGCTGGGACGGAAGTCGCTGTCGTTCTCAAAATCGGTGGAGCTGCATGACAAAGTCATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCGACATTAAAAGTAACTTTTTTGAAATTAGAGATAATATTGGCACAGAAAATATATTCAGTTCGTTTAGTGATTTCTCTTCATTCTTAAACGAAATCATGGACAGTTGTTCATAATAATCAAAAGGAAACAAAAATCATTGTa                                                                                                     <  1:273067/250‑1 (MQ=255)                                                                                                                                                                                       tGCATGACAAAGTCATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCGACATTAAAAGTAACTTTTTTGAAATTAGAGATAATATTGGCACAGAAAATATATTCAGTTCGTTTAGTGATTTCTCTTCATTCTTAAACGAAATCATGGACAGTTGTTCATAATAATCAAAAGGAAACAAAAATCATTGTAAACCATGCCATCTTTTTAGAGATTACCAGCAGAAAATTCCTTGt                                                         >  1:87848/1‑251 (MQ=255)                                                                                                                                                                                                                                                tCATTACCGCCATTAAAAGTAACTTTTTTGAAATTAGAGATAATATTGGCACAGAAAATATATTCAGTTCGTTTAGTGATTTCTCTTCATTCTTAAACGAAATCATGGACAGTTGTTCATAATAATCAAAAGGAAACAAAAATAATTGTAAACCATGCCATCTTTTTAGAGATTACCAGCAGAAAATTCCTTGTTCGAAATTTAGCAGTGCAAGCTATGGAAAAACTTACAACGAAACTACACTCCTTa  >  2:105787/1‑249 (MQ=255)                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                  CCCTTTGACGTGGTGATATGGATGACGGATGGCTGGCCGCTGTATGAATCCCGCCTGAAGGGAAAGCTGCACGTAATCAGCAAGCGATATACGCAGCGAATTGAGCGGCATAACCTGAA‑TCTGAGGCAGCACCTGGCACGGCTGGGACGGAAGTCGCTGTCGTTCTCAAAATCGGTGGAGCTGCATGACAAAGTCATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCGACATTAAAAGTAACTTTTTTGAAATTAGAGATAATATTGGCACAGAAAATATATTCAGTTCGTTTAGTGATTTCTCTTCATTCTTAAACGAAATCATGGACAGTTGTTCATAATAATCAAAAGGAAACAAAAATCATTGTAAACCATGCCATCTTTTTAGAGATTACCAGCAGAAAATTCCTTGTTCGAAATTTAGCAGTGCAAGCTATGGAAAAACTTACAACGGAACTACACTCCTTA  >  NC_000913/3583950‑3584436

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.