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

ATGTTGCCGGGCAACGAGGGAGACCGTCATCCCCGGTTCAAAGCTCTGCTGAACAATTGCGATCTTTTCCTGTGTGGTACGCCGTCTGCGTTTCTCCGGCCCTAAGACATCAATCATCTGTTCTCCAATGACTAGTCTAAAAACTAGTATTAAGACTATCACTTATTTAAGTGATATTGGTTGTCTGGAGATTCAGGGGGCCAGTCTAAACTTGCTCTTTTCTTCTGGTGTTTATCTTTGAAAGATTTGTAGGATAAATAACATTAAGTAAAGCCGCATCGGCTTCTCTGTAATGTTTATCTTTTTCCAGATAACCTGGGGATTCTCCCCCTAAATCAAGCTGTAACGTTATTCTTTCATAACCTTTATCTATGAGGTGGTTCAACCTCGTTGTTGAAAGAGTTTCATCCAGTTCATTGACTTTAGCTTCGCGTAGTAAATCAGAGA  >  NC_000913/2997484‑2997930                                                                                                                                                                                                                |                                                                                                                                                                                                                                                aTGTTGCCGGGCAACGAGGGAGACCGTCATCCCCGGTTCAAAGCTCTTCTGAACAATTGCGATCTTTTCCTGTGTGGTACGCCGTCTGCGTTTCTCCGGCCCTAAGACATCAATCATCTGTTCTCCAATGACTAGTCTAAAAACTAGTATTAAGACTATCACTTATTTAAGTGATATTGGTTGTCTGGAGATTCAGGGGGCCAGTCTAAACTTGCTCTTTTCTTCTGGTGTTTATCTTTGAAAGATTTGTa                                                                                                                                                                                                      <  1:180431/251‑1 (MQ=255)        cGGGCAACGAGGGAGACCGTCATCCCCGGTTCAAAGCTCTGCTGAACAATTGCGATCTTTTCCTGTGTGGTACGCCGTCTGCGTTTCTCCGGCCCTAAGACATCAATCAGCTGTTCTCCAATGACTAGTCGAAAAACTAGTAGTAAGACTATCACTTATGTAAGTGATATTGGTTGTCTGGAGATTCAGGGGGCCAGTCTAAACTTGCTCTTttc                                                                                                                                                                                                                                   <  1:337106/215‑1 (MQ=255)                          gTCATCCCCGGTTCAAAGCTCTGCTGAACAATTGCGATCTTTTCCTGTGTGGTACGCCGTCTGCGTTTCTCCGGCCCTAAGACATCAATCATCTGTTCTCCAATGACTAGTCTAAAAACTAGTATTAAGACTATCACTTATTTAAGTGATATTGGTTGTCTGGAGATTCAGGGGGCCAGTCTAAACtt                                                                                                                                                                                                                                            >  1:167773/1‑188 (MQ=31)                                 ccGGTTCAAAGCTCTGCTGAACAATTGCGATCTTTTCCTGTGTGGTACGCCGTCTGCGTTTCTCCGGCCCTAAGACATCAATCATCTGTTCTCCAATGCCTAGTCTAAAAACTAGTATTAAGACTATCACTTATTTAAGTGATATTGGTTGTCTGGAGATTCAGGGGGCCAGTCTAAACTTGCTCTTTTCTTCTGGTGTTTATCTTTGAAAGATTTGTAGGATAAATAACATTAAGTAAAGCCGCATCGGc                                                                                                                                                                      >  1:432165/1‑251 (MQ=255)                                                   gAACAATTGCGATCTTTTCCTGTGTGGTACGCCGTCTGCGTTTCTCCGGCCCTACGACATCACTCATCTGTTCTCCACTGACTAGTCTAAAAACTAGTATTAAGACTATCACTTATTTAAGTGATATTGGTTGTCTGGAGATTCAGGGGGCCAGTCTCAACTTGCTCTTTTCTTCTGGTGTTTATCTTTGACAGATTTGTa                                                                                                                                                                                                      >  1:211294/1‑201 (MQ=255)                                                             gATCTTTTCCTGTGTGGTACGCCGTCTGCGTTTCTCCGGCCCTAAGACATCAATCATCTGTTCTCCAATGACTAGTCTAAAAACTAGTATTAAGACTATCACTTATTTAAGTGATATTGGTTGTCTGGAGATTCAGGGGGCCAGTCTAAACTTGCTCTTTTCTTCTGGTGTTTATCTTTGAAAGATTTGTAGGATAAATAACATTAAGTAAAGCCGCATCGGCTTCTCTGTAATGTTTATCTTTTTcc                                                                                                                                             >  1:458080/1‑248 (MQ=255)                                                                      cTGTGTGGTACGCCGTCTGCGTTTCTCCGGCCCTAAGAAATCAATCATCTGTTCTCCAATGACTAGTCTAAAAACTAGTATTAAGACTATCACTTATTTAAGTGATATTGGTTGTCTGGAGATTCAGGGGGCCAGTCTAAACTTGCTCTTTTCTTCTGGTGTTTATCTTTGAAAGATTTGTAGGATAAATAACATTAAGTAAAGCCGCATCGGCTTCTCTGTAATGTTTATCTTTTTCCAGATAACCTggg                                                                                                                                 >  1:473204/1‑251 (MQ=255)                                                                          gtgGTACGCCGTCTGCGTTTCTCCGGCCCTAAGACATCAATCATCTGTTCTCCAATGACTAGTATCAAAACTAGTCTTAAGACTATCACTTATTTAAGTGATATTGGTTGTATGGAGATTCAGGGGGCCAGTCTAAACTTGCTATTTTCTTCTGGTGTTTATCTTTGAAAGATTTGTAGGATAAATAACATTAAGTAAAGCCGCATCGGCTTCTCTGTAATGTTTATCTTTTTCCAGATAACCTGGGGAtt                                                                                                                             >  1:382911/1‑251 (MQ=255)                                                                             gTACGCCGTCTGCGTTTCTCCGGCCCTAAGACATCAATCATCTGTTCTCCAATGACTAGTCTAAAAACTAGTATTAAGACTATCACTTATTTAAGTGATATTGGTTGTCTGGAGATTCAGGGGGCCAGTCTAAACTTGCTCGTTTCTTCTGGTGTTTATCTTTGAAAGATTTGTAGGATAAATAACATTAAGTAAAGCCGCATCGGCTTCTCTGTAATGTTTATCTTTTTCCAGATAACCTGGGGATtctc                                                                                                                          <  1:455068/251‑1 (MQ=255)                                                                                           tttCTCCGGCCCTAAGACATCAATCATCTGTTCTCCAATGACTAGTCTAAAAACTAGTATTAAGACTATCACTTATTTAAGTGATATTGGTTGTCTGGAGATTCAGGGGGCCAGTCTAAACTTGCTCTTTTCTTCTGGTGTTTATCTTTGAAAGATTTGTAGGATAAATAACATTAAGTAAAGCCGCATCGGCTTCTCTGTAATGTTTATCTTTTTCCAGATAACCTGGGGATtct                                                                                                                           <  1:208699/236‑1 (MQ=255)                                                                                                ccGGCCCTAAGACATCAATCATCTGTTCTCCAATGACTAGTCTAAAAACTAGTATTAAGACTATCACTTATTTAAGTGATATTGGTTGTCTGGAGATTCAGGGGGCCAGTCTAAACTTGCTCTTTTCTTCTGGTGTTTATCTTTGAAAGATTTGTAGGATAAATAACATTAAGTAAAGc                                                                                                                                                                               >  1:350886/1‑179 (MQ=255)                                                                                                ccGGCCCTAAGACAGCAATCATCTGTTCTCCAATGACTAGTCTAAAAACTAGTATTAAGACTATCACTTATTTAAGTGATATGGGTTGTCNGGANAGTCNGGNGGCCAGTCTAAACTTNCTCTTTTTTTCTGGGGTTTATCTTTGAAAGATTTGTAGGATAAANNANNNGNAGTAAAGc                                                                                                                                                                               <  2:350886/179‑1 (MQ=18)                                                                                                      gtAAGACATCAATCATCTGTTCTCCAATGACTAGTCGAAAAACGAGTATGAAGACTATCACTTATTTAAGTGATATTGGTTGTCTGGAGATTCAGGGGGCCAGTCTAAACTTGCTCTTTTCTTCTGGTGTTTATCttt                                                                                                                                                                                                                  <  1:263294/137‑1 (MQ=255)                                                                                                         aGACATCAATCATCTGTTCTCCAATGACTAGTCTAAAAACTAGTATTAAGACTATCACTTATTTAAGTGATATTGGTTGTCTGGAGATTCAGGGGGCCAGTCTAAACTTGCTCTTTTCTTCTGGTGTTTATCTTTGAAAGATTTGTAGGATAAATAACATTAAGTAAAGCCGCATCGGCTTCTCTGTAATGTTTATCTTTTTCCAGATAACCTGGGGATTCTCCCCCt                                                                                                                     <  1:439225/228‑1 (MQ=255)                                                                                                                               aaTGACTAGTCTAAAAACTAGTATTAAGACTATCACTTATTTAAGTGATATTGGTTGTCTGGAGATTCAGGGGGCCAGTCTAAACTTGCTCTTTTCTTCTGGTGTTTATCTTTGAAAGATTTGTAGGATAAATAACATTAAGTAAAGCCGCATCGGCTTCTCTGTAATGTTTATCTTTTTCCAGATAACCTGGGGATTCTCCCCCTAAATCAAGCTGTAACGGTATTCTTTCATAACCGTTATCTATGAgn                                                                        <  1:150128/251‑2 (MQ=255)                                                                                                                                      aGTCTAAAAACTAGTTTTAAGACTATCACTTATTTAAGTGATATTGGTTGTCTGGAGATTCAGGGGGCCAGTCTAAACTTGCTCTTTTCTTCTGGTGTTTATCTTTGAAAGATTTGTAGGATAAATAACATTAAGTAAAGCCGCATCGGCTTCTCTGTAATGTTTATCTTTTTCCAGATAACCTGGGGATTCTCCCCCTAAATCAAGCTGTAACGTTATTCTTTCATAACCTTTATCTATGAGGTGGTTc                                                                  <  1:296179/250‑1 (MQ=255)                                                                                                                                           aaaaaCTAGTATTAAGACTATCACTTATTTAAGTGATATTGGTTGTCTGGAGATTCAGGGGGCCAGTCTAAACTTGCTCTTTTCTTCTGGTGTTTATCTTTGAAAGATTTGTAGGATAAATAACATTAAGTAAAGCCGCATCGGCTTCTCTGTAATGTTTATCTTTTTCCAGATAACCTGGGGATTCTCCCCCTAAATCAAGCTGTAACGTTATTCTTTCATAACCTTTATCTATTAGGTGGTTCAAcc                                                              <  1:400381/249‑1 (MQ=255)                                                                                                                                            aaaaCTAGTATTAAGACTATCACTTATGTAAGTGATATCTGCTGTCTGGAGCTTCAGGGTGCCAGTCTAAACTTGCTCTTTTCTTCTGGTGTTTATCTTTg                                                                                                                                                                                                                 <  1:445387/101‑1 (MQ=18)                                                                                                                                                              aTCACNTCTTTNANNGACATTGGTTGTCTGGAGANTCCGGGGGCCAGTCTNAACTNGCTCtnttnt                                                                                                                                                                                                                                  >  2:196074/1‑66 (MQ=14)                                                                                                                                                                cACTTATTTAAGTGATATTGGTTGTCTGGAGATTCAGGGGGCCAGTCTAAACTTGCTCTTTTCTTCTGGTGTTTATATTTGCAAGATTTGTAGGATACATAACATTAAGTAAAGCCGCATAGGCTTCTCTGTAATGTTTATAtttt                                                                                                                                                >  1:296777/1‑146 (MQ=37)                                                                                                                                                                                    gTTGTCTGGAGATTCAGGGGGCCAGTCTAAACTTGCTCTTTTCTTCTGGTGTTTATCTTTGAAAGATTTGTAGGATAAATAACATTAAGTAAAGCCGCATCGGcn                                                                                                                                                                     <  1:316593/105‑2 (MQ=255)                                                                                                                                                                                    gTTGTCTGGAGATTCAGGGGGCCAGTCTAAACTTGCTCTTTTCTTCTGGTGTTTATCTTTGAAAGATTTGTAGGATAAATAACATTAAGTAAAGCCGCATCGGCt                                                                                                                                                                     >  2:316593/1‑105 (MQ=255)                                                                                                                                                                                          tNGANATTCANGGGGNCAGTCTAAACTTGCTNTTTGCTTCAGGTGTTTATCNNTNAAAGANTTGTa                                                                                                                                                                                                      <  2:211294/66‑1 (MQ=255)                                                                                                                                                                                            gagaTTCAGGGGGCCAGTCTAAACTTGCTCTTTTCTTCTGgtg                                                                                                                                                                                                                           <  1:492160/43‑1 (MQ=255)                                                                                                                                                                                                     ggggCCAGTCTAATCTTGCTCTTTTTTTCTTGTGTTTAGCTTTGAAAGATTTGTAGGATAAATAACATTAAGTAAAGCCGCATCGGCTTCTCTGTAATGTTTATCTTTTTCCAGATAACCTGGGGAGTCGCCCCGTAAATCAAGCTGTAACGTTATTCTTTCATAACCTTTATCTATGATGTGTTTCAACCTCTTTGTTGAAAGAGTTTCATCCAGTTCATTGCCTTTAGCTTCGCGTAGTAAATCAgaga  <  1:280451/251‑1 (MQ=255)                                                                                                                                                                                                         ccAGTCTAAACTTGCTCTTTTCTTCTGGTGTTTATCTTTGAAAGATTTGTAGGATAAATAACAGTAAGTAAAGCCGCATCGGCTTCTCTGTAATGTTTATCTTTTTCCAGATAACCTGGGGATTCTCCCCCTAAATCAAGCTGTAACGGTAGTCTGTCATAACCTTTATCTATGAGGTGGTTCAACCTCGTTGTTGAAAGAGTTTCATcc                                       <  1:380098/210‑1 (MQ=255)                                                                                                                                                                                                            ntCTAAACTTGCTCTTTTCTTCTGGTGTTTATCTTTGAAAGCTTTGTAGGATAAATAACATTAAACAAAGCCGCATCGGCTTCTCTGTAATATTTATCTTTTTCCAGATAACCTGGGGATTCTCCCCCTAACTCAAGCTGTAACGTTATTCTTTCATAACCTTTATc                                                                               >  1:55328/2‑167 (MQ=255)                                                                                                                                                                                                              cTAAACTTGCTCTTTTCTTCTGGTGTTTATCTTTGAAAGATTTGTAGGATAAATAACATTAAGTAAAGCCGCATCGGCTTCTCTGTAATGTTTATCTTTTTCCAGATAACCTGGGGATTCTCCCCCTAAATCAAGCTGTAACGTTATTCTTTCATAACCTTTATCTATGAGGTGGTTCAACCTc                                                            <  1:190608/184‑1 (MQ=255)                                                                                                                                                                                                                |                                                                                                                                                                                                                                                ATGTTGCCGGGCAACGAGGGAGACCGTCATCCCCGGTTCAAAGCTCTGCTGAACAATTGCGATCTTTTCCTGTGTGGTACGCCGTCTGCGTTTCTCCGGCCCTAAGACATCAATCATCTGTTCTCCAATGACTAGTCTAAAAACTAGTATTAAGACTATCACTTATTTAAGTGATATTGGTTGTCTGGAGATTCAGGGGGCCAGTCTAAACTTGCTCTTTTCTTCTGGTGTTTATCTTTGAAAGATTTGTAGGATAAATAACATTAAGTAAAGCCGCATCGGCTTCTCTGTAATGTTTATCTTTTTCCAGATAACCTGGGGATTCTCCCCCTAAATCAAGCTGTAACGTTATTCTTTCATAACCTTTATCTATGAGGTGGTTCAACCTCGTTGTTGAAAGAGTTTCATCCAGTTCATTGACTTTAGCTTCGCGTAGTAAATCAGAGA  >  NC_000913/2997484‑2997930

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