Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
JC	NC_002947	1,851,535	+TAC	coding (579/2241 nt)	relA →	ATP:GTP 3'‑pyrophosphotransferase

New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_002947	= 1851531	0 (0.000)	17 (1.330) +CTA	18/260	0.1	100%	coding (575/2241 nt)	relA	ATP:GTP 3'‑pyrophosphotransferase
	?	NC_002947	1851532 =	0 (0.000)					coding (576/2241 nt)	relA	ATP:GTP 3'‑pyrophosphotransferase

TAGAGAACCTGCGCAAGATGCTGGTGGCCATGGTCGACGACGTGCGCGTGGCGCTGATCAAGCTGGCCGAACGCACGTGCGCAATCCGGGCGGTCAAGTCCGCCGATGACGAGAAACGCCTGCGTGTCGCGCGTGAAGTGTTCGACATCTACGCGCCGCTCGCGCACCGCCTGGGTATCGATCACATCAAGTGGGAGCTGGAAGACCTGTCCTTCCGCTACCTGGAGCCCGACCAGTACAAGCAGATCGCCAAGCTGTTGCATGAGCGGCGGCTGGACCGCGAGCGCTTCA  >  NC_002947/1851384‑1851674                                                                                                                                                     |                                                                                                                                               agaacctgcgcaagatgctggtggccatggtcgacgacgtgcgcgtggcgctgatcaagctggccgaacgcacgtgcgcaatccgggcggtcaagtccgccgatgacgagaaacgcctgcgtgtcgcgcgtgaagtgttcgacatctact                                                                                                                                               <  1:245214‑M2/2‑1 (MQ=255) agaacctgcgcaagatgctggtggccatggtcgacgacgtgcgcgtggcgctgatcaagctggccgaacgcacgtgcgcaatccgggcggtcaagtccgccgatgacgagaaacgcctgcgtgtcgcgcgtgaagtgttcgacatctact                                                                                                                                               >  2:28055‑M2/149‑150 (MQ=255)                         ccatggtcgacgacgtgcgcgtggcgctgatcaagctggccgaacgcacgtgcgcaatccgggcggtcaagtccgccgatgacgagaaacgcctgcgtgtcgcgcgtgaagtgttcgacatctaCTACGCGCCGCTCGCGCACCGCCTgg                                                                                                                       >  1:287916‑M2/125‑150 (MQ=255)                           atggtcgacgacgtgcgcgtggcgctgatcaagctggccgaacgcacgtgcgcaatccgggcggtcaagtccgccgatgacgagaaacgcctgcgtgtcgcgcgggaagtgttcgacatctaCTACGCGCCGCTCGCGCACCGCCTGGGt                                                                                                                     <  2:117446‑M2/28‑1 (MQ=255)                                 gacgacgtgcgcgtggcgctgatcaagctggccgaacgcacgtgcgcaatccgggcggtcaagtccgccgatgacgagaaacgcctgcgtgtcgcgcgtgaagtgttcgacatctaCTACGCGCCGCTCGCGCACCGCCTGGGTATCGAt                                                                                                               >  2:193371‑M2/117‑150 (MQ=255)                                  acgacgtgcgcgtggcgctgatcaagctggccgaacgcacgtgcgcaatccgggcggtcaagtccgccgatgacgagaaacgcctgcgtgtcgcgcgtgaagtgttcgacatctaCTACGCGCCGCTCGCGCACCGCCTGGGTATCGATc                                                                                                              >  2:257554‑M2/116‑150 (MQ=255)                                   cgacgtgcgcgtggcgctgatcaagctggccgaacgcacgtgcgcaatccgggcggtcaagtccgccgatgacgagaaacgcctgcgtgtcgcgcgtgaagtgttcgacatctaCTACGCGCCGCTCGCGCACCGCCTGGGTATCGATca                                                                                                             <  2:287831‑M2/36‑1 (MQ=255)                                    gacgtgcgcgtggcgctgatcaagctggtcgaacgcacgtgcgcaatccgggcggtcaagtccgccgatgacgagaaacgcctgcgtgtcgcgcgtgacgtgttcgacatctaCTACGCGCCGCTCGCGCCCCGCCTGGGTATCGTTcac                                                                                                            >  2:85327‑M2/114‑150 (MQ=255)                                            cgtggcgctgatcaagctggccgaacgcacgtgcgcaatccgggcggtcaagtccgccgatgacgagaaacgcctgcgtgtcgcgcgtgaagtgttcgacatctaCTACGCGCCGCTCGCGCACCGCCTGGGTATCGATCACATCAAGTg                                                                                                    <  1:257554‑M2/45‑1 (MQ=255)                                                   atgatcaagctggccgaacgcacgtgcgcaatccgggcggtcaagtccgccgatgacgagaaacgcctgcgtgtcgcgcgtgaagtgttcgacatctaCTACGCGCCGCTCGCGCACCGCCTGGGTATCGATCACATCAAGTGGGAGCTg                                                                                             <  1:27851‑M2/52‑1 (MQ=255)                                                                           tgcgcaatccgggcggtcaagtccgccgatgacgagaaacgcctgcgtgtcgcgcgtgaagtgttcgacatctaCTACGCGCCGCTCGCGCACCGCCTGGGTATCGATCACATCAAGTGGGAGCTGGAAGACCTGTCCTTCCGCTACCTg                                                                     >  1:235429‑M2/75‑150 (MQ=255)                                                                             cgcaatccgggcggtcaagtccgccgatgacgagaaacgcctgcgtgtcgcgcgtgaagtgttcgacatctaCTACGCGCCGCTCGCGCACCGCCTGGGTATCGATCACATCAAGTGGGAGCTGGAAGACCTGTCCTTCCGCTACCTGGa                                                                   <  2:235429‑M2/78‑1 (MQ=255)                                                                                    cgggcggtcaagtccgccgatgacgagaaacgcctgcgtgtcgcgcgtgaagtgttcgacatctaCTACGCGCCGCTCGCGCACCGCCTGGGTATCGATCACATCAAGTGGGAGCTGGAAGACCTGTCCTTCCGCTACCTGGAGCCCGAc                                                            <  1:28055‑M2/85‑1 (MQ=255)                                                                                                              gaaacgcctgcgtgtcgcgcgtgaagtgttcgacatctaCTACGCGCCGCTCGCGCACCGCCTGGGTATCGATCACATCAAGTGGGAGCTGGAAGACCTGTCCTTCCGCTACCTGGAGCCCGACCAGTACAAGCAGATCGCCAAGCTGtt                                  >  1:92744‑M2/40‑150 (MQ=255)                                                                                                                  cgcctgcgtgtcgcgcgtgaagtgttcgacatctaCTACGCGCCGCTCGCGCACCGCCTGGGTATCGATCACATCAAGTGGGAGCTGGAAGACCTGTCCTTCCGCTACCTGGAGCCCGACCAGTACAAGCAGATCGCCAAGCTGTTGCAt                              <  1:193371‑M2/115‑1 (MQ=255)                                                                                                                              gcgcgtgaagtgttcgacatctaCTACGCGCCGCTCGCGCACCGc                                                                                                                           >  2:90586‑M2/24‑45 (MQ=255)                                                                                                                              gcgcgtgaagtgttcgacatctaCTACGCGCCGCTCGCGCACCGc                                                                                                                           <  1:90586‑M2/22‑1 (MQ=255)                                                                                                                                              acatctaCTACGCGCCGCTCGCGCACCGCCTGGGTATCGATCACATCAAGTGGGAGCTGGAAGACCTGTCCTTCCGCTACCTGGAGCCCGACCAGTACAAGCAGATCGCCAAGCTGTTGCATGAGCGGCGGCTGGACCGCGAGCGCTtca  <  1:226982‑M2/143‑1 (MQ=255)                                                                                                                                                     |                                                                                                                                               TAGAGAACCTGCGCAAGATGCTGGTGGCCATGGTCGACGACGTGCGCGTGGCGCTGATCAAGCTGGCCGAACGCACGTGCGCAATCCGGGCGGTCAAGTCCGCCGATGACGAGAAACGCCTGCGTGTCGCGCGTGAAGTGTTCGACATCTACGCGCCGCTCGCGCACCGCCTGGGTATCGATCACATCAAGTGGGAGCTGGAAGACCTGTCCTTCCGCTACCTGGAGCCCGACCAGTACAAGCAGATCGCCAAGCTGTTGCATGAGCGGCGGCTGGACCGCGAGCGCTTCA  >  NC_002947/1851384‑1851674

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