Predicted mutation
evidence seq id position mutation annotation gene description
MC JC NC_000913 2,066,152 Δ12,778 bp insH6[yoeF] insH6, yoeA, insD3, insC3, yoeA, isrC, flu, yeeR, yeeS, yeeT, cbeA, cbtA, yeeW, [yoeF]

Missing coverage evidence...
   seq id start end size ←reads reads→ gene description
* * ÷ NC_000913 2066152 2078929 12778 40 [0] [0] 40 insH6–[yoeF] insH6,yoeA,insD3,insC3,yoeA,isrC,flu,yeeR,yeeS,yeeT,cbeA,cbtA,yeeW,[yoeF]

New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913 = 20661510 (0.000)40 (0.610) 37/476 0.3 100% intergenic (‑387/+154) cobU/insH6 cobinamide‑P guanylyltransferase/cobinamide kinase/IS5 family transposase and trans‑activator
?NC_000913 2078930 = 0 (0.000)coding (356/357 nt) yoeF putative uncharacterized protein YoeF

GAAATAAGATGGCATTACGGAATAGCCAGACTTTCTCGTGTTAACGTCCGACAAAATGCGTCATTGACAGCGGTTTCTGTGGATTGACACATTTCATTCAGATAATGAATTAATGCTACTGCCGGAACAACCAGCAAACGGGCATTGTGTTCTGAAATCCAGGCTATTGATTCAACGTCAGCGACAGATCTGCGCTGTAAAAGAACTTGTAAAACAACGTAATAGAATTATCCGGAATGGTGGCGACTATGCACTAGGGAAGGTGCGAATAAGTGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGTTTGAGCCGATTTTTTCTCCCGTAAATGCCTTGAATCAGCCTATTTAGACCGTTTCTTCGCCATTTAAGGCGTTATCCCCAGTTTTTAGTGAGATCTCTCCCACTGACGTATCATTTGGTCCGCCCGAAACAGGTTGGCCAGCGTGAATAACATCGCCAGTTGGTTATCGTT  >  NC_000913/2065902‑2066391
                                                                                                                                                                                                                                                         |                                                                                                                                                                                                                                                
gAAATAAGATGGCATTACGGAATAGCCAGACTTTCTCGTGTTAACGTCCGACAAAATGCGTCATTGACAGCGGTTTCTGTGGATTGACACATTTCATTCAGATAATGAATTAATGCTACTGCCGGAACAACCAGCAAACGGGCATTGTGTTCTGAAATCCAGGCTATTGATTCAACGTCAGCGACAGATCTGCGCTGTAAAAGAACTTGTAAAACAACGTAATAGAATTATCCGGAATGGTGGCGACTATa                                                                                                                                                                                                                                                 <  1:96985‑M1/251‑2 (MQ=255)
       gATGGCATTACGGAATAGCCAGACTTTCTCGTGTTAACGTCCGACAAAATGCGTCATTGACAGCGGTTTCTGTGGATTGACACATTTCATTCAGATAATGAATTAATGCTACTGCCGGAACAACCAGCAAACGGGCATTGTGTTCTGAAATCCAGGCTATTGATTCAACGTCAGCGACAGATCTGCGCTGTAAAAGAACTTGTAAAACAACGTAATAGAATTATCCGGAATGGTGGCGACTATaatttcga                                                                                                                                                                                                                                          <  1:332343‑M1/251‑9 (MQ=255)
                                      ggttAACGTCCGACAAAATGCGTCATTGACAGCGGTTTCTGTGGATTGACACATTTCATTCAGATAATGAATTAATGCTACTGCCGGAACAACCAGCAAACGGGCATTGTGTTCTGAAATCCAGGCTATTGATTCAACGTCAGCGACAGATCTGCGCTGTAAAAGAACTTGTAAAACAACGTAATAGAATTATCCGGAATGGTGGCGACTATaatttcgattggttctcggatgaggatgatgacggtgat                                                                                                                                                                                                           <  2:266884‑M1/250‑40 (MQ=255)
                                          aaCGTCCGACAAAATGCGTCATTGACAGCGGTTTCGGTGGATTGACACATTTCATACAGATAATGAATTAATGCTACTGCCGGAACAACCAGCAAACGGGCATTGTGTTCTGAAATCCAGGCTATTGATTCAACGTCAGCGACAGATCTGCGCTGTAAAAGAACTTGTAAAACAACGTAATAGAATTATCCGGAATGGTGGCGACTATaatttcgattggttctcggatgaggatgatgacggtgatacac                                                                                                                                                                                                       <  2:294391‑M1/251‑44 (MQ=255)
                                                    aaaaTGCGTCATTGACAGAGGTTTCTGTGGATTGACACATTTCATTCAGATAATGAATTAATGCTACTGCCGGAACAACCAGCAAACGGGCATTGTGTTCTGAAATCCAGGCTATTGATTCAACGTCAGCGACAGATCTGCGCTGTAAAAGAACTTGTAAAACAACGTAATAGAATTATCCGGAATGGTGGCGACTATaatttcgattggttctcggatgaggatgatgacggtgatacacagataagaat                                                                                                                                                                                             >  2:415986‑M1/1‑198 (MQ=255)
                                                                          ttCTGTGGATTGACACATTTCATTCAGATAATGAATTAATGCTACTGCCGGAACAACCAGCAAACGGGCATTGTGTTCTGAAATCCAGGCTATTGATTCAACGTCAGCGACAGATCTGCGCTGTAAAAGAACTTGTAAAACAACGTAATAGAATTATCCGGAATGGTGGCGACTATaatttcgattggttctcggatgaggatgatgacggtgatacacagataagaatcaatgagcaatgagggttgcc                                                                                                                                                                        <  2:60933‑M1/250‑75 (MQ=255)
                                                                                   ttGACACATTTCATTCAGATAATGAATTAATGCTACTGCCGGAACAACCAGCAAACGGGCATTGTGTTCTGAAATCCAGGCTATTGATTCAACGTCAGCGACAGATCTGCGCTGTAAAAGAACTTGTAAAACAACGTAATAGAATTATCCGGAATGGTGGCGACTATaatttcgattggttctcggatgaggatgatgacggtgatacacagataagaatcaatgagcaatgagggttgccgggcaaccct                                                                                                                                                              >  1:209088‑M1/1‑167 (MQ=255)
                                                                                    tGACACATTTCATTCAGATAATGAATTAATGCTACTGCCGGAACAACCAGCAAACGGGCATTGTGTTCTGAAATCCAGGCTATTGATTCAACGTCAGCGACAGATCTGCGCTGTAAAAGAACTTGTAAAACAACGTAATAGAATTATCCGGAATGGTGGCGACTATaatttcgattggttctcggatgaggctgatgacggtgatacacagttaagaatcaatgagcactgagggttgccgggcaaccctc                                                                                                                                                             >  2:481786‑M1/1‑166 (MQ=255)
                                                                                           tttCATTCAGATAATGAATTAATGCTACTGCCGGAACAACCAGCAAACGGGCATTGTGTTCTGAAATCCAGGCTATTGATTCAACGTCAGCGACAGATCTGCGCTGTAAAAGAACTTGTAAAACAACGTAATAGAATTATCCGGCATGGTGGCGACTATaatttcgattggttctcggatgaggatgatgacggtgatacacagataagaatcaatgagcaatgagggttgccgggcaaccctcattgaat                                                                                                                                                      >  2:404491‑M1/1‑159 (MQ=255)
                                                                                              cattcaGATAATGAATTAATGCTACTGCCGGAACAACCAGCAAACGGGCATTGTGTTCTGAAATCCAGGCTATTGATTCAACGTCAGCGACAGATCTGCGCTGTAAAAGAACTTGTAAAACAACGTAATAGAATTATCCGGAATGGTGGCGACTATaatttcgattggttctcggatgaggatgatgacggtgatacacagataagaatcaatgagcaatgagggttgccgggcaaccctcattgaataaa                                                                                                                                                   <  1:290473‑M1/251‑96 (MQ=255)
                                                                                                    gATAATGAATTAATGCTACTGCCGGAACAACCAGCAAACGGGCATTGTGTTCTGAAATCCAGGCTATTGATTCAACGTCAGCGACAGATCTGCGCTGTAAAAGAACTTGTAAAACAACGTAATAGAATTATCCGGAATGGTGGCGACTATaatttcgattggttctcggatgaggatgatgacggtgatacacagataagaatcaatgagcaatgagggttgccgggcaaccctcattgaataaaacggga                                                                                                                                             <  1:303642‑M1/251‑102 (MQ=255)
                                                                                                       aaTGAATTAATGCTACTGCCGGAACAACAAGCAAACGGGCATTGTGTTCTGAAATCCAGGCTATTGATGCAACGTCCGCGACAGCTCTGCGCTGTAAAAGAACTTGTAAAACAACGAAATAGAATTATCCGGAATGGTGGCGACTATaatttcgattggttctcggatgaggatgatgacggtgatacacagataagaatcaatgagcaatgagggttgccgggaaaccctcatagaataaatcgggaatt                                                                                                                                          >  1:412464‑M1/1‑147 (MQ=255)
                                                                                                              tAATGCTACTGCCGGAACAACCAGCAAACGGGCATTGTGTTCTGAAATCCAGGCTATTGATTCAACGTCAGCGACAGATCTGCGCTGTAAAAGAACTTGTAAAACAACGTAATAGAATTATCCGGAATGGTGGCGACTATaatttcgattggttctcggatgaggatgatgacggtgatacacagataagaatcaatgagcaatgagggttgccgggcaaccctcattgaataaaacgggaattacttcgc                                                                                                                                   >  1:284750‑M1/1‑140 (MQ=255)
                                                                                                              tAATGCTACTGCCGGAACAACCAGCAAACGGGCATTGTGTTCTGAAATCCAGGCTATTGATTCAACGTCAGCGACAGATCTGCGCTGTAAAAGAACTTGTAAAACAACGTAATAGAATTATCCGGAATGGTGGCGACTATaatttcgattggt                                                                                                                                                                                                                                     <  1:216944‑M1/153‑14 (MQ=255)
                                                                                                              tAATGCTACTGCCGGAACAACCAGCAAACGGGCATTGTGTTCTGAAATCCAGGCTATTGATTCAACGTCAGCGACAGATCTGCGCTGTAAAAGAACTTGTAAAACAACGTAATAGAATTATCCGGAATGGTGGCGACTATaatttcgattggt                                                                                                                                                                                                                                     >  2:216944‑M1/1‑140 (MQ=255)
                                                                                                                  gCTACTGCCGGAACAACCAGCAAACGGGCATTGTGTTCTGAAATCCAGGCTATTGATTCAACGTCAGCGACAGATCTGCGCTGTAAAAGAACTTGTAAAACAACGTAATAGAATTATCCGGAATGGTGGCGACTATaatttcgattggttctcggatgaggatgatgacggtgatacacagataagaatcaatgagcaatgagggttgccgggcaaccctcattgaataaaacgggaattacttcgcttcg                                                                                                                               >  2:154908‑M1/1‑136 (MQ=255)
                                                                                                                     aCTGCCGGAACAACCAGCAAACGGGCATTGTGTTCTGAAATCCAGGCTATTGATTCAACGTCAGCGACAGATCTGCGCTGTAAAAGAACTTGTAAAACAACGTAATAGAATTATCCGGAATGGTGGCGACTATaatttcgattggttctcggatgaggatgatgacggtgatacacagataagaatcaatgagcaatgagggttgccgggcaaccctcattgaataaaacgggaattacttcgcttcgcc                                                                                                                             <  1:144058‑M1/250‑118 (MQ=255)
                                                                                                                     aCTGCCGGAACAACCAGCAAACGGGCATTGTGTTCTGAAATCCAGGCTATTGATTCAACGTCAGCGACAGATCTGCGCTGTAAAAGAACGTGTAAAACAACGTAATAGAATTATCCGGAATGGTGGCGACTATaatttcgattggttctcggatgaggatgatgacggtgatacacagataagaatcaatgagcaatgagggttgccgggcaaccctcattgaataaaacgggaattacttcgcttcgcc                                                                                                                             <  2:64453‑M1/250‑118 (MQ=255)
                                                                                                                          cGGAACAACCAGCAAACGGGCATTGTGTTCTGAAATCCAGGCTATTGATTCAACGTCAGCGACAGATCTGCGCTGTAAAAGAACTTGTAAAACAACGTAATAGAATTATCCGGAATGGTGGCGACTATaatttcgattggttctcggatgaggatgatgacggtgatacacagataagaatcaatgagcaatgagggttgccgggcaaccctcattgaataaaacgggaattacttcgcttcgccgttttt                                                                                                                       >  1:366484‑M1/1‑128 (MQ=255)
                                                                                                                              acaacCAGCAAACGGGCATTGTGTTCTGAAATCCAGGCTATTGATTCAACGTCAGCGACAGATCTGCGCTGTAAAAGAACTTGTAAAACAACGTAATAGAATTATCCGGAATGGTGGCGACTATaatttcgattggttctcggatgaggatgatgacggtgatacacagataagaatcaatgagcaatgagggttgccgggcaaccctcattgaataaancgggaattacttcgcttcgccgtttttcatt                                                                                                                   >  2:231979‑M1/1‑124 (MQ=255)
                                                                                                                                aacCAGCAAACGGGCATTGTGTTCTGAAATCNAGGCTATTGATTCAACGTCAGCGACAGATCTGCGCTGTAAAAGAACTTGTAAAACAACGTAATAGAATTATCCGGAATGGTGGCGACTATaatttcgattggttctcggatgaggatgatgacggtgatacacagataagaatcaatgagcaatgagggttgccgggcaaccctcattgaataaaacgggaattacttcgcttcgccgtttttcatttc                                                                                                                 <  2:238212‑M1/251‑130 (MQ=255)
                                                                                                                                aacCAGCAAACGGGCATTGTGTTCTGAAATCCAGGCTATTGATTCAACGTCAGCGACAGATCTGCGCTGTAAAAGAACTTGTAAAACAACGTAATAGAATTATCCGGAATGGTGGCGACTATaatttcgattggttctcggatgaggatgatgacggtgatacacagataagaatcaatgagcaatgagggttgccgggcaaccctcattgaataaaacgggaattacttcgcttcgccgtttttcatttc                                                                                                                 >  2:108518‑M1/1‑122 (MQ=255)
                                                                                                                                    aGCAAACGGGCATTGTGTTCTGAAATCCAGGCTATTGATTCAACGTCAGCGACAGATCTGCGCTGTAAAAGAACTTGTAAAACAACGTAATAGAATTATCCGGAATGGTGGCGACTATaatttcgattggttctcggatgaggatgatgacggtgatacacagataagaatcaatga                                                                                                                                                                                       <  2:157291‑M1/177‑60 (MQ=255)
                                                                                                                                                          aaaTCCAGGCTATTGATTCAACGTCAGCGACAGATCTGCGCTGTAAAAGAACTTGTAAAACAACGTAATAGAATTATCCGGAATGGTGGCGACTATaatttcgattggttctcggatgaggatgatgacggtgatacacagataagaatcaatgagcaatgagggttgccgggcaaccctcattgaataaaacgggaattacttcgcttcgccgtttttcatttcgcacatagctttcagctttttgg                                                                                          >  1:307473‑M1/1‑96 (MQ=255)
                                                                                                                                                            aTCCAGGCTATTGATTCAACTTCAGCGACAGATCTGCGCTGTAAAAGAACTTGTAAAACAACGTAATAGAATTATCCGGAATGGTGGCGACTATaatttcgattggttctcggatgaggatgatgacggtgata                                                                                                                                                                                                          <  2:213688‑M1/134‑41 (MQ=255)
                                                                                                                                                            aTCCAGGCTATTGATTCAACGTCAGCGACAGATCTGCGCTGTAAAAGAACTTGTAAAACAACGTAATAGAATTATCCGGAATGGTGGCGACTATaatttcgattggttctcggatgaggatgatgacggtgata                                                                                                                                                                                                          >  1:213688‑M1/1‑94 (MQ=255)
                                                                                                                                                                           tCAACGTCAGCGACAGATTTGCGCTGTAAAAGAACTTGTAAAACAACGTAATAGAATTATCCGGAATGGTGGCGACTATaatttcgattggttctcggatgaggatgatgacggtgatacacagataagaatcaatgagcaatgagggttgccgggcaaccctcattgaataaaacgggaattact                                                                                                                                       >  2:416825‑M1/1‑79 (MQ=255)
                                                                                                                                                                           tCAACGTCAGCGACAGATCTGCGCTGTAAAAGAACTTGTAAAACAACGTAATAGAATTATCCGGAATGGTGGCGACTATaatttcgattggttctcggatgaggatgatgacggtgatacacagataagaatcaatgagcaatgagggttgccgggcaaccctcattgaataaaacgggaattact                                                                                                                                       <  1:416825‑M1/186‑108 (MQ=255)
                                                                                                                                                                                   aGCGACAGATCTGCGCTGTAAAAGAACTTGTAAAACAACGTAATAGAATTATCCGGAATGGTGGCGACTATaatttcgattggttctcggatgaggatgatgacggtgatacacagataagaatcaatgagcaatgagggttgccgggcaaccctcattgaataaaacgggaattacttcgcttcgccgtttttcatttcgcccatagctttcagctttttggagatatcgcggcgttctttggagagctc                                                              >  1:347107‑M1/1‑71 (MQ=255)
                                                                                                                                                                                    gCGACAGATCTGCGCTGTAAAAGAACTTGTAAAACAACGTAATAGAATTATCCGGAATGGTGGCGACTATaatttcgattggttctcggatgaggatgatgacggtgatacacagataagaatcaatgagcaatgagggttgccgggcaaccctcattgaataaaacgggaattacttcgcttcgccgtttttcatttcgcccatagctttcagctttttggagatatcgcggcgttctttggagagctcg                                                             <  2:281889‑M1/251‑182 (MQ=255)
                                                                                                                                                                                                       acagaACTTGTAAAACAACGTAATAGAATTATCCGGAATGGTGGCGCCTATaatttcgcttggttctcggatgaggatgatgacggtgatacacagataagaatcaatgagcaatgagggttgccgggcaaccctcattgaataaaacgggaattacttcgcttcgccgtttttcatttcgcccatagctttcagctttttggagatatcgcggcgttctttggcgagctcggcatttttgatgatgtaat                                          >  2:470832‑M1/3‑51 (MQ=255)
                                                                                                                                                                                                          gAACTTGTAAAACAACGTAATAGAATTATACGGAATGATGGCGACTATaatttcgattggttctcggatgaggatgatgacggtgatacacagataagaatcaatgagcaatgagggttgccgggcaaccctcattgaataaaacgggaattacttcgcttcgccgtttttcatttcgcacatagctttcagctttttggagatatcgcggcgttctttggagtgatcggcatttttgatgatgtactcgt                                       >  1:448224‑M1/1‑48 (MQ=255)
                                                                                                                                                                                                                   aaaCAACGTAATAGAATTATCCGGAATGGTGGCGACTATaatttcgattggttctcggatgaggatgatgacggtgatacacagataagaatcaatgagcaatgagggttgccgggcaaccctcattgaataaaacgggaattacttcgcttcgccgtttttcatttcgcccatagctttcagctttttggagatatcgcggcgttctttggagagctcggcatttttgatgatgtaatcgtcaacgcgat                              <  1:12210‑M1/251‑213 (MQ=255)
                                                                                                                                                                                                                               aGAATTATCCGGAATGGTGGCGACTATaatttcgattggttctcggatgaggatgatgacggtgatacacagataagaatcaatgagcaatgagggttgccgggcaaccctcattgaataaaacgggaattacttcgcttcgccgtttttcatttcgcccatagctttcagctttttggagatatcgcggcgttctttggagagctcggcatttttgatgatgtaatcgtcaacgcgatcttcatagtcac                  >  2:487611‑M1/1‑27 (MQ=255)
                                                                                                                                                                                                                                       ccGGAATGGTGGCGACTATaatttcgattggttctcggatgaggatgatgacggtgatacacagataagaatcaatgagcaatgagggttgccgggcaaccctcattgaataaaacgggaattacttcgcttcgccgtttttcatttcgcccatagctttcagctttttggagatatcgcggcgttctttggagagctcggcatttttgatgatgtaatcgtcaacgcgatcttcatagtcacctttcat           <  1:415986‑M1/250‑232 (MQ=255)
                                                                                                                                                                                                                                          gAATGGTGGCGACTATaatttcgattggttctcggctgaggatgatgacggtgatacacagataagaatcaatgagcaatgagggttgccgggcaaccctcattgaataaaacgggaattacttcgcttcgccgtttttcatttcgcccatagctttcagctttttggagatatcgcggcgttctttggagagctcggcatttttgatgatgtaatcgtcaacgcgatcttcatagtcacctttcatgctg       >  2:169618‑M1/1‑16 (MQ=255)
                                                                                                                                                                                                                                           aaTGGTGGCGACTATaatttcgattggttctcggatgaggatgatgacggtgatacacagataagaatcaatgagcaatgagggttgccgggcaaccctcattgaataaaacgggaattacttcgcttcgccgtttttcatttcgcccatagctttcagctttttggagatatcgcggcgttctttggagagctcggcatttttgatgatgtaatcgtcaacgcgatcttcatagtcccctttcatgctgg      >  2:121764‑M1/1‑15 (MQ=255)
                                                                                                                                                                                                                                               gtggCGACTATaatttcgattggttctcggatgaggatgatgacggtgatacacagataagaatcaatgagcaatgagggttgccgggcaaccctcattgaataaaacgggaattacttcgcttcgccgtttttcatttcgcccatagctttcagctttttggagatatcgcggcgttctttggagagctcggcatttttgatgatgtaatcgtcaacgcgatcttcatagtcacctttcatgctggcgat  >  2:96407‑M1/1‑11 (MQ=255)
                                                                                                                                                                                                                                                  gCGACTATaatttcgattggttctcggatgaggatgatgacggtgatacacagataagaatcaatgagcaatgagggttgccgggcaaccctcattgaataaaacgggaattacttcgcttcgccgtttttcatttcgcccatagctttcagctttttggagatatcgcggcgttctttggagagctc                                                              <  1:45876‑M1/188‑181 (MQ=255)
                                                                                                                                                                                                                                                  gCGACTATaatttcgattggttctcggatgaggatgatgacggtgatacacagataagaatcaatgagcaatgagggttgccgggcaaccctcattgaataaaacgggaattacttcgcttcgccgtttttcatttcgcccatagctttcagctttttggagatatcgcggcgttctttggagagctc                                                              >  2:45876‑M1/1‑8 (MQ=255)
                                                                                                                                                                                                                                                         |                                                                                                                                                                                                                                                
GAAATAAGATGGCATTACGGAATAGCCAGACTTTCTCGTGTTAACGTCCGACAAAATGCGTCATTGACAGCGGTTTCTGTGGATTGACACATTTCATTCAGATAATGAATTAATGCTACTGCCGGAACAACCAGCAAACGGGCATTGTGTTCTGAAATCCAGGCTATTGATTCAACGTCAGCGACAGATCTGCGCTGTAAAAGAACTTGTAAAACAACGTAATAGAATTATCCGGAATGGTGGCGACTATGCACTAGGGAAGGTGCGAATAAGTGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGTTTGAGCCGATTTTTTCTCCCGTAAATGCCTTGAATCAGCCTATTTAGACCGTTTCTTCGCCATTTAAGGCGTTATCCCCAGTTTTTAGTGAGATCTCTCCCACTGACGTATCATTTGGTCCGCCCGAAACAGGTTGGCCAGCGTGAATAACATCGCCAGTTGGTTATCGTT  >  NC_000913/2065902‑2066391

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