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)	37 (1.020)	33/486	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)	46 (1.270)	44/486	0.1	100%	coding (256/2235 nt)	relA	GDP/GTP pyrophosphokinase
	?	NC_000913	2996361 =	NA (NA)					noncoding (1331/1331 nt)	IS2	repeat region

CTTACGCCAGAGAAATAACTGGCTGGCTGCTACACCATGTTGCCGGGCAACGAGGGAGACCGTCATCCCCGGTTCAAAGCTCTGCTGAACAATTGCGATCTTTTCCTGTGTGGTACGCCGTCTGCGTTTCTCCGGCCCTAAGACATCAATCATCTGTTCTCCAATGACTAGTCTAAAAACTAGTATTAAGACTATCACTTATTTAAGTGATATTGGTTGTCTGGAGATTCAGGGGGCCAGTCTAAACTTGCTCTTTTCTTCTGGTGTTTATCTTTGAAAGATTTGTAGGATAAATAACATTAAGTAAAGCCGCATCGGCTTCTCTGTAATGTTTATCTTTTTCCAGATAACCTGGGGATTCTCCCCCTAAATCAAGCTGTAACGTTATTCTTTCATAACCTTTATCTATGAGGTGGTTCAACCTCGTTGTTGAAAGAGTTTCATCCAGTTCATTGACTTTAGCTTCGCGTAGTAAATCAGAGACAAAAACAGA  >  NC_000913/2997448‑2997940                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                                          cTTACGCCAGAGAAATAACTGGCTGGCTGCTACACCATGTTGCCGGGCAACGAGGGAGACCGTCATCCCCGGTTCAAAGCTCTGCTGAACAATTGCGATCTTTTCCTGTGTGGTACGCCGTCTGCGTTTCTCCGGCCCTAAGACATCAATCATCTGTTCTCCAATGACTAGTCTAAAAACTAGTATTAAGACTATCACTTATTTAAGTGATATTGGTTGTCTGGAGATTCAGGGGGCCAGTCTAAACTTGc                                                                                                                                                                                                                                                    >  2:245028/1‑251 (MQ=31)       ccAGAGAAATAACTGGCTGGCTGCTACACCATGTTGCCGGGCAACGAGGGAGACCGTCATCCCCGGTTCAAAGCTCTGCTGAACAATTGCGATCTTTTCCTGTGTGGTACGCCGTCTGCGTTTCTCCGGCCCTAAGACATCAATCATCTGTTCTCCAATGACTAGTCTAAAAACTAGTATTAAGACTATCACTTATTTAAGTGATATTGGTTGTCTGGAGATTCAGGGGGCCAGTCTAAACTTGCTCtttt                                                                                                                                                                                                                                              <  2:58173/251‑1 (MQ=33)       caagagacagaacgggctggctGCTACACCATGTTGCCGGGCAACGAGGGAGACCGTCATTCCCGGTTCAAAGCTCTGCTGAACAATTGCGATCTTTTCCTGTGTGGTACGCCGTCTGCGTTTCTCCGGCCCTAAGACATCAATCATCTGTTCTCCAATGACTAGTCTAAAAACTAGTATTAAGACTATCACTGATTTAAGTGATATTGGTTGTCTGGAGACTCAGGGGGACAGTCTAAACTTGCTCttct                                                                                                                                                                                                                                              <  2:16296/237‑3 (MQ=11)          gagaAATAACTGGCTGGCTGCTACACCATGTTGCCGGGCAACGAGGGAGACCGTCATCCCCGGTTCAAAGCTCTGCTGAACAATTGCGATCTTTTCCTGTGTGGTACGCCGTCTGCGTTTCTCCGGCCCTAAGACATCAATCATCTGTTCTCCAATGACTAGTCTAAAAACTAGTATTAAGACTATCACTTATTTAAGTGATATTGGTTGTCTGGAGATTCAGGGGGCCAGTCTTAACTTGCTCTTttctt                                                                                                                                                                                                                                           <  1:327836/251‑1 (MQ=25)                         ggctgCTACACCATGTTGCCGGGCAACGAGGGAGACCGTCATCCCCGGTTCAAAGCTCTGCTGAACAATTGCGATCTTTTCCTGTGTGGTACGCCGTCTGCGTTTCTCCGGCCCTAAGACATCAATCATCTGTTCTCCAATGACTAGTCTAAAAACTAGTATTAAGACTATCACTTATTTAAGTGATATTGGTTGTCTGGAGATTCAGGGGGCCAGTCTAAACTTGCTCTTTTCTTCTGGTGTTTATCttt                                                                                                                                                                                                                            <  1:164246/251‑1 (MQ=255)                                acacCATGTTGCCGGGCAACGAGGGAGACCGTCATCCCCGGTTCAAAGCTCTGCTGAACAATTGCGATCTTTTCCTGTGTGGTACGCCGTCTGCGTTTCTCCGGCCCTAAGACATCAATCATCTGTTCTCCAATGACTAGTCTAAAAACTAGTATTAAGACTATCACTTATTTAAGTGATATTGGTTGTCTGGAGATTCAGGGGGCCAGTCTAAACTTGCTCTTTTCTTCTGGTGTTTATCTTTGAAAGa                                                                                                                                                                                                                      <  1:238925/250‑1 (MQ=255)                                          gCCGGGCAACGAGGGAGACCGTCATCCCCGGTTCAAAGCTCTGCTGAACAATTGCGATCTTTTCCTGTGTGGTACGCCGTCTGCGTTTCTCCGGCCCTAAGACATCAATCATCTGTTCTCCAATGACTAGTCTAAAAACTAGTATTAAGACTATCACTTATTTAAGTGATATTGGTTGTCTGGAGATTCAGGGGGCCAGTCTAAACTTGCTCTTTTCTTCTGGTGTTTATCTTTGAAAGATTTGTAGGat                                                                                                                                                                                                            <  2:10624/250‑1 (MQ=255)                                                cAACGAGGGAGACCGTCATCCCCGGTTCAAAGCTCTGCTGAACAATTGCGATCTTTTCCTGTGTGGTACGCCGTCTGCGTTTCTCCGGCCCTAAGACATCAATCATCTGTTCTCCAATGACTAGTCTAAAAACTAGTATTAAGACTATCACTTATTTAAGTGATATTGGTTGTCTGGAGATTCAGGGGGCCAGTCTAAACTTGCTCTTTTCTTCTGGTGTTTATCTTTGAAAGATTTGTAGGATAAATAAc                                                                                                                                                                                                     >  2:5384/1‑251 (MQ=255)                                                    gAGGGAGACCGTCATCCCCGGTTCAAAGCTCTGCTGAACAATTGCGATCTTTTCCTGTGTGGTACGCCGTCTGCGTTTCTCCGGCCCTAAGACATCAATCATCTGTTCTCCAATGACTAGTCTAAAAACTAGTATTAAGACTATCACTTATTTAAGTGATATTGGTTGTCTGGAGATTCAGGGGGCCAGTCTAAACTTGCTCTTTTCTTCTGGTGTTTATCTTTGAAAGATTTGTAGGATAAATAACAtt                                                                                                                                                                                                  >  1:222803/1‑250 (MQ=255)                                                                        gTTCAAAGCTCTGCTGAACAATTGCGATCTTTTCCTGTGTGGTACGCCGTCTGCGTTTCTCCGGCCCTAAGACATCAATCATCTGTTCTCCAATGACTAGTCTAAAAACTAGTATTAAGACTATCACTTATTTAAGTGATATTGGTTGTCTGGAGATTCAGGGGGCCAGTCTAAACTTGCTCTTTTCTTCTGGTGTTTATCTTTGAAAGATTTGTAGGATAAATATCATTAAGGAAAGCCGCATCGGCTtc                                                                                                                                                                             <  1:315400/251‑1 (MQ=255)                                                                                    gctgAACAATTGCGATCTTTTCCTGTGTGGTACGCCGTCTGCGTTTCTCCGGCCCTAAGACAGCAATCATCTGTGCTCCAATGACTAGTCTAAAAACTAGTATTAAGACTATCACTTAGTTAAGTGATATTGGTTGTCTGGAGATTCAGGGGGCCAGTCGAAACTTGCTCTTTTCTTCTGGTGTTTATCTTTGAAAGATTTGTAGGATAAAGAACATTAAGTAAAGCCGCATCGGCTTCTCTGTAATGttt                                                                                                                                                                 <  1:375163/251‑1 (MQ=255)                                                                                          cAATTGCGATCTTTTCCTGTGTGGTACGCCGTCTGCGTTTCTCCGGCCCTAAGACATCAATCATCTGTTCTCCAATGACTAGTCTAAAAACTAGTATTAAGACTATCACTTATTTAAGTGATATTGGTTGTCTGGAGATTCAGGGGGCCAGTCTAAACTTGCTCTTTTCTTCTGGTGTTTATCTTTGAAAGATTTGTAGGATAAATAACATTAAGTAAAGCCGCATCGGCTTCTCTGTAATGTTTATCttt                                                                                                                                                           <  1:276272/251‑1 (MQ=255)                                                                                          cAATTGCGATCTTTTCCTGTGTGGTACGCCGTCTGCGTTTCTCCGGCCCTAAGACATCAATCATCTGTTCTCCAATGACTAGTCTAAAAACTAGTATTAAGACTATCACTTATTTAAGTGATATTGGTTGTCTGGAGATTCAGGGGGCCAGTCTAAACTTGCTCTTTTCTTCTGGTGTTTATCTTTGAAAGATTTGTAGGATAAATAACATTAAGTAAAGCCGCATCGGCTTCTCTGTAATGTTTATCttt                                                                                                                                                           <  2:222803/251‑1 (MQ=255)                                                                                          cAATTGCGATCTTTTCCTGTGTGGTACGCCGTCTACGTTTCTCCGGCCCTAAGCCATCAATCATCTGTTCTCCAATGACTAGTCTAAAATCTAGTATTAAGACTATCACTTATTTAAGTGATATTGGTTGTCTGGAGATTCAGGGGGCCAGTCTACACTTGCTCTTTTCTTCTTGTGTTTATCTTTGAAAGATTTGTAGGATAAATAACATTAAGTTAAGCCGCATCGGCTTCTGTGTAATGTTTATCttt                                                                                                                                                           <  1:276283/251‑1 (MQ=18)                                                                                                           tgtgtgGTACGCCGTCTGCGTTTCTCCGGCCCTAAGACATCAATCATCTGTTCTCCAATGACTAGTCTAAAAACTAGTATTAAGACTATCACTTATTTAAGTGATATTGGTTGTCTGGAGATTCAGGGGGCCAGTCTAAACTTGCTCTTTTCTTCTGGTGTTTATCTTTGAAAGATTTGTAGGATAAATAACATTAAGTAAAGCCGCATCGGCTTCTCTGTAATGTTTATCTTTTTCCAGATAACCTgggg                                                                                                                                          >  1:300167/1‑251 (MQ=255)                                                                                                                     gCCGTCTGCGTTTCTCCGGCCCTAAGACATCAATCATCTGTTCTCCAATGACTAGTCTAAAAACTAGTATTAAGACTATCACTTATTTAAGTGATATTGGTTGTCTGGAGATTCAGGGGGCCAGTCTAAACTTGCTCTTTTCTTCTGGTGTTTATCTTTGAAAGATTTGTAg                                                                                                                                                                                                               >  2:171869/1‑172 (MQ=255)                                                                                                                     gCCGTCTGCGTTTCGCCGGCCCTAAGACATCAATCATCTGTTCTCCAATGGCTAGTCTAAAAACTAGTAGTAAGACTATCACTTATGTAAGTGATATTGGTGGTCTGGATATTCAGGGGGCCAGTCTAATCGTTCTCTTTTCTTCTGGTGTTTATCTTTGAAAGATTTGTAg                                                                                                                                                                                                               <  1:171869/172‑1 (MQ=17)                                                                                                                            tcgtgtctcCGGCCCTAAGTCATCAAATGTCTGTTCTCCAATGACTAGACTAAAAACTGGTTTTAAGACTATCACTTATTTAAGTGTTATTGGTTGTCTGGAGATTCAGGGGGACAGTCTAAGCTTGCTCTTTTCTTCTGGTGTTTATCTTTGAAAGATTTGTAGGCTATATAACATTAAGTAAAGCAGCATCGGCTTCTCTGTAATGTTTCTCTTTTTCCAGATATCCTGGGGGTTATCCCCCTAAAtg                                                                                                                          <  2:289050/245‑2 (MQ=255)                                                                                                                            gCGTTTCTCCGGCCCTAAGACATCAATCATCTGTTCTCCAATGACTAGTCTAAAAACTAGTATTAAGACTATCACTTATTTAAGTGATAGTGGTTGTCTGGAGATTCAGGGGGCCAGTCTAAACTTGCTCTTTTCTTCTGGTGTTTATCTTTGAAAGATTTGTAGGATAAATAACATTAAGTAAAGCCGCATCGGCTTCTCTGTAATGTTTATCTTTTTCCAGATAACCTGGGGATTCTCCCCCTAAAt                                                                                                                           <  2:296415/249‑1 (MQ=255)                                                                                                                                     cGGCCCTAAGACATCAATCATCTGTTCTCCAATGACTAGTCTAAAAACTAGTATTAAGACTATCACTTATTTAAGTGATATTGGTTGTCTGGAGATTCAGGGGGCCAGTCTAAACTTGCTCTTTTCTTCTGGTGTTTATCTTTGAAAGATTTGTAGGATAAATAACATTAAGTAAAGCCGCATCGGCTTCTCTGTAATGTTTATCTTTTTCCAGATAACCTGGGGATTCTCCCCCTAAATCAAGCTGTAAc                                                                                                                <  2:185249/251‑1 (MQ=255)                                                                                                                                              gACAGCAATCATCTGTTCTCCAATGACTAGTCTAAAAACTAGTATTAAGACTATCACTTATTTATGTGATATTGGTTGTCTGGAGATTCAGGGGGCCAGTCTAAACTTGCTCTTTTCTTCTGGTGTATATCTTTGAAAGATTTGTAGGATAAATAACATTAAGTAAAGCCGCATCGGCTTCTCTGTAATGTTTATCTTTTTCCAGATAACCTGGGGATTCTCCCCCTAAATCAAGCTGTAACGTTATTCtt                                                                                                       <  2:253975/251‑1 (MQ=255)                                                                                                                                                                          aGTCTAAAAACTAGTATTAAGACTATCACTTATTTAAGTGATATTGGTTGTCTGGAGATTCAGGGGGCCAGTCTAAACTTGCTCTTTTCTTCTGGTGTTTATCTTTGAAAGATTTGTAGGATAAATAACATTAAGTAAAGCCGCATCGGCTTCTCTGTAATGTTTATCTTTTTCCAGATAACCTGGGGATTCTCCCCCTAAATCAAGCTGTAACGTTATTCTTTCATAACCTTTATCTATGAGGTGGTTCa                                                                           <  2:6345/251‑1 (MQ=255)                                                                                                                                                                                          ttAAGACTATCACTTATTTAAGTGATATTGGTTGTCTGGAGATTCAGGGGGCCAGTCTAAACTTGCTCTTTTCTTCTGGTGTTTATCTTTGAAAGATTTGTAGGATAAATAACATTAAGTAAAGCCGCATCGGCTTCTCTGTAATGTTTATCTTTTTCCAGATAACCTGGGGATTCTCCCCCTAAATCAAGCTGTAACGTTATTCTTTCATAACCTTTATCTATGAGGTGGTTCAACCTCGTTGTTGaa                                                             >  1:34735/1‑249 (MQ=255)                                                                                                                                                                                          ttAAGACTATAACTTATTTAAGTGAGATTGGTTGTCTGGAGATTAAGGGGGCCAGTATAAACAAGCTCTTTTCTTCTGGTGTTTATCTTTGAAAGATTAGTAGGATAAATAACATTAAGTAAAGCCGCATCGGCTTCTCTGTAATGTTTATCTTTTTCCAGATAACCTGGGGATTCTGCCCCTAAATCAAGCTGTAAAGTTATTCTTTAATAACTTTTATATATGAGGAGGTTCAAGCTCGTTGTGGaa                                                             >  1:34716/1‑249 (MQ=255)                                                                                                                                                                                                    cACTTATTTAAGTGATATTGGTTGTCTGGAGATTCAGGGGGCCAGTCTAAACTTGCTCTTTTCTTCTGGTGTTTATCTTTGAAAGATTTGTAGGATAAATAACATTAAGTAAAGCCGCATCGGCTTCTCTGTAATGTTTATCTTTTTCCAGATAACCTGGGGATTCTCCCCCTAAATCAAGCTGTAACGTTATTCTTTCATAACCTTTATCTATGAGGTGGTTCAACCTCGTTGTTGAAAGAGTTTCATc                                                  <  2:87579/250‑1 (MQ=255)                                                                                                                                                                                                     aCTTATTTAAGTGATATTGGTTGTCTGGAGATTCAGGGGGCCAGTCTAAACTTGCTCTTTTCTTCTGGTGTTTATCTTTGAAAGATTTGTAGGATAAATAACATTAAGTAAAGCCGCATCGGCTTCTCTGTAATGTTTATCTTTTTCCAGATAACCTGGGGATTCTCCCCCTAAATCAAGCTGTAACGTTATTCTTTCATAACCTTTATCTATGAGGTGGTTCAACCTCGTTGTTGAAAGAGTTTCATCCa                                                <  1:245028/251‑1 (MQ=255)                                                                                                                                                                                                                   tatTGGCTGTGTGGAGATTCAGGGGGCCCGTATACACTTGCTCTTTTCTTCTGGTGTTTATCTTTGAAAGATTTGTAGGATAAATAACATTAAGTCAAGCCGCATCGGCTTCTCTGTAATGTTTATCTTTTTCCAGATAACCTGGGGATTATCCCCCTAAATCAAGCTGTCACGTTATTCTTTCATAACCTTTATCTATGAGGTGGTTCAACCTCGTTGTTGAAAGAGTTTCATCAAGTTCATTGCCttt                                   >  2:244258/1‑250 (MQ=255)                                                                                                                                                                                                                                          gggCCAGTCTAAACTTGCTCTTTTCTTCTGGTGTTTATCTTTGAAAGATTTGTAGGATAAATAACATTAAGTAAAGCCGCATCGGCTTCTCTGTAATGTTTATCTTTTTCCAGATAACCTGGGGATTCTCCCCCTAAATCAAGCTGTAACGTTATTCTTTCATAACCTTTATCTATGAGGTGGTTCAACCTCGTTGTTGAAAGAGTTTCATCCAGTTCATTGACTTTAGCTTCGCGTAGTAAATCAgaga            >  2:351277/1‑250 (MQ=255)                                                                                                                                                                                                                                            gCCAGTCTAAACTTGCTCTTTTCTTCTGGTGTTTATCTTTGAAAGATTTGTAGGATAAATAACATTAAGTAAAGCCGCATCGGCTTCTCTGTAATGTTTATCTTTTTCCAGATAACCTGGGGATTCTCCCCCTAAATCAAGCTGTAACGTTATTCTTTCATAACCTTTATCTATGAGGTGGTTCAACCTCGTTGTTGAAAGAGTTTCATCCAGTTCATTGACTTTAGCTTCGCGTAGTAAATCAGAGACa          >  1:184995/1‑250 (MQ=255)                                                                                                                                                                                                                                             ccAGTCTAAACTTGCTCTTTTCTTCTGGTGTTTATCTTTGAAAGATTTGTAGGATAAATAACATTAAGTAAAGCCGCATCGGCTTCTCTGTAATGTTTATCTTTTTCCAGATAACCTGGGGATTCTCCCCCTAAATCAAGCTGTAACGTTATTCTTTCATAACCTTTATCTATGAGGTGGTTCAACCTCGTTGTTGAAAGAGTTTCATCCAGTGCATTGACTTTAGCTTCGCGTAGTAAATCAGAGACaaa        <  2:175459/251‑1 (MQ=255)                                                                                                                                                                                                                                              cAGTCTAAACTTGCTCTTTTCTTCTGGTGTTTATCTTTGAAAGATTTGTAGGATAAATAACATTAAGTAAAGCCGCATCGGCTTCTCTGTAATGTTTATCTTTTTCCAGATAACCTGGGGATTCTCCCCCTAAATCAAGCTGTAACGTTATTCTTTCATAACCTTTATCTATGAGGTGGTTCAACCTCGTTGTTGAAAGAGTTTCATCCAGTTCATTGACTTTAGCTTCGCGTAGTAAATCAGAGACaaa        >  2:79908/1‑250 (MQ=255)                                                                                                                                                                                                                                              cAGTCTAAACTTGCTCTTTTCTTCTGGCGTTTATCTTTGAAAGATTTGTAGGATAAATACCATTAAGTAAAGCCGCATCGGCTTCTCTGTAATGTTTATCTTTTTCCAGATAACCTGGGGATTCTCCCCCTAAATCAAGATGTAACGTTATTCTTTCATAACCTTTATCTATGAGGTGGTTCAACCTCGTTGTTGAACGCGTTTCATCCAGTTCATTGACTTTAGCTTCGCGTAGTAAATCAGAGACaaa        >  2:334273/1‑250 (MQ=255)                                                                                                                                                                                                                                                    aaaCTTGCTCTTTTCTTCTGGTGTTTATCTTTGAAAGATTTGTAGGATAAATAACATTAAGTAAAGCCGCATCGGCTTCTCTGTAATGTTTATCTTTTTCCAGATAACCTGGGGATTCTCCCCCTAAATCAAGCTGTAACGTTATTCTTTCATAACCTTTATCTATGAGGTGGTTCAACCTCGTTGTTGAAAGAGTTTCATCCAGTTCATTGACTTTAGCTTCGCGTAGTAAATCAGAGACAAAAACaga  <  1:189522/250‑1 (MQ=255)                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                                          CTTACGCCAGAGAAATAACTGGCTGGCTGCTACACCATGTTGCCGGGCAACGAGGGAGACCGTCATCCCCGGTTCAAAGCTCTGCTGAACAATTGCGATCTTTTCCTGTGTGGTACGCCGTCTGCGTTTCTCCGGCCCTAAGACATCAATCATCTGTTCTCCAATGACTAGTCTAAAAACTAGTATTAAGACTATCACTTATTTAAGTGATATTGGTTGTCTGGAGATTCAGGGGGCCAGTCTAAACTTGCTCTTTTCTTCTGGTGTTTATCTTTGAAAGATTTGTAGGATAAATAACATTAAGTAAAGCCGCATCGGCTTCTCTGTAATGTTTATCTTTTTCCAGATAACCTGGGGATTCTCCCCCTAAATCAAGCTGTAACGTTATTCTTTCATAACCTTTATCTATGAGGTGGTTCAACCTCGTTGTTGAAAGAGTTTCATCCAGTTCATTGACTTTAGCTTCGCGTAGTAAATCAGAGACAAAAACAGA  >  NC_000913/2997448‑2997940

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