Predicted mutation
evidence seq id position mutation annotation gene description
JC JC NC_000913 2,601,785 Δ1 bp :: IS186 (–) +5 bp coding (585‑589/618 nt) hyfA → hydrogenase 4, 4Fe‑4S subunit

New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913 2514273 =NA (NA)25 (0.750) 23/240 0.6 100% noncoding (1/1345 nt) IS186 repeat region
?NC_000913 2601785 = 0 (0.000)coding (585/618 nt) hyfA hydrogenase 4, 4Fe‑4S subunit
* ? NC_000913 = 2515616NA (NA)23 (0.680) 22/242 0.7 100% noncoding (1344/1345 nt) IS186 repeat region
?NC_000913 = 2601789 0 (0.000)coding (589/618 nt) hyfA hydrogenase 4, 4Fe‑4S subunit

GCGACCTGTGTGATTTCTTGCCAGAAGGTCCGGCCTGCGTTCGCGCTTGCCCGAATCAGGCGTTACGGCTGATCACCGGTGATAGCCTGCAACGTCAGATGAAAGAAAAACAGCGCCTTGCCGCAAGCTGGTTTGCCAATGGCGGGGAGGATCCCCTTTCCCTCACTCAGGAGCAACGCTAATGGATGCCCTGCAATTATTAACCTGGTCGCTGATTCTCTATCTGTTTGCTAGTCTGGCTTCGCTGTTTTTACTCGGTCTGGACAGACTG  >  NC_000913/2601637‑2601907
                                                                                                                                                    |                                                                                                                          
agccccggcattacgtgccgaagtatccagttcttcgggcttaccaatatgggccagaattgctgaccagttatcgtgagagtaattcatcggcacgttaaatcatatcaggcgtaataccacaacccttaagttagcgcttatgggag                                                                                                                            <  1:4771‑M2/1‑1 (MQ=255)
   cccggcattacgtgccgaagtatccagttcttcgggcttaccaatatgggccagaattgctgaccagttatcgtgagagtaattcatcggcacgttaaatcatatcaggcgtaataccacaacccttaagttagcgcttatgggaGGAt                                                                                                                         >  4:67733‑M2/146‑149 (MQ=255)
        tgtgATTTCTTGCCAGAAGGTCCGGCCTGCGTTCGCGCTTGCCCGAATCAGGCGTTACGGCTGATCACCGGTGATAGCCTGCAACGTCAGATGAAAGAAAAACAGCGCCTTGCCGCAAGCTGGTTTGCCAATGGCGGGGAGGATCccca                                                                                                                    <  4:113315‑M2/149‑5 (MQ=255)
              gtgccgaagtatccagttcttcgggcttaccaatatgggccagaattgctgaccagttatcgtgagagtaattcatcggcacgttaaatcatatcaggcgtaataccacaacccttaagttagcgcttatgggaGGATCCCCTTTCCct                                                                                                              >  1:160009‑M2/135‑149 (MQ=255)
                         tccagttcttcgggcttaccaatatgggccagaattgctgaccagttatcgtgagagtaattcatcggcacgttaaatcatatcaggcgtaataccacaacccttaagttagcgcttatgggaGGATccc                                                                                                                      <  4:125663‑M2/7‑1 (MQ=255)
                         tccagttcttcgggcttaccaatatgggccagaattgctgaccagttatcgtgagagtaattcatcggcacgttaaatcatatcaggcgtaataccacaacccttaagttagcgcttatgggaGGATccc                                                                                                                      >  3:125663‑M2/124‑130 (MQ=255)
                          ggTCCGGCCTGCGTTCGCGCTTGCCCGAATCAGGCGTTACGGCTGATCACCGGTGATAGCCTGCAACGTCAGATGAAAGAAAAACAGCGCCTTGCCGCAAGCTGGTTTGCCAATGGCGGGGAGGATCcccataagcgctaacttaaggg                                                                                                  <  4:134105‑M2/149‑23 (MQ=255)
                                gCCTGCGTTCGCGCTTGCCCGAATCAGGCGTTACGGCTGATCACCGGTGATAGCCTGCAACGTCAGATGAAAGAAAAACAGCGCCTTGCCGCAAGCTGGTTTGCCAATGGCGGGGAGGATCcccataagcgctaacttaagggttgaac                                                                                            >  3:134105‑M2/1‑121 (MQ=255)
                                 ttcgggcttaccaatatgggccagaattgctgaccagttatcgtgagagtaattcatcggcacgttaaatcatatcaggcgtaataccacaacccttaagttagcgcttatgggaGGATCCCCTTTCCCTCACTCAGGAGCAACGCTaa                                                                                           >  8:106633‑M2/116‑149 (MQ=255)
                                                  ccGAATCAGGCGTTACGGCTGATCACCGGTGATAGCCTGCAACGTCAGATGAAAGAAAAACAGCGCCTTGCCGCAAGCTGGTTTGCCAATGGCGGGGAGGATCcccataagcgctaacttaagggttgaaccat                                                                                         >  1:21801‑M2/1‑103 (MQ=255)
                                                  ccGAATCAGGCGTTACGGCTGATCACCGGTGATAGCCTGCAACGTCAGATGAAAGAAAAACAGCGCCTTGCCGCAAGCTGGTTTGCCAATGGCGGGGAGGATCcccataagcgctaacttaagggttgaaccat                                                                                         <  2:21801‑M2/134‑32 (MQ=255)
                                                       agaattgctgaccagttatcgtgagagtaattcatcggcacgttaaatcatatcaggcgtaataccacaacccttaagttagcgcttatgggaGGATCCCCTTTCCCTCACTCAGGAGCAACGCTAATGGATGCCCTGCAATTATTAAc                                                                     <  2:160009‑M2/56‑1 (MQ=255)
                                                       agaattgctgaccagttatcgtgagagtaattcatcggcacgttaaatcatatcaggcgtaataccacaacccttaagttagcgcttatgggaGGATCCCCTTTCCCTCACTCAGGAGCAACGCTAATGGATGCCCTGCAATTATTAAc                                                                     <  3:67733‑M2/56‑1 (MQ=255)
                                                          ggCGTTACGGCTGATCACCGGTGATAGCCTGCAACGTCAGATGAAAGAAAAACAGCGCCTTGCCGCAAGCTGGTTTGCCAATGGCGGGGAGGATCcccataagcgctaacttaagggttgaaccatctgaagaatgcgacgcctcggtg                                                                  <  7:124248‑M2/149‑55 (MQ=255)
                                                          ggCGTTACGGCTGATCACCGGTGATAGCCTGCAACGTCAGATGAAAGAAAAACAGCGCCTTGCCGCAAGCTGGTTTGCCAAAGGCGGGGAGGATCcccataagcgctaacttaagggttgacccatctgaagaatgcgacgcctcggtg                                                                  <  8:72640‑M2/149‑55 (MQ=255)
                                                               tgaccagttatcgtgagagtaattcatcggcacgttaaatcatatcaggcgtaataccacaacccttaagttagcgcttatgggaGGATCCCCTTTCCCTCACTCAGGAGCAACGCTAATGGATGCCCTGCAATTATTAACCTGGTCGc                                                             >  3:7727‑M2/86‑149 (MQ=255)
                                                                 cGGCTGATCACCGGTGATAGCCTGCAACGTCAGATGAAAGAAAAACAGCGCCTTGCCGCAAGCTGGTTTGCCAATGGCGGGGAGGATCcccataagcgctaacttaag                                                                                                    <  5:130464‑M2/108‑21 (MQ=255)
                                                                 cGGCTGATCACCGGTGATAGCCTGCAACGTCAGATGAAAGAAAAACAGCGCCTTGCCGCAAGCTGGTTTGCCAATGGCGGGGAGGATCcccataagcgctaacttaag                                                                                                    >  6:130464‑M2/1‑88 (MQ=255)
                                                                  ggCTGATCACCGGTGATAGCCTGCAACGTCAGATGAAAGAAAAACAGCGCCTTGCCGCAAGCTGGTTTGCCAATGGCGGGGAGGATCcccataagcgctaacttaagggttgaaccatctgaagaatgcgacgcctcggtgcctcgtta                                                          >  3:113315‑M2/1‑87 (MQ=255)
                                                                    agttatcgtgagagtaattcatcggcacgttaaatcatatcaggcgtaataccacaacccttaagttagcgcttatgggaGGATCCCCTTTCCCTCACTCAGGAGCAACGCTAATGGATGCCCTGCAATTATTAACCTGGTCGCTGAtt                                                        <  7:106633‑M2/69‑1 (MQ=255)
                                                                          aCCGGTGATAGCCTGCAACGTCAGATGAAAGAAAAACAGCGCCTTGCCGCAAGCTGGTTTGCCAATGGCGGGGAGGATCcccataagcgctaacttaagggttgaaccatctgaagaatgcgacgcctcggtgcctcgttaagacgatg                                                  <  3:19212‑M2/149‑71 (MQ=255)
                                                                            tgagagtaattcatcggcacgttaaatcatatcaggcgtaataccacaacccttaagttagcgcttatgggaGGATCCCCTTTCCCTCACTCAGGAGCAACGCTAATGGATGCCCTGCAATTATTAACCTGGTCGCTGATTCTCTATCt                                                >  5:158188‑M2/73‑149 (MQ=255)
                                                                                       tGCAACGTCAGATGAAAGAAAAACAGCGCCTTGCCGCAAGCTGGTTTGCCAATGGCGGGGAGGATCcccataa                                                                                                                 >  6:126363‑M2/1‑66 (MQ=255)
                                                                                       tGCAACGTCAGATGAAAGAAAAACAGCGCCTTGCCGCAAGCTGGTTTGCCAATGGCGGGGAGGATCcccataa                                                                                                                 <  5:126363‑M2/73‑8 (MQ=255)
                                                                                          aaCGTCAGATGAAAGAAAAACAGCGCCTTGCCGCAAGCTGGTTTGCCAATGGCGGGGAGGATCcccataagcgctaacttaagggttga                                                                                              >  4:16430‑M2/1‑63 (MQ=255)
                                                                                              tCAGATGAAAGAAAATCAGCTCATTGCCGCAAGCTGGTTTGCCAATGGCGGGGAGGATCccaataagcgctaacttaagggttgaaccatctgaagaatgcga                                                                            >  2:64288‑M2/1‑59 (MQ=255)
                                                                                              tCAGATGAAAGAAAAACAGCGCCTTGCCGCAAGCTGGTTTGCCAATGGCGGGGAGGATCcccataagcgctaacttaagggttgaaccatctgaagaatgcga                                                                            <  1:64288‑M2/103‑45 (MQ=255)
                                                                                               cgttaaatcatatcaggcgtaataccacaacccttaagttagcgcttatgggaGGATCCCCTTTCCCTCACTCAGGAGCAACGCTAATGGATGCCCTGCAATTATTAACCTGGTCGCTGATTCTCTATCTGt                                              >  2:94944‑M2/54‑132 (MQ=255)
                                                                                               cgttaaatcatatcaggcgtaataccacaacccttaagttagcgcttatgggaGGATCCCCTTTCCCTCACTCAGGAGCAACGCTAATGGATGCCCTGCAATTATTAACCTGGGCGCTGATTCTCTATCTGt                                              <  1:94944‑M2/79‑1 (MQ=255)
                                                                                                 ttaaatcatatcaggcgtaataccacaacccttaagttagcgcttatgggaGGATCCCCTTTCCCTCACTCAGGAGCAACGCTAATGGATGCCCTGCAATTATTAACCTGGTCGCTGATTCTCTATCTGTTTGCTAGTCTGGCTTCGCt                           >  5:30573‑M2/52‑149 (MQ=255)
                                                                                                  aTGAAAGTAAAACAGCGCCTTGCCGCAAGCTGGTTTGCCAATGGCGGGGAGGATCcccatattcgctaacttaagggtggaaccatctgaagaatgcgccgcctcggtgcctcgttaagacgatgcctcgcgttcttcaattgcgtttt                          >  7:72640‑M2/1‑55 (MQ=255)
                                                                                                    aatcatatcaggcgtaataccacaacccttaagttagcgcttatgggaGGATCCCCTTTCCCTCACTCAGGAGCAACGCTAATGGATGCCCTGCAATTATTAACCTGGTCGCTGATTctct                                                    <  1:116051‑M2/73‑1 (MQ=255)
                                                                                                    aatcatatcaggcgtaataccacaacccttaagttagcgcttatgggaGGATCCCCTTTCCCTCACTCAGGAGCAACGCTAATGGATGCCCTGCAATTATTAACCTGGTCGCTGATTctct                                                    >  2:116051‑M2/49‑121 (MQ=255)
                                                                                                     atcatatcaggcgtaataccacaacccttaagttagcgcttatgggaGGATCCCCTTTCCCTCACTCAGGAGCAACGCTAATGGATGCCCTGCAATTATTAACCTGGTCGCTGATTCTCTATCTGTTTGCTAGTCTGGCTTCGCTGttt                       <  6:158188‑M2/102‑1 (MQ=255)
                                                                                                      tcatatcaggcgtaataccacaacccttaagttagcgcttatgggaGGATCCCCTTTCCCTCACTCAGGAGCAACGCTAATGGATGCCCTGCAATTATTAACCTGGTCGCTGATTCTCTATCTGTTTGCTAGTCTGGCTTCGCTGtttt                      >  5:169450‑M2/47‑149 (MQ=255)
                                                                                                         tatcaggcgtaataccacaacccttaagttagcgcttatgggaGGATCCCCttt                                                                                                                  <  1:70885‑M2/11‑1 (MQ=255)
                                                                                                         tatcaggcgtaataccacaacccttaagttagcgcttatgggaGGATCCCCttt                                                                                                                  >  2:70885‑M2/44‑54 (MQ=255)
                                                                                                           aaaCAGCGCCTTGCCGCAAGCTGGTTTGCCAATGGCGGGGAGGATCcccataagcgctaacttaagggttgaaccatctgaagaatgcgacgcctcggtgcctcgttaagacgatgcctcgcgttcttcaattgcgttttgtaggctgt                 <  2:4710‑M2/149‑104 (MQ=255)
                                                                                                                gcgcCTTGCCGCAAGCTGGTTTGCCAATGGCGGGGAGGATCcccataagcgctaacttaagggttgaaccatctgaagaatgcgacgcctcggtgcctcgttaagacgatgcctcgcgttcttcaattgcgttttgtaggctgtcaggg            >  8:124248‑M2/1‑41 (MQ=255)
                                                                                                                  taataccacaacccttaagttagcgcttatgggaGGATCCCCTTTCCCTCACTCAGGAGCAACGCTAATGGATGCCCTGCAATTATTAACCTGGTCGCTGATTCTCTATCTGTTTGc                                          <  5:169133‑M2/83‑1 (MQ=255)
                                                                                                                  taataccacaacccttaagttagcgcttatgggaGGATCCCCTTTCCCTCACTCAGGAGCAACGCTAATGGATGCCCTGCAATTATTAACCTGGTCGCTGATTCTCTATCTGTTTGc                                          >  6:169133‑M2/35‑117 (MQ=255)
                                                                                                                   ccTTGCCGCAAGCTGGTTTGCCAATGGCGGGGAGGATCcccataagcgctaacttaagggttgaaccatctgaagaatgcgacgcctcggtgcctcgttaagacgatgcctcgcgttcttcaattgcgttttgtaggctgtcagggata         >  4:19212‑M2/1‑38 (MQ=255)
                                                                                                                    ataccacaacccttaagttagcgcttatgggaGGATCACCTTTCCCTCACTCAGGAGCAACGCTAATGGATGCCCTGCAATTATTAACCTGGTCGCTGATTCTCTATCTGTTTGCTAGTCTGGCTTCGCTGTTTTTACTCGGTCTGGAc        <  4:7727‑M2/117‑1 (MQ=255)
                                                                                                                          gCAAGCTGGTTTGCCAATGGCGGGGAGGATCcccataagcgctaacttaagggttgaaccatctgaagaatgcgacgcctcggtgcctcgttaagacgatgcctcgcgttcttcaattgcgttttgtaggctgtcagggatactgtccc  >  1:4710‑M2/1‑31 (MQ=255)
                                                                                                                          gCAAGCTGGTTTGCCAATGGCGGGGAGGATCcccataagcgctaacttaagggttgaaccatctgaagaatgcgac                                                                           >  5:81767‑M2/1‑31 (MQ=255)
                                                                                                                          gCAAGCTGGTTTGCCAATGGCGGGGAGGATCcccataagcgctaacttaagggttgaaccatctgaagaatgcgac                                                                           <  6:81767‑M2/76‑46 (MQ=255)
                                                                                                                                               tgggaGGATCCCCTTTCCCTCACTCAGGAGCAACGCTAATGGATGCCCTGCAATTATTAACCTGGTCGCTGATTCTCTATCTGTTTGCTAGTCTGGCTTCGCTGTTTTTAc                   >  7:10228‑M2/6‑111 (MQ=255)
                                                                                                                                               tgggaGGATCCCCTTTCCCTCACTCAGGAGCAACGCTAATGGATGCCCTGCAATTATTAACCTGGTCGCTGATTCTCTATCTGTTTGCTAGTCTGGCTTCGCTGTTTTTAc                   <  8:10228‑M2/106‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                          
GCGACCTGTGTGATTTCTTGCCAGAAGGTCCGGCCTGCGTTCGCGCTTGCCCGAATCAGGCGTTACGGCTGATCACCGGTGATAGCCTGCAACGTCAGATGAAAGAAAAACAGCGCCTTGCCGCAAGCTGGTTTGCCAATGGCGGGGAGGATCCCCTTTCCCTCACTCAGGAGCAACGCTAATGGATGCCCTGCAATTATTAACCTGGTCGCTGATTCTCTATCTGTTTGCTAGTCTGGCTTCGCTGTTTTTACTCGGTCTGGACAGACTG  >  NC_000913/2601637‑2601907

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