Predicted mutation
evidence seq id position mutation freq annotation gene description
MC JC CP009273 4,303,167 Δ450 bp 100% coding (1‑450/450 nt) rpiB → ribose 5‑phosphate isomerase B/allose 6‑phosphate isomerase

Missing coverage evidence...
   seq id start end size ←reads reads→ gene description
* * ÷ CP009273 4303167 4303616 450 20 [0] [0] 21 rpiB ribose 5‑phosphate isomerase B/allose 6‑phosphate isomerase

New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? CP009273 = 43031660 (0.000)17 (1.650) 17/270 NT 100% intergenic (‑358/‑1) alsR/rpiB d‑allose‑inducible als operon transcriptional repressor; autorepressor; repressor of rpiR/ribose 5‑phosphate isomerase B/allose 6‑phosphate isomerase
?CP009273 4303617 = 0 (0.000)intergenic (+1/‑68) rpiB/yjdP ribose 5‑phosphate isomerase B/allose 6‑phosphate isomerase/putative periplasmic protein

TGAGTGTCTTTAACGCTATTGCGCTGTCAGATGTAACTCTGTAAAACAGATCAGGAAGGCGTAATTCATCGATTTTGGGTGTTTTTAATTCATTAACATCACAAATGTTTTTTGATTGTGAAGTTTTGCACGGACGGGGAAGATGAATGAAAAAGATTGCATTTGGCTGTGATCATGTCGGTTTCATTTTAAAACATGAAATAGTGGCACATTTAGTTGAGCGTGGCGTTGAAGTGATTGATAAAGGAACCTGGTCGTCAGAGCGTACTGATTATCCACATTACGCCAGT  >  CP009273/4303021‑4303310
                                                                                                                                                 |                                                                                                                                                
tGAGTGTCTTTAACGCTATTGCGCTGTCAGATGTAACTCTGTAAAACAGATCAGGAAGGCGTAATTCATCGATTTTGGGTGTTTTTAATTCATTAACATCACAAATGTTTTTTGATTGTGAAGTTTTGCACGGACGGGGAAGATGAgatt                                                                                                                                              <  2:24880‑M1/150‑5 (MQ=255)
 gAGTGTCTTTAACGCTATTGCGCTGTCAGATGTAACTCTGTAAAACAGATCAGGAAGGCGTAATTCATCGATTTTGGGTGTTTTTAATTCATTAACATCACAAATGTTTTTTGATTGTGAAGTTTTGCACGGACGGGGAAGATGAgatt                                                                                                                                              >  1:162005‑M1/1‑145 (MQ=255)
             cGCTATTGCGCTGTCAGATGTAACTCTGTAAAACAGATCAGGAAGGCGTAATTCATCGATTTTGGGTGTTTTTAATTCATTAACATCACAAATGTTTTTTGATTGTGAAGTTTTGCACGGACGGGGAAGATGAgattcatccactacttg                                                                                                                                 >  2:164246‑M1/1‑133 (MQ=255)
                       cTGTCAGATGTAACTCTGTAAAACAGATCAGGAAGGCGTAATTCATCGATTTTGGGTGTTTTTAATTCATTAACATCACAAATGTTTTTTGATTGTGAAGTTTTGCACGGACGGGGAAGATGAgattcatccactacttgcatggatgag                                                                                                                       <  1:16560‑M1/150‑28 (MQ=255)
                                 tAACTCTGTAAAACAGATCAGGAAGGCGTAATTCATCGATTTTGGGTGTTTTTAATTCATTAACATCACAAATGTTTTTTGATTGTGAAGTTTTGCACGGACGGGGAAGATGAgattcatccactacttgcatggatgagtaatgattaa                                                                                                             >  2:61864‑M1/1‑113 (MQ=255)
                                       tGTAAAACAGATCAGGAAGGCGTAATTCATCGATTTTGGGTGTTTTTAATTCATTAACATCACAAATGTTTTTTGATTGTGAAGTTTTGCACGGACGGGGAAGATGAgattcatccactacttgcatggatgagtaatgattaatgtg                                                                                                         <  2:18591‑M1/148‑42 (MQ=255)
                                                    aGGAAGGCGTAATTCATCGATTTTGGGTGTTTTTAATTCATTAACATCACAAATGTTTTTTGATTGTGAAGTTTTGCACGGACGGGGAAGATGAgattcatccactacttgcatggatgagtaatgattaatgtggatagagtttctttt                                                                                          <  2:48395‑M1/150‑57 (MQ=255)
                                                            gTAATTCATCGATTTTGGGTGTTTTTAATTCATTAACATCACAAATGTTTTTTGATTGTGAAGTTTTGCACGGACGGGGAAGATGAgattcatccactacttgcatggatgagtaatgattaatgtggatagagtttctttttgaggttg                                                                                  >  1:61010‑M1/1‑86 (MQ=255)
                                                                      gATTTTGGGTGTTTTTAATTCATTAACATCACAAATGTTTTTTGATTGTGAAGTTTTGCACGGACGGGGAAGATGAgattcatccactacttgcatggatgagtaatgattaatgtggatagagtttctttttgaggttggctaatga                                                                          <  1:104325‑M1/148‑73 (MQ=255)
                                                                          ttGGGTGTTTTTAATTCATTAACATCACAAATGTTTTTTGATTGTGGAGTTTTGCACGGACGGGGAAGATGAgattcatccactacttgcatggatgagtaatgattaatgtggatagagtttctttttgaggttggctaatgaaacgc                                                                     <  1:164246‑M1/149‑78 (MQ=255)
                                                                                 tttttAATTCATTAACATCACAAATGTTTTTTGATTGTGAAGTTTTGCACGGACGGGGAAGATGAgattcatccactacttgcatggatgagtaatgattaatgtggatagagtttctttttgaggttggctaatgaaacgctttccgc                                                              >  1:109843‑M1/1‑65 (MQ=255)
                                                                                   tttAATTCATTAACATCACAAATGTTTTTTGATTGTGAAGTTTTGCACGGACGGGGAAGATGAgattcatccactacttgcatggatgagtaatgattaatgtggatagagtttctttttgaggttggctaatgaaacgctttccgctat                                                           >  2:83051‑M1/1‑63 (MQ=255)
                                                                                       aTTCATTAACATCACAAATGTTTTTTGATTGTGAAGTTTTGCACGGACGGGGAAGATGAgattcatccactacttgcatggatgagtaatgattaatgtggatagagtttctttttgaggttggctaatgaaacgctttccgctatttct                                                       <  1:61864‑M1/150‑92 (MQ=255)
                                                                                         tCATTAACATCACAAATGTTTTTTGATTGTGAAGTTTTGCACGGACGGGGAAGATGAgattcatccactacttgcatggatgagtaatgattaatgtggatagagtttctttttgaggttggctaatgaaacgctttccgctatttcttt                                                     >  1:124141‑M1/1‑57 (MQ=255)
                                                                                                cATCACAAATGTTTTTTGATTGTGAAGTTTTGCACGGACGGGGAAGATGAgattcatccactacttgcatggatgagtaatgattaatgtggatagagtttctttttgaggttggctaatgaaacgctttccgctatttcttttattcac                                              <  2:124141‑M1/150‑101 (MQ=255)
                                                                                                        aTGTTTTTTGATTGTGAAGTTTTGCACGGACGGGGAAGATGAgattcatccactacttgcatggatgagtaatgattaatgtggatagagtttctttttgaggttggctaatgaaacgctttccgctatttcttttattcaccctgctca                                      >  1:89937‑M1/1‑42 (MQ=255)
                                                                                                                                  cggacggGGAAGATGAgattcatccactacttgcatggatgagtaatgattaatgtggatagagtttctttttgaggttggctaattaaacgctttccgctatttcttttattcaccctgctcacgctgtccaccgttccggcccaggcc            >  2:164538‑M1/1‑16 (MQ=255)
                                                                                                                                   ggacggGGAAGATGAgattcatccactacttgcatggatgagtaatgattaatgtggatagagtttctttttgaggttggctaatgaaacgctttccgctatttcttttattcaccctgctcacgctgtccaccgttccggcccaggccg           <  1:165465‑M1/150‑136 (MQ=255)
                                                                                                                                       ggggAAGATGAgattcatccactacttgcatggatgagtaatgattaatgtggatagagtttctttttgaggttggctaatgaaacgctttccgctatttcttttattcaccctgctcacgctgtccaccgttccggcccaggccgatat       >  1:174257‑M1/1‑11 (MQ=255)
                                                                                                                                            aGATGAgattcatccactacttgcatggatgagtaatgattaatgtggatagagtttctttttgaggttggctaatgaaacgctttccgctatttcttttattcaccctgctcacgctgtccaccgttccggcccaggccgatattatcg  <  1:83051‑M1/150‑145 (MQ=255)
                                                                                                                                                 |                                                                                                                                                
TGAGTGTCTTTAACGCTATTGCGCTGTCAGATGTAACTCTGTAAAACAGATCAGGAAGGCGTAATTCATCGATTTTGGGTGTTTTTAATTCATTAACATCACAAATGTTTTTTGATTGTGAAGTTTTGCACGGACGGGGAAGATGAATGAAAAAGATTGCATTTGGCTGTGATCATGTCGGTTTCATTTTAAAACATGAAATAGTGGCACATTTAGTTGAGCGTGGCGTTGAAGTGATTGATAAAGGAACCTGGTCGTCAGAGCGTACTGATTATCCACATTACGCCAGT  >  CP009273/4303021‑4303310

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 17 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: