Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
JC JC	NC_000913	2,913,392	IS2 (–) +5 bp	coding (256‑260/2235 nt)	relA ←	GDP/GTP pyrophosphokinase

New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913	2913392 =	1 (0.050)	28 (1.450)	20/278	0.0	96.4%	coding (260/2235 nt)	relA	GDP/GTP pyrophosphokinase
	?	NC_000913	= 2997691	NA (NA)					noncoding (1/1331 nt)	IS2	repeat region
*	?	NC_000913	= 2913396	1 (0.050)	28 (1.450)	20/278	0.0	96.4%	coding (256/2235 nt)	relA	GDP/GTP pyrophosphokinase
	?	NC_000913	2996361 =	NA (NA)					noncoding (1331/1331 nt)	IS2	repeat region

GTTGCCTTATATAGGGTAACAACGTTGAAAAACCACATGTGATGATATTGTGATAGTTAGCGATGCCGGGTTTGTATCGTGCCTTTTCATTTATTTCCAGCGGTATAGCTGATAAATAATCCCTTAGTACATCAGGAACCTTATAAAAAAACCCGCGTTCTGCAAACCTTTCTCGATAAAAGTGACAACCACTGAGTAAACAGATCAATGGAAGTAAATTCCAATTGTTTATTGGATTTGCCCCTATATTTCCAGACATCTGTTATCACTTAACCCATTACAAGCCCGCTGCCGCAGATATTCCCGTGGCGAGCGATAACCCAGCGCACTATGCGGATGCCATTCGTTATAATGCTCGAACGCCTCTGCAAGGTTCTTTGCTGCCGTTAACCCGTCTGGTTTGGGCATGATACTGATGTAGTCACGCTTTATCGTTTTCACGAAGCTCTCTGCTATTCCGTTACTCTC  >  NC_000913/2996129‑2996596                                                                                                                                                                                                                                         |                                                                                                                                                                                                                                            gttgCCTTATATAGGGTAACAACGTTGAAAAACCACATGTGATGATATTGTGATCGTTAGCGCTCCCGGGTTTGTATCGTGACTTTTCATTTATTTCCAGCGGTCTAGCTGATAAATAATCCCTTAGTACATCAGGAACCTTATAAAAAAACCCGCGTTCTGCAAACCTTTCTCGATAAAAGTGACAACCACTGAGTACACAGATCAATGGAAGTAACTTCCAATTGTTTATTGGATTTGCCCCTATAttt                                                                                                                                                                                                                           >  1:477675/1‑251 (MQ=255)     ccTTATATAGGGTAACAACGTTGAAAAACAACATGTGATGATATTGTGATAGTTAGCGATGCCGGGTTTGTATCGTGCCTTTTCATTTATTTCCAGCGGTATAGCTGATAAATAATCCCTTAGTACATCAGGAACCTTATAAAAAAACCCGCGTTCTGCAAACCTTTCTCGATAAAAGTGACAACCACTGAGTAAACAGATCAATGGAAGTAAATTCCAATTGTTTATTGGATTTGCCCCTATATTTcc                                                                                                                                                                                                                         >  1:159759/1‑249 (MQ=255)                                                                                  ccTTTTCATTTATTTCCAGCGGTATAGCTGATAAATAATCCCTTAGTACATCAGGAACCTTATAAAAAAACCCGCGTTCTGCAAACCTTTCTCGATAAAAGTGACAACCACTGAGTAAACAGATCAATGGAAGTAAATTCCAATTGTTTATTGGATTTGCCCCTATATTTCCAGACATCTGTTATCACTTAACCCATTACAAGCCCGCTGCCGCAGATATTCCCGTGGCGAGCGATAACCCAGCGCACTAt                                                                                                                                          >  1:22794/1‑251 (MQ=255)                                                                                                 ccAGCGGTATAGCTGATAAATAATCCCTTAGTACATCAGGAACCTTATAAAAAAACCCGCGTTCTGCAAACCTTTCTCGATAAAAGTGACAACCACTGAGTAAACAGATCAATGGAAGTAAATTCCAATTGTTTATTGGAGTTGCCCCTATATTTCCAGACATCTGTTATCACTTAACCCATTACAAGCCCGCTGCCGCAGATATTCCCGTGGCGAGCGATAACCCAGCGCACTATGCGGATGCCATTCg                                                                                                                            <  1:264306/250‑1 (MQ=255)                                                                                                           aGCTGATAAATAATCCCTTAGTACATCATGAACCTTATAAAAAAACCCGCGTTCTGCAAACCTTTCTCGATAAAAGTGACAACCACTTAGTAAACAGATCAATGGAAGTAAATTCCAATTGTTTATTGGATTTGcc                                                                                                                                                                                                                                    <  1:273788/136‑1 (MQ=255)                                                                                                             cTGATAAATAATCCCTTAGTACATCAGGAACCTTATAAAAAAACCCGCGTTCTGCAAACCTTTCTCGATAAAAGTGACAACCACTGAGTAAACAGATCAATGGAAGTAAATTCCAATTGTTTATTTGATTTGCCCTTATATTTCCAGACATCTGTTATCACGTAACCCATTACAAGCCCTCTGCCGCAGATATTCCCGTGGCGAGCGATAACCCAgc                                                                                                                                                 <  1:274999/217‑1 (MQ=255)                                                                                                                     taatCCCTTAGTACATCAGGACCCTTATAAAAAAACCCGCGTTCTGCAAACCTTTCTCGATAAAAGTGAAAACCACTGAGTAAACAGATCAATGGAAGTAAATTCCAATTGTTTATTGGATTTGCCCCTATATTTCCAGACATCTGTTATCACTTAACCCATTACAAGCCCGCTGCCGCAGATATTcc                                                                                                                                                                      >  1:154456/1‑188 (MQ=255)                                                                                                                                               ataAAAAAACCCGCGTTCTGCAAACCTTTCTCGATAAAAGTGACAACCACTGAGTAAACAGATCAATGGAAGTAAATTCCAATTGTTTATTGGATTTGCCCCTATATTTCCAGACATCTGTTATCACTTAACCCATTACAAGCCCGCTGCCGCAGATATTCCCGTGGCGAGCGATAACCCAGCGCACTATGCGGATGCCATTCGTTATAATGCTCg                                                                                                                >  1:187643/1‑216 (MQ=255)                                                                                                                                                        nccGCGTTCTGCAAACCTTTCTCGATAAAAGTGACAACCACTGAGTAAACAGATCAATGGAAGTAAATTCCAATTGTTTCTTGGATTTGcc                                                                                                                                                                                                                                    >  1:133878/2‑91 (MQ=255)                                                                                                                                                            cgTTCTGCAAACCTTTCTCGATAAAAGTGACAACCACGGAGTAAACAGATAAATGGTCGTCAATTCCACTTGCTTATTGGACTTGCCCTTATCTTTCCCGACATCTGTTATCACTTAACCCATTACAAGCCCGCTGCCACAGCTATTCCCGTGGCGAGCGATAACCCAGCGCACTATGCGGATGCCATTCGTTATAATGCTCGAACGCCTCTGCAAGGTTCTTTGCTGCCGTTAACCCGTCTGGTTTGGGc                                                                >  1:252015/1‑251 (MQ=17)                                                                                                                                                                cTGCAAACCTTGCTCGATAAAAGGGACAACCACTGAGTAAACAGATCAATGGAATTAAATTCCAATTGTTTATTGGATTTGCCCCTATATTTCCAGACATCTGTTATCACTTAACCCATTACAAGCCCGCTGCCGCAGAGATTCCCGTGGCGAGCGATATCCCAGCGCACTATGCGGATTCCAGGCGTTATAATGCTCGAACGCCTCTGCAAGGTTCTTTGCTGCCGTTAACCCGTCTGGTTTg                                                                   <  1:334082/244‑1 (MQ=18)                                                                                                                                                                              cGATAAAANTNNNANNCACTGAGTAAACAGATCAATGGAAGTAAATTCCAATTGTTTATTNGATTTGCCCCTATATNTCNAGACNTCTNTTATCACTTAACCCATTACAAGCCCGCTGCCGCAGATATTCCCGTGGCGAGCGATAACCCAGCGCACTATGCGGATGCCATTCGTTATACTGCTCGAACGcc                                                                                                          >  2:358918/1‑191 (MQ=255)                                                                                                                                                                              cGATAAAAGTGACAACCACTGAGTAAACAGATCAATGGAAGTAAATTCCAATTGGTTATTGGATTTGCCCCTATATTTCCAGACATCTGTTATCACTTAACCCATTACAAGCCCGCAGCCGCAGATATTCCCTTGTCGAGCGATAACCCAGCGCACTATGCGGATGCCATTCGTTATAATGCTCGAACGcc                                                                                                          <  1:358918/191‑1 (MQ=255)                                                                                                                                                                                      ntGACAAACACTGAGTAAACAGATCAATGGAAGTAAATTCCAATTGTTTATTGGATTTGCCCCTATATTTCCAGACATCTGTTATCACTTAACCCATTACAAGCCCGCTGCCGCAGATATTCCCGTGGCGAGCGATAACCCAGCGCACTATGCGGATGCCATTCGTTATAATGCTCGAACGCCTCTGCAAGGTTCTTTGCTGCCGTTAAcc                                                                              >  1:133411/2‑211 (MQ=255)                                                                                                                                                                                          cAACCACTGAGTAAACAGATCAATGGAAGTAAATTCCAATTGTTTATTGGATTTGCCCCTATATTTCCAGACATCTGTTATCACTTAACCCATTACAAGCCCGCTGCCGCAGATATTCCCGTGGCGAGCGATAAcc                                                                                                                                                     >  1:119159/1‑136 (MQ=255)                                                                                                                                                                                          cAACCACTGAGTAAACAGATCAATGGAAGTAAATTCCAATTGTTTATNGGATTTGCCNNTANATTTCNANNNANNNGTNATCACTTAACCCATTACAANCCNNCTGCCGNAGATATTCCCGTGGCGAGCGATAAcc                                                                                                                                                     <  2:119159/136‑1 (MQ=21)                                                                                                                                                                                                      aaaCANATCAANGNNAGTAAATTCCCATTGTTTANTGGATTTGCCCCTATNTTTCNAGACANCTNt                                                                                                                                                                                                               >  2:218062/1‑66 (MQ=25)                                                                                                                                                                                                      aaaCAGATCAATGGAAGTAAATTCCAATTGTTTATTGGATTTGCCCCTATATTTCCAGACATCTGTTATCACTTAACCCATTACAAGCCCGCTGCCGCAGATATTCCCGTgg                                                                                                                                                                 <  1:218062/112‑1 (MQ=255)                                                                                                                                                                                                          aGATCAATGGAAGTAAATTCCAATTGTTTATTGGATTTGCCACTATATTTCCAGACATCTGTTATCACTTAACCCATTACAAGCCCGCTGCCGCAGATATTCCCGTGGCGAGCGATAACCAAGCGCACTATGCGGATGCCATTCGTTATAATGCTCGAACGCCTCTGCAAGGTTATTTGCTGCCGTTAACCCGCCTGGTTTGGGCATGATACTGATg                                                    >  1:334452/1‑217 (MQ=255)                                                                                                                                                                                                                              cAATTGTTTATTGGATTTGCCCCGAGATTTCCAGACATCTGTTATCACTTAACCCATTACAAGCCCGCTGCCGCAGATATTCCCGAGGCGAGCTATACCCCAGCGCACTAGGCGGATGTCAGTCGTTATAATGCTCGAACGCCTCTGCAAGTTTCTTTGCTGCTGTTAACCCGTCTGGTTTGTGCATGATGCTTATGTATTCACGCTTTATCGTTTTCACGaa                          <  1:153203/223‑1 (MQ=11)                                                                                                                                                                                                                                aTTGTTTATTGGATTTGCCCCTATATTTCCAGACATCTGTTATCACTTAACCCATTACAAGCCCGCTGCCGCAGATATTCCCGTGGCGAGCGATAACCCAGCGCACTATGCGGATGCCATTCGTTATAATGCTCGAACGCCTCTGCAAGGTTCTTtgctgcn                                                                                     <  1:129082/162‑2 (MQ=21)                                                                                                                                                                                                                                  tGTTTATTGGATTTGCCCCTATATTTCCAGACATCTGTTATCACTTAACCCATTACAAGCCCGCTGCCGCAGATATTCCCGTGGCGAGCGATAATCCAGCGCACTATGCGGATGCCATTCGTTATAATGCGCGAACGCCTCTGCAAGTTTCTTTGCTGCCGGTAACCCGTCTGGTTGGGGCATGATACTGATGTAGTCACGCTTTATCGTTTTCACGAAGCTCTCTGCTATTCCGTTACtctc  <  1:320744/243‑1 (MQ=255)                                                                                                                                                                                                                                      tATTGGATTTGCCCCTATATTTCCAGACATCTGTTATCACTTAACCCATTACAAGCCCGCTGCCGCAGATATTCCCGTGGCGAGCGATAACCCAGCGCACTATGCGGATGCCATTCGTTATAATGCTCGAACGCCTCTGCAAGGTTCTTTGCTGCCGTTAACCCGTCTGGTTTGGGCAtn                                                             <  1:150923/180‑2 (MQ=11)                                                                                                                                                                                                                                         |                                                                                                                                                                                                                                            GTTGCCTTATATAGGGTAACAACGTTGAAAAACCACATGTGATGATATTGTGATAGTTAGCGATGCCGGGTTTGTATCGTGCCTTTTCATTTATTTCCAGCGGTATAGCTGATAAATAATCCCTTAGTACATCAGGAACCTTATAAAAAAACCCGCGTTCTGCAAACCTTTCTCGATAAAAGTGACAACCACTGAGTAAACAGATCAATGGAAGTAAATTCCAATTGTTTATTGGATTTGCCCCTATATTTCCAGACATCTGTTATCACTTAACCCATTACAAGCCCGCTGCCGCAGATATTCCCGTGGCGAGCGATAACCCAGCGCACTATGCGGATGCCATTCGTTATAATGCTCGAACGCCTCTGCAAGGTTCTTTGCTGCCGTTAACCCGTCTGGTTTGGGCATGATACTGATGTAGTCACGCTTTATCGTTTTCACGAAGCTCTCTGCTATTCCGTTACTCTC  >  NC_000913/2996129‑2996596

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.