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 =0 (0.000)20 (1.140) 17/292 0.1 100% coding (247/2109 nt) spoT bifunctional (p)ppGpp synthase/hydrolase SpoT
?NC_000913 = 3822651 0 (0.000)coding (252/2109 nt) spoT bifunctional (p)ppGpp synthase/hydrolase SpoT

AAAGCCTGAATCAACTGATTCAAACCTACCTGCCGGAAGACCAAATCAAGCGTCTGCGGCAGGCGTATCTCGTTGCACGTGATGCTCACGAGGGGCAAACACGTTCAAGCGGTGAACCCTATATCACGCACCCGGTAGCGGTTGCCTGCATTCTGGCCGAGATGAAACTCGACTATGAAACGCTGATGGCGGCGCTGCTGCATGACGTGATTGAAGATACTCCCGCCACCTACCAGGATATGGAACAGCTTTTTGGTAAAAGCGTCGCCGAGCTGGTAGAGGGGGTGTCGAAACTTGATAAACTCAAGTTCCGCGATAAGAAAGAGGCGCAGGCCGAAAACTTTCGCAAGATGATTATGGCGATGGTGCAGGATATCCGCGTCATCCTCATCAAACTTGCCGACCGTACCCACAACATGCGCACGCTGGGCTCACTTCGCCCGGACAAACGTCG  >  NC_000913/3822413‑3822866
                                                                                                                                                                                                                                              |                                                                                                                                                                                                                       
tgaatcaactgattcaaacctacctgccggaagaccaaatcaagcgtctgcggcaggcgtatcgcgttgcacgtgatgctcacgaggggcaaacacgttcaagcggtgaaccctatatcacgcacccggtagcgggtgcctgcattctggccgagatgaaactcgactatgaaacgcggatggcggcgctgctgcatgacgtgattgaagatgctcccgccacctaccaggatCAGGATATgg                                                                                                                                                                                                                     >  1:179767‑M1/234‑243 (MQ=255)
      tGAATCAACTGATTCAAACCTACCTGCCGGAAGACCAAATCAAGCGTCTGCGGCAGGCGTATCGCGTTGCACGTGATGCTCACGAGGGGCAAACACGTTCAAGCGGTGAACCCTATATCACGCACCCGGTAGCGGGTGCCTGCATTCTGGCCGAGATGAAACTCGACTATGAAACGCGGATGGCGGCGCTGCTGCATGACGTGATTGAAGATGCTCCCGCCACCTACCAGGATcaggatatgg                                                                                                                                                                                                               >  1:179767‑M2/1‑233 (MQ=255)
          acaACTGATTCAAACCTCCCTGCCGGAAGACCAAATCAAGCGTCTGCGGCAGGCGTATCTCGTTGCACGTGATGCTCACGAGGGGCAAACACGTTCAAGCGGTGAACCCTATATCACGCACCCGGTAGCGGTTGCCTGCATTCTGGCCGAGATGAAACTCGACTATGAAACGCTGATGGCGGCGCTGCTGCATGACGTGATTGAAGATACTCCCGCCACCTACCAGGATcagga                                                                                                                                                                                                                    >  1:32342‑M2/2‑229 (MQ=255)
                         gccggaagaccaaatcaagcgtctgcggcaggcgtatctcgttgcacgtgatgctcacgaggggcaaacacgttcaagcggtgaaccctatatcacgcacccggtagcggttgcctgcattctggccgagatgaaactcgactatgaatcgctgatggcggcgctgctgcatgacgtgattgaagatactcccgccacctaccaggatCAGGATATGGAACAGCTTTTTGGTAAAAGCGTCGCCGAGCTgg                                                                                                                                                                                    <  1:450955‑M1/43‑1 (MQ=255)
                               gCCGGAAGACCAAATCAAGCGTCTGCGGCAGGCGTATCTCGTTGCACGTGATGCTCACGAGGGGCAAACACGTTCAAGCGGTGAACCCTATATCACGCACCCGGTAGCGGTTGCCTGCATTCTGGCCGAGATGAAACTCGACTATGAATCGCTGATGGCGGCGCTGCTGCATGACGTGATTGAAGATACTCCCGCCACCTACCAGGATcaggatatggaacagctttttggtaaaagcgtcgccgagctgg                                                                                                                                                                              <  1:450955‑M2/251‑44 (MQ=255)
                                           gcgtctgcggcatgcgtatctctttgcacgtgatgctcacgaggggcaaacccgttcaagcggtgaaccctatatcacgcacccggtagcggttgcctgcattctggccgagatgaaactcgattatgaaacgctgagggcggcgctgctgcatgacgtgattgaagatactcccgccacctaccaggatCAGGATATGGAACAGCTTTTTGGTAAAAGCGTCGCCGAGCTGGTAGAGTGGGTGTCGAAAc                                                                                                                                                                  >  1:307712‑M1/191‑251 (MQ=255)
                                                 gCGTCTGCGGCATGCGTATCTCTTTGCACGTGATGCTCACGAGGGGCAAACCCGTTCAAGCGGTGAACCCTATATCACGCACCCGGTAGCGGTTGCCTGCATTCTGGCCGAGATGAAACTCGATTATGAAACGCTGAGGGCGGCGCTGCTGCATGACGTGATTGAAGATACTCCCGCCACCTACCAGGATcaggatatggaacagctttttggtaaaagcgtcgccgagctggtagagtgggtgtcgaaac                                                                                                                                                            >  1:307712‑M2/1‑190 (MQ=255)
                                                        acgtatctcgttgcacgtgatgctcacgcggttcaaacacgttcaagcggagaacactatatcacgcacccggtagcggttgcctgcattctggcagagatgacactagactatgaacagctgatggcggcgctgctggatgacgtgattgaagatactcccgccacctaccaggatCAGGCTATGAAACCACTTTTTGGTAAAAGCGTCGCCGCGCTGGTAGAGGGGGTGTCGAAACTTGATAAACTCa                                                                                                                                                      <  1:338839‑M1/73‑1 (MQ=255)
                                                              acgTATCTCGTTGCACGTGATGCTCACGCGGTTCAAACACGTTCAAGCGGAGAACACTATATCACGCACCCGGTAGCGGTTGCCTGCATTCTGGCAGAGATGACACTAGACTATGAACAGCTGATGGCGGCGCTGCTGGATGACGTGATTGAAGATACTCCCGCCACCTACCAGGATcaggctatgaaaccactttttggtaaaagcgtcgccgcgctggtagagggggtgtcgaaacttgataaactca                                                                                                                                                <  1:338839‑M2/249‑74 (MQ=255)
                                                                                                         gtgaaccctatatcacgcacccggtagcggttgcctgctatatggccgagatgaaaatcgagtctgacaagctgatggcggcgctgctgcatgccgtgattgaagatactcccgccacctcacaggatCAGGAAATGGAACAGCTTTTTGGTAAAAGCGTCGCCGAGCTGGAAGAGGGGGTGTCGAAACTTGATAAAc                                                                                                                                                         <  1:241251‑M1/70‑1 (MQ=255)
                                                                                                             accctatatcacgcacccggtagcggttgcctgcattctggccgagatgaaactcgactatgaaacgctgatggcggcgctgctgcatgacgtgattgatgatactcccgccacctaccaggatCAGGATATGGAACAGCTTTTTGGTTAAAGCGTCGCCGAGCTGGTAGAGGGGGTGTCGAAACTTGATAAACTCAAGTTCCGCGATAAGAAAGAGGCGCAGGCCGAAAACTTTCGCAAgat                                                                                                        <  1:439527‑M1/119‑1 (MQ=255)
                                                                                                               gTGAACCCTATATCACGCACCCGGTAGCGGTTGCCTGCTATATGGCCGAGATGAAAATCGAGTCTGACAAGCTGATGGCGGCGCTGCTGCATGCCGTGATTGAAGATACTCCCGCCACCTCACAGGATcaggaaatggaacagctttttggtaaaagcgtcgccgagctggaagagggggtgtcgaaacttgataaac                                                                                                                                                   <  1:241251‑M2/198‑71 (MQ=255)
                                                                                                                   aCCCTATATCACGCACCCGGTAGCGGTTGCCTGCATTCTGGCCGAGATGAAACTCGACTATGAAACGCTGATGGCGGCGCTGCTGCATGACGTGATTGATGATACTCCCGCCACCTACCAGGATcaggatatggaacagctttttggttaaagcgtcgccgagctggtagagggggtgtcgaaacttgataaactcaagttccgcgataagaaagaggcgcaggccgaaaactttcgcaagat                                                                                                  <  1:439527‑M2/243‑120 (MQ=255)
                                                                                                                                        ttcctgcattctggccgagatgaaactcgactatgaagcgctgatggcggagatgccgaatgacgtgattgaagatactcctgccacctaccaggatCAGGATATGGAAAAGCTTTTTGGTAAAAGCGTCg                                                                                                                                                                                             <  1:420781‑M1/34‑1 (MQ=255)
                                                                                                                                         gcctgcattctggccgagatgaaactcgactatgaaacgctgatggcggcgcggctgcatgacgtgatggaagatactcccgccacctaccaggatCAGGATATGGAACAGCTTTTTGGTAAAAGCGTCGCCGAGCTGGTAGAggggg                                                                                                                                                                           >  1:240314‑M1/97‑148 (MQ=255)
                                                                                                                                          cctgcattctggccgagatgaaactcgactatgaaacgctgctggcggcgctgctgcatgacgtgattgaagatactcccgccacctaccaggatCAGGATATGGAACAGCTTTTTGGTAAAAGCGTCGCCGAGCTGGTAGAGGGGGTGTCGAAACTTGATAAACTCAAGTTCCGCGATAAGAAAGAGGCGCAGGCCGAAAACTTTCGCAAGATGATTATGGCGATGGTGCAGGATATc                                                                               <  1:444749‑M1/144‑1 (MQ=255)
                                                                                                                                           ctgcattctggccgagatgaaactcgactatgaaaagctgatggcggcgcagctgcatgacgtgattgacgatactcccgccacctaccaggatCAGGATATGGAACAGCTTTTTGGTAAAAGCGTCGCCGAGCTGGTAGAGGGGGTGTc                                                                                                                                                                       <  1:396718‑M1/56‑1 (MQ=255)
                                                                                                                                              ttcctgCATTCTGGCCGAGATGAAACTCGACTATGAAGCGCTGATGGCGGAGATGCCGAATGACGTGATTGAAGATACTCCTGCCACCTACCAGGATcaggatatggaaaagctttttggtaaaagcgtcg                                                                                                                                                                                       <  1:420781‑M2/129‑35 (MQ=255)
                                                                                                                                               gCCTGCATTCTGGCCGAGATGAAACTCGACTATGAAACGCTGATGGCGGCGCGGCTGCATGACGTGATGGAAGATACTCCCGCCACCTACCAGGATcaggatatggaacagctttttggtaaaagcgtcgccgagctggtagaggggg                                                                                                                                                                     >  1:240314‑M2/1‑96 (MQ=255)
                                                                                                                                                ccTGCATTCTGGCCGAGATGAAACTCGACTATGAAACGCTGCTGGCGGCGCTGCTGCATGACGTGATTGAAGATACTCCCGCCACCTACCAGGATcaggatatggaacagctttttggtaaaagcgtcgccgagctggtagagggggtgtcgaaacttgataaactcaagttccgcgataagaaagaggcgcaggccgaaaactttcgcaagatgattatggcgatggtgcaggatatc                                                                         <  1:444749‑M2/239‑145 (MQ=255)
                                                                                                                                                 cTGCATTCTGGCCGAGATGAAACTCGACTATGAAAAGCTGATGGCGGCGCAGCTGCATGACGTGATTGACGATACTCCCGCCACCTACCAGGATcaggatatggaacagctttttggtaaaagcgtcgccgagctggtagagggggtgtc                                                                                                                                                                 <  1:396718‑M2/150‑57 (MQ=255)
                                                                                                                                                    ggccgagatgaaactcgactatgaaacgctgatggcggcgctgctgcatgacgtgattgaagatactcccgccacctaccaggatCAGGATATGGAACAGCTTTTTGGTAAAAGCGTCGCCGAGCTGGTAGAGGGGGTGTCGAAACTTGATAAACTCAAGTTCcgc                                                                                                                                              >  1:157523‑M1/86‑166 (MQ=255)
                                                                                                                                                          ggCCGAGATGAAACTCGACTATGAAACGCTGATGGCGGCGCTGCTGCATGACGTGATTGAAGATACTCCCGCCACCTACCAGGATcaggatatggaacagctttttggtaaaagcgtcgccgagctggtagagggggtgtcgaaacttgataaactcaagttccgc                                                                                                                                        >  1:157523‑M2/1‑85 (MQ=255)
                                                                                                                                                                                            ctgctgcatgacgtgattgaagatactcccgccacctaccaggatCAGGATATGGAACAGCTTTTTGGTAAAAGCGTCGCCGAGCTGGTAGAGGGGGTGTCGAAACTTGATAAACTCAAGTTCCGCGAtaagan                                                                                                                                      <  2:271130‑M1/89‑2 (MQ=255)
                                                                                                                                                                                             gatgtattgacgtgattgaagatactcccgccacctaccaggatCAGGATATGGAACAGCTTTTTGGTAAAAGCGTCGCCGAGCTGGTAGAGGGGGTGTCGAAACTTGATAAACTCAAGTTCCGCGATAAGAAAGAGGCGCAGGCCGAAAACTTTCGCAAGATGATTATGGCGATGGTGCAGGATATCCGCGTCATCCTCATCAAACTTGc                                                        >  1:222163‑M1/45‑211 (MQ=255)
                                                                                                                                                                                              gctgcatgacgtgattgaagatactcccgccacctaccaggatCAGGATATGGAACAGCTTTTTAGTAAAAGCGTCGCCGAGCTGGTAGAGGGGGTGTCGTAACTTGAGAAACTCAAGTTCCGCTATAAGAAAGAGGCGCAGGACGAAAACTGTCGCAAGATGATGACGGCGGTGGGGCAGGATATCCGCGCCATCCGCATCAAACTTTCCGACCGTACCCACAACATTCGCACGCTGGGCTCACTTCg                 >  1:437796‑M1/44‑249 (MQ=255)
                                                                                                                                                                                              gctgcatgacgtgattgaagatacgccctccacctaccaggatCAGGATATGGAACAGCTTTTTGGTAAAATCGTCGCCGAGCTGGTAGAggg                                                                                                                                                                             >  1:407748‑M1/44‑93 (MQ=255)
                                                                                                                                                                                               ctgcatgacgtgattgaagatactcccgccacctaccaggatCAGGATATGGAACAGCTTTTTGGTAAAAGCGCCGCCGAGCTGGTAGAGGGGGTGTCGAAACTTGATAAACTCAAGTTCCGCGATAAGAAAGAGGCGCAGGCCGAAAACTTTAGCAAGCTGATTATGGCGATGGTTCAGGATATCCGCGTCATCCTCATCAAACTTGCCGACCGGACCCACAACATGCGCACGCTGTGCTCacttccgat              <  1:322296‑M1/209‑5 (MQ=255)
                                                                                                                                                                                                  ctgctgCATGACGTGATTGAAGATACTCCCGCCACCTACCAGGATcaggatatggaacagctttttggtaaaagcgtcgccgagctggtagagggggtgtcgaaacttgataaactcaagttccgcgataagan                                                                                                                                <  2:271130‑M2/134‑90 (MQ=255)
                                                                                                                                                                                                   gatgtattgacgtgaTTGAAGATACTCCCGCCACCTACCAGGATcaggatatggaacagctttttggtaaaagcgtcgccgagctggtagagggggtgtcgaaacttgataaactcaagttccgcgataagaaagaggcgcaggccgaaaactttcgcaagatgattatggcgatggtgcaggatatccgcgtcatcctcatcaaacttgc                                                  >  1:222163‑M2/8‑44 (MQ=255)
                                                                                                                                                                                                    gctgcATGACGTGATTGAAGATACTCCCGCCACCTACCAGGATcaggatatggaacagctttttagtaaaagcgtcgccgagctggtagagggggtgtcgtaacttgagaaactcaagttccgctataagaaagaggcgcaggacgaaaactgtcgcaagatgatgacggcggtggggcaggatatccgcgccatccgcatcaaactttccgaccgtacccacaacattcgcacgctgggctcacttcg           >  1:437796‑M2/1‑43 (MQ=255)
                                                                                                                                                                                                    gctgcATGACGTGATTGAAGATACGCCCTCCACCTACCAGGATcaggatatggaacagctttttggtaaaatcgtcgccgagctggtagaggg                                                                                                                                                                       >  1:407748‑M2/1‑43 (MQ=255)
                                                                                                                                                                                                     ctgcATGACGTGATTGAAGATACTCCCGCCACCTACCAGGATcaggatatggaacagctttttggtaaaagcgccgccgagctggtagagggggtgtcgaaacttgataaactcaagttccgcgataagaaagaggcgcaggccgaaaactttagcaagctgattatggcgatggttcaggatatccgcgtcatcctcatcaaacttgccgaccggacccacaacatgcgcacgctgtgctcacttccgat        <  1:322296‑M2/251‑210 (MQ=255)
                                                                                                                                                                                                      acgtgattgaagatactcctgccacctaccaggatCAGGATATGGAACAGCTTTTTGGTAAAAGCGTCGCCGAGCTGGTAGCGGGAGCGTCGAAACTTGCTAAACTCAAGTTCCGCGATAAGAAAGACGCGCAGCCCGAAAACTTTCGCAAGATGATTATGGCGATGGTGCAGGATATCCGCGTCATCCTCATCAAACTTGCCGACCGTACCCACAACATGAGCACGCTGGGCTCACTTCgcccggccaa        <  1:392160‑M1/215‑5 (MQ=255)
                                                                                                                                                                                                          nattgaagatactcccgccacctaccaggatCAGGATATGGAACAGCTTTTTGGTAAAAGCGTCGCCGAGCTGGTAGAGGGGGTGTCGAAACTTGATAAACTCAAGTTCCGCGATAAGAAAGAGGCGCAGGcc                                                                                                                         <  1:40566‑M1/102‑1 (MQ=255)
                                                                                                                                                                                                            aCGTGATTGAAGATACTCCTGCCACCTACCAGGATcaggatatggaacagctttttggtaaaagcgtcgccgagctggtagcgggagcgtcgaaacttgctaaactcaagttccgcgataagaaagacgcgcagcccgaaaactttcgcaagatgattatggcgatggtgcaggatatccgcgtcatcctcatcaaacttgccgaccgtacccacaacatgagcacgctgggctcacttcgcccggccaa  <  1:392160‑M2/250‑216 (MQ=255)
                                                                                                                                                                                                                natTGAAGATACTCCCGCCACCTACCAGGATcaggatatggaacagctttttggtaaaagcgtcgccgagctggtagagggggtgtcgaaacttgataaactcaagttccgcgataagaaagaggcgcaggcc                                                                                                                   <  1:40566‑M2/132‑103 (MQ=255)
                                                                                                                                                                                                                                     ggatCAGGATATGGAACAGCTTTTTGGTAAAAGCGTCGCCGAGCTGGTCGAgg                                                                                                                                                                              <  1:366724‑M1/49‑1 (MQ=255)
                                                                                                                                                                                                                                           ggatcaggatatggaacagctttttggtaaaagcgtcgccgagctggtcgagg                                                                                                                                                                        <  1:366724‑M2/53‑50 (MQ=255)
                                                                                                                                                                                                                                              |                                                                                                                                                                                                                       
AAAGCCTGAATCAACTGATTCAAACCTACCTGCCGGAAGACCAAATCAAGCGTCTGCGGCAGGCGTATCTCGTTGCACGTGATGCTCACGAGGGGCAAACACGTTCAAGCGGTGAACCCTATATCACGCACCCGGTAGCGGTTGCCTGCATTCTGGCCGAGATGAAACTCGACTATGAAACGCTGATGGCGGCGCTGCTGCATGACGTGATTGAAGATACTCCCGCCACCTACCAGGATATGGAACAGCTTTTTGGTAAAAGCGTCGCCGAGCTGGTAGAGGGGGTGTCGAAACTTGATAAACTCAAGTTCCGCGATAAGAAAGAGGCGCAGGCCGAAAACTTTCGCAAGATGATTATGGCGATGGTGCAGGATATCCGCGTCATCCTCATCAAACTTGCCGACCGTACCCACAACATGCGCACGCTGGGCTCACTTCGCCCGGACAAACGTCG  >  NC_000913/3822413‑3822866

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