Predicted mutation
evidence seq id position mutation annotation gene description
JC NC_002947 4,980,585 +GGC intergenic (‑24/+44) PP_4387 ← / ← flgE hypothetical protein/flagellar hook protein FlgE

New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_002947 = 49805851 (0.030)28 (0.910)
+GGC
27/284 0.2 96.6% intergenic (‑24/+44) PP_4387/flgE hypothetical protein/flagellar hook protein FlgE
?NC_002947 4980586 = 1 (0.030)intergenic (‑25/+43) PP_4387/flgE hypothetical protein/flagellar hook protein FlgE

GCCGGTGACCAGGCTTTGCCATTGATAGAACGGGGCGGGGGCAGCGGGGGCAGGCAGAACCAGCACCAGCAAGGGGAACAAGGCCAGTAACGACAGCCGCTTCATTCGCGTGAACCTTCACATTCTTGG‑‑‑CGCTTTGCGGCCCCAATTACAAACACCAAAGCCAATACAATCATCAGGTCATCTGAATGATGGTCTGCATGATGGTGCTTTCAGTGGAAATGGTCTTGGCATTGGCCTGGTAGTTGCTCTGCGCCTTGATCAGGTTTACCA  >  NC_002947/4980457‑4980726
                                                                                                                                    |                                                                                                                                            
gccggtcaccaggctttgccattgatagaccggggcgggggcagcgggggcaggcagaaccagcaccagcaaggggaacaaggccagtaacgacagccgcttcattcgcgtgaaccttcacattcttggggcCGCAAAGCGGCCCCaa                                                                                                                               <  5:68801‑M2/16‑1 (MQ=255)
  cggggcccaggctctgccattccgacaccggggcgggggcagcgggggcaggcagacccagcaccagcaaggggaacaaggccagtaacgacagccgcttcattcgcgtgaaccttcacattcttggggcCGCTTTGCGGCCCCAAtt                                                                                                                             <  5:164935‑M2/18‑1 (MQ=255)
         caggctttgccattgatagaacggggcgggggcagcgggggcaggcagaaccagcaccagcaaggggaacaaggccagtaacgacagccgcttcattcgcgtgaaccttcacattcttggggcCGCTTTGCGGCCCCAATTACAaacc                                                                                                                      >  1:60127‑M2/124‑147 (MQ=255)
         caggctttgccattgatagaacggggcgggggcagcgggggcaggcagaaccagcaccagcaaggggaacaaggccagtaacgacagccgcttcattcgcgtgaaccttcacattcttggggcCGCTTTGCGGCCCCAATTACAAaca                                                                                                                      >  1:151395‑M2/124‑148 (MQ=255)
              tttgccattgctagaccggggcgggggcagcgggggcaggcagaaccagcaccagcaaggggaacaaggccagtaacgacagccgcttcattcgcgtgaaccttcacattcttggggcCGCTTTGCGGCCCCAATTACAAACACCaaa                                                                                                                 <  2:46627‑M2/30‑1 (MQ=255)
                    attgatagaacggggccggggcagcgggggcaggcagaaccagcaccagcaaggggaacaagcccagtaacgacagccgcttcattcgagtgaaccttcacattcttggcgcCGCTATGCGGCCCCAATTACAAACACCAAACCCAAt                                                                                                           >  4:65188‑M2/113‑148 (MQ=255)
                       gatagaacggggcgggggcagcgggggcaggcagaaccagcaccagcaaggggaacaaggccagtaacgacagccgcttcattcgcgtgaacgtacacattcttggggcCGCTTTGCGGCCCCAAATACAAACACCAAAGCCAATAc                                                                                                         >  5:154595‑M2/110‑147 (MQ=255)
                        atagaacggggcgggggcagcgggggcaggcagaaccagcaccagcaaggggaacaaggccagtaacgacagccgcttcattcgcgtgaaccttcacattcttggggcCGCTTTGCGGCCCCAATTACAAACACCAAAGCCAATACa                                                                                                        >  6:21939‑M2/109‑147 (MQ=255)
                         tagatcggggcggggccagcggggcacggcagaacaagcagaagcaaggggaacaaggccagtaacgacagccgcttcattcgcgtgaagctagacattcttggttcCGCTTTGCGGCCTCAATTACAAACACCAAAGCCAATACAat                                                                                                      <  6:154595‑M2/41‑1 (MQ=255)
                          agaacggggcgggggcagcgggggcaggcagaaccagcaccagcaaggggaacaaggccagtaacgacagccgcttcattcgcgtgaaccttcacattcttggggcCGCTTTGCGGCCCCAATTACAAACACCAAAGCCAATACAatc                                                                                                     >  7:122055‑M2/107‑148 (MQ=255)
                              cggggcgggggcagcgggggcaggcagaaccagcaccagcaaggggaacaaggccagtaacgacagccgcttcattcgcgcgaaccttcacattcttggggcCGCTTTGCGGCCCCAATTACAAACACCAAAGCCAATACAAtcatca                                                                                                 <  3:53466‑M2/46‑1 (MQ=255)
                                                ggcaggcagaaccagcaccagcaaggggaacaaggccagtaacgacagccgcttcattcgcgtgaaccttcacattcttggggcCGCTTTGCGGCCCCAATTACAAACACCAAAGCCAATACAATCATCAGGTCATCTGAATGATGGt                                                                               <  5:21939‑M2/64‑1 (MQ=255)
                                                          accagcaccagcaaggggaacaaggccagtaacgacagccgcttcattcgcgtgaaccttcacattcttggggcCGCTTTGCGGCCCCAATTACAAACACCAAAGCCAATACAATCATCAGGTCATCTGAATGATGGTCTGCatgatg                                                                     <  5:97961‑M2/74‑1 (MQ=255)
                                                            cagcaccagcaAGGGGAACAAGGCCAGTAACGACAGCCGCTTCATTCGCGTGAACCTTCACATTCTTGGGGCCGCGTTGAGTCCCCCAAGACAAACACCCAACCCAATACAAACAACACGTCATCTTAAAGATGGTCTGCattcttgt                                                                   >  1:178953/1‑142 (MQ=255)
                                                               caccagcaaggggaacaaggccagtaacgacagccgcttcattcgcgtgaaccttcacattcttggggcCGCTTTGCGGCCCCAATTACAAACACCAAAGCCAATACAATCATCAGGTCATCTGAATGATGGTCTGCATGATGGTGCt                                                                <  8:27546‑M2/79‑1 (MQ=255)
                                                                accagcaaggggaacaaggccagtaacgacagccgcttcattcgcgtgaaccttcacattcttggggcCGCTTTGCGGCCCCAATTACAAACACCAAAGCCAATACAATCATCAGGTCATCTGAATGATGGTCTGCATGATGGTGCt                                                                >  7:25134‑M2/69‑147 (MQ=255)
                                                                  cagcaaggggaacaaggccagtaacgacagccgcttcattcgcgtgaaccttcacattcttggggcCGCTTTGCGGCCCCAATTACAAACACCAAAGCCAATACAATCATCAGGTCATCTGAATGATGGTCTGCATGATGGTGCTTTc                                                             <  4:395‑M2/82‑1 (MQ=255)
                                                                        ggggaacaaggccagtaacgacagccgcttcattcgcgtgaaccttcacattcttggggcCGCGTTGCGGCCCCAATTACAAACACCAAAGCCAATACAATAATCAGGTCATCTGAATGATGGTCTGCATGATGGTGCTTTCAGTgg                                                        >  5:16647‑M2/61‑147 (MQ=255)
                                                                                 ggccagtaacgacagccgcttcattcgcgtgaaccttcacattcttggggcCGCTTTTCGGCCCCAATTACAAACACCAAAGCCAATACAATCATCAGGTCATCTGAATGATGGTCTGCATGATGGTGCTTTCAGTGGAAATGGTCtt                                              >  8:126642‑M2/52‑148 (MQ=255)
                                                                                   ccagtaacgacagccgcttcattcgcgtgaaccttcacattcttggggcCGCTTTGCGGCCCCAATTACAAACACCAAAGCCAATACAATCATCAGGTCATCTGAATGATGGTCTGCATGATGGTGCTTTCAGTGGAAATGGTCTTgg                                            <  1:51804‑M2/99‑1 (MQ=255)
                                                                                      gtaacgacagccgcttcattcgcgtgaaccttcacattcttggggcCGCTTTGCGGCCCCAATTACAAACACCAAAGCCAATACAATCATCAGGTCATCTGAATGATGGTCTGCATGATGGTGCTTTCAGTGGAAATGGTCTTGGCAt                                         <  2:153141‑M2/102‑1 (MQ=255)
                                                                                       taacgacagccgcttcattcgcgtgaaccttcacattcttggggcCGCTTTGCGGCCCCAATTACAAACACCAAAGCCAATACAATCATCAGGTCATCTGAATGATGGTCTGCATGATGGTGCTTTCAGTGGAAATGGTCTTGGCAtt                                        >  4:14887‑M2/46‑148 (MQ=255)
                                                                                                  gcttcattcgcgtgaaccttcacattcttggggcCGCTTTGCGGCCCCAATTACAAACACCAAAGCCAATACAATCATCAGGTCATCTGAATGATGGTCTGCATGATGGTGCTTTCAGTGGAAATGGTCTTGGCATTGGCCTGGTAGt                             >  2:26808‑M2/35‑148 (MQ=255)
                                                                                                    ttcattcgcgtgaaccttcacattcttggggcCGCTTTGCGGCCCCAATTACAAACACCAAAGCCAATACAATCATCAGGTCATCTGAATGATGGTCTGCATGATGGTGCTTTCAGTGGAAATGGTCTTGGCATTGGCCTGGTAGTTg                           >  3:113585‑M2/33‑148 (MQ=255)
                                                                                                                 accttcacattcttggggcCGCTTTGCGGCCCCAATTACAAACACCAAAGCCAATACAATCATCAGGTCATCTGAATGATGGTCTGCATGATGGTGCTTTCAGTGGAAATGGTCTTGGCATTGGCCTGGTAGTTGCTCTGCGCCTTg               <  2:151395‑M2/128‑1 (MQ=255)
                                                                                                                   cttcacattcttggggcCGCGTTGCGGCCCCAATTACAAACACCAAAGCCAATACAATCATCAGGTCATCTGAATGATGGTCTGCATGATGGTGCTTTCAGTGGAAATGGTCTTGGCATTGGCCTGGTAGTTGCTCTGCGCCTTGa              >  1:18825‑M2/18‑146 (MQ=255)
                                                                                                                    ttcacattcttggggcCGCTTTGCGGCCCCAATTACAAACACCAAAGCCAATACAATCATCAGGTCATCTGAATGATGGTCTGCATGATGGTGCTTTCAGTGGAAATGGTCTTGGCATTGGCCTGGTAGTTGCTCTGCGCCTTGATCa           >  8:7996‑M2/17‑148 (MQ=255)
                                                                                                                             ttggggcCGCTTTGCGGCCCCAATTACAAACACCAAAGCCAATACAATCATCAGGTCATCTGAATGATGGTCTGCATGATGGTGCTTTCAGTGGAAATGGTCTTGGCATTGGCCTGGTAGTTGCTCTGCGCCTTGATCAGGTTTACCa  >  8:16459‑M2/8‑148 (MQ=255)
                                                                                                                                gggcCGCTTTGCGGCCCCAATTACAAACACCAAAGCCAATACAATCATCAGGTCATCTGAATGATGGTCTGCATGATGGTGCTTTCAGTGGAAATGGTCTTGGCATTGGCCTGGTAGTTGCTCTGCGCCTTGATCAGGTTTa     >  6:157563‑M2/5‑142 (MQ=255)
                                                                                                                                 ggcCGCTTTGCGGCCCCAATTACAAACACCAAAGCCAATACAATCATCAGGTCATCTGAATGATGGTCTGCATGATGGTGCTTTCAGTGGAAATGGTCTTGGCATTGGCCTGGTAGTTGCTCTGCGCCTTGATCAGGTTTAc    <  5:157563‑M2/139‑1 (MQ=255)
                                                                                                                                    |                                                                                                                                            
GCCGGTGACCAGGCTTTGCCATTGATAGAACGGGGCGGGGGCAGCGGGGGCAGGCAGAACCAGCACCAGCAAGGGGAACAAGGCCAGTAACGACAGCCGCTTCATTCGCGTGAACCTTCACATTCTTGG‑‑‑CGCTTTGCGGCCCCAATTACAAACACCAAAGCCAATACAATCATCAGGTCATCTGAATGATGGTCTGCATGATGGTGCTTTCAGTGGAAATGGTCTTGGCATTGGCCTGGTAGTTGCTCTGCGCCTTGATCAGGTTTACCA  >  NC_002947/4980457‑4980726

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