New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913_3_bme_pgi 1686709 =40 (0.870)4 (0.090) 4/272 6.3 9.0% intergenic (‑121/+22) fumC/fumA fumarate hydratase (fumarase C),aerobic Class II/fumarate hydratase (fumarase A), aerobic Class I
?NC_000913_3_bme_pgi 1686725 = 43 (0.980)intergenic (‑137/+6) fumC/fumA fumarate hydratase (fumarase C),aerobic Class II/fumarate hydratase (fumarase A), aerobic Class I
Rejected: Coverage evenness skew score above cutoff.
Rejected: Frequency below/above cutoff threshold.

TGTATTCATGACCTGCTCCTCACCTGATTAATTTTTTCTTTCTGTTTTGCTTTCGTTAAGCAACTTTTAGCTCACTTATTATTTACCATTTGATAACAAATGTTTGGTCTTTCGTGCCATGTAAAAAAACCGCCCCGAAGGGCGGCTCTGTTTATTTCACACAGCGGGTGCATTGTGTGAGTTGTATCTGCTGGAAGAAGTCATTTCCTTTATCATCCACAAGGATAAACGCCGGGAAATCTTCCACTTCAATTTTCCAGATGGCTTCCATTCCCAGTTCCGGATATTCAA  >  NC_000913_3_bme_pgi/1686580‑1686870
                                                                                                                                                 |                                                                                                                                                 
tGTATTCATGACCTGCTCCTCACCTGATTAATTTTTTCTTTCTGTTTTGCTTTCGTTAAGCAACTTTTAGCTCACTTATTATTTACCATTTGATAACAAATGTTTGGTCTTTCGTGCCATGTAAAAAAACCGCCCCGAAGGGCGGctct                                                                                                                                                >  2:267163/1‑149 (MQ=255)
     tCATGACCTGCTCCTCACCTGATTAATTTTTTCTTTCTGTTTTGCTTTCGTTAAGCAACTTTTAGCTCACTTATTATTTACCATTTGATAACAAATGTTTGGTCTTTCGTGCCATGTAAAAAAACCGCCCCGAAGGGCGGCTCTGTTTa                                                                                                                                           <  1:366816/149‑1 (MQ=255)
      cATGACCTGCTCCTCACCTGATTAATTTTTTCTTTCTGTTTTGCTTTCGTTAAGCAACTTTTAGCTCACTTATTATTTACCATTTGATAACAAATGTTTGGTCTTTCGTGCCATGTAAAAAAACGGCCCCGAAGGGCGGCTCTGTTTAt                                                                                                                                          <  1:701910/149‑1 (MQ=255)
       aTGACCTGCTCCTCACCTGATTAATTTTTTCTTTCTGTTTTGCTTTCGTTAAGCAACTTTTAGCTCACTTATTATTTACCATTTGATAACAAATGTTTGGTCTTTCGTGCCATGTAAAAAAACCGCCCCGAAGGGCGGCTCTGTTTAtt                                                                                                                                         >  2:101515/1‑149 (MQ=255)
            cTGCTCCTCACCTGATTAATTTTTTCTTTCTGTTTTGCTTTCGTTAAGCAACTTTTAGCTCACTTATTATTTACCATTTGATAACAAATGTTTGGTCTTTCGTGCCATGTAACAAAACCGCCCCGAAGGGCGGCTCTGTTTATTTccca                                                                                                                                    >  1:515653/1‑149 (MQ=255)
            cTGCTCCTCACCTGATTAATTTTTTCTTTCTGTTTTGCTTTCGTTAAGCAACTTTTAGCTCACTTATTATTTACCATTTGATAACAAATGTTTGGTCTTTCGTGCCATGTAAAAAAACCGCCCCGAAGGGCGGCTCTGTTTATTTcaca                                                                                                                                    >  2:251628/1‑149 (MQ=255)
            cTGCTCCTCACCTGATTAATTTTTTCTTTCTGTTTTGCTTTCGTTAAGCAACTTTTAGCTCACTTATTATTTACCATTTGATAACAAATGTTTGGTCTTTCGTGCCATGTAAAAAAACCGCCCCGAAGGGCGGCTCTGTTTATTTcaca                                                                                                                                    <  1:251628/149‑1 (MQ=255)
                   tcACCTGATTAATTTTTTCTTTCTGTTTTGCTTTCGTTAAGCAACTGTTAGCTCTTTTATTATTTACCATTTGATAACAAATGTTTGGTCTTTCGTTCGATCTAAAAAACCCGCCCCGAAGGGCGGCTCTGTCTAATTCCAACAGCggg                                                                                                                             >  2:87950/1‑149 (MQ=255)
                         gATTAATTTTTTCTTTCTGTTTTGCTTTCGTTAAGCAACTTTTAGCTCACTTATTATTTACCATTTGATAACAAATGTTTGGTCTTTCGTGCCATGTAAAAAAACCGCCCCGAAGGGCGGCTCTGTTTATTTCACACAGCGGGTGCAtt                                                                                                                       <  1:647635/149‑1 (MQ=255)
                               ttttttCTTTCTGTTTTGCTTTCGTTAAGCAACTTTTAGCTCACTTATTATTTACCATTTGATAACAAATGTTTGGTCTTTCCTGCTTGTTATAAAAACCGCCCCGAAGGGCGGCTCTGTTTATTTCACACAGCGGGTGCATTGTGTGa                                                                                                                 >  1:594421/1‑149 (MQ=255)
                                          tGTTTTGCTTTCGTTAAGCAACTTTTAGCTCACTTATTATTTACCATTTGATAACAAATGTTTGGTCTTTCGTGCCATGTAAAAAAACCGCCCCGAAGGGGGGCTCTGTTTCTTTCACACAGCGGGTGCATGGGGTGAGGTGTAGctgc                                                                                                      >  1:248010/1‑149 (MQ=255)
                                            ttttGCTTTCGTTAAGCAACTTTTAGCTCACTTATTATTTACCATTTGATAACAAATGTTTGGTCTTTCGTGCCATGTAAAAAAACCGCCCCGAAGGGCGGCTCTGTTTATTTCACACAGCGGGTGCATTGTGTGAGTTGTATctgctg                                                                                                    >  1:637172/1‑149 (MQ=255)
                                            ttttGCTTTCGTTAAGCAACTTTTAGCTCACTTATTATTTACCATTTGATAACAAATGTTTGGTATTTCGTGCCATGTAAAAAAACCGCCCCGAAGGGCGGCTATGTTTATTTCACACAGCGGGTGCATTGGTTGAGTTGTAGctgctg                                                                                                    >  1:695224/1‑149 (MQ=255)
                                            ttttGCTTGTGTTAAGCAACTTTTAGCTCAGTTCTTATTTACCATTTGATAACAAATGTTTGGTCTTTCGTGCCATGTACCAAAACCGCCCCGAAGGGCGGCTCTGTTTATTTCACACAGCGGGTGCATTGTGTGAGTTGTATctgctg                                                                                                    <  2:220959/149‑1 (MQ=255)
                                                  tttCGTTAAGCAACTTTTAGCTCACTTATTATTTACCATTTGATAACAAATGTTTGGTCTTTCGTGCCATTTAATAAAACCGCCCCGAAGGGCGGCTCTGTTTATTTCACACAGCGGGTGCATTGTGTGAGTTGTATCTGCTGgaagaa                                                                                              >  2:725710/1‑149 (MQ=255)
                                                     cGTTAAGCAACTTTTAGCTCACTTATTATTTACCATTTGATAACAAATGTTTGGTCTTTCGTGCCATGTAAAAAAACCGCCCCGAAGGGCGGCTCTGTTTATTTCACACAGCGGGTGCATTGTGTGAGTTGTATCTGCTGGAAGAAGTc                                                                                           <  2:111052/149‑1 (MQ=255)
                                                          aGCAACTTTTAGCTCACTTATTATTTACCATTTGATAACAAATGTTTGGTCTTTCGTGCCATGTAAAAAAACCGCCCCGAAGGGCGGCTCTGTTTATTTCACACAGCGGGTGCATTGTGTGAGTTGTATCTGCTGGAAGAAGTCATTTc                                                                                      <  1:101515/149‑1 (MQ=255)
                                                             aaCTTTTAGCTCACTTATTAATTACCATTTGATAACAAATGTTTGGTCTTTCGTGCCATGTAAAAAAACCGCCCCGAAGGGCGGCTCTGTTTATTTCACACAGCGGGTGCATTGTGTGAGTTGTATCTGCTGGAAGAAGTCATTTCCtt                                                                                   <  2:586279/149‑1 (MQ=255)
                                                                  ttAGCTCACTTATTATTTACCATTTGATAACAAATGTTTGGTCTTTCGTGCCATGTAAAAAAACCGCCCCGAAGGGCGGCTCTGTTTATTTCACACAGCGGGTTCATTGTGTGAGTTGTATCTGCTGGAAGAAGTCATTTCCTTTatca                                                                              >  1:169445/1‑149 (MQ=255)
                                                                        cACTTATTATTTACCATTTGATAACAAATGTTTGGTCTTTCGTGCCATGTAAAAAAACCGCCCCGAAGGGCGGCTCTGTTTATTTCACACAGCGGGTGCATTGTGTGAGtt                                                                                                              >  2:196026/1‑111 (MQ=255)
                                                                        cACTTATTATTTACCATTTGATAACAAATGTTTGGTCTTTCGTGCCATGTAAAAAAACCGCCCCGAAGGGCGGCTCTGTTTATTTCACACAGCGGGTGCATTGTGTGAGtt                                                                                                              <  1:196026/111‑1 (MQ=255)
                                                                           ttattaTTTACCATTTGATAACAAATGTTTGGTCTTTCGTGCCATGTAAAAAAACCGCCCCGAAGGGCGGCTCTGTTTATTTCACACAGCGGGTGCATTGTGTGAGTTGTATCTGCTGGAAGAAGTCATTTCCTTTATCATCCACAAgg                                                                     >  1:416194/1‑149 (MQ=255)
                                                                               tattTACCATTTGATAACAAATGTTTGGTCTTTCGTGCCATGTAAAAAAACCGCCCCGAAGGGCGGCTCTGTTTATTTCACACAGCGGGTGCATTGTGTGAGTTGTATCTGCTGGAAGAAGTCATTTCCTTTATCATCCACAAGGATaa                                                                 >  1:18047/1‑149 (MQ=255)
                                                                               tattTACCATTTGATAACAAATGTTTGGTCTTTCGTGCCATGTAAAAAAACCGCCCCGAAGGGCGGCTCTGTTTATTTCACACAGCGGGTGCATTGTGTGAGTTGTATCTGCTGGAAGAAGTCATTTCCTTTATCATCCACAAGGATaa                                                                 >  2:324305/1‑149 (MQ=255)
                                                                                     ccATTTGATAACAAATGTTTGGTCTTTCGTGCCATGTAAAAAAACCGCCCCGAAGGGCGGCTCTGTTTATTTCACACAGCGTGTGCATTGTGTGAGTTGTATCTGCTGGAAGAAGTCATTTCCTTTATCATCCACAAGGATAAACGCCg                                                           <  1:267163/149‑1 (MQ=255)
                                                                                                 aaaTGTTTGGTCTTTCGTGCCATGTAAAAAAACCGCCCCGAAGGGCGGCTCTGTTTATTTCACACAGCGGGTGCATTGTGTGAGTTGTATCTGCTGGAAGAAGTCATTTCCTTTATCATCCACAAGGATAAACGCCGGGAAATCTTCCa                                               <  2:637172/149‑1 (MQ=255)
                                                                                                  aaTGTTTGGTCTTTCGTGCCATGTAAAAAAACCGCCCCGAAGGGCGGCTCTGTTTATTTCACACAGCGGGTGCATTGTGTGAGTTGTATCTGCTGGAAGAAGTCATTTCCTTTATCATCCACAAGGATAAACGCCGGGAAATCTTCCAc                                              <  1:324305/149‑1 (MQ=255)
                                                                                                       ttGGTCTTTCGTGCCATGTAAAAAAACCGCCCCGAAGGGCGGCTCTGTTTATTTCACACAGCGGGTGCATTGTGTGATTGGTATctgct                                                                                                     <  2:300660/89‑1 (MQ=255)
                                                                                                       ttGGTCTTTCGTGCCATGTAAAAAAACCGCCCCGAAGGGCGGCTCTGTTTATTTCACACAGCGGGTGCATTGTGTGAGTTGTATctgct                                                                                                     >  1:300660/1‑89 (MQ=255)
                                                                                                         gctCTTTCGTGCCATGTAGAAAGACAGTCCCGAAGGGCGGCTTCGTTTATTTCACACGGCGGGTGCATTGTGTGAGTTGTATCTGCTGCAAGAACTCATTTCCTTTCTCATCCACAAGGATAAACGCCGCGAAATCTTCCACTTCAAtt                                       <  2:348085/147‑1 (MQ=255)
                                                                                                           tCTTTCGTGCCATGTAAAAAAACCGCCCCGAAGGGCGGCTCTGTTTATTTCACACAGCGGGTGCATTGTGTGAGTTGTATCTGCTGGAAGAAGTCATTTCCTTTATCATCCACAAGGATAAACGCCGGGAAATCTTCCACTTCAAtttt                                     >  1:173585/1‑149 (MQ=255)
                                                                                                              ttCGTGCCATGTAAAAAAACCGCCCCGAAGGGCGGCTCTGTTTATTTCACACAGCGGGTGCATTGTGTGAGTTGTATCTGCTGGAAGAAGTCATTTCCTTTATCATCCACAAGGATAAACGCCGGGAAATCTTCCACTTCAATTTTCCa                                  <  1:773101/149‑1 (MQ=255)
                                                                                                               tCGTGCCATGTAAAAAAACCGCCCCGAAGGGCGGCTCTGTTTATTTCACACAGCGGGTGCATTGTGTGAGTTGTATCTGCTGGAAGAAGTCATTTCCTTTATCATCCACAAGGATAAACGCCGGGAAATCTTCCACt                                             <  2:709024/137‑1 (MQ=255)
                                                                                                               tCGTGCCATGTAAAAAAACCGCCCCGAAGGGCGGCTCTGTTTATTTCACACAGCGGGTGCATTGTGTGAGTTGTATCTGCTGGAAGAAGTCATTTCCTTTATCATCCACAAGGATAAACGCCGGGAAATCTTCCACt                                             >  1:709024/1‑137 (MQ=255)
                                                                                                                 gtaccATGTAACAAAACGGCCCTGACGGGCGGCTCTGTTTATTTCGCACCGCGGGTGTATTGTGTGAGTTGTATCTTCTGGAAGAAGTCATTTCCTTTATCATCCACAAGGATAAACGCCGGGAGATCTTCCACTTCAATTTTCCAGAt                               <  2:484836/146‑1 (MQ=255)
                                                                                                                             aaaaCCGCCCCGAAGGGCGGCTCTGTTTATTTCACACAGCGGGTGCATTGTGTGAGTTGTATCTGCTGGAAGAAGTCATTTCCTTTATCATCCACAAGGATAAACGCCGGGAAAACTTCCCCTTCAATTTTCCAGATGGCTTCCATTcc                   >  2:263112/1‑149 (MQ=255)
                                                                                                                                   gCCCCGACGGCCGGCTCTGTTTATTTCACACAGCGGGTGCATTGTGTGAGTTGTATCTGCTGGAAGAAGTCATTTCCTTTATCATCCACAAGGATAAACGCCGGGAAATCTTCCACTCCAATTTTCCAGATGGCTTCCATACCCactgc             >  2:345847/1‑145 (MQ=37)
                                                                                                                                   gCCCCGAAGTGCGGGTCTGTTTATTTCACACGGCGGGTGCATTGTGTGAGTTGTATCTGCTGGAAGAAGTCATTGCCTTTATCATCCACAAGGATAAACGCCGGGAAATCTTCCACTTCAATTTTCCAGATGGCTTCCATTCCCAGTTc             <  2:787489/149‑1 (MQ=255)
                                                                                                                                     cccGAAGGGCGGCTCTGTTTATTTCACACAGCGGGTGCATTGTGTGAGTTGTATCTGCTGGAAGAAGTCATTTCCTTTATCATCCACAAGGATAAACGCCGGGAAATCTTCCACTTCAAtt                                       <  2:181513/121‑1 (MQ=255)
                                                                                                                                     cccGAAGGGCGGCTCTGTTTATTTCACACAGCGGGTGCATTGTGTGAGTTGTATCTGCTGGAAGAAGTCATTTCCTTTATCATCCACAAGGATAAACGCCGGGAAATCTTCCACTTCAAtt                                       >  1:181513/1‑121 (MQ=255)
                                                                                                                                     cccGAAGGGCGGCTCTGTTTATTTCACACAGCGGGTGCATTGTGTGAGTTGTATCTGCTGGAAGAAGTCATTTCCTTTATCATCCACAAGGATAAACGCCGGGAAATCTTCCACTTCAATTTTCCAGATGGCTTACATTACCactgcag           >  1:349837/1‑143 (MQ=255)
                                                                                                                                       cGAAGGGCGGCTCTGTTTATTTCACACAGCGGGTGCATTGTGTGAGTTGTATCTGCTGGCAGAAGTCATTTCCTTTATCATCCACAAGGATAAACGCCGGGAAATCTTCCCCTTCAATTTTCCGGCTGTCTTCCATTCCCAGTGCCGGa         >  2:337436/1‑149 (MQ=38)
                                                                                                                                       cGAAGGGCGGCTCTGTTTATTTCACACAGCGGGTGCATTGTGTGAGTTGTATCTGCTGGAAGAAGTCATTTCCTTTATCATCCACAAGGATAAACGACGGGAACTTTTCCACTTCAATTTTCCAGATGGCTTCCATTCCCAGTTCCGGa         >  1:247470/1‑149 (MQ=255)
                                                                                                                                          aggcctgcTCTGTTCATTCCACACAGCGGGTGCATTGTGTGAGTTGGATCTGCTGGAAGAAGGCATTTCCTTTATCATCCACAAGGATAAACGCCGGGAAATCTTCCACTTCAATTTTCCAGATGGCTTCCATTCCCAGTTCCGGatat      <  2:12108/143‑1 (MQ=37)
                                                                                                                                              cggcTCTGTTTATTTCACACTGCGGGTGCATTGTGTGAGTTGTATCTGCTGGAACAAGTCATTTCCTTTATCATCCACAAGGATAAACGCCGGGAAATCTTCCACTTCTATTCTCCAGATGGCTTCCATTCCCAGTTCCGGATATTCaa  <  2:615203/149‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                                 
TGTATTCATGACCTGCTCCTCACCTGATTAATTTTTTCTTTCTGTTTTGCTTTCGTTAAGCAACTTTTAGCTCACTTATTATTTACCATTTGATAACAAATGTTTGGTCTTTCGTGCCATGTAAAAAAACCGCCCCGAAGGGCGGCTCTGTTTATTTCACACAGCGGGTGCATTGTGTGAGTTGTATCTGCTGGAAGAAGTCATTTCCTTTATCATCCACAAGGATAAACGCCGGGAAATCTTCCACTTCAATTTTCCAGATGGCTTCCATTCCCAGTTCCGGATATTCAA  >  NC_000913_3_bme_pgi/1686580‑1686870

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