New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913 1958447 =42 (0.650)3 (0.050) 3/242 NT 7.2% intergenic (‑315/+73) torY/cutC TMAO reductase III (TorYZ), cytochrome c‑type subunit/copper homeostasis protein
?NC_000913 1958474 = 38 (0.630)intergenic (‑342/+46) torY/cutC TMAO reductase III (TorYZ), cytochrome c‑type subunit/copper homeostasis protein

ATTTAAATCCCATAATTAATGTGAATATATACAATATTTTTAGCACAACGCTATGCGCACTGTCACCCGCTCAACAGATCACGGTTATCGTTCGTTTTTTATACTGTTCAGGGATAAAAAAAGGCCCCTGTTGAAATTGCAGGGGCCTGGTACGAGCAAGCATCATATTGGGCGACATGATGCAACGGTAAAAATCATTTGGCCTGATGGCGTTCAATGATTCCTTTCATTTCAGCAACCGCCGCCCCGTCTACGATATAACGCGAATACTCGTCCGCGTGTTCATCTGAT  >  NC_000913/1958326‑1958616
                                                                                                                                                    |                                                                                                                                              
aTTTAAATCCCATAATTAATGTGAATATATACAATATTTTTAGCACAACGCTATGCGCACTGTCACCCGCTCAACAGATCACGGTTATCGTTCGTTTTTTATACTGTTCAGGGATAAAAAAAGGCCCCTGTTGAAATTGCAGGGGCCTg                                                                                                                                                >  1:18548/1‑149 (MQ=255)
aTTTAAATCCCATAATTAATGTGAATATATACAATATTTTTAGCACAACGCTATGCGCACTGTCACCCGCTCAACAGATCACGGTTATCGTTCGTTTTTTATACTGTTCAGGGATAAAAAAAGGCCCCTGTTGAAATTGCAGGGGCCTg                                                                                                                                                >  3:38549/1‑149 (MQ=255)
 tttAAATCCCATAATTAATGTGAATATATACAATATTTTTAGCACAACGCTATGCGCACTGTCACCCGCTCAACAGATCACGGTTATCGTTCGTTTTTTATACTGTTCAGGGATAAAAAAAGGCCCCTGTTGAAATTGCAGGGGCCTgg                                                                                                                                               >  3:348003/1‑149 (MQ=255)
      aTCCCATAATTAATGTGAATATATACAATATTTTTAGCACAACGCTATGCGCACTGTCACCCGCTCAACAGATCACGGTTATCGTTCGTTTTTTATACTGTTCAGGGATAAAAAAAGGCCCCTGTTGAAATTGCAGGGGCCTGGTACGa                                                                                                                                          <  2:401978/149‑1 (MQ=255)
           aTAATTAATGTGAATATATACAATATTTTTAGCACAACGCTATGCGCACTGTCACCCGCTCAACAGATCACGGTTATCGTTCGTTTTTTATACTGTTCAGGGATAAAAAAAGGCCCCTGTTGAAATTGCAGGGGCCTGGTACGagcaag                                                                                                                                     >  4:515784/1‑149 (MQ=255)
               ttaatGTGAATATATACAATATTTTTAGCACAACGCTATGCGCACTGTCACCCGCTCAACAGATCACGGTTATCGTTCGTTTTTTATACTGTTCAGGGATAAAAAAAGCCCCCTGTTGAAATTGCAGGGGCCTGGTACGAGCAAGcatc                                                                                                                                 <  4:340109/149‑1 (MQ=255)
                        atatatACAATATTTTTAGCACAACGCTATGCGCACTGTCACCCGCTCAACAGATCACGGTTATCGTTCGTTTTTTATACTGTTCAGGGATAAAAAAAGGCCCCTGTTGAAATTGCAGGGGCCTGGTACGAGCAAGCATCATATTGGGc                                                                                                                        <  1:413794/149‑1 (MQ=255)
                                aaTATTTTTAGCACAACGCTATGCGCACTGTCACCCGCTCAACAGATCACGGTTATCGTTCGTTTTTTATACTGTTCAGGGATAAAAAAAGGCCCCTGTTGAAATTGCAGGGGCCTGGTACGAGCAAGCATCATATTGGGCGACatgat                                                                                                                >  4:551796/1‑149 (MQ=255)
                                aaTATTTTTAGCACAACGCTATGCGCACTGTCACCCGCTCAACAGATCACGGTTATCGTTCGTTTTTTATACTGTTCAGGGATAAAAAAAGGCCCCTGTTGAAATTGCAGGGGCCTGGTACGAGCAAGCATCATATTGGGCGACatgat                                                                                                                >  4:311266/1‑149 (MQ=255)
                                 atatTTTTAGCACAACGCTATGCGCACTGTCACCCGCTCAACAGATCACGGTTATCGTTCGTTTTTTATACTGTTCAGGGATAAAAAAAGGCCCCTGTTGAAATTGCAGGGGCCTGGTACGAGCAAGCATCATATTGGGCGACatgatg                                                                                                               >  3:447050/1‑149 (MQ=255)
                                 atatTTTTAGCACAACGCTATGCGCACTGTCACCCGCTCAACAGATCACGGTTATCGTTCGTTTTTTATACTGTTCAGGGATAAAAAAAGGCCCCTGTTGAAATTGCAGGGGCCTGGTACGAGCAAGCATCATATTGGGCGACatgatg                                                                                                               <  4:121319/149‑1 (MQ=255)
                                  tatTTTTAGCACAACGCTATGCGCACTGTCACCCGCTCAACAGATCACGGTTATCGTTCGTTTTTTATACTGTTCAGGGATAAAAAAAGGCCCCTGTTGAAATTGCAGGGGCCTGGTACGAGCAAGCATCATATTGGGCGACATGATGc                                                                                                              <  3:515784/149‑1 (MQ=255)
                                   atTTTTAGCACAACGCTATGCGCACTGTCACCCGCTCAACAGATCACGGTTATCGTTCGTTTTTTATACTGTTCAGGGATAAAAAAAGGCCCCTGTTGAAATTGCAGGGGCCTGGTACGAGCAAGCATCATATTGGGCGACATGATGCa                                                                                                             >  1:43670/1‑149 (MQ=255)
                                              aaCGCTATGCGCACTGTCACCCGCTCAACAGATCACGGTTATCGTTCGTTTTTTATACTGTTCAGGGATAAAAAAAGGCCCCTGTTGAAATTGCAGGGGCCTGGTACGAGCAAGCATCATATTGGGCGACATGATGCAACGGTAAAAAt                                                                                                  >  2:282092/1‑149 (MQ=255)
                                                      gcgcACTGTCACCCGCTCAACAGATCACGGTTATCGTTCGTTTTTTATACTGTTCAGGGATAAAAAAAGGCCCCTGTTGAAATTGCAGGGGCCTGGTACGAGCAAGCATCATATTGGGCGACATGATGCAACGGTAAAAATCATTTGGc                                                                                          <  4:38549/149‑1 (MQ=255)
                                                      gcgcACTGTCACCCGCTCAACAGATCACGGTTATCGTTCGTTTTTTATACTGTTCAGGGATAAAAAAAGGCCCCTGTTGAAATTGCAGGGGCCTGGTACGAGCAAGCATCATATTGGGCGACATGATGCAACGGTAAAAATCATTTGGc                                                                                          <  2:18548/149‑1 (MQ=255)
                                                      gcgcACTGTCACCCGCTCAACAGATCACGGTTATCGTTCGTTTTTTATACTGTTCAGGGATAAAAAAAGGCCCCTGTTGAAATTGCAGGGGCCTGGTACGAGCAAGCATCATATTGGGCGACATGATGCAACGGTAAAAATCATTTGGc                                                                                          <  4:100337/149‑1 (MQ=255)
                                                           cTGTCACCCGCTCAACAGATCACGGTTATCGTTCGTTTTTTATACTGTTCAGGGATAAAAAAAGGCCCCTGTTGAAATTGCAGGGGCCTGGTACGAGCAAGCATCATATTGGGCGACATGATGCAACGGTAAAAATCATTTGGCCTGAt                                                                                     <  2:274878/149‑1 (MQ=255)
                                                                aCCCGCTCAACAGATCACGGTTATCGTTCGTTTTTTATACTGTTCAGGGATAAAAAAAGGCCCCTGTTGAAATTGCAGGGGCCTGGTACGAGCAAGCATCATATTGGGCGACATGATGCAACGGTAAAAATCATTTGGCCTGATGGCGt                                                                                >  3:229547/1‑149 (MQ=255)
                                                                  ccGCTCAACAGATCACGGTTATCGTTCGTTTTTTATACTGTTCAGGGATAAAAAAAGGCCCCTGTTGAAATTGCAGGGGCCTGGTACGAGCAAGCATCATATTGGGCGACATGATGCAACGGTAAAAATCATTTGGCCTGATGGCGTTc                                                                              >  4:102088/1‑149 (MQ=255)
                                                                  ccGCTCAACAGATCACGGTTATCGTTCGTTTTTTATACTGTTCAGGGATAAAAAAAGGCCCCTGTTGAAATTGCAGGGGCCTGGTACGAGCAAGCATCATATTGGGCGACATGATGCAACGGTAAAAATCATTTGGCCTGATGGCGTTc                                                                              >  1:282802/1‑149 (MQ=255)
                                                                  ccGCTCAACAGATCACGGTTATCGTTCGTTTTTTATACTGTTCAGGGATAAAAAAAGGCCCCTGTTGAAATTGCAGGGGCCTGGTACGAGCAAGCATCATATTGGGCGACATGATGCAACGGTAAAAATCATTTGGCCTGATGGCGTTc                                                                              >  1:137712/1‑149 (MQ=255)
                                                                    gCTCAACAGATCACGGTTATCGTTCGTTTTTTATACTGTTCAGGGATAAAAAAAGGCCCCTGTTGAAATTGCAGGGGCCTGGTACGAGCAAGCATCATATTGGGCGACATGATGCAACGGTAAAAATCATTTGGCCTGATGGCGTTCaa                                                                            <  1:315250/149‑1 (MQ=255)
                                                                           aGATCACGGTTATCGTTCGTTTTTTATACTGTTCAGGGATAAAAAAAGGCCCCTGTTGAAATTGCAGGGGCCTGGTACGAGCAAGCATCATATTGGGCGACATGATGCAACGGTAAAAATCATTTGGCCTGATGGCGTTCAATGATTcc                                                                     <  3:99976/149‑1 (MQ=255)
                                                                            gATCACGGTTATCGTTCGTTTTTTATACTGTTCAGGGATAAAAAAAGGCCCCTGTTGAAATTGCAGGGGCCTGGTACGAGCAAGCATCATATTGGGCGACATGATGCAACGGTAAAAATCATTTGGCCTGATGGCGTTCAATGATTCCt                                                                    >  2:179323/1‑149 (MQ=255)
                                                                               cACGGTTATCGTTCGTTTTTTATACTGTTCAGGGATAAAAAAAGGCCCCTGTTGAAATTGCAGGGGCCTGGTACGAGCAAGCATCATATTGGGCGACATGATGCAACGGTAAAAATCATTTGGCCTGATGGCGTTCAATGATTCCtttc                                                                 <  3:2163/149‑1 (MQ=255)
                                                                                 cGGTTATCGTTCGTTTTTTATACTGTTCAGGGATAAAAAAAGGCCCCTGTTGAAATTGCAGGGGCCTGGTACGAGCAAGCATCATATTGGGCGACATGATGCAACGGTAAAAATCATTTGGCCTGATGGCGTTCAATGATTCCtttcat                                                               <  3:311266/149‑1 (MQ=255)
                                                                                   gTTATCGTTCGTTTTTTATACTGTTCAGGGATAAAAAAAGGCCCCTGTTGAAATTGCAGGGGCCTGGTACGAGCAAGCATCATATTGGGCGACATGATGCAACGGTAAAAATCATTTGGCCTGATGGCGTTCAATGATTCCtttcattt                                                             >  3:266584/1‑149 (MQ=255)
                                                                                     tATCGTTCGTTTTTTATACTGTTCAGGGATAAAAAAAGGCCCCTGTTGAAATTGCAGGGGCCTGGTACGAGCAAGCATCATATTGGGCGACATGATGCAACGGTAAAAATCATTTGGCCTGATGGCGTTCAATGATTCCtttcatttca                                                           <  1:27002/149‑1 (MQ=255)
                                                                                      aTCGTTCGTTTTTTATACTGTTCAGGGATAAAAAAAGGCCCCTGTTGAAATTGCAGGGGCCTGGTACGAGCAAGCATCATATTGGGCGACATGATGCAACGGTAAAAATCATTTGGCCTGATGGCGTTCAATGATTCCTTTCATTTCAg                                                          <  1:521002/149‑1 (MQ=255)
                                                                                         gttcgttTTTTATACTGTTCAGGGATAAAAAAAGGCCCCTGTTGAAATTGCAGGGGCCTGGTACGAGCAAGCATCATATTGGGCGACATGATGCAACGGTAAAAATCATTTGGCCTGATGGCGTTCAATGATTCCTTTCATTTCAGCaa                                                       >  1:565452/1‑149 (MQ=255)
                                                                                          ttcgttTTTTATACTGTTCAGGGATAAAAAAAGGCCCCTGTTGAAATTGCAGGGGCCTGGTACGAGCAAGCATCATATTGGGCGACATGATGCAACGGTAAAAATCATTTGGCCTGATGGCGTTCAATGATTCCTTTCATTTCAGCAAc                                                      >  4:116869/1‑149 (MQ=255)
                                                                                                   tataCTGTTCAGGGATAAAAAAAGGCCCCTGTTGAAATTGCAGGGGCCTGGTACGAGCAAGCATCATATTGGGCGACATGATGCAACGGTAAAAATCATTTGGCCTGATGGCGTTCAATGATTCCtttc                                                                 <  4:222690/129‑1 (MQ=255)
                                                                                                   tataCTGTTCAGGGATAAAAAAAGGCCCCTGTTGAAATTGCAGGGGCCTGGTACGAGCAAGCATCATATTGGGCGACATGATGCAACGGTAAAAATCATTTGGCCTGATGGCGTTCAATGATTCCtttc                                                                 >  3:222690/1‑129 (MQ=255)
                                                                                                           tCAGGGATAAAAAAAGGCCCCTGTTGAAATTGCAGGGGCCTGGTACGAGCAAGCATCATATTGGGCGACATGATGCAACGGTAAAAATCATTTGGCCTGATGGCGTTCAATGATTCCTTTCATTTCAGCAACCGCCGCCCCGTCTACGa                                     <  2:137712/149‑1 (MQ=255)
                                                                                                           tCAGGGATAAAAAAAGGCCCCTGTTGAAATTGCAGGGGCCTGGTACGAGCAAGCATCATATTGGGCGACATGATGCAACGGTAAAAATCATTTGGCCTGATGGCGTTCAATGATTCCTTTCATTTCAGCAACCGCCGCCCCGTCTACGa                                     <  2:282802/149‑1 (MQ=255)
                                                                                                               ggATAAAAAAAGGCCCCTGTTGAAATTTCAGGGGCCTGGTACGAGCAAGCATCATATTGGGCGACATGATGCAACGGTAAAAATCATTTGGCCTGATGGCGTTCAATGATTCCTTTCATTTCAGCAACCGCCGCCCCGTCTACGAtata                                 >  2:254122/1‑149 (MQ=255)
                                                                                                                     aaaaaGGCCCCTGTTGAAATTGCAGGGGCCTGGTACGAGCAAGCATCATATTGGGCGACATGATGCAACGGTAAAAATCATTTGGCCTGATGGCGTTCAATGAtt                                                                       <  4:159640/105‑1 (MQ=255)
                                                                                                                     aaaaaGGCCCCTGTTGAAATTGCAGGGGCCTGGTACGAGCAAGCATCATATTGGGCGACATGATGCAACGGTAAAAATCATTTGGCCTGATGGCGTTCAATGAtt                                                                       >  3:159640/1‑105 (MQ=255)
                                                                                                                          GGCCCCTGTTGAAATTGCAGGGGCCTGGTACGAGCAAGCATCATATTGGGCGACATGATGCAACGGTAAAAATCATTTGGCCTGATGGCGTTCAATGATTCCTTTCATTTCAGCAACCGCCGCCCCGTCTACGATATAACGCGAATACt                      <  3:417773‑M1/149‑1 (MQ=255)
                                                                                                                             ccctgcaatttcaacaggggcctGGTACGAGCAAGCATCATATTGGGCGACATGATGCAACGGTaaa                                                                                                     <  2:469302‑M2/44‑1 (MQ=255)
                                                                                                                             ccctgcaatttcaacaggggcctGGTACGAGCAAGCATCATATTGGGCGACATGATGCAACGGTaaa                                                                                                     >  1:469302‑M2/24‑67 (MQ=255)
                                                                                                                             CCCTGTTGAGATTGCAGGGGCCTGGTACGAGCAAGCATCATATTGGGCGACATGATGCAACGGTAAAAATCATTTGGCCTGATGGCGTTCAATGATTCCTTTCATTTCAGCAACCGCCGCCCCGTCGACGATATAACGCGAATACTCGt                   <  1:512185‑M1/149‑1 (MQ=255)
                                                                                                                              cctgcaatttcaacaggggcctGGGACGAGCAAGCATCATATTGGGCGACATGATGCAACGGTAAAAATCATTTGGCCTGATGGCGTTCAATGATTCCTTTCATTTCAGCAACCGCCGCCCCGTCTACGATATAACGCGAATACTCGTc                  >  2:335443‑M2/23‑149 (MQ=255)
                                                                                                                               CTGTTGAAATTGCAGGGGCCTGGTACGAGCAAGCATCATATTGGGCGACATGATGCAACGGTAAAAATCATTTGGCCTGATGGCGTTCAATGATTCCTTTCATTTCAGCAACCGCCGCCCCGTCTACGATATAACGCGAATACTCGTcc                 >  4:229547‑M1/1‑149 (MQ=255)
                                                                                                                                 GTTGAAATTGCAGGGGCCTGGTACGAGCAAGCATCATATTGGGCGACATGATGCAACGGTAAAAATCATTTGGCCTGATGGCGTTCAATGATTCCTTTCATTTCAGCAACCGCCGCCCCGTCTACGATATAACGCGAATACTCGTCcgc               <  4:108659‑M1/149‑1 (MQ=255)
                                                                                                                                 GTTGAAATTGCAGGGGCCTGGTACGAGCAAGCATCATATTGGGCGACATGATGCAACGGTAAAAATCATTTGGCCTGATGGCGTTCAATGATTCCTTTCATTTCAGCAACCGCCGCCCCGTCTACGATATAACGCGAATACTCGTCcgc               <  1:422924‑M1/149‑1 (MQ=255)
                                                                                                                                  ttGAAATTGCAGGGGCCTGGTACGAGCAAGCATCATATTGGGCGACATGATGCAACGGTAAAAATCATTTGGCCTGATGGCGTTCAATGATTCCtt                                                                   >  2:478118/1‑96 (MQ=255)
                                                                                                                                  ttGAAATTGCAGGGGCCTGGTACGAGCAAGCATCATATTGGGCGACATGATGCAACGGTAAAAATCATTTGGCCTGATGGCGTTCAATGATTCCtt                                                                   <  1:478118/96‑1 (MQ=255)
                                                                                                                                              gggcctGGTACGAGCAAGCATCATATTGGGCGACATGATGCAACGGTAAAAATCATTTGGCCTGATGGCGTTCAATGATTCCTTTCATTTCAGCAACCGCCGCCCCGTCTACGATATAACGCGAATACTCGTCCGCGTGTTCATCtgat  <  1:254122‑M2/143‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                              
ATTTAAATCCCATAATTAATGTGAATATATACAATATTTTTAGCACAACGCTATGCGCACTGTCACCCGCTCAACAGATCACGGTTATCGTTCGTTTTTTATACTGTTCAGGGATAAAAAAAGGCCCCTGTTGAAATTGCAGGGGCCTGGTACGAGCAAGCATCATATTGGGCGACATGATGCAACGGTAAAAATCATTTGGCCTGATGGCGTTCAATGATTCCTTTCATTTCAGCAACCGCCGCCCCGTCTACGATATAACGCGAATACTCGTCCGCGTGTTCATCTGAT  >  NC_000913/1958326‑1958616

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