Predicted mutation
evidence seq id position mutation annotation gene description
MC JC NC_000913 4,091,572 3461 bp→GCATGCGGATCCTCTAGAGTCGACATATG [rhaD]rhaB [rhaD], rhaA, rhaB

Missing coverage evidence...
   seq id start end size ←reads reads→ gene description
* * ÷ NC_000913 4091572 4095032 3461 149 [0] [0] 173 [rhaD]–rhaB [rhaD],rhaA,rhaB

New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913 = 40915710 (0.000)103 (1.070)
+GCATGCGGATCCTCTAGAGTCGACATATG
58/92 0.1 100% coding (283/825 nt) rhaD rhamnulose‑phosphate aldolase
?NC_000913 4095033 = 0 (0.000)intergenic (‑4/‑284) rhaB/rhaS rhamnulokinase/RhaS

Only 100 of 173 total aligned reads displayed.

CTCGTAACGCGCCAGCATCACGCGCCCACTGGATGCGCCGAGATCGACGGCGACACAATTGCGAAAGGTCATAATGTGATCCTGCTGAATTTCATTACGACCAGTCTAAAAAGCGCCTGAATTCGCGACCTTCTCGTTACTGACAGGAAAA  >  NC_000913/4094958‑4095108
                                                                           |                                                                           
tggttaacccccaaagaatgtggtagcccgcgccgtcgctgtcgacgcatgcggatcctctagagtcgacatatgg                                                                             >  1:217859‑M2/76‑76 (MQ=255)
tggttaacccccaaagaatgtggtagcccgcgccgtcgctgtcgacgcatgcggatcctctagagtcgacatatgg                                                                             <  1:2420282‑M2/1‑1 (MQ=255)
 ggttaacccccaaagaatgtggtagcccgcgccgtcgctgtcgacgcatgcggatcctctagagtcgacatatgGt                                                                            <  2:4021102‑M2/2‑1 (MQ=255)
  gttaacccccaaagaatgtggtagcccgcgccgtcgctgtcgacgcatgcggatcctctagagtcgacatatgGTg                                                                           >  1:2456423‑M2/74‑76 (MQ=255)
    taacccccaaagaatgtggtagcccgcgccgtcgctgtcgacgcatgcggatcctctagagtcgacatatgGTGAt                                                                         >  1:2175305‑M2/72‑76 (MQ=255)
       cccccaaagaatgtggtagcccgcgccgtcgctgtcgacgcatgcggatcctctagagtcgacatatgGTGATCct                                                                      >  2:2934266‑M2/69‑76 (MQ=255)
       cccccaaagaatgtggtagcccgcgccgtcgctgtcgacgcatgcggatcctctagagtcgacatatgGTGATCct                                                                      <  1:2795550‑M2/8‑1 (MQ=255)
        ccccaaagaatgtggtagcccgcgccgtcgctgtcgacgcatgcggatcctctagagtcgacatatgGTGATCctg                                                                     <  2:1560582‑M2/9‑1 (MQ=255)
          ccaaagaatgtggtagcccgcgccgtcgctgtcgacgcatgcggatcctctagagtcgacatatgGTGATCctgct                                                                   <  2:5461003‑M2/11‑1 (MQ=255)
            aaagaatgtggtagcccgcgccgtcgctgtcgacgcatgcggatcctctagagtcgacatatgGTGATCCTGCTGa                                                                 <  1:2879493‑M2/13‑1 (MQ=255)
               gaatgtggtagcccgcgccgtcgctgtcgacgcatgcggatcctctagagtcgacatatgGTGATCCTGCTGAAtt                                                              <  2:3219961‑M2/16‑1 (MQ=255)
               gaatgtggtagcccgcgccgtcgctgtcgacgcatgcggatcctctagagtcgacatatgGTGATCCTCCTGAAtt                                                              >  2:924878‑M2/61‑76 (MQ=255)
                aatgtggtagcccgcgccgtcgctgtcgacgcatgcggatcctctagagtcgacatatgGTGATCCTGCTGAAttt                                                             >  1:1027434‑M2/60‑76 (MQ=255)
                aatgtggtagcccgcgccgtcgctgtcgacgcatgcggatcctctagagtcgacatatgGTGATCCTGCTGAAttt                                                             <  2:2127728‑M2/17‑1 (MQ=255)
                aatgtggtagcccgcgccgtcgctgtcgacgcatgcggatcctctagagtcgacatatgGTGATCCTGCTGAAttt                                                             <  2:3488026‑M2/17‑1 (MQ=255)
                 atgtggtagcccgcgccgtcgctgtcgacgcatgcggatcctctagagtcgacatatgGTGATCCTGCTGAATTTc                                                            >  2:4444058‑M2/59‑76 (MQ=255)
                  tgtgttcgcccgcgccgtcgctgtcgccgcatgcggatcctctcgagtcgacatatgGTGATCCTGCTGAATTTCa                                                           <  1:1879092‑M2/19‑1 (MQ=255)
                  tgtggtagcccgcgccgtcgctgtcgacgcatgcggatcctctagagtcgacatatgGTGATCCTGCTGAATTTCa                                                           >  1:2232209‑M2/58‑76 (MQ=255)
                   gtggtagcccgcgccgtcgctgtcgacgcatgcggatcctctagagtcgacatatgGTGATCCTGCTGAATTTCAt                                                          <  2:4974029‑M2/20‑1 (MQ=255)
                    tggtagcccgcgccgtcgctgtcgacgcatgcggatcctctagagtcgacatatgGTGATCCTGCTGAATTTCAtt                                                         <  1:4213132‑M2/21‑1 (MQ=255)
                    gggtggccggcgccgtcgctgtcgacgcatgcggatcctctagagtcgacatatgGTGATCCTGCTGAATTTCAtt                                                         <  2:5780500‑M2/21‑1 (MQ=255)
                    gggtagcccgcgccgtcgctgtcgacgcatgcggatcctctagagtcgacatatgGTGATCCTGCTGAATTTCAtt                                                         <  2:2535665‑M2/21‑1 (MQ=255)
                     ggtagcccgcgccgtcgctgtcgacgcatgcggatcctctagagtcgacatatgGTGATCCTGCTGAATTTCATTa                                                        <  2:2226531‑M2/22‑1 (MQ=255)
                        agcccgcgccgtcgctgtcgacgcatgcggatcctctagagtcgacatatgGTGATCCTGCTGAATTTCATTACGa                                                     >  1:3340758‑M2/52‑76 (MQ=255)
                        agcccgcgccgtcgctgtcgacgcatgcggatcctctagagtcgacatatgGTGATCCTGCTGAATTTCATTACGa                                                     <  2:2857818‑M2/25‑1 (MQ=255)
                        agcccgcgccgtcgctgtcgacgcatgcggatcctctagagtcgacatatgGTGATCCTGCTGAATTTCATTACGa                                                     >  1:4668427‑M2/52‑76 (MQ=255)
                           ccgcgccgtcgctgtcgacgcatgctgatcctctagagtcgacatatgGTGATCCTGATGCATTTCATTACGACCa                                                  >  2:4699455‑M2/49‑76 (MQ=255)
                           ccgcgccgtcgctgtcgacgcatgcggatcctctagagtcgacatatgGTGATCCTGCTGAATTTCATTACGACCa                                                  >  2:4523366‑M2/49‑76 (MQ=255)
                           ccgcgccgtcgctgtcgacgcatgcggatcctctagagtcgacatatgGTGATCCTGCTGAATTTCATTACGACCa                                                  <  2:3541924‑M2/28‑1 (MQ=255)
                           ccgcgccgtcgctgtcgacgcatgcggatcctctagagtcgacatatgGTGATCCTGCTGAATTTCATTACGACCa                                                  >  2:971653‑M2/49‑76 (MQ=255)
                           ccgcgccgtcgctgtcgacgcatgcggatcctctagagtcgacatatgGTGATCCTGCTGAATTTCATTACGACCa                                                  >  1:1116567‑M2/49‑76 (MQ=255)
                            cgcgccgtcgctgtcgacgcatgcggatcctctagagtcgacatatgGTGATCCTGCTGAATTTCATTACGACCAg                                                 >  2:2365005‑M2/48‑76 (MQ=255)
                            cgcgccgtcgctgtcgacgcatgcggatcctctagagtcgacatatgGTGATCCTGCTGAATTTCATTACGACCAg                                                 >  2:3336171‑M2/48‑76 (MQ=255)
                             gcgccgtcgctgtcgacgcatgcggatcctctagagtcgacatatgGTGATCCTGCTGAATTTCATTACGACCAGt                                                >  2:4315745‑M2/47‑76 (MQ=255)
                              cgccgtcgctgtcgacgcatgcggatcctctagagtcgacatatgGTGATCCTGCTGAATTTCATTACGACCAGTc                                               >  2:809503‑M2/46‑76 (MQ=255)
                               gccgtcgctgtcgacgcatgcggatcctctagagtcgacatatgGTGATCCTGCTGAATTTCATTACGACCAGTCt                                              >  1:301030‑M2/45‑76 (MQ=255)
                               gccgtcgctgtcgacgcatgcggatcctctagagtcgacatatgGTGATCCTGCTGAATTTCATTACGACCAGTCt                                              >  1:3082268‑M2/45‑76 (MQ=255)
                                ccgtcgctgtcgacgcatgcggatcctctagagtcgacatatgGTGATCCTGCTGAATTTCATTACGACCAGTCTa                                             >  1:165277‑M2/44‑76 (MQ=255)
                                ccgtcgctgtcgacgcatgcggatcctctagagtcgacatatgGTGATCCTGCTGAATTTCATTACGACCAGTCTa                                             <  2:560567‑M2/33‑1 (MQ=255)
                                 cgtcgctgtcgacgcatgcggatcctctagagtcgacatatgGTGATCCTGCTGAATTTCATTACGCCCAGTCTaa                                            >  2:1376‑M2/43‑76 (MQ=255)
                                 cgtcgctgtcgacgcatgcggatcctctagagtcgacatatgGAGATCCTACTGAATTTCATTACGACTAGTCTaa                                            >  1:5728809‑M2/43‑76 (MQ=255)
                                  gtcgctgtcgacgcatgcggatcctctagagtcgacatatgGTGATCCTGCTGAATTTCATTGCGACCAGTCTaaa                                           <  2:1529495‑M2/35‑1 (MQ=255)
                                  gtcgctgtcgacgcatgcggatcctctagagtcgacatatgGTGATCCTGCTGAATTTCATTACGACCAGTCTaaa                                           >  2:3052237‑M2/42‑76 (MQ=255)
                                   tcgcggtcgacgcatgcggatcctctagagtcgacatatgGTGATCCTGCTGAATTTCATTACGACCAGTCTaaaa                                          <  2:5249769‑M2/36‑1 (MQ=255)
                                    cgctgtcgacgcatgcggatcctctagagtcgacatatgGTGATCCTGCTGAATTTCATTACGACCAGTCTaaaaa                                         >  2:5384111‑M2/40‑76 (MQ=255)
                                     gctgtcgacgcatgcggatcctctagagtcgacatatgGTGATCCTGCTGAATTTCATTACGACCAGTCTAAAAAg                                        >  2:227830‑M2/39‑76 (MQ=255)
                                       tgtcgacgcatgcggatcctctagagtcgacatatgGTGATCCTGCTGAATTTCATTACGACCAGTCTAAAAAgcg                                      >  1:4681556‑M2/37‑76 (MQ=255)
                                       tgtcgacgcatgcggatcctctagagtcgacatatgGTGATCCTGCTGAATTTCATTACGACCAGTCTAAAAAgcg                                      <  2:4254327‑M2/40‑1 (MQ=255)
                                       tgtcgacgcatgcggatcctctagagtcgacatatgGTGATCCTGCTGAATTTCATTACGACCAGTCTAAAAAgcg                                      <  1:1894946‑M2/40‑1 (MQ=255)
                                        gtcgacgcatgcggatcctctagagtcgacatatgGTTATCCTGCTAAATTTCATTACGACCAGTCTAAAAAgcgc                                     <  1:3898188‑M2/41‑1 (MQ=255)
                                          cgacgcatgcggatcctctagagtcgacatatgGTGATCCTGCTGAATTTCATTACGACCAGTCTAAAAAGCGCCt                                   <  1:26767‑M2/43‑1 (MQ=255)
                                           gacgcatgcggatcctctagagtcgacatatgGTGATCCTGCTGAATTTCATTACGACCAGTCTAAAAAGCGCCTg                                  <  2:4668426‑M2/44‑1 (MQ=255)
                                           gacgcatgcggatcctctagagtcgacatatgGTGATCCTGCTGAATTTCATTACGACCAGTCTAAAAAGCGCCTg                                  <  2:2289902‑M2/44‑1 (MQ=255)
                                            acgcatgcggatcctctagagtcgacatatgGTGATCCTGCTGAATTTCATTACGACCAGTCTAAAAAGCGCCTGa                                 >  1:3567655‑M2/32‑76 (MQ=255)
                                            acgcatgcggatcctctagagtcgacatatgGTGATCCTGCTGAATTTCATTACGACCAGTCGAAAAAGCGCCTGa                                 >  1:5528990‑M2/32‑76 (MQ=255)
                                             cgcatgcggatcctctagagtcgacatatgGTGATCCTGCTGAATTTCATTACGACCAGTCTAAAAAGCGCCTGaa                                >  2:3654695‑M2/31‑76 (MQ=255)
                                             cgcatgcggatcctctagagtcgacatatgGTGATCCTGCTGAATTTCATTACGACCAGTCTAAAAAGCGCCTGaa                                <  2:4039961‑M2/46‑1 (MQ=255)
                                              gcatgcggatcctctagagtcgacatatgGTGATCCTGCTGAATTTCATTACGACCAGTCTAAAAAGCGCCTGAAt                               >  1:831052‑M2/30‑76 (MQ=255)
                                              gcatgcggatcctctagagtcgacatatgGTGATCCTGCTGAATTTCATTACGACCAGTCTAAAAAGCGCCTGAAt                               >  1:2023637‑M2/30‑76 (MQ=255)
                                               catgcggatcctctagagtcgacatatgGTGATCCTGCTGAATTTCATTACGACCAGTCTAAAAAGCGCCTGAAtt                              <  1:5119873‑M2/48‑1 (MQ=255)
                                               catgcggatcctctagagtcgacatatgGTGATCCTGCTGAATTTCATTACGACCAGTCTAAAAAGCGCCTGAAtt                              <  1:4776226‑M2/48‑1 (MQ=255)
                                               catgcggatcctctagagtcgacatatgGTGATCCTGCTGAATTTCATTACGACCAGTCTAAAAAGCGCCTGAAtt                              >  1:4509160‑M2/29‑76 (MQ=255)
                                               catgcggatcctctaaagtcgacatatgGTGATCCTGCTGAATTTCATTACGACCAGTCTAAAAAGCGCCTGAAtt                              <  1:4374932‑M2/48‑1 (MQ=255)
                                                atgcggatcctctagagtcgacatatgGTGATCCTGCTGAATTTCATTACGACCAGTCTAAAAAGCGCCTGAATTc                             <  2:276113‑M2/49‑1 (MQ=255)
                                                atgcggatcctctagagtcgacatatgGTGATCCTGCTGAATTTCATTACGACCAGTCTAAAAAGCGCCTGAATTc                             >  1:3443668‑M2/28‑76 (MQ=255)
                                                atgcggatcctctagagtcgacatatgGTGATCCTGCTGAATTTCATTACGACCAGTCTAAAAAGCGCCTGAATTc                             >  1:5245150‑M2/28‑76 (MQ=255)
                                                    ggatcctctagagtcgacatatgGTGATCCTGCTGAATTTCATTACGACCAGTCTAAAAAGCGCCTGAATTCGCGa                         >  2:839639‑M2/24‑76 (MQ=255)
                                                       tcctctagagtcgacatatgGTGATCCTGCTGAATTTCATTACGACCAGTCTAAAAAGCGCCTGAATTCGCGACCt                      >  2:3605811‑M2/21‑76 (MQ=255)
                                                        cctctagagtcgacatatgGTGATCCTGCTGAATTTCATTACGACCAGTCTAAAAAGCGCCTGAATTCGCGACCtt                     <  2:4478067‑M2/57‑1 (MQ=255)
                                                        cctctagagtcgacatatgGTGATCCTGCTGAATTTCATTACGACCAGTCTAAAAAGCGCCTGAATTCGCGACCtt                     <  2:1167684‑M2/57‑1 (MQ=255)
                                                         ctctagagtcgacatatgGTGATCCTGCTGAATTTCATTACGACCAGTCTAAAAAGCGCCTGAATTCGCGACCTtc                    <  2:4561778‑M2/58‑1 (MQ=255)
                                                         ctctagagtcgacatatgGTGATCCTGCTGAATTTCATTACGACCAGTCTAAAAAGCGCCTGAATTCGCGACCTtc                    <  2:2706223‑M2/58‑1 (MQ=255)
                                                         ctctagagtcgacatatgGTGATCCTGCTGAATTTCATTACGACCAGTCTAAAAAGCGCCTGAATTCGCGACCTtc                    <  2:2093042‑M2/58‑1 (MQ=255)
                                                          tctagagtcgacatatgGTGATCCTGCTGAATTTCATTACGACCAGTCTAAAAAGCGCCTGAATTCGCGAcc                       <  1:709381‑M2/55‑1 (MQ=255)
                                                          tctagagtcgacatatgGTGATCCTGCTGAATTTCATTACGACCAGTCTAAAAAGCGCCTGAATTCGCGAcc                       >  2:709381‑M2/18‑72 (MQ=255)
                                                           ctatagtcgacatatgGTGATCCTGCTGAATTTCATTACGACCAGTCTAAAAAGCGCCTGAATTCGCGCCCTtctc                  <  1:3612417‑M2/60‑1 (MQ=255)
                                                           ctagagtcgacatatgGTGATCCTGCTGAATTTCATTACGACCAGTCTAAAAAGCGCCTGAATTCGCGACCTtctc                  <  2:753422‑M2/60‑1 (MQ=255)
                                                           ctagagtcgacatatgGTGATCCTGCTGAATTTCATTACGACCAGTCTAAAAAGCGCCTGAATTCGCGACCTtctc                  <  1:3655953‑M2/60‑1 (MQ=255)
                                                            tagagtcgacatatgGTGATCCTGCTGAATTTCATTACGACCAGTCTAAAAAGCGCCTGAATTCGCGACCTTCTCg                 >  2:46286‑M2/16‑76 (MQ=255)
                                                             agagtcgacatatgGTGATCCTGCTGAATTTCATTACGACCAGTCTAAAAAGCGCCTGAATTCGCGACCTTCTCGt                >  1:4536721‑M2/15‑76 (MQ=255)
                                                              gagtcgacatatgGTGATCCTGCTGAATTTCATTACGACCAGTCTAAAAAGCGCCTGAATTCGCGATCTTCTCGtt               >  1:2914478‑M2/14‑76 (MQ=255)
                                                               agtcgacatatgGTGATCCTGCTGAATTTCATTACGACCAGTCTAAAAAGCGCCTGAATTCGCGACCTTCTCGTTa              <  2:4849387‑M2/64‑1 (MQ=255)
                                                                ntcgacatatgGTGATCCTGCTGAATTTCATTACGACCCGTCTAAAAAGCGCCTGAATTCGCGACCTTCTCGTTAc             >  1:1947990‑M2/12‑76 (MQ=255)
                                                                gtcgacatatgGTGATCCTGCTGAATTTCATTACGACCAGTCTAAAAAGCGCCTGAATTCGCGACCTTCTCGTTAc             >  1:3953757‑M2/12‑76 (MQ=255)
                                                                gtcgacataggGTGATCCTGCTGAATTTCATTACGACCAGTCTAAAAAGCGCCTGAATTCGCGACCTTCTCGTTAc             <  2:1947990‑M2/65‑1 (MQ=255)
                                                                 tctacatatgTTGATCCTGCTGAATTTCATTACGACCAGTCTAAAAAGCGCCTGAATTCGCGACCTTCTCGTTACt            <  2:5113119‑M2/66‑1 (MQ=255)
                                                                  cgacatatgGTGATCCTGCTGAATTTCATTACGACCAGTCTAAAAAGCGCCTGAATTCGCGACCTTCTCGTTACTg           <  1:4509747‑M2/67‑1 (MQ=255)
                                                                  cgacataggGTGATCCTGCTGAATTTCATTACGACCAGTCTAAAAAGCGCCTGAATTCGCGACCTTCTCGTTACTg           <  1:1352685‑M2/67‑1 (MQ=255)
                                                                   gacatatgGTGATCCTGCTGAATTTCATTACGACCAGTCTAAAAAGCGCCTGAATTCGCGACCTTCTCGTTACTGa          >  2:5782329‑M2/9‑76 (MQ=255)
                                                                   gacatatgGTGATCCTGCTGAATTTCATTACGACCAGTCTAAAAAGAGCCTGAATTCGCGACCTTCTCGGTactag          >  2:2159728‑M2/9‑74 (MQ=255)
                                                                     catatgGTGATCCTGCTGAATTTCATTACGACCAGTCTAAAAAGCGCCTGAATTCGCGACCTTCTCGTTACTGACa        >  1:4569264‑M2/7‑76 (MQ=255)
                                                                     catatgGTGATCCTGCTGAATTTCATTACGACCAGGCTAAAAAGCGACTGAATTCGCGACCGGCTCGTTACTGACa        <  1:1425999‑M2/70‑1 (MQ=255)
                                                                       tatgGTGATCCTGCTGAATTTCATTACGACCAGTCTAAAAAGCGCCTGAATTCGCGACCTTCTCGTTACTGACAgg      >  1:4092644‑M2/5‑76 (MQ=255)
                                                                        atgGTGATCCTGCTGAATTTCATTACGACCAGTCTAAAAAGCGCCTGAATTCGCGACCTTCTCGTTACTGACAGGa     >  1:1756010‑M2/4‑76 (MQ=255)
                                                                        atgGTGATCCTGCTGAATTTCATTACGACCAGTCTAAAAAGCGCCTGAATTCGCGACCTTCTCGTTACTGACAGGa     >  1:1444976‑M2/4‑76 (MQ=255)
                                                                        atgGTGATCCTGCTGAATTTCATTACGACCAGTCTAAAAAGCGCCTGAATTCGCGACCTTCTCGTTACTGACAGGa     <  2:1159670‑M2/73‑1 (MQ=255)
                                                                          gGTGATCCTGCTGAATTTCATTACGACCAGTCTAAAAAGCGCCTGAATTCGCGACCTTCTCGTTACTGACAGGaaa   >  1:4065914‑M2/2‑76 (MQ=255)
                                                                          gGTGATCCTGCTGAATTTCATTACGACCAGTCTAAAAAGCGCCTGAATTCGCGACCTTCTCGTTACTGACAGGaaa   >  1:775443‑M2/2‑76 (MQ=255)
                                                                           gtgATCCTGCTGAATTTCATTACGACCAGTCTAAAAAGCGCCTGAATTCGCGACCTTCTCGTTACTGACAGGaaaa  >  1:5005205/1‑76 (MQ=255)
                                                                           gtgATCCTGCTGAATTTCATTACGACCAGTCTAAAAAGCGCCTGAATTCGCGACCTTCTCGTTACTGACAGGaaaa  <  2:2129355/76‑1 (MQ=255)
                                                                           |                                                                           
CTCGTAACGCGCCAGCATCACGCGCCCACTGGATGCGCCGAGATCGACGGCGACACAATTGCGAAAGGTCATAATGTGATCCTGCTGAATTTCATTACGACCAGTCTAAAAAGCGCCTGAATTCGCGACCTTCTCGTTACTGACAGGAAAA  >  NC_000913/4094958‑4095108

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