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 =	0 (0.000)	40 (0.990)	37/482	0.2	100%	coding (260/2235 nt)	relA	GDP/GTP pyrophosphokinase
	?	NC_000913	= 2997691	NA (NA)					noncoding (1/1331 nt)	IS2	repeat region
*	?	NC_000913	= 2913396	0 (0.000)	53 (1.320)	47/482	0.1	100%	coding (256/2235 nt)	relA	GDP/GTP pyrophosphokinase
	?	NC_000913	2996361 =	NA (NA)					noncoding (1331/1331 nt)	IS2	repeat region

CATTGCCAATGGTTGTTGCCTTATATAGGGTAACAACGTTGAAAAACCACATGTGATGATATTGTGATAGTTAGCGATGCCGGGTTTGTATCGTGCCTTTTCATTTATTTCCAGCGGTATAGCTGATAAATAATCCCTTAGTACATCAGGAACCTTATAAAAAAACCCGCGTTCTGCAAACCTTTCTCGATAAAAGTGACAACCACTGAGTAAACAGATCAATGGAAGTAAATTCCAATTGTTTATTGGATTTGCCCCTATATTTCCAGACATCTGTTATCACTTAACCCATTACAAGCCCGCTGCCGCAGATATTCCCGTGGCGAGCGATAACCCAGCGCACTATGCGGATGCCATTCGTTATAATGCTCGAACGCCTCTGCAAGGTTCTTTGCTGCCGTTAACCCGTCTGGTTTGGGCATGATACTGATGTAGTCACGCTTTATCGTTTTCACGAAGCTCTCTGCTATTCCGTTACTCT  >  NC_000913/2996115‑2996595                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                           cATTGCCAATGGTTGTTGCCTTATATAGGGTAACAACGTTGAAAAACCACATGTGATGATATTGTGATAGTTAGCGATGCCGGGTTTGTATCGTGCCTTTTCATTTATTTCCAGCGGTATAGCTGATAAATAATCCCTTAGTACATCAGGAACCTTATAAAAAAACCCGCGTTCTGCAAACCTTTCTCGATAAAAGTGACAACCACTGAGTAAACAGATCAATGGAAGTAAATTCCAATTGTTTATTGGAt                                                                                                                                                                                                                                        <  2:17580/251‑1 (MQ=255)    tGCCAATGGTTGTTGCCTTATATAGGGTAACAACGTTGAAAAACCACATGTGATGATATTGTGATAGTTAGCGATGCCGGGTTTGTATCGTGCCTTTTCATTTATTTCCAGCGGTATAGCTGATAAATAATCCCTTAGTACATCAGGAACCTTATAAAAAAACCCGCGGTCTGCAAACCTTTCTCGATAAAGGTGACAACCACTGAGTAAACAGATCAATGGAAGTAAATTCCAATTGGTTATGGGagtt                                                                                                                                                                                                                                      <  2:52427/250‑4 (MQ=255)                 tgCCTTATATAGGGTAACAACGTTGAAAAACCACATGTGATGATATTGTGATAGTTAGCGATGCCGGGTTTGTATCGTGCCTTTTCATTTATTTCCAGCGGTATAGCTGATAAATAATCCCTTAGTACATCAGGAACCTTATAAAAAAACCCGCGTTCTGCAAACCTTTCTCGATAAAAGTGACAACCACTGAGTAAACAGATCAATGGAAGTAAATTCCAATTGTTTATTGGATTTGCCCCTATATTTcc                                                                                                                                                                                                                        >  1:313019/1‑251 (MQ=255)                   ccTTATAGAGGGTAACAACGTTGAAAAACCACATGTGATGATATTGTGATAGTTAGCGATGCCGGGTTTGTATCGTGCCTTTTCATTTATTTCCAGCGGTATAGCTGATAAATAATCCCTTAGTACATCAGGAACCTTATAAAAAAACCCGCGTTCTGCAAACCTTTCTCGATAAAAGTGACAACCACTGAGTAAACAGATCAATGGAAGTAAATTCCAATTGTTTATTGGATTTGCCCCTATATTTCCAg                                                                                                                                                                                                                      <  1:382384/251‑1 (MQ=255)                             ggTAACAACGTTGAAAAACCACATGTGATGATATTGTGATAGTTAGCGATGCCGGGTTTGTATCGTGCCTTTTCATTTATTTCCAGCGGTATAGCTGATAAATAATCCCTTAGTACATCAGGAACCTTATAAAAAAACCCGCGTTCTGCAAACCTTTCTCGATAAAAGTGACAACCACTGAGTAAACAGATCAATGGAAGTAAATTCCAATTGTTTATTGGATTTGCCCCTATATTTCCAGACATCTGtt                                                                                                                                                                                                             >  2:357713/1‑250 (MQ=255)                                aacaacGTTGAAAAACCACATGTGATGATATTGTGATAGTTAGCGATGCCGGGTTTGTATCGTGCCTTTTCATTTATTTCCAGCGGTATAGCTGATAAATAATCCCTTAGTACATCAGGAACCTTATAAAAAAACCCGCGTTCTGCAAACCTTTCTCGATAAAAGTGACAACCACTGAGTAAACAGATCAATGGAAGTAAATTCCAATTGTTTATTGGATTTGCCCCTATATTTCCAGACATCTGTTATc                                                                                                                                                                                                          <  1:352311/250‑1 (MQ=255)                                    acGTTGAAAAACCACATGTGATGATATTGTGATAGTTAGCGATGCCGGGTTTGTATCGTGCCTTTTCATTTATTTCCAGCGGTATAGCTGATAAATAATCCCTTAGTACATCAGGAACCTTATAAAAAAACCCGCGTTCTGCAAACCTTTCTCGATAAAAGTGACAACCACTGAGTAAACAGATCAATGGAAGTAAATTCCAATTGTTTATTGGATTTGCCCCTATATTTCCAGACATCTGTTATCACtt                                                                                                                                                                                                      >  1:267137/1‑250 (MQ=255)                                         gAAAAACCACATGTGATGATATTGTGATAGTTAGCGATGCCGGGTTTGTATCGTGCCTTTTCATTTATTTCCAGCGGTATAGCTGATAAATAATCCCTTAGTACATCAGGAACCTTATAAAAAAACCCGCGTTCTGCAAACCTTTCTCGATAAAAGTGACAACCACTGAGTAAACAGATCAATGGAAGTAAATTCCAATTGTTTATTGGATTTGCCCCTATATTTCCAGACATCTGTTATCACTTAAccc                                                                                                                                                                                                 <  2:36381/250‑1 (MQ=255)                                               ccACATGTGATGATATTGTGATAGTTAGCGATGCCGGGTTTGTATCGTGCCTTTTCATTTATTTCCAGCGGTATAGCTGATAAATAATCCCTTAGTACATCAGGAACCTTATAAAAAAACCCGCGTTCTGCAAACCTTTCTCGATAAAAGTGACAACCACTGAGTAAACAGATCAATGGAAGTAAATTCCAATTGTTTATTGGATTTGCCCCTATATTTCCAGACATCTGTTATCACTTAACCCATTACa                                                                                                                                                                                           >  2:386960/1‑250 (MQ=255)                                                      tgatgaTATTGTGATAGTTAGCGATGCCGGGTTTGTATCGTGCCTTTTCATTTATTTCCAGCGGTATAGCTGATAAATAATCCCTTAGTACATCAGGAACCTTATAAAAAAACCCGCGTTCTGCACACCTTTCTCGATAAAAGTGACAACCACTGAGTAAACAGATCAATGGAAGTAAATTCCAATTGTTTATTGGATTTGCCCCTATATTTCCAGACATCTGTTATCACTTAACCCATTACAAGCCCGCt                                                                                                                                                                                   >  1:277222/1‑251 (MQ=255)                                                       gatgatATTGTGATAGTTAGCGATGCCGGGTTTGTATCGTGCCTTTTCATTTATTTCCAGCGGTATAGCTGATAAATAATCCCTTAGTACATCAGGAACCTTATAAAAAAACCCGCGTTCTGCAAACCTTTCTCGATAAAAGTGACAACCACTGAGTAAACAGATCAATGGAAGTAAATTCCAATTGTTTATTGGATTTGCCCCTATATTTCCAGACATCTGTTATCACTTAACCCATTACAAGCCCGCt                                                                                                                                                                                   <  2:352212/250‑1 (MQ=255)                                                                     aGTTAGCGATGCCGGGTTTGTATCGTGCCTTTTCATTTATTTCCAGCGGTATAGCTGATAAATAATCCCTTAGTACATCAGGAACCTTATAAAAAAACCCGCGTTCTGCAAACCTTTCTCGATAAAAGTGACAACCACTGAGTAAACAGATCAATGGAAGTAAATTCCAATTGTTTATTGGATTTGCCCCTATATTTCCAGACATCTGTTATCACTTAACCCATTACAAGCCCGCTGCCGCAGATATTccc                                                                                                                                                                    >  2:38143/1‑251 (MQ=255)                                                                          gCGATGCCGGGTTTGGATCGTGCCTTTTCATTTGTTTCCAGCGGTATAGCTGATAAATAATCCCTTAGTACATCAGGAACCTTATAACAAAACCCGCGTTCTGCAAACCTTTCTCGATAAAAGTGACAACCACTGAGTAAACAGATCAATGGAAGTAAATTCCAATTGTTTATTGGATTTGCCCCTATATTTCCAGACATCTGTTATCACTTAACCCATTACAAGCCCGCTGCCGCAGATATTCCCGTGgc                                                                                                                                                               >  2:168178/1‑251 (MQ=255)                                                                           cGATGCCGGGTTTGTATCGTGCCTTTTCATTTATTTCCAGCGGTATAGCTGATAAATAATCCCTTAGTACATCAGGAACCTTATAAAAAAACCCGCGTTCTGCAAACCTTTCTCGATAAAAGTGACAACCACTGAGTAAACAGATCAATGGAAGTAAATTCCAATTGTTTATTGGATTTGCCCCTATATTTCCAGACATCTGTTATCACTTAACCCATTACAAGCCCGCTGCCGCAGATATTCCCGTGgcg                                                                                                                                                              >  1:134658/1‑251 (MQ=255)                                                                                 cGGGTTTGTATCGTGCCTTTTCATTTATTTCCAGCGGTATAGCTGATAAATAATCCCTTAGTACATCAGGAACCTTATAAAAAAACCCGCGTTCTGCAAACCTTTCTCGATAAAAGTGACAACCACTGAGTAAACAGATCAATGGAAGTAAATTCCAATTGTTTATTGGATTTGCCCCTATATTTCCAGACATCTGTTATCACTTAACCCATTACAAGCCCGCTGCCGCAGATATTCCCGTGGCGAGCGAt                                                                                                                                                        >  1:337655/1‑251 (MQ=255)                                                                                            cGTGCCTTTTCATTTATTTCCAGCGGTATAGCTGATAAATAATCCCTTAGTACATCAGGAACCTTATAAAAAAACCCGCGTTCTGCAAACCTTTCTCGATAAAAGTGACAACCACTGAGTAAACAGATCAATGGAAGTAAATTCCAATTGTTTATTGGATTTGCCCCTATATTGCCAGACATCTGTGATCACGTAACCCATTACAAGCCCGCTGCCGCAGATATTCCCGTGGCGAGCGATAACCCAGCGCa                                                                                                                                             <  2:179705/251‑1 (MQ=255)                                                                                                ccTTTTCATTGATTGCCAGCGGTATAGCTGATAAATAATCCCTTAGTACATCAGGAACCTGATAAAAAAACCCGCGTTCTGCAAACCTTTCGCGATAAAAGTGACAACCACTGAGTAAACAGATCAATGGAAGTAAATGCCAAGTGTGTATTGGATTTGCCCCTATAGTTCTAGACATCTGTGATCACGTAACCCATTACAAGCCCGCTGCCGCAGATAGTCCCCTTTCGAGCGGGAACGCAGCGCACTAt                                                                                                                                         <  2:44331/251‑1 (MQ=37)                                                                                                                 aGCGGTATAGCTGATAAATAATCCCTTAGTACATCAGGAACCTTATAAAAAAACCCGCGTTCTGCAAACCTTTCTCGATAAAAGTGACAACCACTGAGTAAACAGATCAATGGAAGTAAATTCCAATTGTTTATTGGATTTGCCCCTATATTTCCAGACATCTGTTATCACTTAACCCATTACAAGCCCGCTGCCGCAGATATTCCCGTGGCGAGCGATAACCCAGCGCACTATGCGGATGCCATTCGTta                                                                                                                        <  2:89858/251‑1 (MQ=255)                                                                                                                   cGGTATAGCTGATAAATAATCCCTTAGTACATCAGGAACCTTATAAAAAAACCCGCGTTCTGCAAACCTTTCTCGATAAAAGTGACAACCACTGAGTAAACAGATCAATGGAAGTAAATTCCAATTGTTTATTGGATTTGCCCCTATATTTCCAGACATCTGTTATCACTTAACCCATTACAAGCCCGCTGCCGCAGATATTCCCGTGGCGAGCGATAACCCAGCGCACTATGCGGATGCCATTCGTtat                                                                                                                       >  2:188722/1‑250 (MQ=255)                                                                                                                         agcggATAAATAATCCCTTAGTACATCAGGAACCTTATAAAACAACCCGCGTTCTGCACACCTTTCTCGATAAACATGACAACAACTGAGTAAACGGATCAATGGAAGTAAATTCCAATTGTTAATTGGATTTGCCCCTCTATTTCCAGACATCTGTTCTCACTTAACCCATTACAAGCCCGCTGCCGCAGATATTCCCGTGGCGAGCGCTAACCCAGCGCACTATTCGCATCCCTTTCGTTATAATGccc                                                                                                                >  2:5174/5‑249 (MQ=21)                                                                                                                            tGATAAATAATCCCTTAGTACATCAGGAACCTTATAAAAAAACCCGCGTTCTGCAAACCTTTCTCGATAAAAGTGACAACCACTGAGTAAACAGATCAATGGAAGTAAATTCCAATTGTTTATTGGATTTGCCCCTATATTTCCAGACATCTGTTATCACTTAACCCATTACAAGCCCGCTGCCGCAGATATTCCCGTGGCGAGCGATAACCCAGCGCACTATGCGGATGCCATTCGTTATAATGCTCGac                                                                                                             >  2:215209/1‑250 (MQ=255)                                                                                                                                            aGTACATCAGGAACCTTATAAAAAAACCCGCGTTCTGCAAACCTTTCTCGATAAAAGTGACAACCACTGAGTAAACAGATCAATGGAAGTAAATTCCAATTGTTTATTGGATTTGCCCCTATATTTCCAGACATCTGTTATCACTTAACCCATTACAAGCCCGCTGCCGCAGATATTCCCGTGGCGAGCGATAACCCAGCGCACTATGCGGATGCCATTCGTTATAATGCTCGAACGCCTCTGCAAGGTTc                                                                                             >  1:218604/1‑251 (MQ=255)                                                                                                                                               aCATCAGGAACCTTATAAAAAAACCCGCGTTCTGCAAACCTTTCTCGATAAAAGTGACAACCACTGAGTAAACAGATCAATGGAAGTAAATTCCAATTGTTTATTGGATTTGCCCCTATATTTCCAGACATCTGTTATCACTTAACCCATTACAAGCCCGCTGCCGCAGATATTCCCGTGGCGAGCGATAACCCAGCGCACTATGCGGATGCCATTCGTTATAATGCTCGAACGCCTCTGCAAGGTTCttt                                                                                          >  1:121339/1‑251 (MQ=255)                                                                                                                                                               aaaaaacCCCGCGTTCTGCAAACCTTTCTCGATAAAAGTGACAACCACTGAGTAAACAGATCAATGGAAGTAAATTCCAATTGTTTATTGGATTTGCCCCTATATTTCCAGACATCTGTTATCACTTAACCCATTACAAGCCCGCTGCCGCAGATATTCCCGTGGCGAGCGATAACCCAGCGCACTATGCGGATGCCATTCGTTATAATGCTCGAACGCCTCTGCAAGGTTCTTTGCTGCCGTTAACCCGt                                                                          >  1:331112/1‑251 (MQ=255)                                                                                                                                                                               tGCAAACCTTTCTCGATAAAAGTGACAACCACTGAGTAAACAGATCAATGGAAGTAAATTCCAATTGTTTATTGGATTTGCCCCTATATTTCCAGACATCTGTTATCACTTAACCCATTACAAGCCCGCTGCCGCAGATATTCCCGTGGCGAGCGATAACCCAGCGCACTATGCGGATGCCATTCGTTATAATGCTCGAACGCCTCTGCAAGGTTCTTTGCTGCCGTTAACCCGTCTGGTTTGGGCATGAt                                                          >  2:183712/1‑251 (MQ=255)                                                                                                                                                                                  aaaCCTTTCTCGATAAAAGTGACAACCACTGAGTAAACAGATCAATGGAAGAAAATTCCAATTGTTTATTGGATTTGCCCCTATATTTCCAGACATCTGTTATCACTTAACCCATTACAAGCCCGCTGCCGCAGATATTCCCGTGGCGAGCGATAACCCAGCGCACTATGCGGATGCCATTCGTTATAATGCTCGAACGCCTCTGCAAGGTTCTTTGCTGCCGTTAACCCTTCTGGTTTGGTCATGAtcct                                                       >  1:109113/1‑248 (MQ=255)                                                                                                                                                                                    aCCTTTCTCGATAAAAGTGACAACCACTGAGTAAACAGATCAATGGAAGTAAATTCCAATTGTTTATTGGATTTGCCCCTATATTTCCAGACATCTGTTATCACTTAACCCATTACAAGCCCGCTGCCGCAGATATTCCCGTGGCGAGCGATAACCCAGCGCACTATGCGGATGCCATTCGTTATAATGCTCGAACGCCTCTGCAAGGTTCTTTGCTGCCGTTAACCCGTCTGGTTTGGGCATGATACTg                                                      >  1:69280/1‑250 (MQ=255)                                                                                                                                                                                                               tGAGTAAACAGATCAATGGAAGTAAATTCCAATTGTTTATTGGATTTGCCCCTATATTTCCAGACATCTGTTATCACTTAACCCATTACAAGCCCGCTGCCGCAGATATTCCCGTGGCGAGCGATAACCCAGCGCACTATGCGGATGCCATTCGTTATAATGCTCGAACGCCTCTGCAAGGTTCTTTGCTGCCGTTAACCCGTCTGGTTTGGGCATGATACTGATGTAGTCACGCTTTATCGTTTTCACg                           >  1:166955/1‑250 (MQ=255)                                                                                                                                                                                                                     aaCAGATCAATGGAAGTAAATTCCAATTGTTTATTGGATTTGCCCCTATATTTCCAGACATCTGTTATCACTTAACCCATTACAAGCCCGCTGCCGCAGATATTCCCGTGGCGAGCGATAACCCAGCGCACTATGCGGATGCCATTCGTTATAATGCTCGAACGCCTCTGCAAGGTTCTTTGCTGCCGTTAACCCGTCTGGTTTGGGCATGATACTGATGTAGTCACGCTTTATCGTTTTCACGAAGctct                    >  2:114173/1‑251 (MQ=255)                                                                                                                                                                                                                        aGATCAATGGAAGTAAATTCCAATTGTTTATTGGATTTGCCCCTATATTTCCAGACATCTGTTATCACTTAACCCATTACAAGCCCGCTGCCGCAGATATTCCCGTGGCGAGCGATAACCCAGCGCACTATGCGGATGCCATTCGTTATAATGCTCGAACGCCTCTGCAAGGTTCTTTGCTGCCGTTAACCCGTCTGGTTTGGGCATGATACTGATGTAGTCACGCTTTATCGTTTTCACGAAGCTCTCTg                 >  1:339465/1‑251 (MQ=255)                                                                                                                                                                                                                                       aaTTCCAATTGTTTATTGGATTTGCCCCTATATTTCCAGACATCTGTTATCACTTAACCCATTACAAGCCCGCTGCCGCAGATATTCCCGTGGCGAGCGATAACCCAGCGCACTATGCGGATGCCATTCGTTATAATGCTCGAACGCCTCTGCAAGGTTCTTTGCTGCCGTTAACCCGTCTGGTTTGGGCATGATACTGCTGTAGTCACGCTTTATCGTTTTCACGAAGCTCTCTGCTATTCCGTTActct  >  2:16539/1‑251 (MQ=25)                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                           CATTGCCAATGGTTGTTGCCTTATATAGGGTAACAACGTTGAAAAACCACATGTGATGATATTGTGATAGTTAGCGATGCCGGGTTTGTATCGTGCCTTTTCATTTATTTCCAGCGGTATAGCTGATAAATAATCCCTTAGTACATCAGGAACCTTATAAAAAAACCCGCGTTCTGCAAACCTTTCTCGATAAAAGTGACAACCACTGAGTAAACAGATCAATGGAAGTAAATTCCAATTGTTTATTGGATTTGCCCCTATATTTCCAGACATCTGTTATCACTTAACCCATTACAAGCCCGCTGCCGCAGATATTCCCGTGGCGAGCGATAACCCAGCGCACTATGCGGATGCCATTCGTTATAATGCTCGAACGCCTCTGCAAGGTTCTTTGCTGCCGTTAACCCGTCTGGTTTGGGCATGATACTGATGTAGTCACGCTTTATCGTTTTCACGAAGCTCTCTGCTATTCCGTTACTCT  >  NC_000913/2996115‑2996595

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