New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913 = 408127339 (1.060)6 (0.180) 6/224 2.8 13.2% coding (588/636 nt) fdoI formate dehydrogenase‑O, cytochrome b556 subunit
?NC_000913 = 4081294 43 (1.260)coding (567/636 nt) fdoI formate dehydrogenase‑O, cytochrome b556 subunit
Rejected: Frequency below/above cutoff threshold.

TATAAATTCTTGAGCCGAGGGAACAATAACGGCGGAATCATATCCGCCGTACGTTTCTCGCTTGAACCCAGCTCATCTTGCGGGATTATGCGAATACTCATTCAGCTTTCTTTTCCGTTGTCTTGCGGACCTCACGGTACCAGCGCGGGTGATGTTTCTTCGCCCAGGCGCTGGTAACCCATCCTTCCACCATGGCGGTAATCGTGCCTTTCACCCAAAGGGCGGCGTAGATATGCACCATGATAACCACAATTAACGCCACTGCGGCAAATGAATGCAGCATTAACGC  >  NC_000913/4081128‑4081416
                                                                                                                                                 |                                                                                                                                               
tataAATTCTTGAGCCGAGGGAACAATAACGGCGGAATCATATCCGCCGTACGTTTCTCGCTTGAACCCAGCTCATCTTGCGGGGTTATGCGAATACTCATTCAGCTTTCTTTTCCGTTGTCTTGCGGACCTCACGGTACCAGCGCgg                                                                                                                                               <  4:32104/148‑1 (MQ=255)
   aaaTTCTTGAGCCGAGGGAACAATAACGGCGGAATCATATCCGCCGTACGTTTCTCGCTTGAACCCAGCTCATCTTGCGGGATTATGCGAATACTCATTCAGCTTTCTTTTCCGTTGTCTTGCGGACCTCACGGTACCAGCGCGGGTGa                                                                                                                                           >  8:116056/1‑149 (MQ=255)
           gAGCCGAGGGAACAATAACGGCGGAATCATATCCGCCGTACGTTTCTCGCTTGAACCCAGCTCATCTTGCGGGATTATGCGAATACTCATTCAGCTTTCTTTTCCGTTGTCTTGCGGACCTCACGGTACCAGCGCGGGTGATGTttctt                                                                                                                                   >  7:163242/1‑149 (MQ=255)
              ccGAGGGAACAATAACGGCGGAATCATATCCGCCGTACGTTTCTCGCTTGAACCCAGCTCATCTTGCGGGATTATGCGAATACTCATTCAGCTTTCTTTTCCGTTGTCTTGCGGACCTCACGGTACCAGCGCGGGTGATGTTTCTTCGc                                                                                                                                >  6:179961/1‑149 (MQ=255)
                  gggAACAATAACGGCGGAATCATATCCGCCGTACGTTTCTCGCTTGAACCCAGCTCATCTTGCGGGATTATGCGAATACTCATTCAGCTTTCTTTTCCGTTGTCTTGCGGACCTCACGGTACCAGCGCGGGTGATGTTTCTTCGCCCAg                                                                                                                            <  8:163242/149‑1 (MQ=255)
                        aaTAACGGCGGAATCATATCCGCCGTACGTTTCTCGCTTGAACCCAGCTCATCTTGCGGGATTATGCGAATACTCATTCAGCTTTCTTTTCCGTTGTCTTGCGGAACTCACGGTACCAGCGCGGGTGATGTTTCTTCGCCCAGgcgctg                                                                                                                      <  6:176488‑M2/149‑7 (MQ=255)
                         aTAACGGCGGAATCATATCCGCCGTACGTTTCTCGCTTGAACCCAGCTCATCTTGCGGGATTATGCGAATACTCATTCAGCTTTCTTTTCCGTTGTCTTGCGGACCTCACGGTACCAGCGCGGGTGATGTTTCTTCGCCCAGgcgctgg                                                                                                                     >  7:116056‑M2/1‑142 (MQ=255)
                                   aaTCATATCCGCCGTACGTTTCTCGCTTGAACCCAGCTCATCTTGCGGGATTATGCGAATACTCATTCAGCTTTCTTTTCCGTTGTCTTGCGGACCTCACGGTACCAGCGCGGGTGATGTTTCTTCGCCCAGGCGCTGGTAACCCATcc                                                                                                           >  5:181294/1‑149 (MQ=255)
                                       atatCCGCCGTACGTTTCTCGCTTGAACCCAGCTCATCTTGCGGGATTATGCGAATACTCATTCAGCTTTCTTTTCCGTTGTCTTGCGGACCTCACGGTACCAGCGCctgggcgaagaaacatcacccgcgctggtaccgtgag                                                                                                            <  7:136766‑M1/144‑38 (MQ=255)
                                       atatCCGCCGTACGTTTCTCGCTTGAACCCAGCTCATCTTGCGGGATTATGCGAATACTCATTCAGCTTTCTTTTCCGTTGTCTTGCGGACCTCACGGTACCAGCGCctgggcgaagaaacatcacccgcgctggtaccgtgag                                                                                                            >  8:136766‑M1/1‑107 (MQ=255)
                                       atatCCGCCGTACGTTTCTCGCTTGAACCCAGCTCATCTTGCGGGATTATGCGAATACTCATTCAGCTTTCTTTTCCGTTGTCTTGCGGACCTCACGGTACCAGCGCGGGTGATGTTTCTTCGCCCAGGCGCTGGTAACCCATCCTTcc                                                                                                       >  3:100806/1‑149 (MQ=255)
                                          tCCGCCGTACGTTTCTCGCTTGAACCCAGCTCATCTTGCGGGATTATGCGAATACTCATTCAGCTTTCTTTTCCGTTGTCTTGCGGACCTCACGGTACCAGCGCGGGTGATGTTTCTTCGCCCAGGCGCTGGTAACCCATCCTTccac                                                                                                     >  2:7838/1‑148 (MQ=255)
                                                 tACGTTTCTCGCTTGAACCCAGCTCATCTTGCGGGATTATGCGAATACTCATTCAGCTTTCTTTTCCGTTGTCTTGCGGACCTCACGGTACCAGCGCctgggcgaagaaacatcacccgcgctggtacc                                                                                                                 >  5:75455‑M1/1‑97 (MQ=255)
                                                 tACGTTTCTCGCTTGAACCCAGCTCATCTTGCGGGATTATGCGAATACTCATTCAGCTTTCTTTTCCGTTGTCTTGCGGACCTCACGGTACCAGCGCctgggcgaagaaacatcacccgcgctggtacc                                                                                                                 <  6:75455‑M1/129‑33 (MQ=255)
                                                  aCGTTTCTCGCTTGAACCCAGCTCATCTTGCGGGATTATGCGAATACTCATTCAGCTTTCTTTTCCGTTGTCTTGCGGACCTCACGGTACCAGCGCctgggcgaagaaacatcacccgcgctgg                                                                                                                     >  2:38364‑M1/1‑96 (MQ=255)
                                                  aCGTTTCTCGCTTGAACCCAGCTCATCTTGCGGGATTATGCGAATACTCATTCAGCTTTCTTTTCCGTTGTCTTGCGGACCTCACGGTACCAGCGCctgggcgaagaaacatcacccgcgctgg                                                                                                                     <  1:38364‑M1/124‑29 (MQ=255)
                                                           gCTTGAACCCAGCTCATCTTGCGGGATTATGCGAATACTCATTCAGCTTTCTTTTCCGTTGTCTTGCGGACCTCACGGTACCAGCGCGGGTGATGTTTCTTCGCCCAGGCGCTGGTAACCCATcc                                                                                                           >  8:93973/1‑125 (MQ=255)
                                                           gCTTGAACCCAGCTCATCTTGCGGGATTATGCGAATACTCATTCAGCTTTCTTTTCCGTTGTCTTGCGGACCTCACGGTACCAGCGCGGGTGATGTTTCTTCGCCCAGGCGCTGGTAACCCATcc                                                                                                           <  7:93973/125‑1 (MQ=255)
                                                           gCTTGAACCCAGCTCATCTTGCGGGATTATGCGAATACTCATTCAGCTTTCTTTTCCGTTGTCTTGCGGACCTCACGGTACCAGCGCGGGTGATGTTTCTTCGCCCAGGCGCTGGTAACCCATCCTTCCACCATGGCGGTAATCGt                                                                                      <  4:100806/146‑1 (MQ=255)
                                                                 aCCCAGCTCATCTTGCGGGATTATGCGAATACTCATTCAGCTTTCTTTTCCGTTGTCTTGCGGACCTCACGGTACCAGCGCGGGTGATGTTTCTTCGCCCAGGCGCTGGTAACCCATCCTTCCACCATGGCGGTAATCGTGCCtt                                                                                 >  2:132769/1‑145 (MQ=255)
                                                                 aCCCAGCTCATCTTGCGGGATTATGCGAATACTCATTCAGCTTTCTTTTCCGTTGTCTTGCGGACCTCACGGTACCAGCGCGGGTGATGTTTCTTCGCCCAGGCGCTGGTAACCCATCCTTCCACCATGGCGGTAATCGTGCCtt                                                                                 <  1:132769/145‑1 (MQ=255)
                                                                   ccAGCTCATCTTGCGGGATTATGCGAATACTCATTCAGCTTTCTTTTCCGTTGTCTTGCGGACCTCACGGTACCAGCGCGGGTGATGTTTCTTCGCCCAGGCGCTGGTAACCCATCCTTCCACCATGGCGGTAATCGTGcc                                                                                   >  4:6423/1‑141 (MQ=255)
                                                                   ccAGCTCATCTTGCGGGATTATGCGAATACTCATTCAGCTTTCTTTTCCGTTGTCTTGCGGACCTCACGGTACCAGCGCGGGTGATGTTTCTTCGCCCAGGCGCTGGTAACCCATCCTTCCACCATGGCGGTAATCGTGcc                                                                                   <  3:6423/141‑1 (MQ=255)
                                                                          aTCTTGCGGGATTATGCGAATACTCATTCAGCTTTCTTTTCCGTTGTCTTGCGGACCTCACGGTACCAGCGCggg                                                                                                                                              <  1:39857/75‑1 (MQ=255)
                                                                          aTCTTGCGGGATTATGCGAATACTCATTCAGCTTTCTTTTCCGTTGTCTTGCGGACCTCACGGTACCAGCGCggg                                                                                                                                              >  2:39857/1‑75 (MQ=255)
                                                                          aTCTTGCGGGATTATGCGAATACTCATTCAGCTTTCTTTTCCGTTGTCTTGCGGACCTCACGGTACCAGCGCGGGTGATGTTTCTTCGCCCAGGCGCTGGTAACCCATCCTTCCACCATGGCGGTAATCGTGCCTTTCACCCAAAgg                                                                      <  1:123195/147‑1 (MQ=255)
                                                                             ttGCGGGATTATGCGAATACTCATTCAGCTTTCTTTTCCGTTGTCTTGCGGACCTCACGGTACCAGCGCGGGTGATGTTTCTTCGCCCAGGCGCTGGTAACCCATCCTTCCACCATGGCGGTAATCGTGCCTTTCACCCAAAGggcggc                                                                 >  6:12674/1‑149 (MQ=255)
                                                                               gCGGGATTATGCGAATACTCATTCAGCTTTCTTTTCCGTTGTCTTGCGGACCTCACGGTACCAGCGCGGGTGATGtt                                                                                                                                       <  2:99894/77‑1 (MQ=255)
                                                                               gCGGGATTATGCGAATACTCATTCAGCTTTCTTTTCCGTTGTCTTGCGGACCTCACGGTACCAGCGCGGGTGATGtt                                                                                                                                       >  1:99894/1‑77 (MQ=255)
                                                                                   gATTATGCGAATACTCATTCAGCTTTCTTTTCCGTTGTCTTGCGGACCTCACGGTACCAGCGCGGGTGATGTTTCTTCGCCCAGGCGCTGGTAACCCATCCTTCCACCATGGCGGTAATCGTGCCTTTCACCCAAAGGGCGGCGTAGa                                                            >  4:36333/1‑148 (MQ=255)
                                                                                                cTCATTCAGCTTTCTTTTCCGTTGTCTTGCGGACCTCACGGTACCAGCGCGGGTGATGTTTCTTCGCCCAGGCGCTGGTAACCCATCCTTCCACCATGGCGGTAATCGTGCCTTTCACCCAAAGGGCGGCGTAGATATGCACCATGAt                                               <  5:179961/148‑1 (MQ=255)
                                                                                                 tCATTCAGCTTTCTTTTCCGTTGTCTTGCGGACCTCACGGTACCAGCGCGGGTGATGTTTCTTCGCCCAGGCGCTGGTAACCCATCCTTCCACCATGGCGGTAATCGTGCCTTTCACCCAAAGGGCGGCGTAGATATGCACCATGATaa                                             >  6:74425/1‑149 (MQ=255)
                                                                                                    ttCAGCTTTCTTTTCCGTTGTCTTGCGGACCTCACGGTACCAGCGCGGGTGATGTTTCTTCGCCCAGGCGCTGGTAACCCATCCTTCCACCATGGCGGTAATCGTGCCTTTCACCCAAAGGGCGGCGTAGATATGCACCATGATAACca                                          <  1:7838/149‑1 (MQ=255)
                                                                                                      cAGCTTTCTTTTCCGTTGTCTTGCGGACCTCACGGTACCAGCGCGGGTGATGTTTCTTCGCCCAGGCGCTGGTAACCCATCCTTCCACCATGGCGGTAATCGTGCCTTTCACCCAAAGGGCGGCGTAGATATGCACCATGATAACcac                                         <  5:12674/148‑1 (MQ=255)
                                                                                                        gCTTTCTTTTCCGTTGTCTTGCGGACCTCACGGTACCAGCGCGGGTGATGTTTCTTCGCCCAGGCGCTGGTAACCCATCCTTCCACCATGGCGGTAATCGTGCCTTTCACCCAAAGGGCGGCGTAGATATGCACCATGATAACCACAAt                                      >  6:32256/1‑149 (MQ=255)
                                                                                                         ctttctttTCCGTTGTCTTGCGGACCTCACGGTACCAGCGCGGGTGATGTTTCTTCGCCCAGGCGCTGGTAACCCATCCTTCCACCATGGCGGTAATCGTGc                                                                                    >  3:64332/1‑102 (MQ=255)
                                                                                                         ctttctttTCCGTTGTCTTGCGGACCTCACGGTACCAGCGCGGGTGATGTTTCTTCGCCCAGGCGCTGGTAACCCATCCTTCCACCATGGCGGTAATCGTGc                                                                                    <  4:64332/102‑1 (MQ=255)
                                                                                                         ctttctttTCCGTTGTCTTGCGGACCTCACGGTACCAGCGCGGGTGATGTTTCTTCGCCCAGGCGCTGGTAACCCATCCTTCCACCATGGCGGTAATCGTGCCTTTCACCCAAAGGGCGGCGTAGATATGCACCATGATAACCACAAtt                                     <  5:32256/149‑1 (MQ=255)
                                                                                                                    gTTGTCTTGCGGACCTCACGGTACCAGCGCGGGTGATGTTTCTTCGCCCAGGCGCTGGTAACCCATCCTTCCACCATGGCGGTAATCGTGCCTTTCACCCAAAGGGCGGCGTAGATATGCACCATGATAACCACAATTAACGCCACTGc                          >  8:160276/1‑149 (MQ=255)
                                                                                                                        tCTTGCGGACCTCACGGTACCAGCGCGGGTGATGTTTCTTCGCCCAGGCGCTGGTAACCCATCCTTCCACCATGGCGGTAATCGTGCCTTTCACCCAAAGGGCGGCGTAGATATGCACCATGATAACCACAATTAACGCCACTGCGGc                       <  5:74425/148‑1 (MQ=255)
                                                                                                                         cTTGCGGACCTCACGGTACCAGCGCGGGTGATGTTTCTTCGCCCAGGCGCTGGTAACCCATCCTTCCACCATGGCGGTAATCGTGCCTTTCACCCAAAGGGCGGCGTAGATATGCACCATGATaa                                             <  1:196336/125‑1 (MQ=255)
                                                                                                                         cTTGCGGACCTCACGGTACCAGCGCGGGTGATGTTTCTTCGCCCAGGCGCTGGTAACCCATCCTTCCACCATGGCGGTAATCGTGCCTTTCACCCAAAGGGCGGCGTAGATATGCACCATGATaa                                             >  2:196336/1‑125 (MQ=255)
                                                                                                                          ttGCGGACCTCACGGTACCAGCGCGGGTGATGTTTCTTCGCCCAGgcgctg                                                                                                                      <  1:10848‑M2/51‑7 (MQ=255)
                                                                                                                          ttGCGGACCTCACGGTACCAGCGCGGGTGATGTTTCTTCGCCCAGgcgctg                                                                                                                      >  2:10848‑M2/1‑45 (MQ=255)
                                                                                                                                  cTCACGGTACCAGCGCGGGTGATGTTTCTTCGCCCAGgcgctggtaccgtgaggtccgcaagacaacggaaaagaaagctgaatgagtattcgcataatcccgcaagatgagctgggttcaagcgagaaacgtacggcggatat                 <  7:136766‑M2/144‑108 (MQ=255)
                                                                                                                                  cTCACGGTACCAGCGCGGGTGATGTTTCTTCGCCCAGgcgctggtaccgtgaggtccgcaagacaacggaaaagaaagctgaatgagtattcgcataatcccgcaagatgagctgggttcaagcgagaaacgtacggcggatat                 >  8:136766‑M2/1‑37 (MQ=255)
                                                                                                                                  cTCACGGTACCAGCGCGGGTGATGTTTCTTCGCCCAGGCGCTGGTAACCCATCCTTCCACCATGGCGGTAATCGTGCCTTTCACCCAAAGGGCGGCGTAGATATGCACCATGATAACCACAATTAACGCCACTGCGGCAAATGAATgca            <  3:36333/149‑1 (MQ=255)
                                                                                                                                     aCGGTACCAGCGCGGGTGATGTTTCTTCGCCCAGGCGCTGGTAACCCATCCTTCCACCATGGCGGTAATCGTGCCTTTCACCCAAAGGGCGGCGTAGATATGCACCATGATAACCACAATTAACGCCACTGCGGCAAATGAATgcagc          <  5:176488/148‑1 (MQ=255)
                                                                                                                                       ggTACCAGCGCGGGTGATGTTTCTTCGCCCAGgcgctggtaccgtgaggtccgcaagacaacggaaaagaaagctgaatgagtattcgcataatcccgcaagatgagctgggttcaagcgagaaacgta                           <  6:75455‑M2/129‑98 (MQ=255)
                                                                                                                                       ggTACCAGCGCGGGTGATGTTTCTTCGCCCAGgcgctggtaccgtgaggtccgcaagacaacggaaaagaaagctgaatgagtattcgcataatcccgcaagatgagctgggttcaagcgagaaacgta                           >  5:75455‑M2/1‑32 (MQ=255)
                                                                                                                                       ggTACCAGCGCGGGTGATGTTTCTTCGCCCAGGCGCTGGTAACCCATCCTTCCACCATGGCGGTAATCGTGCCTTTCACCCAAAGGGCGGCGTAGATATGCACCATGATAACCACAATTAACGCCACTGCGGCAAATGAATGCAGCAtt       <  6:181294/149‑1 (MQ=255)
                                                                                                                                        gTACCAGCGCGGGTGATGTTTCTTCGCCCAGGCGCTGGTAACCCATCCTTCCACCATGGCGGTAATCGTGCCTTTCACCCAAAGGGCGGCGTAGATATGCACCATGATAACCACAATTAACGCCACTGCGGCAAATGAATGCAGCAtt       >  6:127533/1‑148 (MQ=255)
                                                                                                                                           ccAGCGCctgggcgaagaaacatcacccgcgctggtaccgtgaggtccgcaagacaacggaaaagaaagctgaatgagtattcgcataatcccgcaagatgagctgggttcaagcgagaaacgtacggcggatatgattccgccgttat   >  7:116056‑M1/1‑7 (MQ=255)
                                                                                                                                           ccAGCGCGGGTGATGTTTCTTCGCCCAGgcgctggtaccgtgaggtccgcaagacaacggaaaagaaagctgaatgagtattcgcataatcccgcaagatgagctgggttcaagcgagaaacgt                            >  2:38364‑M2/1‑28 (MQ=255)
                                                                                                                                           ccAGCGCGGGTGATGTTTCTTCGCCCAGgcgctggtaccgtgaggtccgcaagacaacggaaaagaaagctgaatgagtattcgcataatcccgcaagatgagctgggttcaagcgagaaacgt                            <  1:38364‑M2/124‑97 (MQ=255)
                                                                                                                                           ccAGCGCGGGTGATGTTTCTTCGCCCAGGCGCTGGTAACCCATCCTTCCACCATGGCGGTAATCGTGCCTTTCACCCAAAGGGCGGCGTAGATATGCACCATGATAACCACAATTAACGCCACTGCGGCAAATGAATGCAGCATTAAc    <  7:160276/148‑1 (MQ=255)
                                                                                                                                            cAGCGCctgggcgaagaaacatcacccgcgctggtaccgtgagttccgcaagacaacggaaaagaaagctgaatgagtattcgcataatcccgcaagatgagctgggttcaagcgagaaacgtacggcggatatgattccgccgttatt  <  6:176488‑M1/149‑144 (MQ=255)
                                                                                                                                            cAGCGCctgggcgaagaaacatcacccgcgctggtaccgtgaggtccgcaa                                                                                                    >  2:10848‑M1/1‑6 (MQ=255)
                                                                                                                                            cAGCGCctgggcgaagaaacatcacccgcgctggtaccgtgaggtccgcaa                                                                                                    <  1:10848‑M1/51‑46 (MQ=255)
                                                                                                                                                 |                                                                                                                                               
TATAAATTCTTGAGCCGAGGGAACAATAACGGCGGAATCATATCCGCCGTACGTTTCTCGCTTGAACCCAGCTCATCTTGCGGGATTATGCGAATACTCATTCAGCTTTCTTTTCCGTTGTCTTGCGGACCTCACGGTACCAGCGCGGGTGATGTTTCTTCGCCCAGGCGCTGGTAACCCATCCTTCCACCATGGCGGTAATCGTGCCTTTCACCCAAAGGGCGGCGTAGATATGCACCATGATAACCACAATTAACGCCACTGCGGCAAATGAATGCAGCATTAACGC  >  NC_000913/4081128‑4081416

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.