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‑M1/148‑20 (MQ=255)
  cggggcccaggctcTGCCATTCCGACACCGGGGCGGGGGCAGCGGGGGCAGGCAGACCCAGCACCAGCAAGGGGAACAAGGCCAGTAACGACAGCCGCTTCATTCGCGTGAACCTTCACATTCTTGGggccgctttgcggccccaatt                                                                                                                             <  5:164935‑M1/142‑22 (MQ=255)
         cAGGCTTTGCCATTGATAGAACGGGGCGGGGGCAGCGGGGGCAGGCAGAACCAGCACCAGCAAGGGGAACAAGGCCAGTAACGACAGCCGCTTCATTCGCGTGAACCTTCACATTCTTGGggccgctttgcggccccaattacaaacc                                                                                                                      >  1:60127‑M1/1‑120 (MQ=255)
         cAGGCTTTGCCATTGATAGAACGGGGCGGGGGCAGCGGGGGCAGGCAGAACCAGCACCAGCAAGGGGAACAAGGCCAGTAACGACAGCCGCTTCATTCGCGTGAACCTTCACATTCTTGGggccgctttgcggccccaattacaaaca                                                                                                                      >  1:151395‑M1/1‑120 (MQ=255)
              tttGCCATTGCTAGACCGGGGCGGGGGCAGCGGGGGCAGGCAGAACCAGCACCAGCAAGGGGAACAAGGCCAGTAACGACAGCCGCTTCATTCGCGTGAACCTTCACATTCTTGGggccgctttgcggccccaattacaaacaccaaa                                                                                                                 <  2:46627‑M1/148‑34 (MQ=255)
                    aTTGATAGAACGGGGCCGGGGCAGCGGGGGCAGGCAGAACCAGCACCAGCAAGGGGAACAAGCCCAGTAACGACAGCCGCTTCATTCGAGTGAACCTTCACATTCTTGGcgccgctatgcggccccaattacaaacaccaaacccaat                                                                                                           >  4:65188‑M1/1‑109 (MQ=255)
                       gATAGAACGGGGCGGGGGCAGCGGGGGCAGGCAGAACCAGCACCAGCAAGGGGAACAAGGCCAGTAACGACAGCCGCTTCATTCGCGTGAACGTACACATTCTTGGggccgctttgcggccccaaatacaaacaccaaagccaatac                                                                                                         >  5:154595‑M1/1‑106 (MQ=255)
                        aTAGAACGGGGCGGGGGCAGCGGGGGCAGGCAGAACCAGCACCAGCAAGGGGAACAAGGCCAGTAACGACAGCCGCTTCATTCGCGTGAACCTTCACATTCTTGGggccgctttgcggccccaattacaaacaccaaagccaataca                                                                                                        >  6:21939‑M1/1‑105 (MQ=255)
                         tAGATCGGGGCGGGGCCAGCGGGGCACGGCAGAACAAGCAGAAGCAAGGGGAACAAGGCCAGTAACGACAGCCGCTTCATTCGCGTGAAGCTAGACATTCTTGGttccgctttgcggcctcaattacaaacaccaaagccaatacaat                                                                                                      <  6:154595‑M1/148‑45 (MQ=255)
                          aGAACGGGGCGGGGGCAGCGGGGGCAGGCAGAACCAGCACCAGCAAGGGGAACAAGGCCAGTAACGACAGCCGCTTCATTCGCGTGAACCTTCACATTCTTGGggccgctttgcggccccaattacaaacaccaaagccaatacaatc                                                                                                     >  7:122055‑M1/1‑103 (MQ=255)
                              cggggcggggGCAGCGGGGGCAGGCAGAACCAGCACCAGCAAGGGGAACAAGGCCAGTAACGACAGCCGCTTCATTCGCGCGAACCTTCACATTCTTGGggccgctttgcggccccaattacaaacaccaaagccaatacaatcatca                                                                                                 <  3:53466‑M1/148‑50 (MQ=255)
                                                ggcaggcaGAACCAGCACCAGCAAGGGGAACAAGGCCAGTAACGACAGCCGCTTCATTCGCGTGAACCTTCACATTCTTGGggccgctttgcggccccaattacaaacaccaaagccaatacaatcatcaggtcatctgaatgatggt                                                                               <  5:21939‑M1/148‑68 (MQ=255)
                                                          accagcaccagcAAGGGGAACAAGGCCAGTAACGACAGCCGCTTCATTCGCGTGAACCTTCACATTCTTGGggccgctttgcggccccaattacaaacaccaaagccaatacaatcatcaggtcatctgaatgatggtctgcatgatg                                                                     <  5:97961‑M1/148‑78 (MQ=255)
                                                            cagcaccagcaAGGGGAACAAGGCCAGTAACGACAGCCGCTTCATTCGCGTGAACCTTCACATTCTTGGGGCCGCGTTGAGTCCCCCAAGACAAACACCCAACCCAATACAAACAACACGTCATCTTAAAGATGGTCTGCattcttgt                                                                   >  1:178953/1‑142 (MQ=255)
                                                               caccagcaAGGGGAACAAGGCCAGTAACGACAGCCGCTTCATTCGCGTGAACCTTCACATTCTTGGggccgctttgcggccccaattacaaacaccaaagccaatacaatcatcaggtcatctgaatgatggtctgcatgatggtgct                                                                <  8:27546‑M1/148‑83 (MQ=255)
                                                                accagcaAGGGGAACAAGGCCAGTAACGACAGCCGCTTCATTCGCGTGAACCTTCACATTCTTGGggccgctttgcggccccaattacaaacaccaaagccaatacaatcatcaggtcatctgaatgatggtctgcatgatggtgct                                                                >  7:25134‑M1/1‑65 (MQ=255)
                                                                  cagcaAGGGGAACAAGGCCAGTAACGACAGCCGCTTCATTCGCGTGAACCTTCACATTCTTGGggccgctttgcggccccaattacaaacaccaaagccaatacaatcatcaggtcatctgaatgatggtctgcatgatggtgctttc                                                             <  4:395‑M1/148‑86 (MQ=255)
                                                                        ggggAACAAGGCCAGTAACGACAGCCGCTTCATTCGCGTGAACCTTCACATTCTTGGggccgcgttgcggccccaattacaaacaccaaagccaatacaataatcaggtcatctgaatgatggtctgcatgatggtgctttcagtgg                                                        >  5:16647‑M1/1‑57 (MQ=255)
                                                                                 ggCCAGTAACGACAGCCGCTTCATTCGCGTGAACCTTCACATTCTTGGggccgcttttcggccccaattacaaacaccaaagccaatacaatcatcaggtcatctgaatgatggtctgcatgatggtgctttcagtggaaatggtctt                                              >  8:126642‑M1/1‑48 (MQ=255)
                                                                                   ccAGTAACGACAGCCGCTTCATTCGCGTGAACCTTCACATTCTTGGggccgctttgcggccccaattacaaacaccaaagccaatacaatcatcaggtcatctgaatgatggtctgcatgatggtgctttcagtggaaatggtcttgg                                            <  1:51804‑M1/148‑103 (MQ=255)
                                                                                      gTAACGACAGCCGCTTCATTCGCGTGAACCTTCACATTCTTGGggccgctttgcggccccaattacaaacaccaaagccaatacaatcatcaggtcatctgaatgatggtctgcatgatggtgctttcagtggaaatggtcttggcat                                         <  2:153141‑M1/148‑106 (MQ=255)
                                                                                       tAACGACAGCCGCTTCATTCGCGTGAACCTTCACATTCTTGGggccgctttgcggccccaattacaaacaccaaagccaatacaatcatcaggtcatctgaatgatggtctgcatgatggtgctttcagtggaaatggtcttggcatt                                        >  4:14887‑M1/1‑42 (MQ=255)
                                                                                                  gCTTCATTCGCGTGAACCTTCACATTCTTGGggccgctttgcggccccaattacaaacaccaaagccaatacaatcatcaggtcatctgaatgatggtctgcatgatggtgctttcagtggaaatggtcttggcattggcctggtagt                             >  2:26808‑M1/1‑31 (MQ=255)
                                                                                                    ttCATTCGCGTGAACCTTCACATTCTTGGggccgctttgcggccccaattacaaacaccaaagccaatacaatcatcaggtcatctgaatgatggtctgcatgatggtgctttcagtggaaatggtcttggcattggcctggtagttg                           >  3:113585‑M1/1‑29 (MQ=255)
                                                                                                                 aCCTTCACATTCTTGGggccgctttgcggccccaattacaaacaccaaagccaatacaatcatcaggtcatctgaatgatggtctgcatgatggtgctttcagtggaaatggtcttggcattggcctggtagttgctctgcgccttg               <  2:151395‑M1/147‑132 (MQ=255)
                                                                                                                   cTTCACATTCTTGGggccgcgttgcggccccaattacaaacaccaaagccaatacaatcatcaggtcatctgaatgatggtctgcatgatggtgctttcagtggaaatggtcttggcattggcctggtagttgctctgcgccttga              >  1:18825‑M1/1‑14 (MQ=255)
                                                                                                                    ttCACATTCTTGGggccgctttgcggccccaattacaaacaccaaagccaatacaatcatcaggtcatctgaatgatggtctgcatgatggtgctttcagtggaaatggtcttggcattggcctggtagttgctctgcgccttgatca           >  8:7996‑M1/1‑13 (MQ=255)
                                                                                                                             ttGGggccgctttgcggccccaattacaaacaccaaagccaatacaatcatcaggtcatctgaatgatggtctgcatgatggtgctttcagtggaaatggtcttggcattggcctggtagttgctctgcgccttgatcaggtttacca  >  8:16459‑M1/1‑4 (MQ=255)
                                                                                                                                gggccgctttgcggccccaattacaaacaccaaagccaatacaatcatcaggtcatctgaatgatggtctgcatgatggtgctttcagtggaaatggtcttggcattggcctggtagttgctctgcgccttgatcaggttta     >  6:157563‑M1/1‑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.