Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
JC	NC_000913	3,822,651	(CAGGAT)1→2	coding (252/2109 nt)	spoT →	bifunctional (p)ppGpp synthase/hydrolase SpoT

New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913	3822646 =	1 (0.020)	27 (0.670)	25/482	0.6	94.7%	coding (247/2109 nt)	spoT	bifunctional (p)ppGpp synthase/hydrolase SpoT
	?	NC_000913	= 3822651	2 (0.050)					coding (252/2109 nt)	spoT	bifunctional (p)ppGpp synthase/hydrolase SpoT

ATTTCGAGAGTTTCACGGGCGATGCGGCGACGTTTGTCCGGGCGAAGTGAGCCCAGCGTGCGCATGTTGTGGGTACGGTCGGCAAGTTTGATGAGGATGACGCGGATATCCTGCACCATCGCCATAATCATCTTGCGAAAGTTTTCGGCCTGCGCCTCTTTCTTATCGCGGAACTTGAGTTTATCAAGTTTCGACACCCCCTCTACCAGCTCGGCGACGCTTTTACCAAAAAGCTGTTCCATATCCTG‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  <  NC_000913/3822893‑3822646 ‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑ATCCTGGTAGGTGGCGGGAGTATCTTCAATCACGTCATGCAGCAGCGCCGCCATCAGCGTTTCATAGTCGAGTTTCATCTCGGCCAGAATGCAGGCAACCGCTACCGGGTGCGTGATATAGGGTTCACCGCTTGAACGTGTTTGCCCCTCGTGAGCATCACGTGCAACGAGATACGCCTGCCGCAGACGCTTGATTTGGTCTTCCGGCAGGTAGGTTTGAATCAGTTGATTCAGGCTTTCAAAC  <  NC_000913/3822651‑3822408                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                              ATTTCGAGAGTTTCACGGGCGATGCGGCGACGTTTGTCCGGGCGAAGTGAGCCCAGCGTGCGCATGTTGTGGGTACGGTCGGCAAGTTTGATGAGGATGACGCGGATATCCTGCACCATCGCCATAATCATCTTGCGAAAGTTTTCGGCCTGCGCCTCTTTCTTATCGCGGAACGTGAGTTTATCAAGTTTCGACACCCCCTCTACCAGCTCGGCGACGCTTTTACCAAAAAGCTGTTCCATATCCTGATC                                                                                                                                                                                                                                                   <  2:47077/251‑1         AGTTTCACGGGCGATGCGGCGACGTTTGTCCGGGCGAAGTGAGCCCAGCGTGCGCATGTTGTGGGTACGGTCGGCAAGTTTGATGAGGATGACGCGGATATCCTGCACCATCGCCATAATCATCTTGCGAAAGTTTTCGGCCTGCGCCTCTTTCTTATCGCGGAACTTGAGTTTATCAAGTTTCGACACCCCCTCTACCAGCTCGGCGACGCTTTTACCAAAAAGCTGTTCCATATCCTGATCCTGGTAGG                                                                                                                                                                                                                                           >  2:86400/1‑251             TCACGGGCGATGCGGCGACGTTTGTCCGGGCGAAGTGAGCCCAGCGTGCGCATGTTGTGGGTACGGTCGGCAAGTTTGATGAGGATGACGCGGATATCCTGCACCATCGCCATAATCATCTTGCGAAAGTTTTCGGCCTGCGCCTCTTTCTTATCGCGGAACTTGAGTTTATCAAGTTTCGACACCCCCTCTACCAGCTCGGCGACGCTTTTACCAAAAAGCTGTTCCATATCCTGATCCTGGTAGGTGG                                                                                                                                                                                                                                        <  1:282102/250‑1               ACGGGCGATGCGGCGACGTTTGTCCGGGCGAAGTGAGCCCAGCGTGCGCATGTTGTGGGTACGGTCGGCAAGTTTGATGAGGATGACGCGGATATCCTGCACCATCGCCATAATCATCTTGCGAAAGTTTTCGGCCTGCGCCTCTTTCTTATCGCGGAACTTGAGTTTATCAAGTTTCGACACCCCCTCTACCAGCTCGGCGACGCTTTTACCAAAAAGCTGTTCCATATCCTGATCCTGGTAGGTGGCGG                                                                                                                                                                                                                                     >  1:88617/1‑251                 gtgcgtTGCGGCGACGTTTGTCCGGGCTAAGTGGGCCCAGCGTGCGCATGTTGTGGGTACGGTCGGCAAGTTTGATGAGGAGGACGCGGATATCCTGCACCATCGCCATAATCATCTTGCGAAAGTTTTCGGTCTGCGCCTCTTTCTTATCGCGGAACTTGAGTTTATCAAGTTTCGACACCCCCTCTACCAGCTCGGCGACGCTTTTACCAAAAAGCTGTTCCAAATCCGGAGCCTGGTAGGTGTCGGGA                                                                                                                                                                                                                                   <  2:383172/245‑1                                           gGAAGTGAGCCCAGCGTGCGCATGTTGTGGGTACGGTCGGCAAGTTTGATGAGGATGACGCGGATATCCTGCACCATCGCCATAATCATCTTGCGAAAGTTTTCGGCCTGCGCCTCTTTCTTATCGCGGAACTTGAGTTTATCAAGTTTCGACACCCCCTCTACCAGCTCGGCGACGCTTTTACCAAAAAGCTGTTCCATATCCTGATCCTGGTAGGTGGCGGGAGTATCTTCAATCACGTCATGCAGCAG                                                                                                                                                                                                         <  2:293387/250‑1                                                               CATGTTGTGGGTACGGTCGGCAAGTTTGATGAGGATGACGCGGATATCCTGCACCATCGCCATAATCATCTTGCGAAAGTTTTCGGCCTGCGCCTCTTTCTTATCGCGGAACTTGAGTTTATCAAGTTTCGACACCCCCTCTACCAGCTCGGCGACGCTTTTACCAAAAAGCTGTTCCATATCCTGATCCTGGTAGGTGGCGGGAGTATCTTCAATCACGTCATGCAGCAGCGCCGCCATCAGCGTTTCAT                                                                                                                                                                                     <  2:251917/251‑1                                                                                             GAGGATGACGCGGATATCCTGCACCATCGCCATAATCATCTTGCGAAAGTTTTCGGCCTGCGCCTCTTTCTTATCGCGGAACTTGAGTTTATCAAGTTTCGACACCCCCTCTACCAGCTCGGCGACGCTTTTACCAAAAAGCTGTTCCATATCCTGATCCTGGTAGGTGGCGGGAGTATCTTCAATCACGTCATGCAGCAGCGCCGCCATCAGCGTTTCATAGTCGAGTTTCATCTCGGCCAGAATGCAGG                                                                                                                                                       <  2:340150/251‑1                                                                                                        GGATATCCTGCACCATCGCCATAATCATCTTGCGAAAGTTTTCGGCCTGCGCCTCTTTCTTATCGCGGAACTTGAGTTTATCAAGTTTCGACACCCCCTCTACCAGCTCGGCGACGCTTTTACCTAAAAGCTGTTCCATATCCTGATCCTGGTAGGTGGCGGGAGTATCTTCAATCACGTCATGCAGCAGCGCCGCCATCAGCGTTTCATAGTCGAGTTTCATCTCGGCCAGAATGCAGGCAACCGCTACC                                                                                                                                            >  2:281029/1‑251                                                                                                         GATATCCTGCACCATCGCCATAATCATCTTGCGAAAGTTTTCGGCCTGCGCCTCTTTCTTATCGCGGAACTTGAGTTTATCAAGTTTCGACACCCCCTCTACCAGCTCGGCGACGCTTTTACCAAAAAGCTGTTCCATATCCTGATCCTGGTAGGTGGCGGGAGTATCTTCAATCACGTCATGCAGCAGCGCCGCCATCAGCGTTTCATAGTCGAGTTTCATCTCGGCCAGAATGCAGGCAACCGCTACCG                                                                                                                                           <  1:121125/251‑1                                                                                                         GATATCCTGCACCATCGCCATAATCATCTTGCGAAAGTTTTCGGCCTGCGCCTCTTTCTTATCGCGGAACTTGAGTTTATCAAGTTTCGACACCCCCTCTACCAGCTCGGCGACGCTTTTACCAAAAAGCTGTTCCATATCCTGATCCTGGTAGGTGGCGGGAGTATCTTCAATCACGTCATGCAGCAGCGCCGCCATCAGCGTTTCATAGTCGAGTTTCATCTCGGCCAGAATGCAGGCAACCGCTACCG                                                                                                                                           <  1:121140/251‑1                                                                                                                     CATCGCCATAATCATCTTGCGAAAGTTTTCGGCCTGCGCCTCTTTCTTATCGCGGAACTTGAGTTTATCAAGTTTCGACACCCCCTCTACCAGCTCGGCGACGCTTTTACCAAAAAGCTGTTCCATATCCTGATCCTGGTAGGTGGCGGGAGTATCTTCAATCACGTCATGCAGCAGCGCCGCCATCAGCGTTTCATAGTCGAGTTTCATCTCGGCCAGAATGCAGGCAACCGCTACCGGGTGCGTGATAT                                                                                                                               <  2:160666/251‑1                                                                                                                            ATAATCATCTTGCGAAAGTTTTCGGCCTGCGCCTCTTTCTTATCGCGGAACTTGAGTTTATCAAGTTTCGACACCCCCTCTACCAGCTCGGCGACGCTTTTACCAAAAAGCTGTTCCATATCCTGATCCTGGTAGGTGGCGGGAGTATCTTCAATCACGTCATGCAGCAGCGCCGCCATCAGCGTTTCATAGTCGAGTTTCATCTCGGCCAGAATGCAGGCAACCGCTACCGGGTGCGTGATATAGGGTTC                                                                                                                        >  1:119344/1‑251                                                                                                                                                        GCGCCTCGTTCTTAGCGCGGAACTTGAGTTTATCAAGTTTCGACACCCCCGCTACCAGCTCGGCGACGCTTTTACCAAAAAGCTGTTCCATATCCTGCTCCTGGTAGGTGGCGGGAGTATCTTCAATCACGTCATGCAGCAGCGCCGCCATCAGCGCTTCATAGTCGAGTTTCATCTCGGCCAGAATGCAGGCAACCGCTACCGGGTGCGTGATATAGGGTTCACCGCTTGAACGTTTTTGCCCCTCGTGA                                                                                            <  2:246473/251‑1                                                                                                                                                                       CGCGGAACTTGAGTTTATCAAGTTTCGACACCCCCTCTACCAGCTCGGCGACGCTTTTACCAAAAAGCTGTTCCATATCCTGATCCTGGTAGGTGGCGGGAGTATCTTCAATCACGTCATGCAGCAGCGCCGCCATCAGCGTTTCATAGTCGAGTTTCATCTCGGCCAGAATGCAGGCAACCGCTACCGGGTGCGTGATATAGGGTTCACCGCTTGAACGTGTTTGCCCCTCGTG                                                                                             <  1:319160/235‑1                                                                                                                                                                       CGCGGAACTTGAGTTTATCAAGTTTCGACACCCCCTCTACCAGCTCGGCGACGCTTTTACCAAAAAGCTGTTCCATATCCTGATCCTGGTAGGTGGCGGGAGTATCTTCAATCACGTCATGCAGCAGCGCCGCCATCAGCGTTTCATAGTCGAGTTTCATCTCGGCCAGAATGCAGGCAACCGCTACCGGGTGCGTGATATAGGGTTCACCGCTTGAACGTGTTTGCCCCTCGTG                                                                                             >  2:319160/1‑235                                                                                                                                                                          GGAACTTGAGTTTATCAAGTTTCGACACCCCCTCTACCAGCTCGGCGACGCTTTTACCAAAAAGCTGTTCCATATCCTGATCCTGGTAGGTGGCGGGAGTATCTTCAATCACGTCATGCAGCAGCGCCGCCATCAGCGTTTCATAGTCGAGTTTCATCTCGGCCAGAATGCAGGCAACCGCTACCGGGTGCGTGATATAGGGTTCACCGCTTGAACGTGTTTGCCCCTCGTGAGCATCACGTGCAACGAGA                                                                          <  2:327345/251‑1                                                                                                                                                                            AACTTGAGTTTATCAAGTGTCGACACCCCCTCTACCAGCTCGGCGACGCTTTTACCAAAAAGCTGTTCCATATCCTGATCCTGGTAGGTGGCGGGAGTATCTTCAATCACGTCATGCAGCAGCGCCGCCATCAGCGTTTCATAGTCGAGTTTCATCTCGGCCAGAATGCAGGCAACCGCTACCGGGTGCGTGATATAGGGTTCACCGCTTGAACGTGTTTGCCCCTCGTGAGCATCACGTGCAACGAGATA                                                                        <  2:162050/251‑1                                                                                                                                                                                                  ACACCCCCTCTACCAGCTCGGCGACGCTTTTACCAAAAAGCTGTTCCATATCCTGATCCTGGTAGGTGGCGGGAGTATCTTCAATCACGTCATGCAGCAGCGCCGCCATCAGCGTTTCATAGTCGAGTTTCATCTCGGCCAGAATGCAGGCAACCGCTACCGGGTGCGTGATATAGGGTTCACCGCTTGAACGTGTTTGCCCCTCGTGAGCATCACGTGCAACGAGATACGCCTGCCGCAGACGCTTGATT                                                  >  1:261933/1‑251                                                                                                                                                                                                      CCCCTCTACCAGCTCGGCGACGCTTTTACCAAAACGCTGTTCCATATCCTGATCCTGGTAGGTGGCGGGAGTATCTTCAATCACGTCATGCAGCAGCGCCGCCATCCGCGTTTCATAGTCGAGTTTCATCTCGGCCAGAATGCAGGCAACCGCTACCGGGTGCGTGATATAGGGTTCACCGCTTGAACGTGTTTGCCCCTCGTGAGCATCACGTGCAACGAGATACGCCTGCCGCAGACGCTTGATTTGGT                                              >  2:149334/1‑251                                                                                                                                                                                                             ACCAGCTCGGCGACGCTTTTACCAAAAAGCTGTACCATATCCTGATCCTGGTAGGTGGCGGGAGTATCTTCAATCACGTCATGCAGCAGCGCCGCCATCAGCGTTTCATAGTCGAGTTTCATCTCGGCCAGAATGCAGGCAACCGCTACCGGGTGCGTGATATAGGGTTCACCGCTTGAACGTGTTTGCCCCTCGTGAGCATCACGTGCAACGAGATACGCCTGCCGCAGACGCTTGATTTGGTCTTCCGG                                       <  2:369819/251‑1                                                                                                                                                                                                                                  CCAAAAAGCTGTTCCATATCCTGATCCTGGTAGGTGGCGGGAGTATCTTCAATCACGTCATGCAGCAGCGCCGCCATCAGCGTTTCATAGTCGAGTTTCATCTCGGCCAGAATGCAGGCAACCGCTACCGGGTGCGTGATATAGGGTTCACCGCTTGAACGTGTTTGCCCCTCGTGAGCATCACGTGCAACGAGATACGCCTGCCGCAGACGCTTGATTTGGTCTTCCGGCAGGTAGGTTTGAATCAGTTG                  <  1:360094/251‑1                                                                                                                                                                                                                                       AAGCTGTTCCATATCCTGATCCTGGTAGGTGGCGGGAGTATCTTCAATCACGTCATGCAGCAGCGCCGCCATCAGCGTTTCATAGTCGAGTTTCATCTCGGCCAGAATGCAGGCAACCGCTACCGGGTGCGTGATATAGGGTTCACCGCTTGAACGTGTTTGCCCCTCGTGAGCATCACGTGCAACGAGATACGCCTGCCGCAGACGCTTGATTTGGTCTTCCGGCAGGTAGGTTTGAATCAGTTGATTC              >  1:348642/1‑250                                                                                                                                                                                                                                            GTTCCATATCCTGATCCTGGTAGGTGGCGGGAGTATCTTCAATCACGTCATGCAGCAGCGCCGCCATCAGCGTTTCATAGTCGAGTTTCATCTCGGCCAGAATGCAGGCAACCGCTACCGGGTGCGTGATATAGGGTTCACCGCTTGAACGTGTTTGCCCCTCGTGAGCATCACGTGCAACGAGATACGCCTGCCGCAGACGCTTGATTTGGTCTTCCGGCAGGTAGGTTTTAATCAGTTGATTCAGGCTT        >  2:30948/1‑251                                                                                                                                                                                                                                            GTTCCATATCCTGATCCTGGTAGGTGGCGGGAGTATCTTCAATCACGTCATGCAGCAGCGCCGCCATCAGCGTTTCATAGTCGAGTTTCATCTCGGCCAGAATGCAGGCAACCGCTACCGGGTGCGTGATATAGGGTTCACCGCTTGAAC                                                                                                             >  2:318115/1‑150                                                                                                                                                                                                                                            GTTCCATATCCTGATCCTGGTAGGTGGCGGGAGTATCTTCAATCACGTCATGCAGCAGCGCCGCCATCAGCGTTTCATAGTCGAGTTTCATCTCGGCCAGAATGCAGGCAACCGCTACCGGGTGCGTGATATAGGGTTCACCGCTTGAAC                                                                                                             <  1:318115/150‑1                                                                                                                                                                                                                                                   ATCCTGATCCTGGTAGGTGGCGGGAGTATCTTCAATCACGTCATGCAGCAGCGCCGCCATCAGCGTTTCATAGTCGAGTTTCATCTCGGCCAGAATGCAGGCAACCGCTACCGGGTGCGTGATATAGGGTTCACCGCTTGAACGTGTTTGCCCCTCGTGAGCATCACGTGCAACGAGATACGCCTGCCGCAGACGCTTGATTTGGTCTTCCGGCAGGTAGGTTTGAATCAGTTGATTCAGGCTTTCAAAC  >  2:256161/1‑250                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                              ATTTCGAGAGTTTCACGGGCGATGCGGCGACGTTTGTCCGGGCGAAGTGAGCCCAGCGTGCGCATGTTGTGGGTACGGTCGGCAAGTTTGATGAGGATGACGCGGATATCCTGCACCATCGCCATAATCATCTTGCGAAAGTTTTCGGCCTGCGCCTCTTTCTTATCGCGGAACTTGAGTTTATCAAGTTTCGACACCCCCTCTACCAGCTCGGCGACGCTTTTACCAAAAAGCTGTTCCATATCCTG‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  <  NC_000913/3822893‑3822646 ‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑ATCCTGGTAGGTGGCGGGAGTATCTTCAATCACGTCATGCAGCAGCGCCGCCATCAGCGTTTCATAGTCGAGTTTCATCTCGGCCAGAATGCAGGCAACCGCTACCGGGTGCGTGATATAGGGTTCACCGCTTGAACGTGTTTGCCCCTCGTGAGCATCACGTGCAACGAGATACGCCTGCCGCAGACGCTTGATTTGGTCTTCCGGCAGGTAGGTTTGAATCAGTTGATTCAGGCTTTCAAAC  <  NC_000913/3822651‑3822408

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