Predicted mutation
evidence seq id position mutation annotation gene description
MC JC NC_000913 4,518,060 Δ7 bp coding (170‑176/522 nt) fecI ← RNA polymerase sigma factor FecI

Missing coverage evidence...
   seq id start end size ←reads reads→ gene description
* * ÷ NC_000913 4518060 4518066 7 38 [0] [0] 38 fecI RNA polymerase sigma factor FecI

New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913 = 45180590 (0.000)35 (0.680) 25/264 0.8 100% coding (177/522 nt) fecI RNA polymerase sigma factor FecI
?NC_000913 4518067 = 0 (0.000)coding (169/522 nt) fecI RNA polymerase sigma factor FecI

AAGGCGCTCCCCCCTCCGGCATAAGCGCCAGCATCTCCAGATACGCTTTTTCCAGCGCGTTTCGGCGAAACAGGTCCACCATCACGCGTTTGGCGATAGTGCAGAGGAAGGAGCGAGGATCGCGGATCGTCGAGAGCGTTTCGCTGACCATTACCCGCAAAAAAGTGTCCTGGGCAATGTCATCTGCATCAAAAGCAGACTGGAGTTTGCGCGTCAGCCAGCTTTTCAACCAGCCGTGATGTGTGCCATAAAGCGACTCGA  >  NC_000913/4517938‑4518198
                                                                                                                         |                                                                                                                                           
aaGGCGCTCCCCCCTCCGGCATAAGCGCCAGCATCTCCAGATACGCTTTTTCCAGCGCGTTTCGGCGAAACAGGTCCACCATCACGCGTTTGGCGATAGTGCAGAGGAAGGAGCGAGGATGGtcgagagcgtttcgctga                                                                                                                           >  1:195787‑M1/1‑122 (MQ=255)
aaGACGCTCCCCCCTCCGGCATAAGCGCCAGCATCTCCATATACGCTTTTTCGAGCGCGTTTCGGCGAAACAGGTCCACCATCACGCGTTTGGCGATAGTGCAGAGGAAGGAGCGAGGATCGtcgagagcgtttcgctga                                                                                                                           >  2:401656‑M1/1‑122 (MQ=255)
  ggCGCTCCCCCCTCCGGCATAAGCGCCAGCATCTCCAGATACGCTTTTTCCAGCGCGTTTCGGCGAAACAGGTCCACCATCACGCGTTTGGCGATAGTGCAGAGGAAGGAGCGAGGATCGtcgagagcgtttcgctgacc                                                                                                                         >  2:835698‑M1/1‑120 (MQ=255)
  ggCGCTCCCCCCTCCGGCATAAGCGCCAGCATCTCCAGATACGCTTTTTCCAGCGCGTTTCGGCGAAACAGGTCCACCATCACGCGTTTGGCGATAGTGCAGAGGAAGGAGCGAGGATCGtcgagagcg                                                                                                                                    <  1:483993‑M1/129‑10 (MQ=255)
  ggCGCTCCCCCCTCCGGCATAAGCGCCAGCATCTCCAGATACGCTTTTTCCAGCGCGTTTCGGCGAAACAGGTCCACCATCACGCGTTTGGCGATAGTGCAGAGGAAGGAGCGAGGATCGtcgagagcg                                                                                                                                    >  2:483993‑M1/1‑120 (MQ=255)
      cTCCCCCCTCCGGCATAAGCGCCAGCATCTCCAGATACGCTTTTTCCAGCGCGTTTCGGCGAAACAGGTCCACCATCAGGCGTTTGGCGATAGTGCAGAGGAAGGAGCGAGGATCGtcgagagcgtttcgctgaccatta                                                                                                                     >  1:324227‑M1/1‑116 (MQ=255)
       tCCCCCCTCCGGCATAAGCGCCAGCATCTCCAGATACGCTTTTTCCAGCGCGTTTCGGCGAAACAGGTCCACCATCACGCGTTTGGCGATAGTGCAGAGGAAGGAGCGAGGATCGtcgagagcgtttcgctgaccattac                                                                                                                    >  2:508178‑M1/1‑115 (MQ=255)
       tCCCCCCTCCGGCATAAGCGCCAGCATCTCCAGATACGCTTTTTCCAGCGCGTTTCGGCGAAACAGGTCCACCATCACGCGTTTGGCGATAGTGCAGAGGAAGGAGCGAGGATCGtcgagagcgtttcgctgaccattac                                                                                                                    >  1:19160‑M1/1‑115 (MQ=255)
       tCCCCCCTCCGGCATAAGCGCCAGCATCTCCAGATACGCTTTTTCCAGCGCGTTTCGGCGAAACAGGTCCACCATCACGCGTTTGGCGATAGTGCAGAGGAAGGAGCGAGGATCGtcgagagcgtttcgctgaccattac                                                                                                                    >  2:726340‑M1/1‑115 (MQ=255)
       tCCCCCCTCCGGCATAAGCGCCAGCATCTCCAGATACGCTTTTTCCAGCGCGTTTCGGCGAAACAGGTCCACCATCACGCGTTTGGCGATAGTGCAGAGGAAGGAGCGAGGATCGtcgagagcgtttcgctgaccattac                                                                                                                    >  1:886001‑M1/1‑115 (MQ=255)
              tCCGGCATAAGCGCCAGCATCTCCAGATACGCTTTTTCCAGCGCGTTTCGGCGAAACAGGTCCACCATCACGCGTTTGGCGATAGTGCAGAGGAAGGAGCGAGGATCGtcgagagcgtttcgctgaccattacccgcaaa                                                                                                             >  1:879500‑M1/1‑108 (MQ=255)
                cGGCATAAGCGCCAGCATCTCCAGATACGCTTTTTCCAGCGCGTTTCGGCGAAACAGGTCCACCATCACGCGTTTGGCGATAGTGCAGAGGAAGGAGCGAGGATCGtcgagagcgtttcgctgaccattacccgcaaaaa                                                                                                           >  1:452747‑M1/1‑106 (MQ=255)
                 ggCATAAGCGCCAGCATCTCCAGATACGCTTTTTCCAGCGCGTTTCGGCGAAACAGGTCCACCATCACGCGTTTGGCGATAGTGCAGAGGAAGGAGCGAGGATCGtcgagagcgtttcgctgaccattacccgcaaaaaa                                                                                                          <  2:509755‑M1/140‑36 (MQ=255)
                 ggCATAAGCGCCAGCATCTCCAGATACGCTTTTTCCAGCGCGTTTCGGCGAAACAGGTCCACCATCACGCGTTTGGCGATAGTGCAGAGGAAGGAGCGAGGATCGtcgagagcgtttcgctgaccattacccgcaaaaaa                                                                                                          <  1:208933‑M1/140‑36 (MQ=255)
                      aaGCGCCAGCATCTCCAGATACGCTTTTTCCAGCGCGTTTCGGCGAAACAGGTCCACCATCAGGCGTTTGGCGATAGTGCAGAGGAAGGAGCGAGGATCGtcgagagcgtttcgctgaccattacccgcaaaaaagtgtc                                                                                                     <  2:324227‑M1/140‑41 (MQ=255)
                      aaGCGCCAGCATCTCCAGATACGCTTTTTCCAGCGCGTTTCGGCGAAACAGGTCCACCATCACGCGTTTGGCGATAGTGCAGAGGAAGGAGCGAGGATCGtcgagagcgtttcgctgaccattacccgcaaaaaagtgtc                                                                                                     <  2:783948‑M1/140‑41 (MQ=255)
                      aaGCGCCAGCATCTCCAGATACGCTTTTTCCAGCGCGTTTCGGCGAAACAGGTCCACCATCACGCGTTTGGCGATAGTGCAGAGGAAGGAGCGAGGATCGtcgagagcgtttcgctgaccattacccgcaaaaaagtgtc                                                                                                     <  2:367966‑M1/140‑41 (MQ=255)
                          gcCAGCATCTCCAGATACGCTTTTTCCAGCGCGTTTCGGCGAAACAGGTCCACCATCACGCGTTTGGCGATAGTGCAGAGGAAGGAGCGAGGATCGtcgagagcgtttcgctgaccattacccgcaaaaaagtgtcctgg                                                                                                 <  1:726340‑M1/140‑45 (MQ=255)
                          gcCAGCATCTCCAGATACGCTTTTTCCAGCGCGTTTCGGCGAAACAGGTCCACCATCACGCGTTTGGCGATAGTGCAGAGGAAGGAGCGAGGATCGtcgagagcgtttcgctgaccattacccgcaaaaaagtgtcctgg                                                                                                 <  1:906469‑M1/140‑45 (MQ=255)
                                  ctcCAGATACGCTTTTTCCAGCGCGTTTCGGCGAAACAGGTCCACCATCACGCGTTTGGCGATAGTGCAGAGGAAGGAGCGAGGATCGtcgagagcgtttcgctgaccattacccgcaaaaaagtgtcctgggcaatgtc                                                                                         >  1:699150‑M1/1‑88 (MQ=255)
                                       gATACGCTTTTTCCAGCGCGTTTCGGCGAAACAGGTCCACCATCACGCGTTTGGCGATAGTGCAGAGGAAGGAGCGAGGATCGtcgagagcgtttcgctgaccattacccgcaaaaaagtgtcctgggcaatgtcatctg                                                                                    <  2:257452‑M1/140‑58 (MQ=255)
                                          aCGCTTTTTCCAGCGCGTTTCGGCGAAACAGGTCCACCATCACGCGTTTGGCGATAGTGCAGAGGAAGGAGCGAGGATCGtcgagagcgtttcgctgaccattacccgcaaaaaagtgtcctgggcaatgtcatctgcat                                                                                 >  2:292066‑M1/1‑80 (MQ=255)
                                                tttCCAGCGCGTTTCGGCGAAACAGGTCCACCATCACGCGTTTGGCGATAGTGCAGAGGAAGGAGCGAGGATCGtcgagagcgtttcgctgaccattacccgcaaaaaagtgtcctgggcaatgtcatctgcatcaaaag                                                                           <  2:169349‑M1/140‑67 (MQ=255)
                                                tttCCAGCGCGTTTCGGCGAAACAGGTCCACCATCACGCGTTTGGCGATAGTGCAGAGGAAGGAGCGAGGATCGtcgagagcgtttcgctgaccattacccgcaaaaaagtgtcctgggcaatgtcatctgcatcaaaag                                                                           <  1:232799‑M1/140‑67 (MQ=255)
                                                                         gTCCACCATCACGCGTTTGGCGATAGTGCAGAGGAAGGAGCGAGGATCGtcgagagcgtttcgctgaccattacccgcaaaaaagtgtcctgggcaatgtcatctgcatcaaaagcagactggagtttgcgcgtcagcca                                                  <  1:292066‑M1/140‑92 (MQ=255)
                                                                                   aCGCGTTTGGCGATAGTGCAGAGGAAGGAGCGAGGATCGtcgagagcgtttcgctgaccattacccgcaaaaaagtgtcctgggcaatgtcatctgcatcaaaagcagactggagtttgcgcgtcagccagcttttcaac                                        <  2:19160‑M1/140‑102 (MQ=255)
                                                                                   aCGCGTTTGGCGATAGTGCAGAGGAAGGAGCGAGGATCGtcgagagcgtttcgctgaccattacccgcaaaaaagtgtcctgggcaatgtcatctgcatcaaaagcagactggagtttgcgcgtcagccagcttttcaac                                        >  1:462589‑M1/1‑39 (MQ=255)
                                                                                    cgcgTTTGGCGATAGTGCAGAGGAAGGAGCGAGGATCGtcgagagcgtttcgctgaccattacccgcaaaaaagtgtcctgggcaatgtcatctgcatcaaaagcagactggagtttgcgcgtcagccagcttttcaacc                                       >  1:329312‑M1/1‑38 (MQ=255)
                                                                                         ttGGCGATAGTGCAGAGGAAGGAGCGAGGATCGtcgagagcgtttcgctgaccattacccgcaaaaaagtgtcctgggcaatgtcatctgcatcaaaagcagactggagtttgcgcgtcagccagcttttcaaccagccg                                  >  2:877526‑M1/1‑33 (MQ=255)
                                                                                             cGATAGTGCAGAGGAAGGAGCGAGGATCGtcgagagcgtttcgctgaccattacccgcaaaaaagtgtcctgggcaatgtcatctgcatcaaaagcagactggagtttgcgcgtcagccagcttttcaaccagccgtgat                              <  2:462589‑M1/140‑112 (MQ=255)
                                                                                                tAGTGCAGAGGAAGGAGCGAGGATCGtcgagagcgtttcgctgaccattacccgcaaaaaagtgtcctgggcaatgtcatctgcatcaaaagcagactggagtttgcgcgtcagccagcttttcaaccagccgtgatgtg                           <  2:699150‑M1/140‑115 (MQ=255)
                                                                                                   tGCAGAGGAAGGAGCGAGGATCGtcgagagcgtttcgctgaccattacccgcaaaaaagtgtcctgggcaatgtcatctgcatcaaaagcagactggagtttgcgcgtcagccagcttttcaaccagccgtgatgtgtgc                        >  1:275457‑M1/1‑23 (MQ=255)
                                                                                                             ggaGCGAGGATCGtcgagagcgtttcgctgaccattacccgcaaaaaagtgtcctgggcaatgtcatctgcatcaaaagcagactggagtttgcgcgtcagccagcttttcaaccagccgtgatgtgtgccataaagcga              >  1:616943‑M1/1‑13 (MQ=255)
                                                                                                               aGCGAGGATCGtcgagagcgtttcgctgaccattacccgcaaaaaagtgtcctgggcaatgtcatctgcatcaaaagcagactggagtttgcgcgtcagccagcttttcaaccagccgtgatgtgtgccataaagcgact            >  2:103173‑M1/1‑11 (MQ=255)
                                                                                                                  gAGGATCGtcgagagcgtttcgctgaccattacccgcaaaaaagtgtcctgggcaatgtcatctgcatcaaaagcagactggagtttgcgcgtcagccagcttttcaaccagccgtgatgtgtgccataaagcgactcga         >  2:50824‑M1/1‑8 (MQ=255)
                                                                                                                      aTCGtcgagagcgtttcgctgaccattacccgcaaaaaagtgtcctgggcaatgtcatctgcatcaaaagcagactggagtttgcgcgtcagccagcttttcaaccagccgtgatgtgtgccataaagcgactcgaacgt     >  1:476767‑M1/1‑4 (MQ=255)
                                                                                                                      aTCGtcgagagcgtttcgctgaccattacccgcaaaaaagtgtcctgggcaatgtcatctgcatcaaaagcagactggagtttgcgcgtcagccagcttttcaaccagccgtgatgtgtgccataaagcgactcgaacgt     >  2:760628‑M1/1‑4 (MQ=255)
                                                                                                                         gtcgagagcgtttcgctgaccattacccgcaaaaaagtgtcctgggcaatgtcatctgcatcaaaagcagactggagtttgcgcgtcagccagcttttcaaccagccgtgatgtgtgccataaagcgactcgaacgttaa  <  2:616943‑M1/140‑140 (MQ=255)
                                                                                                                         |                                                                                                                                           
AAGGCGCTCCCCCCTCCGGCATAAGCGCCAGCATCTCCAGATACGCTTTTTCCAGCGCGTTTCGGCGAAACAGGTCCACCATCACGCGTTTGGCGATAGTGCAGAGGAAGGAGCGAGGATCGCGGATCGTCGAGAGCGTTTCGCTGACCATTACCCGCAAAAAAGTGTCCTGGGCAATGTCATCTGCATCAAAAGCAGACTGGAGTTTGCGCGTCAGCCAGCTTTTCAACCAGCCGTGATGTGTGCCATAAAGCGACTCGA  >  NC_000913/4517938‑4518198

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 40 ≤ 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.