Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
JC JC	NC_000913	4,124,498	IS2 (+) +5 bp	intergenic (‑42/+110)	cytR ← / ← priA	DNA‑binding transcriptional repressor CytR/primosomal protein N'

New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913	4124498 =	2 (0.040)	48 (1.010)	46/486	0.1	95.9%	intergenic (‑42/+114)	cytR/priA	DNA‑binding transcriptional repressor CytR/primosomal protein N'
	?	NC_000913	= 4499511	NA (NA)					noncoding (1331/1331 nt)	IS2	repeat region
*	?	NC_000913	= 4124502	2 (0.040)	41 (0.870)	37/484	0.3	95.2%	intergenic (‑46/+110)	cytR/priA	DNA‑binding transcriptional repressor CytR/primosomal protein N'
	?	NC_000913	4498181 =	NA (NA)					noncoding (1/1331 nt)	IS2	repeat region

TATTTTCTTTGTTAAAAAAGACAAACGGTATTAACTGATGTATTTACTATTTACCGCTCCCTGCAGGGGGGATTATCAGAATCCGCGGTTATTCATTATTTATATACATTCAATAAAAAAGTAACCTCATAAAAAATCCTGTACTCACAGGGTGATGCGGGATAACTGACGTAAGGAGTGCTGCAGGGATTTGATGGTCCCTTGTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGTTTTTAGACTAGTCATTGGAGAACAGATGATTGATGTCTTAGGGCCGGAGAAACGCAGACGGCGTACCACACAGGAAAAGATCGCAATTGTTCAGCAGAGCTTTGAACCGGGGATGACGGTCTCCCTCGTTGCCCGGCAACATGGTGTAGCAGCCAGCCAGTTATTTCTCTGGC  >  NC_000913/4497977‑4498419                                                                                                                                                                                                             |                                                                                                                                                                                                                                               tattTTCTTTGTTAAAAAAGACAAACTGTATTAACTGATGTATTTACTATTTACCGCTCCCTGCAGGGGGGATTATCAGAATCCGCGGTTATTCATTATTTATATACATTCAATAAAAAAGTAACCTCATAAAAAATCCTGTACTCACAGGGTGATGCGGGATAACTGACGTAAGGAGTGCTGCAGGGATTTGATGGTCCCTTGTgctgatatgaatacggatattttatacagaatacttttcagtaca                                                                                                                                                                                                   <  2:348040‑M1/250‑46 (MQ=255)      tCTTTGTTAAAAAAGACAAACGGTATTAACTGATGTATTTACTATTTACCGCTCACTGCAGGGGGGATTATCAGAATCCGCGGTTATTCATTATTTATATACATTCAATAAAAAAGTAACCTCATAAAAAATCCTGTTCTCACAGGGTGATGCGGGATAACTGACGTAAGGAGTGCTGCAGTGATTTGATTGTCCCTTGTgctgatatgaatacggatattttatacagaatacttttcagtacaagcga                                                                                                                                                                                              <  1:499831‑M1/250‑51 (MQ=255)         ttGTTAAAAAAGACAAACGGTATTAACTGATGTATTTACTATTTACCGCTCCCTGCAGGGGGGATTATCAGAATCCGCGGTTATTCATTATTTATATACATTCAATAAAAAAGTAACCTCATAAAAAATCCTGTACTCACAGGGTGATGCGGGATAACTGACGTAAGGAGTGCTGCAGGGATTTGATGGTCCCTTGTgctgatatgaatacggatattttatacagaatacttttcagtacaagcgactgt                                                                                                                                                                                          >  1:83299‑M1/1‑197 (MQ=255)                                              actatttacCGCTCCCTGCAGGGGGGATTATCAGAATCCGCGGTTATTCATTATTTATATACATTCAATAAAAAAGTAACCTCATAAAAAATCCTGTACTCACAGGGTGATGCGGGATAACTGACGTAAGGAGTGCTGCAGGGATTTGATGGTCCCTTGTgctgatatgaatacggatattttatacagaatacttttcagtacaagcgactgtattacagtctgccatatatgattacctgtgatcagg                                                                                                                                                      <  2:459007‑M1/250‑91 (MQ=255)                                               ctatttacCGCTCACTGCAGGGGGGATTATCAGAATCCGCGGTTATTCATTATTTATATACATTCAATAAAAAAGTAACCTCATAAAAAATCCTGTACTCACAGGGTGATGCGGGATAACTGACGTAAGGAGTGCTGCAGGGATTTGATGGTCCCTTGTgctgatatgaatacggatattttatacagaatacttttcagtacaagcgactgtattacagtctgccatatatgattacctgtgatcagga                                                                                                                                                     >  1:155596‑M1/1‑159 (MQ=255)                                                    tacCGCTCCCTGCAGGGGGGATTATCAGAATCCGCGGTTATTCATTATTTATATACATTCAATAAAAAAGTAACCTCATAAAAAATCCTGTACTCACAGGGTGATGCGGGATAACTGACGTAAGGAGTGCTGCAGGGCTTTGATGGTCCCTTGTgctgatatgaatacggatattttatacagaatacttttcagtacaagcgactgtattacagtctgccatatatgattacctgtgatcaggaagaat                                                                                                                                                >  1:494462‑M1/1‑154 (MQ=255)                                                     acCGCTCCCTGCAGGGGGGATTATCAGAATCCGCGGTTATTCATTATTTATATACATTCAATAAAAAAGTAACCTCATAAAAAATCCTGTACTCACAGGGTGATGCGGGATAACTGACGTAAGGAGTGCTGCAGGGATTTGATGGTCCCTTGTgctgatatgaatacggatattttatacagaatacttttcagtacaagcgactgtattacagtctgccatatatgattacctgtgatcaggaagaata                                                                                                                                               >  1:258169‑M1/1‑153 (MQ=255)                                                     acCGCTCCCTGCAGGGGGGATTATAAGAATCCGCGGTTATTCATTATTTATATACATTCAATAAAAAAGTAACCTCATAAAAAATCCTGTACTCAAAGGGTGATGCGGGATAACTGACGTAAGGAGTGCTGCAGGGATTTGATGGTCCCTTGTgctgatatgaatacggatattttatacagaatacttttcagtacaagagcctgtatgacagtctgcaatatatgattacctgtgatcaggaaaaata                                                                                                                                               >  2:357391‑M1/1‑153 (MQ=255)                                                                 aGGGTGGATTATCAGAATCCGCGGTTATTCATTATTTATATACATTCAATAAAAAAGTAACCTCATAAAAAATCCTGTACTCACAGGGTGATGCGGGATAACTGACGTAAGGAGTGCGGCAGGGATTTGATGGTCCCTTGTgctgatatgaatacggatattttatacagaatacttttcagtacaagcgactgtattacagtctgccatatatgattacctgtgatcaggaagaataagtggctgaggtt                                                                                                                                  <  2:346‑M1/251‑111 (MQ=255)                                                                                 aTCCGCGGTTATTCATTATTTATATACATTCAATTAAAAAGTAACCTCATAAAAAATCCTGTACTCACAGGGTGATGCGGGATAAATTACGTAAGGAGTGCTGCCGGGATTTGATGGTCCCTTGTgctgatatgaatacgggtattttatacaacatacttttcagtacaaacgactgtattaaagtctgccatatatgattacctttgatcaggaagaataaggggctggggtttcaaactttccgccgg                                                                                                                  >  2:210388‑M1/1‑125 (MQ=255)                                                                                    cgcgGTTATTCATTATTTATATACATTCAATAAAAAAGTAACCTCATAAAAAATCCTGTACTCACAGGGTGATGCGGGATAACTGACGTAAGGAGTGCTGCAGGGATTTGATGGTCCCTTGTgctgatatgaatacggatattttatacagaatacttttcagtacaagcgactgtattacagtctgccatatatgattacctgtgatcaggaagaataagtggctgaggtttcaaactgtctgccggtat                                                                                                               >  1:406801‑M1/1‑122 (MQ=255)                                                                                        gTTATTCATTATTTATATACATTCAATAAAAAAGTAACCTCATAAAAAATCCTGTACTCACAGGGTGATGCGGGATAACTGACGTAAGGAGTGCTGCAGGGATTTGATGGTCCCTTGTgctgatatgaatacggatattttatacagaatacttttcagtacaagcgactgtattacagtctgccaaatatgattacctgtgatcaggaagaataagtggctgaggtttcaaactgtctgccggtatattc                                                                                                           >  2:21742‑M1/1‑118 (MQ=255)                                                                                                  atttatATACATTCAATAAAAAAGTAACCTCATAAAAAATCCTGTACTCACAGGGTGATGCGGGATAACTGACGTAAGGAGTGCTGCAGGGATTTGATGGTCCCTTGTgctgatatgaatacggatattttatacagaatacttttcagtacaagcgactgtattacagtctgccatatatgattacctgtgatcaggaagaataagtggctgaggtttcaaactgtctgccggtatattcctctctccct                                                                                                 >  2:366612‑M1/1‑108 (MQ=255)                                                                                                              tCAATAAAAAAGTAACCTCATAAAAAATCCTGTACTCACAGGGTGATGCGGGATAACTGACGTAAGGAGTGCTGCAGGGATTTGATGGTCCCTTGTgctgatatgaatacggatattttatacagaatacttttcagtacaagcgactgtattacagtctgccatatatgattacctgtgatcaggaagaataagtggctgaggtttcaaactgtctgccggtatattcctctctcccttaaaaaccata                                                                                      <  2:149726‑M1/250‑155 (MQ=255)                                                                                                                    aaaaaGTAACCTCATAAAAAATCCTGTACTCACAGGGTGATGCGGGATAACTGACGTAAGGAGTGCTGCAGGGATTTGATGGTCCCTTGTgctgatatgactacggatattttatacagaatacttttcagtacaagcgactgtattacagtctgccatatatgattacctgtgatcaggaagaataagtggctgaggtttcaaactgtctgccggtatattcctctctcccttaaaaaccataacaggtt                                                                               >  1:195338‑M1/1‑90 (MQ=255)                                                                                                                     aaaaGTAACCTCATAAAAAATCCTGTACTCACAGGGTGATGCGGGATAACTGACGTAAGGAGTGCTGCAGGGATTTGATGGTCCCTTGTgctgatatgaatacgggtattttatacagaatacttttcagtacaagcgactgtattacatactgccatatatgattacctgtgatcaggaagaataagtgactgaggtttcaaactgtctgccggtatattcctctctcccttaaaaaccataacaggtt                                                                               >  1:488457‑M1/1‑89 (MQ=255)                                                                                                                            accacataaAAAATCCTGTACTCACAGGGTGATGCGGGATAACTGACGTAAGGAGTGCTGCAGGGATTTGATGGTCCCTTGTgctgatatgaatacggatattttatacagaatacttttcagtacaagcgactgtattacagtctgccatatatgattacctgtgatcaggaagaataagtggctgaggtttcaaactgtctgccggtatattcctctctcccttaaaaaccataacaggttagttatct                                                                       <  2:494462‑M1/247‑170 (MQ=255)                                                                                                                                 aTAAAAAATCCTGTACTCACAGGGCGATGCGGGATAACTGACGTAAGGAGTGCTGCAGGGATTTGATGGTCCCTTGTgctgatatgaatacggatattttatacagaatacttttcagtacaagcgactgtattacagtctgccatatatgattacctgtgatcaggaagaataagtggctgaggtttcaaactgtctgccggtatattcctctctcacttaaaaaccataacaggttagttatcttcatc                                                                  >  1:165564‑M1/1‑77 (MQ=255)                                                                                                                                  tAAAAAATCCTGTACTCACAGGGAGATGCGGGATAACTGACGTAAGGAGTGCTGCAGGGATTTGATGGTCCCTTGTgctgatatgaatacggatattttatacagaatacttttcagtacaagcgactgtattacagtctgccatatatgattacctgtgatcaggaagaataagtggctgaggtttcaaactgtctgccggtatattcctctctcccttaaaaaccataacaggttagttatcttagtct                                                                 >  2:415714‑M1/1‑76 (MQ=255)                                                                                                                                   aaaaaaTCCTGTACTCACAGGGTGATGCGGGATAACTGACGTAAGGAGTGCTGCAGGGATTTGATGGTCCCTTGTgctgatatgaatacggatattttatacagaatacttttcagtacaagcgactgtattacagtctgccatatatgattacctgtgatcaggaagaataagtggctgaggtttcaaactgtctgccggtatattcctctctcccttaaaaaccataacaggttagttatcttcgtctg                                                                <  1:415714‑M1/251‑177 (MQ=255)                                                                                                                                                   acaGGGTGATGCGGGATAACTGACGTAAGGAGTGCTGCAGGGATTTGATGGTCCCTTGTgctgatatgaatacggatattttatacagaatacttttcagtacaagcgactgtattacagtctgccatatatgattacctgtgatcaggaagaataagtggctgaggtttcaaactgtctgccggtatattcctctctcccttaaaaaccataacaggttagttatcttcgtctgtcagcggattgaatgg                                                <  1:473295‑M1/251‑193 (MQ=255)                                                                                                                                                           aTGCGGGATAACTTACGTAAGGAGTGCTGCAGGTATTTGATGGTCCCTTGTgctgctatgaatacggatattttatacagaatacttttcagtacaagcgactgtattacagtctgccatatatgattacctgtgatcaggaagaataagtggctgaggtttcaaaatgtctgccggtatattcctctctcccttaaaaaccataacaggttagttatcttcgtctgtcagcggattgaatggcgga                                            >  2:178006‑M1/1‑51 (MQ=255)                                                                                                                                                               gggATAACTGACGTAAGGAGTGCTGCAGGGATTTGATGGTCCCTTGTgctgatatgaatacggatattttatacagaatacttttcagtacaagcgactgtattacagtctgccatatatgattacctgtgatcaggaagaataagtggctgaggtttcaaactgtctgccggtatattcctctctcccttaaaaaccataacaggttagttatcttcgtctgtcagcggattgaatggcggtatattttc                                    >  1:277211‑M1/1‑47 (MQ=255)                                                                                                                                                                         aCGTAAGGAGTGCTGCAGGGATTTGATGGTCCCTTGTgctgatatgaatacggatattttatacagaatacttttcagtacaagcgactgtattacagtctgccatatatgattacctgtgatcaggaagaataagtggctgaggtttcaaactgtctgccggtatattcctctctcccttaaaaaccataacaggttagttatcttcgtctgtcagcggattgaatggcggtatattttcagtacaaata                          <  1:210388‑M1/251‑215 (MQ=255)                                                                                                                                                                             aaGGAGTGCTGCAGGGATTTGATGGTCCCTTGTgctgatatgaatacggatattttatacagaatacttttcagtacaagcgactgtattacagtctgccatatatgattacctgtgatcaggaagagtaagtggctgaggtttcaaactgtctgccggtatattcctctctcccttaaaaaccataacaggttagttatcttcgtctgtcagcggattgaatggcggtatattttcagtacaaataccgg                      <  1:158464‑M1/251‑219 (MQ=255)                                                                                                                                                                                                 gATGGTCCCTTGTgctgatatgaatacggatattttatacagaatacttttcagtacaagcgactgtattacagtctgccaaatatgattacctgtgatcaggaagaataagtggctgaggtttcaaactgtctgccggtatattcctctctcccttaaaaaccataacaggttagttatcttcgtctgtcagcggattgaatggcggtatattttcagtacaaataccggtcatccactgaaaaatacc   <  1:21742‑M1/250‑238 (MQ=255)                                                                                                                                                                                                 gATGGTCCCTTGTgctgatatgaatacggatatgttattcagaatacttttcagtacaagcgactttattacagtctgccatatatgattacctgtgatcaggaataataagtggctgaggtttcaaactgtctgccggtatattcctctctcccttaaaaaccataacaggttagttatcttcgtctgtcagcggattgaatggcggtatattttcagtacaaataccggtcatccactgaaaaatacct  <  2:112832‑M1/251‑239 (MQ=255)                                                                                                                                                                                                             |                                                                                                                                                                                                                                               TATTTTCTTTGTTAAAAAAGACAAACGGTATTAACTGATGTATTTACTATTTACCGCTCCCTGCAGGGGGGATTATCAGAATCCGCGGTTATTCATTATTTATATACATTCAATAAAAAAGTAACCTCATAAAAAATCCTGTACTCACAGGGTGATGCGGGATAACTGACGTAAGGAGTGCTGCAGGGATTTGATGGTCCCTTGTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGTTTTTAGACTAGTCATTGGAGAACAGATGATTGATGTCTTAGGGCCGGAGAAACGCAGACGGCGTACCACACAGGAAAAGATCGCAATTGTTCAGCAGAGCTTTGAACCGGGGATGACGGTCTCCCTCGTTGCCCGGCAACATGGTGTAGCAGCCAGCCAGTTATTTCTCTGGC  >  NC_000913/4497977‑4498419

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.