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)23 (0.670)
+GGC
20/284 0.9 95.9% 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

CGGCCGGTGACCAGGCTTTGCCATTGATAGAACGGGGCGGGGGCAGCGGGGGCAGGCAGAACCAGCACCAGCAAGGGGAACAAGGCCAGTAACGACAGCCGCTTCATTCGCGTGAACCTTCACATTCTTGG‑‑‑CGCTTTGCGGCCCCAATTACAAACACCAAAGCCAATACAATCATCAGGTCATCTGAATGATGGTCTGCATGATGGTGCTTTCAGTGGAAATGGTCTTGGCATTGGCCTGGTAGTTGCTCTGCGCCTTGATCAGG  >  NC_002947/4980455‑4980719
                                                                                                                                      |                                                                                                                                     
cggccggtgaccaggctttgccattgatagaccggggcgggggcagcgggggcaggcagaaccagcaccagcaaggggaacaaggccagtaacgacagccgcttcattcgcgtgaaccttcacattcttggggcCGCTTTGCGGccc                                                                                                                           <  7:133383‑M2/13‑1 (MQ=255)
  gccggtgaccaggctttgccatggatagaacggggcgggggcagcgggggcaggcagaaccagcaccagcaaggggaacaaggccagtaacgacagccgcttcattcgcgtgaaccttcacattcttggggcCGCAAAGCGGCCCCaa                                                                                                                        <  4:95802‑M2/16‑1 (MQ=255)
                      attgatagaacggggcgggggcagcgggggcaggcagaaccagcaccagcaaggggaacaaggccagtaacgacagccgcttcattcgcgtgaaccttcacattcttggggcCGCTTTGCGGCCCCAATTACAAACACCAAAGCCaa                                                                                                     >  8:27732‑M2/113‑147 (MQ=255)
                      attgatagaacggggcgggggcagcgggggcaggcagaaccagcaccagcaaggggaacaaggccagtaacgacagccgcttcattcgcgtgaaccttcacattcttggggcCGCTTTGCGGCCCCAATTACAAACACCAAAGCCAAt                                                                                                    >  3:175493‑M2/113‑148 (MQ=255)
                      attgatagaacggggcgggggcagcgggggcaggcagaaccagcaccagcaaggggaacaaggccagtaacgacagccgcttcattcgcgtgaaccttcacattcttggggcCGCTTTGCGGCCCCAATTACAAACACCAAAGCCAAt                                                                                                    >  4:109978‑M2/113‑148 (MQ=255)
                           tagacgggggcgggggcagcgggggcaggcagaaccagcaccagcaaggggaacaaggccagtaacgacagccgcttcattcgcgtgaaccttcacattcttggggcCGCTTCGCGGCCCCAATTTCAAACACCAAAGCCAATACaa                                                                                                <  5:120865‑M2/40‑1 (MQ=255)
                               acggggcgggggcagcgggggcaggcagaaccagcaccagcaaggggaacaaggccagtaacgacagccgcttcattcgcgtgaaccttcacattcttgcggcCGCTTTGCGGCCCCAATTACAAACACCAAAGCCAAAACAatca                                                                                             >  4:68726‑M2/104‑146 (MQ=255)
                                 ggggcgggggcagcgggggcaggcagaaccagcaccagcaaggggaaccaggccagtaacgacagccgcttcattcgcgtgccccttcacattcttggggcCGCTTTGCGGCCCCAATTACAAACACCAAAGCCAATACAATCATCAg                                                                                         >  5:66356‑M2/102‑148 (MQ=255)
                                 ggggcgggggcagcgggggcaggcagaaccagcaccagcaaggggaacaaggccagtaacgacagccgcttcattcgcgtgaaccttcacattcttggggcCGCTTTGCGGCCCCAATTACAAACACCAAAGCCAATACCATCATCAg                                                                                         >  5:69435‑M2/102‑148 (MQ=255)
                                 ggggcgggggcagcgggggcaggcagaaccagcaccagcaaggggaacaaggccagtaacgacagccgcttcattcgcgtgaaccttcacattcttggggcCGCTTTGCGGCCCCAATTACAAACACCAAAGCCAATACAATCATCAg                                                                                         <  2:167047‑M2/47‑1 (MQ=255)
                                                ggggcaggcagaaccagcaccagcaaggggaacaaggccagtaacgacagccgcttcattcgcgtgaaccttcacattcttggggcCGCTTTGCGGCCCCAATTACAAACACCAAAGCCAATACAATCATCAGGTCATCTGAatgatg                                                                          <  4:94449‑M2/62‑1 (MQ=255)
                                                            accagcaccagcaaggggaacaaggccagtaacgacagccgcttcattcgcgtgaaccttcacattcttggggcCGCTTTGCGGCCCCAATTACAAACACCAAAGCCAATACAATCATCAGGTCATCTGAATGATGGTCTGCatgatg                                                              <  3:145795‑M2/74‑1 (MQ=255)
                                                              cagcaccagcaaggggaacaaggccagtaacgacagccgcttcattcgcgtgaaccttcacattcttggggcCGCTTTGCGGCCCCAATTACAAACACCAAAGCCAATACAATCATCAGGTCATCTGAATGATGGTCTGCATGATGGt                                                            <  1:178117‑M2/76‑1 (MQ=255)
                                                                 cagcagcaAGGCGACCAAGGCCTGTAACGCCACCCGCATCATAGGCGTGAACCTTCACATTGATGGGGCCGCTTTGCGGCCCCAATTAAAAACACCAAAGCCAATACATTCATCAGGTCATCTGAATGATTGTCTTCATGATGGGGCt                                                         <  4:7744/145‑1 (MQ=255)
                                                                     agcaaggggaacaaggccagtaacgacagccgcttcattcgcgtgaaccttcacattcttggggcCGCTTTGCGGCCCCAATTACAAACACCAAAGCCAATACAATCATCAGGTCATCTGAATGATGGTCTGCATGATGGTGCTTTCa                                                     <  7:27732‑M2/83‑1 (MQ=255)
                                                                              aacaaggccagtaacgacagccgcttcattcgcgtgaaccttcacattcttggggcCGCTTTGCGGCCCCAATTACAAACACCAAAGCCAATACAATCATCAGCTCATCTGAATGATGGTCTGCAAGATGGTGCTTTCAGTGGAAatc                                            >  3:82571‑M2/57‑147 (MQ=255)
                                                                                      cagtaacgacagccgcttcattcgcgtgaaccttcacattcttggggcCGCTTTGCGGCCCCAATTACAAACACCAAAGCCAATACAATCATCAGGTCATCTGAATGATGGTCTGCATGATGGTGCTTTCAGTGGAAATGGTCTTGGc                                    >  8:94647‑M2/49‑148 (MQ=255)
                                                                                         taacgacagccgcttcattcgcgtgaaccttcacattcttggggcCGCTTTGCGGCCCCAATTACAAACACCAAAGCCAATACAATCATCAGGTCATCTGAATGATGGTCTGCATGATGGTGCTTTCAGTGGAAATGGTCTTGGCAtt                                 <  4:54919‑M2/103‑1 (MQ=255)
                                                                                                   cgcttcattcgcgtgaaccttcacattcttggggcCGCTTTGCGGCCCCAATTACAAACACCAAAGCCAATACAATCATCAGGTCATCTGAATGATGGTCTGCATGATGGTGCTTTCAGTGGAAATGGTCTTGGCATTGGCCTGGTAg                       >  5:96805‑M2/36‑148 (MQ=255)
                                                                                                                tgaaccttcacattcttggggcCGCTTTGCGGCCCCAATTACAAACACCAAAGCCAATACAATAATCAGGTCATCTGAATGATGGTCTGCATGATGGTGCTTTCAGTGGAAATGGTCTTGGCATTGGCCTGGTAGTTGCTCTGCGCCt          >  2:55150‑M2/23‑148 (MQ=255)
                                                                                                                 gaaccttcacattcttggggcCGCTTTGCGGCCCCAATTACAAACACCAAAGCCAATACAATCATCAGGTCATCTGAATGATGGTCTGCATGATGGTGCTTTCAGTGGAAATGGTCTTGGCATTGGCCTGGTAGTTGCTCTGCGCCtt         <  3:117047‑M2/127‑1 (MQ=255)
                                                                                                                   accttcacattcttggggcCGCTTTGCGGCCCCAATTACAAACACCAAAGCCAATACAATCATCAGGTCATCTGAATGATGGTCTGCATGATGGTGCTTTCAGTGGAAATGGTCTTGGCATTGGCCTGGTAGTTGCTCTGCGCCTTGa       <  3:154799‑M2/129‑1 (MQ=255)
                                                                                                                    ccttcacattcttggggcCGCTTTGCGGCCCCAATTACAAACACCAAAGCCAATACAATCATCAGGTCATCTGAATGATGGTCTGCATGATGGTGCTTTCAGTGGAAATGGTCTTGGCATTGGCCTGGTAGTTGCTCTGCGCCTTGAt      <  4:42657‑M2/130‑1 (MQ=255)
                                                                                                                         acattcttggggcCGCTTTGCGGCCCCAATTACAAACACCAAAGCCAATACAAGCATCAGGGCATCTGAATGATGGTCTGCATGATGGTGCTTTCAGTGGAAATGGGCTTGGCATTGGCCTGGTAGTTGCTCTGCGCCTTGATCAgg  <  4:28272‑M2/134‑1 (MQ=255)
                                                                                                                                      |                                                                                                                                     
CGGCCGGTGACCAGGCTTTGCCATTGATAGAACGGGGCGGGGGCAGCGGGGGCAGGCAGAACCAGCACCAGCAAGGGGAACAAGGCCAGTAACGACAGCCGCTTCATTCGCGTGAACCTTCACATTCTTGG‑‑‑CGCTTTGCGGCCCCAATTACAAACACCAAAGCCAATACAATCATCAGGTCATCTGAATGATGGTCTGCATGATGGTGCTTTCAGTGGAAATGGTCTTGGCATTGGCCTGGTAGTTGCTCTGCGCCTTGATCAGG  >  NC_002947/4980455‑4980719

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.