Predicted mutation
evidence seq id position mutation annotation gene description
MC JC NC_002947 1,648,886 Δ12 bp intergenic (‑123/+162) oprB‑II ← / ← PP_1446 carbohydrate‑selective porin/TonB‑dependent receptor

Missing coverage evidence...
   seq id start end size ←reads reads→ gene description
* * ÷ NC_002947 1648886 1648897 12 31 [0] [0] 31 oprB‑II/PP_1446 carbohydrate‑selective porin/TonB‑dependent receptor

New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_002947 = 16488850 (0.000)31 (0.770) 29/262 0.5 100% intergenic (‑122/+174) oprB‑II/PP_1446 carbohydrate‑selective porin/TonB‑dependent receptor
?NC_002947 1648898 = 0 (0.000)intergenic (‑135/+161) oprB‑II/PP_1446 carbohydrate‑selective porin/TonB‑dependent receptor

TTTAGGCAGTTGGAACATGAGATAGCGCTATCTTTTGATTATTGAAAAAAACGGCAAGCGCATCGAAGGGCGGGCCCCTCAATGGATATGGCGACAAGATTGTTCGATGAAGCGTTTCAGCTTGCGGCGGGGAGGATATAGGGGAAGTTATAAGAAAAAAGACAAATTGTCGCAGAGGATTGTTGTTGAACGGGCAACAGTCGTCATTAATCGAGGCAGGTGTAGGAGCGGGTTCACCCGCGAATGCGCCGGTTGCCCTACCATCGTATTCGCGGGTGAACCCGT  >  NC_002947/1648746‑1649030
                                                                                                                                           |                                                                                                                                                 
tttAGGCAGTTGGAACATGAGATAGCGCTATCTTTTGATTATTGAAAAAAACGGCAAGCGCATCGAAGGGCGGGCCCCTCAATGGATATGGCGACAAGATTGTTCGATGAAGCGTTTCAGCTTGCGGCGGGGAGGATATAagaaaaaaga                                                                                                                                         >  2:657187‑M1/1‑140 (MQ=255)
       aGTTGGAACATGAGATAGCGCTATCTTTTGATTATTGAAAAAAACGGCAAGCGCATCGAAGGGCGGGCCCCTCAATGGATATGGCGACAAGATTGTTCGATGAAGCGTTTCAGCTTGCGGCGGGGAGGATATAagaaaaaagacaaattg                                                                                                                                  >  1:472011‑M1/1‑133 (MQ=255)
            gAACATGAGATAGCGCTATCTTTTGATTATTGAAAAAAACGGCAAGCGCATCGAAGGGCGGGCCCCTCAATGGATATGGCGACAAGATTGTTCGATGAAGCGTTTCAGATTGCGGTGGGGAGGAGATAagaaaaaagacaaattgtcgca                                                                                                                             >  2:779231‑M1/1‑128 (MQ=255)
                    gaTAGCGCTATCTTTTGATTATTGAAAAAAACGGCAAGCGCATCGAAGGGCGGGCCCCTCAATGGATATGGCGACAAGATTGTTCGATGAAGCGTTTCAGCTTGCGGCGGGGAGGATATAagaaaaaagacaaattgtcgcagaggattg                                                                                                                     <  1:657187‑M1/150‑31 (MQ=255)
                         cgcTATCTTTTGATTATTGAAAAAAACGGCAAGCGCATCGAAGGGCGGGCCCCTCAATGGATATGGCGACAAGATTGTTCGATGAAGCGTTTCAGCTTGCGGCGGGGAGGATATAagaaaaaagacaaattgtcgcagaggattgttgtt                                                                                                                <  2:472011‑M1/150‑36 (MQ=255)
                                            aaaaaaaCGGCAAGCGCATCGAAGGGCGGGCCCCTCAATGGATATGGCGACAAGATTGTTCGATGAAGCGTTTCAGCTTGCGGCGGGGAGGATATAagaaaaaagacaaatt                                                                                                                                   <  1:566962‑M1/112‑17 (MQ=255)
                                            aaaaaaaCGGCAAGCGCATCGAAGGGCGGGCCCCTCAATGGATATGGCGACAAGATTGTTCGATGAAGCGTTTCAGCTTGCGGCGGGGAGGATATAagaaaaaagacaaatt                                                                                                                                   >  2:566962‑M1/1‑96 (MQ=255)
                                                 aaCGGCAAGCGCATCGAAGGACGGGCCCCTCAATGGATATGGCGACAAGATTGTTCGATGAAGCGTTTCAGCTTGCGGCGGGGAGGATATAagaaaaaagacaaattgtcgcagaggattgttgttgaacgggcaacagtcgtcattaat                                                                                        >  2:487551‑M1/1‑91 (MQ=255)
                                                    ggCAAGCGCATCGAAGGGCGGGCCCCTCAATGGATATGGCGACAAGATTGTTCGATGAAGCGTTTCAGCTTGCGGCGGGGAGGATATAagaaaaaagacaaattgtcgcagaggattgttgttgaacgggcaacagtcgtcattaatcga                                                                                     <  1:529760‑M1/150‑63 (MQ=255)
                                                    ggCAAGCGCATCGAAGGGCGGGCCCCTCAATGGATATGGCGACAAGATTGTTCGATGAAGCGTTTCAGCTTGCGGCGGGGAGGATATAagaaaaaagacaaattgtcgcagaggattgttgttgaacgggcaacagtcgtcattaatcga                                                                                     <  2:451015‑M1/150‑63 (MQ=255)
                                                           gcATCGAAGGGCGGGCCCCTCAATGGATATGGCGACAAGATTGTTCGATGAAGCGTTTCAGCTTGCGGCGGGGAGGATATAagaaaaaagacaaattgtcgcaga                                                                                                                           >  1:861923‑M1/1‑81 (MQ=255)
                                                           gcATCGAAGGGCGGGCCCCTCAATGGATATGGCGACAAGATTGTTCGATGAAGCGTTTCAGCTTGCGGCGGGGAGGATATAagaaaaaagacaaattgtcgcaga                                                                                                                           <  2:861923‑M1/105‑25 (MQ=255)
                                                                gAAGGGCGGGCCCCTCAATGGATATGGCGACAAGATTGTTCGATGAAGCGTTTCAGCTTGCGGCGGGGAGGATATAagaaaaaagacaaattgtcgcagaggattgttgttgaacgggcaaca                                                                                                    >  2:515380‑M1/1‑76 (MQ=255)
                                                                gAAGGGCGGGCCCCTCAATGGATATGGCGACAAGATTGTTCGATGAAGCGTTTCAGCTTGCGGCGGGGAGGATATAagaaaaaagacaaattgtcgcagaggattgttgttgaacgggcaaca                                                                                                    <  1:515380‑M1/123‑48 (MQ=255)
                                                                gAAGGGCGGGCCCCTCAATGGATATGGCGACAAGATTGTTCGATGAAGCGTTTCAGCTTGCGGCGGGGAGGATATAagaaaaaagacaaat                                                                                                                                    >  1:213032‑M1/1‑76 (MQ=255)
                                                                gAAGGGCGGGCCCCTCAATGGATATGGCGACAAGATTGTTCGATGAAGCGTTTCAGCTTGCGGCGGGGAGGATATAagaaaaaagacaaat                                                                                                                                    <  2:213032‑M1/91‑16 (MQ=255)
                                                                           cccTCAATGGATATGGCGACAAGATTGTTCGATGAAGCGTTTCAGCTTGCGGCGGGGAGGATATAagaaaaaagacaaattgtcgcagaggattgttgtttaacgggcaacagtcgtcattaatcgaggcaggtgtaggagcgggttcac                                                              >  1:835287‑M1/1‑65 (MQ=255)
                                                                             cTCAATGGATATGGCGACAAGATTGTTCGATGAAGCGTTTCAGCTTGCGGCGGGGAGGATATAagaaaaaagacaaattgtcgcagaggattgttgttgaacgggcaacagtcgtcattaatcgaggcaggtgtaggagcgggttcaccc                                                            >  1:479204‑M1/1‑63 (MQ=255)
                                                                              tCAATGGATATGGCGACAAGATTGTTCGATGAAGCGTTTCAGCTTGCGGCGGGGAGGATATAagaaaaaagacaaattgtcgcagaggattgttgttgaacgggcaacagtcgtcattaatcgaggcaggtgtaggagcgggttcacccg                                                           >  2:416124‑M1/1‑62 (MQ=255)
                                                                                            gACAAGATTGTTCGATGAAGCGTTTCAGCTTGCGGCGGGGAGGATATAagaaaaaagacaaattgtcgcagaggattgttgttgaacgggcaacagtcgtcattaatcgaggcaggtgtaggagcgggttcacccgcgaatgcgccggtt                                             >  2:1017428‑M1/1‑48 (MQ=255)
                                                                                             aCAAGATTGTTCGATGAAGCGTTTCAGCTTGCGGCGGGGAGGATATAagaaaaaagacaaattgtcgcagaggattgttgttgaacgggcaacagtcgtcattaatcgaggcaggtgtaggagcgggttcacccgcgaatgcgccggttg                                            >  1:544872‑M1/1‑47 (MQ=255)
                                                                                               aaGATTGTTCGATGAAGCGTTTCAGCTTGCGGCGGGGAGGATATAagaaaaaagacaaattgtcgcagaggattgttgttgaacgggcaacagtcgtcattaatcgaggcaggtgtaggagcgggttcacccgcgaatgcgccggttgcc                                          <  1:779231‑M1/150‑106 (MQ=255)
                                                                                                        cGATGAAGCGTTTCAGCTTGCGGCGGGGAGGATATAagaaaa                                                                                                                                             >  1:992553‑M1/1‑36 (MQ=255)
                                                                                                        cGATGAAGCGTTTCAGCTTGCGGCGGGGAGGATATAagaaaa                                                                                                                                             <  2:992553‑M1/42‑7 (MQ=255)
                                                                                                             aaGCGTTTCAGCTTGCGGCGGGGAGGATATAagaaaaaagacaaattgtcgcagaggattgttgttgaacgggcaacagtcgtcattaatcgaggcaggtgtaggagcgggttcacccgcgaatgcgccggttgccctaccatcgtattc                            >  1:239488‑M1/1‑31 (MQ=255)
                                                                                                                 gTTTCAGCTTGCGGCGGGGAGGATATAagaaaaaagacaaattgtcgcagaggattgttgttgaacgggcaacagtcgtcattaatcgaggcaggtgtaggagcgggttcacccgcgaatgcgccggttgccctaccatcgtattcgcgg                        <  2:544872‑M1/150‑124 (MQ=255)
                                                                                                                         ttGCGGCGGGGAGGATATAagaaaaaagacaaattgtcgcagaggattgttgttgaacgggcaacagtcgtcattaatcgaggcaggtgtaggagcgggttcacccgc                                                          <  2:318375‑M1/108‑90 (MQ=255)
                                                                                                                         ttGCGGCGGGGAGGATATAagaaaaaagacaaattgtcgcagaggattgttgttgaacgggcaacagtcgtcattaatcgaggcaggtgtaggagcgggttcacccgc                                                          >  1:318375‑M1/1‑19 (MQ=255)
                                                                                                                                   gaggaTATAagaaaaaagacaaattgtcgcagaggattgttgt                                                                                                                 >  1:85336‑M1/1‑9 (MQ=255)
                                                                                                                                   gaggaTATAagaaaaaagacaaattgtcgcagaggattgttgt                                                                                                                 <  2:85336‑M1/43‑35 (MQ=255)
                                                                                                                                       atatAagaaaaaagacaaattgtcgcagaggattgttgttgaacgggcaacagtcgtcattaatcgaggcaggtgtaggagcgggttcacccgcgaatgcgccggttgccctaccatcgtattcgcgggtgaacccgttcctacaccgat  <  1:416124‑M1/150‑146 (MQ=255)
                                                                                                                                           |                                                                                                                                                 
TTTAGGCAGTTGGAACATGAGATAGCGCTATCTTTTGATTATTGAAAAAAACGGCAAGCGCATCGAAGGGCGGGCCCCTCAATGGATATGGCGACAAGATTGTTCGATGAAGCGTTTCAGCTTGCGGCGGGGAGGATATAGGGGAAGTTATAAGAAAAAAGACAAATTGTCGCAGAGGATTGTTGTTGAACGGGCAACAGTCGTCATTAATCGAGGCAGGTGTAGGAGCGGGTTCACCCGCGAATGCGCCGGTTGCCCTACCATCGTATTCGCGGGTGAACCCGT  >  NC_002947/1648746‑1649030

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