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.

GGGAACAATAACGGCGGAATCATATCCGCCGTACGTTTCTCGCTTGAACCCAGCTCATCTTGCGGGATTATGCGAATACTCATTCAGCTTTCTTTTCCGTTGTCTTGCGGACCTCACGGTACCAGCGCGGGTGATGTTTCTTCGCCCAGGCGCTGGTAACCCATCCTTCCACCATGGCGGTAATCGTGCCTTTCACCCAAAGGGCGGCGTAGATATGCACCATGATAACCACAATTAACGCCACTGCGGCAAATGAATGCAGCATTAACGCGAATCGGATCACCGGGATTGAG  >  NC_000913/4081146‑4081438
                                                                                                                                                    |                                                                                                                                                
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)
                     atatCCGCCGTACGTTTCTCGCTTGAACCCAGCTCATCTTGCGGGATTATGCGAATACTCATTCAGCTTTCTTTTCCGTTGTCTTGCGGACCTCACGGTACCAGCGCGGGTGATGTTTCTTCGCCCAGGCGCTGGTAACCCATCCTTcc                                                                                                                             >  3:100806/1‑149 (MQ=255)
                        tCCGCCGTACGTTTCTCGCTTGAACCCAGCTCATCTTGCGGGATTATGCGAATACTCATTCAGCTTTCTTTTCCGTTGTCTTGCGGACCTCACGGTACCAGCGCGGGTGATGTTTCTTCGCCCAGGCGCTGGTAACCCATCCTTccac                                                                                                                           >  2:7838/1‑148 (MQ=255)
                                         gCTTGAACCCAGCTCATCTTGCGGGATTATGCGAATACTCATTCAGCTTTCTTTTCCGTTGTCTTGCGGACCTCACGGTACCAGCGCGGGTGATGTTTCTTCGCCCAGGCGCTGGTAACCCATcc                                                                                                                                 <  7:93973/125‑1 (MQ=255)
                                         gCTTGAACCCAGCTCATCTTGCGGGATTATGCGAATACTCATTCAGCTTTCTTTTCCGTTGTCTTGCGGACCTCACGGTACCAGCGCGGGTGATGTTTCTTCGCCCAGGCGCTGGTAACCCATcc                                                                                                                                 >  8:93973/1‑125 (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)
                                                        aTCTTGCGGGATTATGCGAATACTCATTCAGCTTTCTTTTCCGTTGTCTTGCGGACCTCACGGTACCAGCGCGGGTGATGTTTCTTCGCCCAGGCGCTGGTAACCCATCCTTCCACCATGGCGGTAATCGTGCCTTTCACCCAAAgg                                                                                            <  1:123195/147‑1 (MQ=255)
                                                           ttGCGGGATTATGCGAATACTCATTCAGCTTTCTTTTCCGTTGTCTTGCGGACCTCACGGTACCAGCGCGGGTGATGTTTCTTCGCCCAGGCGCTGGTAACCCATCCTTCCACCATGGCGGTAATCGTGCCTTTCACCCAAAGggcggc                                                                                       >  6:12674/1‑149 (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                                       >  8:136766‑M2/1‑37 (MQ=255)
                                                                                                                cTCACGGTACCAGCGCGGGTGATGTTTCTTCGCCCAGgcgctggtaccgtgaggtccgcaagacaacggaaaagaaagctgaatgagtattcgcataatcccgcaagatgagctgggttcaagcgagaaacgtacggcggatat                                       <  7:136766‑M2/144‑108 (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)
                                                                                                                         ccAGCGCGGGTGATGTTTCTTCGCCCAGgcgctggtaccgtgaggtccgcaagacaacggaaaagaaagctgaatgagtattcgcataatcccgcaagatgagctgggttcaagcgagaaacgt                                                  >  2:38364‑M2/1‑28 (MQ=255)
                                                                                                                         ccAGCGCGGGTGATGTTTCTTCGCCCAGgcgctggtaccgtgaggtccgcaagacaacggaaaagaaagctgaatgagtattcgcataatcccgcaagatgagctgggttcaagcgagaaacgt                                                  <  1:38364‑M2/124‑97 (MQ=255)
                                                                                                                         ccAGCGCGGGTGATGTTTCTTCGCCCAGGCGCTGGTAACCCATCCTTCCACCATGGCGGTAATCGTGCCTTTCACCCAAAGGGCGGCGTAGATATGCACCATGATAACCACAATTAACGCCACTGCGGCAAATGAATGCAGCATTAAc                          <  7:160276/148‑1 (MQ=255)
                                                                                                                                 ggTGATGTTTCTTCGCCCAGGCGCTGGTAACCCATCCTTCCACCATGGCGGTAATCGTGCCTTTCACCCAAAGGGCGGCGTAGATATGCACCATGATAACCACAATTAACGCCACTGCGGCAAATGAATGCAGCATTAAc                          >  6:49031/1‑140 (MQ=255)
                                                                                                                                 ggTGATGTTTCTTCGCCCAGGCGCTGGTAACCCATCCTTCCACCATGGCGGTAATCGTGCCTTTCACCCAAAGGGCGGCGTAGATATGCACCATGATAACCACAATTAACGCCACTGCGGCAAATGAATGCAGCATTAAc                          <  5:49031/140‑1 (MQ=255)
                                                                                                                                       gTTTCTTCGCCCAGGCGCTGGTAACCCATCCTTCCACCATg                                                                                                                       >  3:102671/1‑41 (MQ=255)
                                                                                                                                       gTTTCTTCGCCCAGGCGCTGGTAACCCATCCTTCCACCATg                                                                                                                       <  4:102671/41‑1 (MQ=255)
                                                                                                                                       gTTTCTTCGCCCAGGCGCTGGTAACCCATCCTTCCACCATGGCGGTAATCGTGCCTTTCACCCAAAGGGCGGCGTa                                                                                    <  6:33805/76‑1 (MQ=255)
                                                                                                                                       gTTTCTTCGCCCAGGCGCTGGTAACCCATCCTTCCACCATGGCGGTAATCGTGCCTTTCACCCAAAGGGCGGCGTa                                                                                    >  5:33805/1‑76 (MQ=255)
                                                                                                                                        tttCTTCGCCCAGGCGCTGGTAACCCATCCTTCCACCATGGCGGTAATCGTGCCTTTCACCCAAAGGGCGGCGTAGATATGCACCATGATAACCACAATTAACGCCACTGCGGCAAATGAATgcag                                 <  2:120982/126‑1 (MQ=255)
                                                                                                                                        tttCTTCGCCCAGGCGCTGGTAACCCATCCTTCCACCATGGCGGTAATCGTGCCTTTCACCCAAAGGGCGGCGTAGATATGCACCATGATAACCACAATTAACGCCACTGCGGCAAATGAATgcag                                 >  1:120982/1‑126 (MQ=255)
                                                                                                                                        tttCTTCGCCCAGGCGCTGGTAACCCATCCTTCCACCATGGCGGTAATCGTGCCTTTCACCCAAAGGGCGGCGTAGATATGCACCATGATAACCACAATTAACGCCACTGCGGCAAATGAATGCAGCATTAACGCGAATCGGATCACCg          >  3:25819/1‑149 (MQ=255)
                                                                                                                                            ttcGCCCAGGCGCTGGTAACCCATCCTTCCACCATGGCGGTAATCGTGCCTTTCACCCAAAGGGCGGCGTAGATATGCACCATGATAACCACAATTAACGCCACTGCGGCAAATGAATGCAGCATTAACGCGAATCGGATCACCGGGAt      >  1:120704/1‑149 (MQ=255)
                                                                                                                                              cGCCCAGGCGCTGGTAACCCATCCTTCCACCATGGCGGTAATCGTGCCTTTCACCCAAAGGGCGGCGTAGATATGCACCATGATAACCACAATTAACGCCACTGCGGCAAATGAATGCAGCATTAACGCGAATCGGATCACCGGGAtt     <  5:127533/148‑1 (MQ=255)
                                                                                                                                                 ccAGGCGCTGGTAACCCATCCTTCCACCATGGCGGTAATCGTGCCTTTCACCCAAAGGGCGGCGTAGATATGCACCATGATAACCACAATTAACGCCACTGCGGCAAATGAATGCAGCATTAACGCGAATCGGATCACCGGGATTgag  >  6:170574/1‑148 (MQ=255)
                                                                                                                                                 ccAGGCGCTGGTAACCCATCCTTCCACCATGGCGGTAATCGTGCCTTTCACCCAAAGGGCGGCGTAGATATGCACCATGATAACCACAATTAACGCCACTGCGGCAAATGAATGCAGCATTAACGCGAATCGGATCACCGGGATTgag  >  2:70223/1‑148 (MQ=255)
                                                                                                                                                    |                                                                                                                                                
GGGAACAATAACGGCGGAATCATATCCGCCGTACGTTTCTCGCTTGAACCCAGCTCATCTTGCGGGATTATGCGAATACTCATTCAGCTTTCTTTTCCGTTGTCTTGCGGACCTCACGGTACCAGCGCGGGTGATGTTTCTTCGCCCAGGCGCTGGTAACCCATCCTTCCACCATGGCGGTAATCGTGCCTTTCACCCAAAGGGCGGCGTAGATATGCACCATGATAACCACAATTAACGCCACTGCGGCAAATGAATGCAGCATTAACGCGAATCGGATCACCGGGATTGAG  >  NC_000913/4081146‑4081438

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.