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

GCATTTCGCCGGAATCAGCCATCTCTTCCTGCCACGCAATCAGTTTACGCAGCCAGGCAATCCGGTCTTCATGTCCCGAACGTGCGCCGCCAGCAGCCGCGCCCTCTTTATATTTCCAGTGCGCAGCAACACCCAACTCTGCATCTTCATGCATCTGTTTGGTGCGGATTTGGATCTCAACGGTTTTTCCACCCGGCCCCAGAACCACGGTATGAATAGACTGATAACCGTTTGGTTTCGGGTTAGCGACGTAATCGTCAAACTCATCCGGCAGGTGGCGATAGTG  >  NC_000913/2912472‑2912757                                                                                                                                                  |                                                                                                                                             gCATTTCGCCGGAATCAGCCATCTCTTCCTGCCACGCAATCAGTTTACGCAGCCAGGCAATCCGGTCTTCATGTCCCGAACGTGCGCCGCCAGCAGCCGCGCCCTCTTTATATTTCCAGTGCGCAGCAACACCCAACTCTGCATCTgtttg                                                                                                                                         >  2:118220‑M1/1‑146 (MQ=255)         ccGGAATCAGCCATCTCTTCCTGCCACGCAATCAGTTTACGCAGCCAGGCAATCCGGTCTTCATGTCCCGAACGTGCGCCGCCAGCAGCCGCGCCCTCTTTATATTTCCAGTGCGCAGCAACACCCAACTCTGCATCTgtttggtgcggat                                                                                                                                 <  1:701054‑M1/151‑14 (MQ=255)          cGGAATCAGCCATCTCTTCCTGCCACGCAATCAGTTTACGCAGCCAGGCAATCCGGTCTTCATGTCCCGAACGTGCGCCGCCAGCAGCCGCGCCCTCTTTATATTTCCAGTGCGCAGCAACACCCAACTCTGCATCTgtttggtgcggatt                                                                                                                                >  1:593137‑M1/1‑137 (MQ=255)            gAATCAGCCATCTCTTCCTGCCACGCAATCAGTTTACGCAGCCAGGCAATCCGGTCTTCATGTCCCGAACGTGCGCCGCCAGCAGCCGCGCCCTCTTTATATTTCCAGTGCGCAGCAACACCCAACTCTGCATCTgtttggtgcggatttg                                                                                                                              <  1:210779‑M1/151‑17 (MQ=255)              aTCAGCCATCTCTTCCTGCCACGCAATCAGTTTACGCAGCCAGGCAATCCGGTCTTCATGTCCCGAACGTGCGCCGCCAGCAGCCGCGCCCTCTTTATATTTCCAGTGCGCAGCAACACCCAACTCTGCATCTgtttggtgcggatttgga                                                                                                                            <  1:749813‑M1/151‑19 (MQ=255)              aTCAGCCATCTCTTCCTGCCACGCAATCAGTTTACGCAGCCAGGCAATCCGGTCTTCATGTCCCGAACGTGCGCCGCCAGCAGCCGCGCCCTCTTTATATTTCCAGTGCGCAGCAACACCCAACTCTGCATCTgtttggtgcggatttgg                                                                                                                             >  2:459027‑M1/1‑133 (MQ=255)              aTCAGCCATCTCTTCCTGCCACGCAATCAGTTTACGCAGCCAGGCAATCCGGTCTTCATGTCCCGAACGTGCGCCGCCAGCAGCCGCGCCCTCTTTATATTTCCAGTGCGCAGCAACACCCAACTCTGCATCTgtttggtgcggatttcg                                                                                                                             >  1:710164‑M1/1‑133 (MQ=255)                   ccATCTCTTCCTGCCACGCAATCAGTTTACGCAGCCAGGCAATCCGGTCTTCATGTCCCGAACGTGCGCCGCCAGCAGCCGCGCCCTCTTTATATTTCCAGTGCGCAGCAACACCCAACTCTGCATCTgtttggtgcggatttggatctca                                                                                                                       <  1:566232‑M1/151‑24 (MQ=255)                      tctcTTCCTGCCACGCAATCAGTTTACGCAGCCAGGCAATCCGGTCTTCATGTCCCGAACGTGCGCCGCCAGCAGCCGCGCCCTCTTTATATTTCCAGTGCGCAGCAACACCCAACTCTGCATCTgtttggtgcggatttggatctcaacg                                                                                                                    >  2:576179‑M1/1‑125 (MQ=255)                           tCCTGCCACGCAATCAGTTTACGCAGCCAGGCAATCCGGTCTTCATGTCCCGAACGGGCGCCGCCAGCAGCCGCGCCCTCTTTATATTTCCAGTGCGCAGCAACACCCAACTCTGCATCTgtttggtgcggatttggatctcaacggtt                                                                                                                 <  2:810394‑M1/149‑30 (MQ=255)                               gCCACGCAATCAGTTTACGCAGCCAGGCAATCCGGTCTTCATGTCCCGAACGTGCGCCGCCAGCAGCCGCGCCCTCTTTATATTTCCAGTGCGCAGCAACACCCAACTCTGCATCTgtttggtgcggatttggatctcaacggtttttcca                                                                                                           <  2:468328‑M1/151‑36 (MQ=255)                               gCCACGCAATCAGTTTACGCAGCCAGGCAATCCGGTCTTCATGTCCCGAACGTGCGCCGCCAGCAGCCGCGCCCTCTTTATATTTCCAGTGCGCAGCAACACCCAACTCTGCATCTgtttggtgcggatttggatctcaacggtttttc                                                                                                             <  2:131700‑M1/149‑34 (MQ=255)                               gCCACGCAATCAGTTTACGCAGCCAGGCAATCCGGTCTTCATGTCCCGAACGTGCGCCGCCAGCAGCCGCGCCCTCTTTATATTTCCAGTGCGCAGCAACACCCAACTCTGCATCTgtttggtgcggatttggatctcaacggtttttc                                                                                                             <  2:593171‑M1/149‑34 (MQ=255)                               gCCACGCAATCAGTTTACGCAGCCAGGCAATCCGGTCTTCATGTCCCGAACGTGCGCCGCCAGCAGCCGCGCCCTCTTTATATTTCCAGTGCGCAGCAACACCCAACTCTGCATCTgtttggtgcggatttggatctcaacggtttttc                                                                                                             <  1:258373‑M1/149‑34 (MQ=255)                                 cACGCAATCAGTTTACGCAGCCAGGCAATCCGGTCTTCATGTCCCGAACGTGCGCCGCCAGCAGCCGCGCCCTCTTTATATTTCCAGTGCGCAGCAACACCCAACTCTGCATCTgtttggtgcggatttggatctcaacggtttttccacc                                                                                                         <  1:728258‑M1/151‑38 (MQ=255)                                     cAATCAGTTTACGCAGCCAGGCAATCCGGTCTTCATGTCCCGAACGTGCGCCGCCAGCAGCCGCGCCCTCTTTATATTTCCAGTGCGCAGCAACACCCAACTCTGCATCTgtttggtgcggatttggatctcaacggtttttccacccggc                                                                                                     >  2:35053‑M1/1‑110 (MQ=255)                                      aaTCAGTTTACGCAGCCAGGCAATCCGGTCTTCATGTCCCGAACGTGCGCCGCCAGCAGCCGCGCCCTCTTTATATTTCCAGTGCGCAGCAACACCCAACTCTGCATCTgtttggtgcggatttggatctcaacggtttttccacccggcc                                                                                                    >  2:358878‑M1/1‑109 (MQ=255)                                           gTTTACGCAGCCAGGCAATCCGGTCTTCATGTCCCGAACGTGCGCCGCCAGCAGCCGCGCCCTCTTTATATTTCCAGTGCGCAGCAACACCCAACTCTGCATCTgtttggtgcggatttggatctcaacggtttttccacccggcccc                                                                                                  >  1:471126‑M1/1‑104 (MQ=255)                                               aCGCAGCCAGGCAATCCGGTCTTCATGTCCCGAACGTGCGCCGCCAGCAGCCGCGCCCTCTTTATATTTCCAGTGCGCAGCAACACCCAACTCTGCATCTgtttggtgcggatttggatctcaacggtttttccacccggccccagaacc                                                                                            >  2:205521‑M1/1‑100 (MQ=255)                                                cGCAGCCAGGCAATCCGGTCTTCATGTCCCGAACGGGCGCCGCCAGCAGCCGCGCCCTCTTTATATTTCCAGTGCGCAGCAACACCCAACTCTGCATCTgtttggtgcggatttggatctcaacggtgtttccacccggccccagaacc                                                                                            <  1:16747‑M1/149‑51 (MQ=255)                                                  cAGCCAGGCAATCCGGTCTTCATGTCCCGAACGTGCGCCGCCAGCAGCCGCGCCCTCTTTATATTTCCAGTGCGCAGCAACACCCAACTCTGCATCTgtttggtgcggatttggatctcaacggtttttccacccggccccagaaccacgg                                                                                        <  2:555276‑M1/151‑55 (MQ=255)                                                   aGCCAGGCAATCCGGTCTTCATGTCCCGAACGTGCGCCGCCAGCAGCCGCGCCCTCTTTATATTTCCAGTGCGCAGCAACACCCAACTCTGCATCTgtttggtgcggatttggatctcaacggtttttccacccggccccagaaccacggt                                                                                       >  1:438834‑M1/1‑96 (MQ=255)                                                     ccAGGCAATCCGGTCTTCATGTCCCGAACGTGCGCCGCCAGCAGCCGCGCCCTCTTTATATTTCCAGTGCGCAGCAACACCCAACTCTGCATCTgtttggtgcggatttggatctcaacggtttttccacccggccccagaaccacggta                                                                                      >  2:512275‑M1/1‑94 (MQ=255)                                                      cAGGCAATCCGGTCTTCATGTCCCGAACGTGCGCCGCCAGCAGCCGCGCCCTCTTTATATTTCCAGTGCGCAGCAACACCCAACTCTGCATCTgtttggtgcggatttggatctcaacggtttttccacccggccccagaaccacggtatg                                                                                    <  1:2725‑M1/151‑59 (MQ=255)                                                           aaTCCGGTCTTCATGTCCCGAACGTGCGCCGCCAGCAGCCGCGCCCTCTTTATATTTCCAGTGCGCAGCAACACCCAACTCTGCATCTgtttggtgcggatttggatctcaacggtttttccacccggccccagaaccacggtatgaat                                                                                 >  1:11483‑M1/1‑88 (MQ=255)                                                            aTCCGGTCTTCATGTCCCGAACGTGCGCCGCCAGCAGCCGCGCCCTCTTTATATTTCCAGTGCGCAGCAACACCCAACTCTGCATCTgtttggtgcggatttggatctcaacggtttttccacccggccccagaaccacggtatgaat                                                                                 >  2:390505‑M1/1‑87 (MQ=255)                                                             tctggtctTCATGTCCCGAACGTGCGCCGCCAGCAGCCGCGCCCTCTTTATATTTCCAGTGCGCAGCAACACCCAACTCTGCATCTgtttggtgcggatttggatctcaacggtttttccacccggccccagaaccacggtatgaatagac                                                                             >  2:312023‑M1/4‑86 (MQ=255)                                                                ggTCTTCATGTCCCGAACGTGCGCCGCCAGCAGCCGCGCCCTCTTTATATTTCCAGTGCGCAGCAACACCCAACTCTGCATCTgtttggtgcggatttggatctcaacggtttttccacccggccccagaaccacggtatgaatagactga                                                                          <  1:390549‑M1/151‑69 (MQ=255)                                                                ggTCTTCATGTCCCGAACGTGCGCCGCCAGCAGCCGCGCCCTCTTTATATTTCCAGTGCGCAGCAACACCCAACTCTGCATCTgtttggtgcggatttggatctcaacggtttttccacccggccccagaaccacggtatgaatagactg                                                                           >  2:138759‑M1/1‑83 (MQ=255)                                                                      cATGTCCCGAACGTGCGCCGCCAGCAGCCGCGCCCTCTTTATATTTCCAGTGCGCAGCAACACCCAACTCTGCATCTgtttggtgcggatttggatctcaacggtttttccacccggccccagaaccacggtatgaatagactgataaccg                                                                    <  1:27990‑M1/151‑75 (MQ=255)                                                                      cATGTCCCGAACGTGCGCCGCCAGCAGCCGCGCCCTCTTTATATTTCCAGTGCGCAGCAACACCCAACTCTGCATCTgtttggtgcggatttggatctcaacggtttttccacccggccccagaaccacggtatgaatagactgataacc                                                                     <  2:93657‑M1/150‑74 (MQ=255)                                                                         gTCCCGAACGTGCGCCGCCAGCAGCCGCGCCCTCTTTATATTTCCAGTGCGCAGCAACACCCAACTCTGCATCTgtttggtgcggatttggatctcaacggtttttccacccggccccagaaccacggtatgaatagactgataaccgtt                                                                  <  2:397316‑M1/150‑77 (MQ=255)                                                                             cGAACGTGCGCCGCCAGCAGCCGCGCCCTCTTTATATTTCCAGTGCGCAGCAACACCCAACTCTGCATCTgtttggtgcggatttggatctcaacggtttttccacccggccccagaaccacggtatgaatagactgataaccgtttggtt                                                             >  1:401689‑M1/1‑70 (MQ=255)                                                                              gAACGTGCGCCGCCAGCAGCCGCGCCCTCTTTATATTTCCAGTGCGCAGCAACACCCAACTCTGCATCTgtttggtgcggatttggatctcaacgttttttccccccggccccagaaccacggtatgaatagactgataaccgtttggtt                                                             <  1:80149‑M1/150‑82 (MQ=255)                                                                                    gcgcCGCCAGCAGCCGCGCCCTCTTTATATTTCCAGTGCGCAGCAACACCCAACTCTGCATCTgtttggtgcggatttggatctcaacggtttttccacccggccccagaaccacggtatgaatagactgataaccgtttggtttcgggt                                                       <  1:568389‑M1/150‑88 (MQ=255)                                                                                     cgccgcCAGCAGCCGCGCCCTCTTTATATTTCCAGTGCGCAGCAACACCCAACTCTGCATCTgtttggtgcggatttggatctcaacggtttttcaacccggcaccagaaccacggtatgaatagactgataaccgtttggtttcgggtt                                                      >  2:811695‑M1/1‑62 (MQ=255)                                                                                          ccAGCAGCCGCGCCCTCTTTATATTTCCAGTGCGCAGCAACACCCAACTCTGCATCTgtttggtgcggatttggatctcaacggtttttccacccggccccagaaccacggtatgaatagactgataaccgtttggtttcgggttagcgac                                                <  2:489154‑M1/151‑95 (MQ=255)                                                                                                gcCGCGCCCTCTTTATATTTCCAGTGCGCAGCAACACCCAACTCTGCATCTgtttggtgcggatttggatctcaacggtttttccacccggccccagaaccacggtatgaatagactgataaccgtttggtttcgggttagcgacgtaatc                                          <  1:610465‑M1/151‑101 (MQ=255)                                                                                                gcCGCGCCCTCTTTATATTTCCAGTGCGCAGCAACACCCAACTCTGCATCTgtttggtgcggatttggatctcaacggtttttccacccggccccagaaccacggtatgaatagactgataaccgtttggtttcgggttagcgacgtaatc                                          <  2:296533‑M1/151‑101 (MQ=255)                                                                                                     gcCCTCTTTATATTTCCAGTGCGCAGCAACACCCAACTCTGCATCTgtttggtgcggatttggatctcaacggtttttccacccggccccagaaccacggtatgaatagactgataaccgtttggtttcgggttagcgacgtaatcgtca                                      >  2:53039‑M1/1‑46 (MQ=255)                                                                                                        ctctTTATATTTCCAGTGCGCAGCAACACCCAACTCTGCATCTgtttggtgcggatttggatctcaacggtttttccacccggccccagaaccacggtatgaatagactgataaccgtttggtttcgggttagcgacgtaatcgtcaaact                                  >  1:400719‑M1/1‑43 (MQ=255)                                                                                                           tttATATTTCCAGTGCGCAGCAACACCCAACTCTGCATCTgtttggtgcggatttggatctcaacggtttttccacccggccccagaaccacggtatgaatagactgataaccgtttggtttcgggttagcgacgtaatcgtcaaactcat                               >  2:569291‑M1/1‑40 (MQ=255)                                                                                                            ttATATTTCCAGTGCGCAGCAACACCCAACTCTGCATCTgtttggtgcggatttggatctcaacggtttttccacccggccccagaaccacggtatgaatagactgataaccgtttggtttcgggttagcgacgtaatcgtcaaactca                                <  2:57371‑M1/149‑111 (MQ=255)                                                                                                              atatTTCCAGTGCGCAGCAACACCCAACTCTGCATCTgtttggtgcggatttggatctcaacggtttttccacccggccccagaaccacggtatgaatagactgataaccgtttggtttcgggttagcgacgtaatcgtcaaactcatcc                             <  1:31458‑M1/150‑114 (MQ=255)                                                                                                              atatTTCCAGTGCGCAGCAACACCCAACTCTGCATCTgtttggtgcggatttggatctcaacggtttttccacccggccccagaaccacggtatgaatagactgataaccgtttggtttcgggttagcgacgtaatcgtcaaacgcatcc                             <  1:435029‑M1/150‑114 (MQ=255)                                                                                                                  ttCCAGTGCGCAGCAACACCCAACTCTGCATCTgtttggtgcggatttggatctcaacggtttttccacccggccccagaaccacggtatgaatagactgataaccgtttggtttcgggttagcgacgtaatcgtcaaactcatccggca                         <  1:88120‑M1/150‑118 (MQ=255)                                                                                                                     cAGTGCGCAGCAACACCCAACTCTGCATCTgtttggtgcggatttggatctcaacggtttttccacccggccccagaaccacggtatgaatagactgataaccgtttggtttcgggttagcgacgtaatcgtcaaactcatccggcaggtg                     >  1:521526‑M1/1‑30 (MQ=255)                                                                                                                     cAGTGCGCAGCAACACCCAACTCTGCATCTgtttggtgcggatttggatctcaacggtttttccacccggccccagaaccacggtatgaatagactgataaccgtttggtttcgggttagcgacgtaatcgtcaaactcatccggcaggtg                     <  2:656444‑M1/151‑122 (MQ=255)                                                                                                                        tGCGCAGCAACACCCAACTCTGCATCTgtttggtgcggatttggatctcaacggtttttccacccggccccagaaccacggtatgaatagactgataaccgtttggtttcgggttagcgacgtaatcgtcaaactcatccggcaggtggcg                  >  2:345158‑M1/1‑27 (MQ=255)                                                                                                                          cgcAGCAACACCCAACTCTGCATCTgtttggtgcggatttggatctcaacggtttttccacccggccccagaaccacggtatgaatagactgataaccgtttggtttcgggttagcgacgtaatcgtcaaactcatccagcaggtggcgat                >  1:389812‑M1/1‑25 (MQ=255)                                                                                                                                  cacCCAACTCTGCATCTgtttggtgcggatttggatctcaacggtttttccacccggccccagaaccacggtatgaatagactgataaccgtttggtttcgggttagcgacgtaatcgtcaaactcatccggcaggtggcgatagtgagtg        >  2:459825‑M1/1‑17 (MQ=255)                                                                                                                                  cacCCAACTCTGCATCTgtttggtgcggatttggatctcaacggtttttccacccggccccagaaccacggtatgaatagactgataaccgtttggtttcgggttagcgacgtaatcgtcaaactcatccggcaggtggcgatagtgagtg        >  2:342660‑M1/1‑17 (MQ=255)                                                                                                                                        aCTCTGCATCTgtttggggcggatttggagctcaacggtttttccacccggccccagaaccacggtatgaatagactgataaccgtttggtttcgggttagcgacgtaatcgtcaaactcatccggcaggtggcgatagtgagtgtgcact  >  1:376620‑M1/1‑11 (MQ=255)                                                                                                                                                  |                                                                                                                                             GCATTTCGCCGGAATCAGCCATCTCTTCCTGCCACGCAATCAGTTTACGCAGCCAGGCAATCCGGTCTTCATGTCCCGAACGTGCGCCGCCAGCAGCCGCGCCCTCTTTATATTTCCAGTGCGCAGCAACACCCAACTCTGCATCTTCATGCATCTGTTTGGTGCGGATTTGGATCTCAACGGTTTTTCCACCCGGCCCCAGAACCACGGTATGAATAGACTGATAACCGTTTGGTTTCGGGTTAGCGACGTAATCGTCAAACTCATCCGGCAGGTGGCGATAGTG  >  NC_000913/2912472‑2912757

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: ‑