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)	34 (1.070)	30/492	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)	55 (1.730)	52/492	0.0	100%	coding (256/2235 nt)	relA	GDP/GTP pyrophosphokinase
	?	NC_000913	2996361 =	NA (NA)					noncoding (1331/1331 nt)	IS2	repeat region

ATTGCCAATGGTTGTTGCCTTATATAGGGTAACAACGTTGAAAAACCACATGTGATGATATTGTGATAGTTAGCGATGCCGGGTTTGTATCGTGCCTTTTCATTTATTTCCAGCGGTATAGCTGATAAATAATCCCTTAGTACATCAGGAACCTTATAAAAAAACCCGC‑GTTCTGCAAACCTTTCTCGATAAAAGTGACAACCACTGAGTAAAC‑AGATCAATGGAAGTAAATTCCAA‑TTGTTTATTGGATTT‑GCCCCTATATTTCCAGACATCTGTTATCACTTAACCCATTACAAGCCCGCTGCCGCAGATATTCCCGTGGCGAGCGATAACCCAGCGCACTATGCGGATGCCATTCGTTATAATGCTCGAACGCCTCTGCAAGGTTCTTTGCTGCCGTTAACCCGTCTGGTTTGGGCATGATACTGATGTAGTCACGCTTTATCGTTTTCACGAAGCTCTCTGCTATTCCGTTACTCTCCGGACTCCGCAC  >  NC_000913/2996116‑2996608                                                                                                                                                                                                                                                         |                                                                                                                                                                                                                                                         aTTGCCAATGGTTGTTGCCTTATATAGGGTAACAACGTTGAAAAACCACATGTGATGATATTGTGATAGTTAGCGATGCCGGGTTTGTATCGTGCCCTTTCATTTATTTCCAGCGGTATAGCTGATAAATAATCCCTTAGTACATCAGGAACCTTATAAAAAAACCCGC‑GTTCTGCAAACCTTTCTCGATAAAAGTGACAACCACTGAGTAAAC‑AGATCAATGGAAGTAAATTCCAA‑TTGTTTATTGGAtt                                                                                                                                                                                                                                                     >  1:129371/1‑251 (MQ=255)                     tatataGGGTAACAACGTTGAAAAACCACATGTGATGATATTGTGATAGTTAGCGATGCCGGGTTTGTATCGTGCCTTTTCATTTATTTCCAGCGGTATAGCTGATAAATAATCCCTTAGTACATCAGGAACCTTATAAAAAAACCCGC‑GTTCTGCAAACCTTTCTCGATAAAAGTGACAACCACTGAGTAAAC‑AGATCAATGGAAGTAAATTCCAA‑TTGTTTATTGGATTT‑GCCCCTATATTTCCAGAc                                                                                                                                                                                                                                 >  2:134368/1‑250 (MQ=255)                     tatataGGGTAACAACGTTGAAAAACCACATGTGATGATATTGTGATAGTTAGCGATGCCGGGATTGTATCGTGCCTTTTCATTTATTTCCAGCGGTATAGCTGATAAATAATCCCTTAGTACATCAGGAACCTTATAAAAAAACCCGC‑GTTCTGCAAACCTTTCTCGATAAAAGTGACAACCACTGAGTAAAC‑AGATCAATGGAAGTAAATTCCAA‑TTGTTTATTGGATTT‑GCCCCTATATTTCCAGAc                                                                                                                                                                                                                                 >  2:175052/1‑250 (MQ=255)                                     gTTGAAAAACCACATGTGATGATATTGTGATAGTTAGCGATGCCGGGTTTGTATCGTGCCTTTTCATTTATTTCCAGCGGTATAGCTGATAAATAATCCCTTAGTACATCAGGAACCTTATAAAAAAACCCGC‑GTTCTGCAAACCTTTCTCGATAAAAGTGACAACCACTGAGTAAAC‑AGATCAATGGAAGTAAATTCCAA‑TTGTTTATTGGATTT‑GCCCCTATATTTCCAGACATCTGTTATCACTTAAc                                                                                                                                                                                                                >  2:281360/1‑251 (MQ=255)                                                      gatgatATTGTGATAGTTAGCGATGCCGGGTTTGTATCGTGCCTTTTCATTTATTTCCAGCGGTATAGCTGATAAATAATCCCTTAGTACATCAGGAACCTTATAAAAAAACCCGC‑GTTCTGCAAACCTTTCTCGATAAAAGTGACAACCACTGAGTAAAC‑AGATCAATGGAAGTAAATTCCAA‑TTGTTTATTGGATTT‑GCCCCTATATTTCCAGACATCTGTTATCACTTAACCCATTACAAGCCCGCt                                                                                                                                                                                                <  1:302776/250‑1 (MQ=255)                                                      gatgatATTGTGATAGTTAGCGATGCCGGGTTTGTATCGTGCCTTTTCATTTATTTCCAGCGGTATAGCTGATAAATAATCCCTTAGTACATCAGGAACCTTATAAAAAAACCCGC‑GTTCTGCAAACCTTTCTCGATAAAAGTGACAACCACTGAGTAAAC‑AGATCAATGGAAGTAAATTCCAA‑TTGTTTATTGGATTT‑GCCCCTATATTTCCAGACATCTGTTATCACTTAACCCATTACAAGCCCGCTg                                                                                                                                                                                               >  2:239934/1‑251 (MQ=255)                                                       atgatATTGTGATAGTTAGCGATGCCGGGTTTGTATCGTGCCTTTTCATTTATTTCCAGCGGTATAGCTGATAAATAATCCCTTAGTACATCAGGAACCTTATAAAAAAACCCGC‑GTTCTGCAAACCTTTCTCGATAAAAGTGACAACCACTGAGTAAAC‑AGATCAATGGAAGTAAATTCCAA‑TTGTTTATTGGATTT‑GCCCCTATATTTCCAGACATCTGTTATCACTTAACCCATTACAAGCCCGCTGc                                                                                                                                                                                              >  1:277389/1‑251 (MQ=255)                                                                    aGTTAGCGATGCCGGGTTTGTATCGGGCCGTTTCATTTATTTCCAGCGGTATAGCTGATAAATAATCCCTTAGTACATCAGGAACCTTATAAAAAAACCCGC‑GTTCTGCAAACCTTTCTCGATAAAAGTGACAACCACTGAGTAAAC‑AGATCAATGGAAGTAAATTCCAA‑TTGTTTATTGGATTT‑GCCCCTATATTTCCAGACATCTGTTATCACTTAACCCATTACAAGCCCGCTGCCGCAGATATTccc                                                                                                                                                                                 <  2:314870/251‑1 (MQ=255)                                                                      ttggcgaggccggtttgtttccttGCCTTTCCATTTATTCCCACCGGTTTAGTTATTAAATAACCCCTAAGT‑CAACCAGAACCCTATAAAAAAAACCGCCTTTCTGCAAACCTTTTTGGATAAAA‑TGGCAAACAATTAATAAACCAGATCAATG‑AAATAAAATCCAAATTGTT‑ATTTGATTTTGCCCCTA‑AATTCCAGACAACCGGTATCACTTAACCCATTACAAGCCCGCTGCCGCAGATATTCCCGt                                                                                                                                                                               <  2:313619/227‑1 (MQ=255)                                                                                                cTTTTCATTTATTTCCAGCGGTATAGCTGATAAATAATCCCTTAGTACATCAGGAACCTTATAAAAAAACCCGC‑GTTCTGCAAACCTTTCTCGATAAAAGTGACAACCACTGAGTAAAC‑AGATCAATGGAAGTAAATTCCAA‑TTGTTTATTGGATTT‑GCCCCTATATTTCCAGACATCTGTTATCACTTAACCCATTACAAGCCCGCTGCCGCAGATATTCCCGTGGCGAGCGATAACCCAGCGCACTATg                                                                                                                                                     >  2:93876/1‑251 (MQ=255)                                                                                                          atttCCAGCGGTATAGCTGATAAATAATCCCTTAGTACATCAGGAACCTTATAAAAAAACCCGC‑GTTCTGCAAACCTTTCTCGATAAAAGTGACAACCACTGAGTAAAC‑AGATCAATGGAAGTAAATTCCAA‑TTGTTTATTGGATTT‑GCCCCTATATTTCCAGACATCTGTTATCACTTAACCCATTACAAGCCCGCTGCCGCAGATATTCCCGTGGCGAGCGATAACCCAGCGCACTATGCGGATGCCa                                                                                                                                            <  1:93876/250‑1 (MQ=255)                                                                                                                      ataGCTGATAAATAATCCCTTAGTACATCAGGAACCTTATAAAAAAACCCGC‑GTTCTGCAAACCTTTCTCGATAAAAGTGACAACCACTGAGTAAAC‑AGATCAATGGAAGTAAATTCCAA‑TTGTTTATTGGATTT‑GCCCCTATATTTCCAGACATCTGTTATCACTTAACCCATTACAAGCCCGCTGCCGCAGATATTCCCGTGGCGAGCGATAACCCAGCGCACTATGCGGATGCCATTCGTTATAATGc                                                                                                                               >  1:80073/1‑251 (MQ=255)                                                                                                                                          tAGTACATCAGGAACCTTATAAAAAAACCCGC‑GTTCTGCAAACCTTTCTCGATAAAAGTGACAACCACTGAGTAAAC‑AGATCAATGGAAGTAAATTCCAA‑TTGTTTATTGGATTT‑GCCCCTATATTTCCAGACATCTGTTATCACTTAACCCATTACAAGCCCGCTGCCGCAGATATTCCCGTGGCGAGCGATAACCCAGCGCACTATGCGGATGCCATTCGTTATAATGCTCGAACGCCTCTGCAAGGtt                                                                                                           >  1:106312/1‑251 (MQ=255)                                                                                                                                                  cAGGAACCTTATAAAAAAACCCGC‑GTTCTGCAAACCTTTCTCGATAAAAGTGACAACCACTGAGTAAAC‑AGATCAATGGAAGTAAATTCCAA‑TTGTTTATTGGATTT‑GCCCCTATATTTCCAGACATCTGTTATCACTTAACCCATTACAAGCCCGCTGCCGCAGATATTCCCGTGGCGAGCGATAACCCAGCGCACTATGCGGATGCCATTCGTTATAATGCTCGAACGCCTCTGCACGGTTCTTtgctg                                                                                                   >  2:45702/1‑251 (MQ=255)                                                                                                                                                    ggAACCTTATAAAAAAACCCGC‑GTTCTGCAAACCTTTCTCGATAAAAGTGACAACCACTGAGTAAAC‑AGATCAATGGAAGTAAATTCCAA‑TTGTTTATTGGATTT‑GCCCCTATATTTCCAGACATCTGTTATCACTTAACCCATTACAAGCCCGCTGCCGCAGATATTCCCGTGGCGAGCGATAACCCAGCGCACTATGCGGATGCCATTCGTTATAATGCTCGAACGCCTCTGCAAGGTTCTTTGCTGcc                                                                                                 <  1:238591/251‑1 (MQ=255)                                                                                                                                                           tataAAAAAACCCGC‑GTTCTGCAAACCTTTCTCGATAAAAGTGACAACCACTGAGTAAAC‑AGATCAATGGAAGTAAATTCCAA‑TTGTTTATTGGATTT‑GCCCCTATATTTCCAGACATCTGTTATCACTTAACCCATTACAAGCCCGCTGCCGCAGATATTCCCGTGGCGAGCGATAACCCAGCGCACTATGCGGATGCCATTCGTTATAATGCTCGAACTCCTCTGCAAGGTTCTTTGCTGCCGTTAAc                                                                                           <  1:268552/250‑1 (MQ=255)                                                                                                                                                           tataAAAAAACCCGC‑GTTCTGCAAACCTTTCTCGATAAAAGTGACAACCACTGAGTAAAC‑AGATCAATGGAAGTAAATTCCAA‑TTGTTTATTGGATTT‑GCCCCTATATTTCCAGACATCTGTTATCACTTAACCCATTACAAGCCCGCTGCCGCAGATATTCCCGTGGCGAGCGATAACCCAGCGCACTATGCGGATGCCATTCGTTATAATGCTCGAACGCCTCTGCAAGGTTCTTTGCTGCCGTTAAcc                                                                                          >  1:68688/1‑251 (MQ=255)                                                                                                                                                                    aCCCGC‑GTTCTGCAAACCTTTCTCGATAAAAGTGACAACCACTGAGTAAAC‑AGATCAATGGAAGTAAATTCCAA‑TTGTTTATTGGATTT‑GCCCCTATATTTCCAGACATCTGTTATCACTTAACCCATTACAAGCCCGCTGCCGCAGATATTCCCGTGGCGAGCGATAACCCAGCGCACTATGCGGATGCCATTCGTTATAATGCTCGAACGCCTCTGCAAGGTTCTTTGCTGCCGTTAACCCGTCTGGtt                                                                                 >  1:269810/1‑251 (MQ=255)                                                                                                                                                                    aCCCGC‑GTTCTGCAAACCTTTCTCGATAAAAGTGACAACCACTGAGTAAAC‑AGATCAATGGAAGTAAATTCCAA‑TTGTTTATTGGATTT‑GCCCCTATATTTCCAGACATCTGTTATCACTTAACCCATTACAAGCCCGCTGCCGCAGATATTCCCGTGGCGAGCGATAACCCAGCGCACTATGCGGATGCCATTCGTTATAATGCTCGAACGCCTCTGCAAGGTTCTTTGCTGCCGTTAACCCGTCTGGtt                                                                                 >  1:203307/1‑251 (MQ=255)                                                                                                                                                                           gTTCTGCAAACCTTTCTCGATAAAAGTGACAACCACTGAGTAAAC‑AGATCAATGGAAGTAAATTCCAA‑TTGTTTATTGGATTT‑GCCCCTATATTTCCAGACATCTGTTATCACTTAACCCATTACAAGCCCGCTGCCGCAGATATTccc                                                                                                                                                                                 <  1:319633/149‑1 (MQ=255)                                                                                                                                                                           gTTCTGCAAACCTTTCTCGATAAAAGTGACAACCACTGAGTAAAC‑AGATCAATGGAAGTAAATTCCAA‑TTGTTTATTGGATTT‑GCCCCTATATTTCCAGACATCTGTTATCACTTAACCCATTACAAGCCCGCTGCCGCAGATATTccc                                                                                                                                                                                 >  2:319633/1‑149 (MQ=255)                                                                                                                                                                                      cTTTCTCGATAAACGTGACAACCACTGAGTAAAC‑AGCTCAATGGAAGTAAATTCCAA‑TTGTTTATTGGATTT‑GCCCCTATATTTCCAGACATCTGTTATCACTTAACCCATTACAAGCCCGCTGCCGCAGATATTCCCGTGGCGAGCGATAACCCAGCGCACTATGCGGATGCCATTCGTTATAATGCTCGAACGCCTCTGCAAGGTTCTTTGCTGCCGTTAACCCGTCTGGTTTGGGCATGATACTGATg                                                                >  1:260384/1‑251 (MQ=255)                                                                                                                                                                                                 aaaGTGACAACCACTGATTAAAC‑AGATCAATGGAAGTAAATTCCAA‑TTGTTTATTGGATTT‑GCCCCTATATTTCCAGACATCTGTTATCACTTAACCCATTACAAGCCCGCTGCCGCAGATATTCCCGTGGCGAGCGATAACCCAGCGCACTATGCGGATGCCATTCGTTATAATGCTCGAACGCCTCTGCAAGGTTCTTTGCTGCCGTTAACCCGTCTGGTTTGGGCATGATGCTGATTTAGGCACGCtt                                                     <  2:310330/251‑1 (MQ=255)                                                                                                                                                                                                  aaGTGACAACCACTGAGTAAAC‑AGATCAATGGAAGTAAATTCCAA‑TTGTTTATTGGATTT‑GCCCCTATATTTCCAGACATCTGTTATCACTTAACCCATTACAAGCCCGCTGCCGCAGATATTCCCGTGGCGAGCGATAACCCAGCGCACTATGCGGATGCCATTCGTTATAATGCTCGAACGCCTCTGCAAGGTTCTTTGCTGCCGTTAACCCGTCTGGTTTGGGCATGATACTGATGTAGTCACGCttt                                                    >  1:10964/1‑251 (MQ=255)                                                                                                                                                                                                             aCTGAGTAAAC‑AGATCACTGGAAGTAAATTCCAA‑TTGTTTATTGGATTT‑GCCCCTATATTTCCAGACATCTGTTATCACTTAACCCATTACAAGCCCGCTGCCGCAGATATTCCCGTGGCGAGCGATAACCCAGCGCACTATGCGGATGCCATTCGTTATAATGCTCGAACGCCTCTGCAAGGTTCTTTGCTGCCGTTAACCCGTCTGGTTTGGGCATGATACTGATGTAGTCACGCTTTATCGTTTTCAc                                         >  2:214612/1‑251 (MQ=255)                                                                                                                                                                                                                gAGTAAAC‑AGATCAATGGAAGTAAATTCCAA‑TTGTTTATTGGATTT‑GCCCCTATATTTCCAGACATCTGTTATCACTTAACCCATTACAAGCCCGCTGCCGCAGATATTCCCGTGGCGAGCGATAACCCAGCGCACTATGCGGATGCCATTCGTTATAATGCTCGAACGCCTCTGCAAGGTTCTTTGCTGCCGTTAACCCGTCTGGTTTGGGCATGATACTGATGTAGTCACGCTTTATCGTTTTCACGaa                                      <  1:71007/251‑1 (MQ=255)                                                                                                                                                                                                                         aGATCAATGGAAGTAAATTCCAA‑TTGTTTATTGGATTT‑ACCCCTATATTTCCAGACATTTGTTATCACTTAACCCATTACACGCCCGATGCCGCAGATATTCCCGTGGCGAGCGATAACACAGCGCACTATGCGGATGCCATTCGTTATAATGCTCGAACGCCTCTGCAAGGTTCTTTGCTGCCGTTAACCCGTCTGGTTTCGGCATGATACTGATGTAGTCACGCTTTATCGTTTTCACGAAGCTCTCTg                              >  1:18089/1‑251 (MQ=255)                                                                                                                                                                                                                                     gTAAATTCCAA‑TTGTTTATTGGATTT‑GCCCCTATATTTCCAGACATCTGTTATCACTTAACCCATTACAAGCCCGCTGCCGCAGATATTCCCGTGGCGAGCGATAACCCAGCGCACTATGCGGATGCCATTCGTTATAATGCTCGAACGCCTCTGCAAGGTTCTTTGCTGCCGTTAACCCGTCTGGTTTGGGCATGATACTGATGTAGTCACGCTTTATCGTTTTCACGAAGCTCTCTGCTATTCCGTTAc                  >  2:311675/1‑251 (MQ=37)                                                                                                                                                                                                                                                      tATTGGATTT‑GCCCCTATATTTCCAGACATCTGTTATCACTTAACCCATTACTAGCCCGCTGCCGCAGATATTCCCGTGGCGAGCGATAACCAAGAGCACTATGCGGATGCCATTCGTTATCATGCTCGCACGCCTCTGAAAGGTTCTTTGCTGCCGTTAACCCGTCTGGTTTGGGCATGATACTGATGTAGCCACGCTTTATCGTTTTCACGCAGCTCTCTGCTCTTCCGTTACTCTCCGGACTCCGCAc  >  2:18502/1‑251 (MQ=11)                                                                                                                                                                                                                                                         |                                                                                                                                                                                                                                                         ATTGCCAATGGTTGTTGCCTTATATAGGGTAACAACGTTGAAAAACCACATGTGATGATATTGTGATAGTTAGCGATGCCGGGTTTGTATCGTGCCTTTTCATTTATTTCCAGCGGTATAGCTGATAAATAATCCCTTAGTACATCAGGAACCTTATAAAAAAACCCGC‑GTTCTGCAAACCTTTCTCGATAAAAGTGACAACCACTGAGTAAAC‑AGATCAATGGAAGTAAATTCCAA‑TTGTTTATTGGATTT‑GCCCCTATATTTCCAGACATCTGTTATCACTTAACCCATTACAAGCCCGCTGCCGCAGATATTCCCGTGGCGAGCGATAACCCAGCGCACTATGCGGATGCCATTCGTTATAATGCTCGAACGCCTCTGCAAGGTTCTTTGCTGCCGTTAACCCGTCTGGTTTGGGCATGATACTGATGTAGTCACGCTTTATCGTTTTCACGAAGCTCTCTGCTATTCCGTTACTCTCCGGACTCCGCAC  >  NC_000913/2996116‑2996608

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