Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
MC JC	NC_000913	2,912,618	Δ10 bp	coding (1025‑1034/2235 nt)	relA ←	(p)ppGpp synthetase I/GTP pyrophosphokinase

Missing coverage evidence...
			seq id	start	end	size	←reads	reads→	gene	description
*	*	÷	NC_000913	2912618	2912627	10	53 [0]	[0] 54	relA	(p)ppGpp synthetase I/GTP pyrophosphokinase

New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913	= 2912617	0 (0.000)	53 (1.080)	46/284	0.2	100%	coding (1035/2235 nt)	relA	(p)ppGpp synthetase I/GTP pyrophosphokinase
	?	NC_000913	2912628 =	0 (0.000)					coding (1024/2235 nt)	relA	(p)ppGpp synthetase I/GTP pyrophosphokinase

GGAATCAGCCATCTCTTCCTGCCACGCAATCAGTTTACGCAGCCAGGCAATCCGGTCTTCATGTCCCGAACGTGCGCCGCCAGCAGCCGCGCCCTCTTTATATTTCCAGTGCGCAGCAACACCCAACTCTGCATCTTCATGCATCTGTTTGGTGCGGATTTGGATCTCAACGGTTTTTCCACCCGGCCCCAGAACCACGGTATGAATAGACTGATAACCGTTTGGTTTCGGGTTAGCGACGTAATCGTCAAACTCATCCGGCAGGTGGCGATAGTGAGTGTGCACTATCCCC  >  NC_000913/2912482‑2912773                                                                                                                                                   |                                                                                                                                                  gcatttcgccggaatcagccatctcttcctgccacgcaatcagtttacgcagccaggcaatccggtcttcatgtcccgaacgtgcgccgccagcagccgcgccctctttatatttccagtgcgcagcaacacccaactctgcatctGTTTg                                                                                                                                               >  2:118220‑M2/147‑151 (MQ=255)         ccggaatcagccatctcttcctgccacgcaatcagtttacgcagccaggcaatccggtcttcatgtcccgaacgtgcgccgccagcagccgcgccctctttatatttccagtgcgcagcaacacccaactctgcatctGTTTGGTGCGGAt                                                                                                                                       <  1:701054‑M2/13‑1 (MQ=255)          cggaatcagccatctcttcctgccacgcaatcagtttacgcagccaggcaatccggtcttcatgtcccgaacgtgcgccgccagcagccgcgccctctttatatttccagtgcgcagcaacacccaactctgcatctGTTTGGTGCGGAtt                                                                                                                                      >  1:593137‑M2/138‑151 (MQ=255)            gaatcagccatctcttcctgccacgcaatcagtttacgcagccaggcaatccggtcttcatgtcccgaacgtgcgccgccagcagccgcgccctctttatatttccagtgcgcagcaacacccaactctgcatctGTTTGGTGCGGATTTg                                                                                                                                    <  1:210779‑M2/16‑1 (MQ=255)              atcagccatctcttcctgccacgcaatcagtttacgcagccaggcaatccggtcttcatgtcccgaacgtgcgccgccagcagccgcgccctctttatatttccagtgcgcagcaacacccaactctgcatctGTTTGGTGCGgatttcg                                                                                                                                   >  1:710164‑M2/134‑148 (MQ=255)              atcagccatctcttcctgccacgcaatcagtttacgcagccaggcaatccggtcttcatgtcccgaacgtgcgccgccagcagccgcgccctctttatatttccagtgcgcagcaacacccaactctgcatctGTTTGGTGCGGATTTgg                                                                                                                                   >  2:459027‑M2/134‑150 (MQ=255)              atcagccatctcttcctgccacgcaatcagtttacgcagccaggcaatccggtcttcatgtcccgaacgtgcgccgccagcagccgcgccctctttatatttccagtgcgcagcaacacccaactctgcatctGTTTGGTGCGGATTTGGa                                                                                                                                  <  1:749813‑M2/18‑1 (MQ=255)                   ccatctcttcctgccacgcaatcagtttacgcagccaggcaatccggtcttcatgtcccgaacgtgcgccgccagcagccgcgccctctttatatttccagtgcgcagcaacacccaactctgcatctGTTTGGTGCGGATTTGGATCTCa                                                                                                                             <  1:566232‑M2/23‑1 (MQ=255)                      tctcttcctgccacgcaatcagtttacgcagccaggcaatccggtcttcatgtcccgaacgtgcgccgccagcagccgcgccctctttatatttccagtgcgcagcaacacccaactctgcatctGTTTGGTGCGGATTTGGATCTCAACg                                                                                                                          >  2:576179‑M2/126‑151 (MQ=255)                           tcctgccacgcaatcagtttacgcagccaggcaatccggtcttcatgtcccgaacgggcgccgccagcagccgcgccctctttatatttccagtgcgcagcaacacccaactctgcatctGTTTGGTGCGGATTTGGATCTCAACGGtt                                                                                                                       <  2:810394‑M2/29‑1 (MQ=255)                               gccacgcaatcagtttacgcagccaggcaatccggtcttcatgtcccgaacgtgcgccgccagcagccgcgccctctttatatttccagtgcgcagcaacacccaactctgcatctGTTTGGTGCGGATTTGGATCTCAACGGTTTTTc                                                                                                                   <  2:131700‑M2/33‑1 (MQ=255)                               gccacgcaatcagtttacgcagccaggcaatccggtcttcatgtcccgaacgtgcgccgccagcagccgcgccctctttatatttccagtgcgcagcaacacccaactctgcatctGTTTGGTGCGGATTTGGATCTCAACGGTTTTTc                                                                                                                   <  1:258373‑M2/33‑1 (MQ=255)                               gccacgcaatcagtttacgcagccaggcaatccggtcttcatgtcccgaacgtgcgccgccagcagccgcgccctctttatatttccagtgcgcagcaacacccaactctgcatctGTTTGGTGCGGATTTGGATCTCAACGGTTTTTc                                                                                                                   <  2:593171‑M2/33‑1 (MQ=255)                               gccacgcaatcagtttacgcagccaggcaatccggtcttcatgtcccgaacgtgcgccgccagcagccgcgccctctttatatttccagtgcgcagcaacacccaactctgcatctGTTTGGTGCGGATTTGGATCTCAACGGTTTTTCCa                                                                                                                 <  2:468328‑M2/35‑1 (MQ=255)                                 cacgcaatcagtttacgcagccaggcaatccggtcttcatgtcccgaacgtgcgccgccagcagccgcgccctctttatatttccagtgcgcagcaacacccaactctgcatctGTTTGGTGCGGATTTGGATCTCAACGGTTTTTCCAcc                                                                                                               <  1:728258‑M2/37‑1 (MQ=255)                                     caatcagtttacgcagccaggcaatccggtcttcatgtcccgaacgtgcgccgccagcagccgcgccctctttatatttccagtgcgcagcaacacccaactctgcatctGTTTGGTGCGGATTTGGATCTCAACGGTTTTTCCACCCGGc                                                                                                           >  2:35053‑M2/111‑151 (MQ=255)                                      aatcagtttacgcagccaggcaatccggtcttcatgtcccgaacgtgcgccgccagcagccgcgccctctttatatttccagtgcgcagcaacacccaactctgcatctGTTTGGTGCGGATTTGGATCTCAACGGTTTTTCCACCCGGcc                                                                                                          >  2:358878‑M2/110‑151 (MQ=255)                                           gtttacgcagccaggcaatccggtcttcatgtcccgaacgtgcgccgccagcagccgcgccctctttatatttccagtgcgcagcaacacccaactctgcatctGTTTGGTGCGGATTTGGATCTCAACGGTTTTTCCACCCGGcccc                                                                                                        >  1:471126‑M2/105‑148 (MQ=255)                                               acgcagccaggcaatccggtcttcatgtcccgaacgtgcgccgccagcagccgcgccctctttatatttccagtgcgcagcaacacccaactctgcatctGTTTGGTGCGGATTTGGATCTCAACGGTTTTTCCACCCGGCCCCAGAAcc                                                                                                  >  2:205521‑M2/101‑150 (MQ=255)                                                cgcagccaggcaatccggtcttcatgtcccgaacgggcgccgccagcagccgcgccctctttatatttccagtgcgcagcaacacccaactctgcatctGTTTGGTGCGGATTTGGATCTCAACGGTGTTTCCACCCGGCCCCAGAAcc                                                                                                  <  1:16747‑M2/50‑1 (MQ=255)                                                  cagccaggcaatccggtcttcatgtcccgaacgtgcgccgccagcagccgcgccctctttatatttccagtgcgcagcaacacccaactctgcatctGTTTGGTGCGGATTTGGATCTCAACGGTTTTTCCACCCGGCCCCAGAACCACgg                                                                                              <  2:555276‑M2/54‑1 (MQ=255)                                                   agccaggcaatccggtcttcatgtcccgaacgtgcgccgccagcagccgcgccctctttatatttccagtgcgcagcaacacccaactctgcatctGTTTGGTGCGGATTTGGATCTCAACGGTTTTTCCACCCGGCCCCAGAACCACGGt                                                                                             >  1:438834‑M2/97‑151 (MQ=255)                                                     ccaggcaatccggtcttcatgtcccgaacgtgcgccgccagcagccgcgccctctttatatttccagtgcgcagcaacacccaactctgcatctGTTTGGTGCGGATTTGGATCTCAACGGTTTTTCCACCCGGCCCCAGAACCACGGTa                                                                                            >  2:512275‑M2/95‑150 (MQ=255)                                                      caggcaatccggtcttcatgtcccgaacgtgcgccgccagcagccgcgccctctttatatttccagtgcgcagcaacacccaactctgcatctGTTTGGTGCGGATTTGGATCTCAACGGTTTTTCCACCCGGCCCCAGAACCACGGTATg                                                                                          <  1:2725‑M2/58‑1 (MQ=255)                                                           aatccggtcttcatgtcccgaacgtgcgccgccagcagccgcgccctctttatatttccagtgcgcagcaacacccaactctgcatctGTTTGGTGCGGATTTGGATCTCAACGGTTTTTCCACCCGGCCCCAGAACCACGGTATGAAt                                                                                       >  1:11483‑M2/89‑149 (MQ=255)                                                            atccggtcttcatgtcccgaacgtgcgccgccagcagccgcgccctctttatatttccagtgcgcagcaacacccaactctgcatctGTTTGGTGCGGATTTGGATCTCAACGGTTTTTCCACCCGGCCCCAGAACCACGGTATGAAt                                                                                       >  2:390505‑M2/88‑148 (MQ=255)                                                             tctggtcttcatgtcccgaacgtgcgccgccagcagccgcgccctctttatatttccagtgcgcagcaacacccaactctgcatctGTTTGGTGCGGATTTGGATCTCAACGGTTTTTCCACCCGGCCCCAGAACCACGGTATGAATAGAc                                                                                   >  2:312023‑M2/87‑151 (MQ=255)                                                                ggtcttcatgtcccgaacgtgcgccgccagcagccgcgccctctttatatttccagtgcgcagcaacacccaactctgcatctGTTTGGTGCGGATTTGGATCTCAACGGTTTTTCCACCCGGCCCCAGAACCACGGTATGAATAGACTg                                                                                 >  2:138759‑M2/84‑150 (MQ=255)                                                                ggtcttcatgtcccgaacgtgcgccgccagcagccgcgccctctttatatttccagtgcgcagcaacacccaactctgcatctGTTTGGTGCGGATTTGGATCTCAACGGTTTTTCCACCCGGCCCCAGAACCACGGTATGAATAGACTGa                                                                                <  1:390549‑M2/68‑1 (MQ=255)                                                                      catgtcccgaacgtgcgccgccagcagccgcgccctctttatatttccagtgcgcagcaacacccaactctgcatctGTTTGGTGCGGATTTGGATCTCAACGGTTTTTCCACCCGGCCCCAGAACCACGGTATGAATAGACTGATAAcc                                                                           <  2:93657‑M2/73‑1 (MQ=255)                                                                      catgtcccgaacgtgcgccgccagcagccgcgccctctttatatttccagtgcgcagcaacacccaactctgcatctGTTTGGTGCGGATTTGGATCTCAACGGTTTTTCCACCCGGCCCCAGAACCACGGTATGAATAGACTGATAACCg                                                                          <  1:27990‑M2/74‑1 (MQ=255)                                                                         gtcccgaacgtgcgccgccagcagccgcgccctctttatatttccagtgcgcagcaacacccaactctgcatctGTTTGGTGCGGATTTGGATCTCAACGGTTTTTCCACCCGGCCCCAGAACCACGGTATGAATAGACTGATAACCGtt                                                                        <  2:397316‑M2/76‑1 (MQ=255)                                                                             cgaacgtgcgccgccagcagccgcgccctctttatatttccagtgcgcagcaacacccaactctgcatctGTTTGGTGCGGATTTGGATCTCAACGGTTTTTCCACCCGGCCCCAGAACCACGGTATGAATAGACTGATAACCGTTTGGtt                                                                   >  1:401689‑M2/71‑151 (MQ=255)                                                                              gaacgtgcgccgccagcagccgcgccctctttatatttccagtgcgcagcaacacccaactctgcatctGTTTGGTGCGGATTTGGATCTCAACGTTTTTTCCCCCCGGCCCCAGAACCACGGTATGAATAGACTGATAACCGTTTGGtt                                                                   <  1:80149‑M2/81‑1 (MQ=255)                                                                                    gcgccgccagcagccgcgccctctttatatttccagtgcgcagcaacacccaactctgcatctGTTTGGTGCGGATTTGGATCTCAACGGTTTTTCCACCCGGCCCCAGAACCACGGTATGAATAGACTGATAACCGTTTGGTTTCGGGt                                                             <  1:568389‑M2/87‑1 (MQ=255)                                                                                     cgccgccagcagccgcgccctctttatatttccagtgcgcagcaacacccaactctgcatctGTTTGGTGCGGATTTGGATCTCAACGGTTTTTCAACCCGGCACCAGAACCACGGTATGAATAGACTGATAACCGTTTGGTTTCGGGtt                                                            >  2:811695‑M2/63‑150 (MQ=255)                                                                                          ccagcagccgcgccctctttatatttccagtgcgcagcaacacccaactctgcatctGTTTGGTGCGGATTTGGATCTCAACGGTTTTTCCACCCGGCCCCAGAACCACGGTATGAATAGACTGATAACCGTTTGGTTTCGGGTTAGCGAc                                                      <  2:489154‑M2/94‑1 (MQ=255)                                                                                                gccgcgccctctttatatttccagtgcgcagcaacacccaactctgcatctGTTTGGTGCGGATTTGGATCTCAACGGTTTTTCCACCCGGCCCCAGAACCACGGTATGAATAGACTGATAACCGTTTGGTTTCGGGTTAGCGACGTAATc                                                <  1:610465‑M2/100‑1 (MQ=255)                                                                                                gccgcgccctctttatatttccagtgcgcagcaacacccaactctgcatctGTTTGGTGCGGATTTGGATCTCAACGGTTTTTCCACCCGGCCCCAGAACCACGGTATGAATAGACTGATAACCGTTTGGTTTCGGGTTAGCGACGTAATc                                                <  2:296533‑M2/100‑1 (MQ=255)                                                                                                     gccctctttatatttccagtgcgcagcaacacccaactctgcatctGTTTGGTGCGGATTTGGATCTCAACGGTTTTTCCACCCGGCCCCAGAACCACGGTATGAATAGACTGATAACCGTTTGGTTTCGGGTTAGCGACGTAATCGTCa                                            >  2:53039‑M2/47‑150 (MQ=255)                                                                                                        ctctttatatttccagtgcgcagcaacacccaactctgcatctGTTTGGTGCGGATTTGGATCTCAACGGTTTTTCCACCCGGCCCCAGAACCACGGTATGAATAGACTGATAACCGTTTGGTTTCGGGTTAGCGACGTAATCGTCAAACt                                        >  1:400719‑M2/44‑151 (MQ=255)                                                                                                           tttatatttccagtgcgcagcaacacccaactctgcatctGTTTGGTGCGGATTTGGATCTCAACGGTTTTTCCACCCGGCCCCAGAACCACGGTATGAATAGACTGATAACCGTTTGGTTTCGGGTTAGCGACGTAATCGTCAAACTCAt                                     >  2:569291‑M2/41‑151 (MQ=255)                                                                                                            ttatatttccagtgcgcagcaacacccaactctgcatctGTTTGGTGCGGATTTGGATCTCAACGGTTTTTCCACCCGGCCCCAGAACCACGGTATGAATAGACTGATAACCGTTTGGTTTCGGGTTAGCGACGTAATCGTCAAACTCa                                      <  2:57371‑M2/110‑1 (MQ=255)                                                                                                              atatttccagtgcgcagcaacacccaactctgcatctGTTTGGTGCGGATTTGGATCTCAACGGTTTTTCCACCCGGCCCCAGAACCACGGTATGAATAGACTGATAACCGTTTGGTTTCGGGTTAGCGACGTAATCGTCAAACTCATcc                                   <  1:31458‑M2/113‑1 (MQ=255)                                                                                                              atatttccagtgcgcagcaacacccaactctgcatctGTTTGGTGCGGATTTGGATCTCAACGGTTTTTCCACCCGGCCCCAGAACCACGGTATGAATAGACTGATAACCGTTTGGTTTCGGGTTAGCGACGTAATCGTCAAACGCATcc                                   <  1:435029‑M2/113‑1 (MQ=255)                                                                                                                  ttccagtgcgcagcaacacccaactctgcatctGTTTGGTGCGGATTTGGATCTCAACGGTTTTTCCACCCGGCCCCAGAACCACGGTATGAATAGACTGATAACCGTTTGGTTTCGGGTTAGCGACGTAATCGTCAAACTCATCCGGCa                               <  1:88120‑M2/117‑1 (MQ=255)                                                                                                                     cagtgcgcagcaacacccaactctgcatctGTTTGGTGCGGATTTGGATCTCAACGGTTTTTCCACCCGGCCCCAGAACCACGGTATGAATAGACTGATAACCGTTTGGTTTCGGGTTAGCGACGTAATCGTCAAACTCATCCGGCAGGTg                           >  1:521526‑M2/31‑151 (MQ=255)                                                                                                                     cagtgcgcagcaacacccaactctgcatctGTTTGGTGCGGATTTGGATCTCAACGGTTTTTCCACCCGGCCCCAGAACCACGGTATGAATAGACTGATAACCGTTTGGTTTCGGGTTAGCGACGTAATCGTCAAACTCATCCGGCAGGTg                           <  2:656444‑M2/121‑1 (MQ=255)                                                                                                                        tgcgcagcaacacccaactctgcatctGTTTGGTGCGGATTTGGATCTCAACGGTTTTTCCACCCGGCCCCAGAACCACGGTATGAATAGACTGATAACCGTTTGGTTTCGGGTTAGCGACGTAATCGTCAAACTCATCCGGCAGGTGGCg                        >  2:345158‑M2/28‑151 (MQ=255)                                                                                                                          cgcagcaacacccaactctgcatctGTTTGGTGCGGATTTGGATCTCAACGGTTTTTCCACCCGGCCCCAGAACCACGGTATGAATAGACTGATAACCGTTTGGTTTCGGGTTAGCGACGTAATCGTCAAACTCATCCAGCAGGTGGCGAt                      >  1:389812‑M2/26‑151 (MQ=255)                                                                                                                                  cacccaactctgcatctGTTTGGTGCGGATTTGGATCTCAACGGTTTTTCCACCCGGCCCCAGAACCACGGTATGAATAGACTGATAACCGTTTGGTTTCGGGTTAGCGACGTAATCGTCAAACTCATCCGGCAGGTGGCGATagtgagtg              >  2:459825‑M2/18‑151 (MQ=255)                                                                                                                                  cacccaactctgcatctGTTTGGTGCGGATTTGGATCTCAACGGTTTTTCCACCCGGCCCCAGAACCACGGTATGAATAGACTGATAACCGTTTGGTTTCGGGTTAGCGACGTAATCGTCAAACTCATCCGGCAGGTGGCGATagtgagtg              >  2:342660‑M2/18‑151 (MQ=255)                                                                                                                                        actctgcatctGTTTGGGGCGGATTTGGAGCTCAACGGTTTTTCCACCCGGCCCCAGAACCACGGTATGAATAGACTGATAACCGTTTGGTTTCGGGTTAGCGACGTAATCGTCAAACTCATCCGGCAGGTGGCGATAGTGAGTGTGCACt        >  1:376620‑M2/12‑151 (MQ=255)                                                                                                                                                tccgTTTGGTGCGGATTTGGATCTCAACGGTTTTTCCACCCGGCCCCAGAACCACGGTATGAATAGACTGATAACCGTTTGGTTTCGGGTTAGCGACGTAATCGTCAAACTCATCCGGCAGGTGGCGATAGTGAGTGTGCACTATcccc  >  1:251719/4‑149 (MQ=255)                                                                                                                                                   |                                                                                                                                                  GGAATCAGCCATCTCTTCCTGCCACGCAATCAGTTTACGCAGCCAGGCAATCCGGTCTTCATGTCCCGAACGTGCGCCGCCAGCAGCCGCGCCCTCTTTATATTTCCAGTGCGCAGCAACACCCAACTCTGCATCTTCATGCATCTGTTTGGTGCGGATTTGGATCTCAACGGTTTTTCCACCCGGCCCCAGAACCACGGTATGAATAGACTGATAACCGTTTGGTTTCGGGTTAGCGACGTAATCGTCAAACTCATCCGGCAGGTGGCGATAGTGAGTGTGCACTATCCCC  >  NC_000913/2912482‑2912773

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 24 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑