New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913 1972772 =42 (1.220)6 (0.190) 6/238 2.7 13.8% intergenic (‑81/+64) tar/cheW methyl‑accepting chemotaxis protein II/purine‑binding chemotaxis protein
?NC_000913 1972820 = 36 (1.120)intergenic (‑129/+16) tar/cheW methyl‑accepting chemotaxis protein II/purine‑binding chemotaxis protein
Rejected: Frequency below/above cutoff threshold.

TGTAACAGTGCGAATACCCCCAGCACCATTACCAACAGCGTGACTACGCGGATACGGTTAATCATAAGGCACCTTCCTGAAAACAAGTTGATCTCGTTATCGGCAAGGAGGGGGGAAACTTTATTGCTGATGCCACCCGCCGCGAAATTGAAATAAAAAACCCGATGCGCAGATCATCGGGTTCATTTCAATTGAGGAAATCGGGAGAATTACGCCACTTCTGACGCCGCGCTATCTAACAGCGCCATCTCTTCGCTGTTCAGCAGTTTTTCGATGTTCACCAGAATC  >  NC_000913/1972627‑1972914
                                                                                                                                                 |                                                                                                                                              
tGTAACAGTGCGAATACCCCCAGCACCATTACCAACAGCGTGACTACGCGGATACGGTTAATCATAAGGCACCTTCCTGAAAACAAGTTGATCTCGTTATCGGCAAGGAGGGGGGAAACTTTATTGCTGATGCCACCCGCCGCGAAAtt                                                                                                                                             >  1:132994/1‑149 (MQ=255)
ctttcaacagaatatctcactggactgggcgcactgggcgaccggatgttgattctggtgaacatcgaaaaactgctgaacagcgaagagatggcgctgttagatagcgcggcgtcagaagtggcgtaattctcccgatttcctcAAtt                                                                                                                                             >  1:74163‑M1/146‑149 (MQ=255)
   tcaacagaatatctcactggactgggcgcactgggcgaccggatgttgattctggtgaacatcgaaaaactgctgaacagcgaagagatggcgctgttagatagcgcggcgtcagaagtggcgtaattctcccgatttcctcAATTg                                                                                                                                            >  1:30687‑M1/143‑147 (MQ=255)
    caacagaatatctcactggactgggcgcactgggcgaccggatgttgattctggtgaacatcgaaaaactgctgaacagcgaagagatggcgctgttagatagcgcggcgtcagaagtggcgtaattctcccgatttcctcAATTGaa                                                                                                                                          >  7:19638‑M1/142‑148 (MQ=255)
     aacagaatatctcactggactgggcgcactgggcgaccggatgttgattctggtgaacatcgaaaaactgctgaacagcgaagagatggcgctgttagatagcgcggcgtcagaagtggcgtaattctcccgatttcctcAAtt                                                                                                                                             >  3:24262‑M1/141‑144 (MQ=255)
     aacagaatatctcactggactgggcgcactgggcgaccggatgttgattctggtgaacatcgaaaaactgctgaacagcgaagagatggcgctgttagatagcgcggcgtcagaagtggcgtaattctcccgatttcctcAAtt                                                                                                                                             <  4:24262‑M1/4‑1 (MQ=255)
           gAATACCCCCAGCACCATTACCAACAGCGTGACTACGCGGATACGGTTAATCATAAGGCACCTTCCTGAAAACAAGTTGATCTCGTTATCGGCAAGGAGGGGGGAAACTTTATTGCTGATGCCACCCGCCGCGAAATTGAAATa                                                                                                                                       <  6:154873/144‑1 (MQ=255)
           gAATACCCCCAGCACCATTACCAACAGCGTGACTACGCGGATACGGTTAATCATAAGGCACCTTCCTGAAAACAAGTTGATCTCGTTATCGGCAAGGAGGGGGGAAACTTTATTGCTGATGCCACCCGCCGCGAAATTGAAATa                                                                                                                                       >  5:154873/1‑144 (MQ=255)
           atatctcactggactgggcgcactgggcgaccggatgttgattctggtgaacatcgaaaaactgctgaacagcgaagagatggcgctgttagatagcgcggcgtcagaagtggcgtaattctcccgatttcctcAATTGAAAt                                                                                                                                        >  4:48120‑M1/135‑143 (MQ=255)
           atatctcactggactgggcgcactgggcgaccggatgtagattctggtgaacatcgaaaaactgctgaacagcgaagagatggcgctgttagatagcgcggcgtcagaagtggcgtaattctcccgatttcctcAATTGAAAt                                                                                                                                        <  3:48120‑M1/9‑1 (MQ=255)
            tatctcactggactgggcgcactgggcgaccggatgttgattctggtgaacatcgaaaaactgctgaacagcgaagagatggcgctgttagatagcgcggcgtcagaagtggcgtaattctcccgatttcctcAATTGAAAt                                                                                                                                        >  8:150826‑M1/134‑142 (MQ=255)
            tatctcactggactgggcgcactgggcgaccggatgttgattctggtgaacatcgaaaaactgctgaacagcgaagagatggcgctgttagatagcgcggcgtcagaagtggcgtaattctcccgatttcctcAATTGAAAt                                                                                                                                        <  7:150826‑M1/9‑1 (MQ=255)
              tctcactggactgggcgcactgggcgaccggatgttgattctggtgaacatcgaaaaactgctgaacagcgaagagatggcgctgttagatagcgcggcgtcagaagtggcgtaattctcccgatttcctcAATTGAAATAAAAAAccc                                                                                                                               <  5:3767‑M1/18‑1 (MQ=255)
                cccccAGCACCATTACCAACAGCGTGACTACGCGGATACGGTTAATCATAAGGCACCTTCCTGAAAACAAGTTGATCTCGTTATCGGCAAGGAGGGGGGAAACTTTATTGCTGATGCCACCCGCCGCGAAATTGAAATAAAAAACCCGa                                                                                                                             >  5:118953/1‑149 (MQ=255)
                   ccAGCACCATTACCAACAGCGTGACTACGCGGATACGGTTAATCATAAGGCACCTTCCTGAAAACAAGTTGATCTCGTTATCGGCAAGGAGGGGGGAAACTTTATTGCTGATGCCACCCGCCGCGAAATTGaa                                                                                                                                          >  4:122774/1‑133 (MQ=255)
                   ccAGCACCATTACCAACAGCGTGACTACGCGGATACGGTTAATCATAAGGCACCTTCCTGAAAACAAGTTGATCTCGTTATCGGCAAGGAGGGGGGAAACTTTATTGCTGATGCCACCCGCCGCGAAATTGaa                                                                                                                                          <  3:122774/133‑1 (MQ=255)
                      gCACCATTACCAACAGCGTGACTACGCGGATACGGTTAATCATAAGGCACCTTCCTGAAAACAAGTTGATCTCGTTATCGGCAAGGAGGGGGGAAACTTTATTGCTGATGCCACCCGCCGCGAAATTGAAATAAAAAACCCGATGCGCa                                                                                                                       <  2:132994/149‑1 (MQ=255)
                          cATTACCAACAGCGTGACTACGCGGATACGGTTAATCATAAGGCACCTTCCTGAAAACAAGTTGATCTCGTTATCGGCAAGGAGGGGGGAAACTTTATTGCTGATGCCACCCGCCGCGAAATTGAAATAAAAAACCCGATgcgc                                                                                                                        >  4:144908/1‑144 (MQ=255)
                          cATTACCAACAGCGTGACTACGCGGATACGGTTAATCATAAGGCACCTTCCTGAAAACAAGTTGATCTCGTTATCGGCAAGGAGGGGGGAAACTTTATTGCTGATGCCACCCGCCGCGAAATTGAAATAAAAAACCCGATgcgc                                                                                                                        <  3:144908/144‑1 (MQ=255)
                           aTTACCAACAGCGTGACTACGCGGATACGGTTAATCATAAGGCACCTTCCTGAAAACAAGTTGATCTCGTTATCGGCAAGGAGGGGGGAAACTTTATTGCTGATGCCACCCGCCGCGAAATTGAAATAAAAAACCCGCTGCGCAGatca                                                                                                                  >  8:1669/1‑149 (MQ=255)
                                              cgcgGATACGGTTAATCCTCAGGCACCTTCCTGACAACAAGTTGATCTCGTTATCGGCAAGGAGGGGGGAAACTTTATTGCTGATGCCACCCGCCGCGAAATTGAAATAAAAAACCCGATGCGCAGATCATCGGGTTCATTTCAATTGa                                                                                               <  3:22707/149‑1 (MQ=255)
                                                ttgattctggtgaacatcgaaaaactgctgaacagcgaagagatggcgctgttagatagcgcggcgtcagaagtggcgtaattctcccgatttcctcAATTGAAATAAAAAACCCGATGCGCAGATCATCGGGTTCATTTCAATTGAg                                                                                              >  6:3767‑M1/98‑148 (MQ=255)
                                                  gATACGGTTAATCATAAGGCACCTTCCTGAAAACAAGTTGATCTCGTTATCGGCAAGGAGGGGGGAAACTTTATTGCTGATGCCACCCGCCGCGAAATTGAAATAAAAAACCCGATGCGCAGATCATCGGGTTCATTTCAATTGAGGaa                                                                                           <  7:1669/149‑1 (MQ=255)
                                                       ggTTAATCATAAGGCACCTTCCTGAAAACAAGTTGATCTCGTTATCGGCAAGGAGGGGGGAAACTTTATTGCTGATGCCACCCGCCGCGAAATTGAAATAAAAAACCCGATGCGCAGATCATCGGGTTCATTTCaa                                                                                                   <  6:24013/136‑1 (MQ=255)
                                                       ggTTAATCATAAGGCACCTTCCTGAAAACAAGTTGATCTCGTTATCGGCAAGGAGGGGGGAAACTTTATTGCTGATGCCACCCGCCGCGAAATTGAAATAAAAAACCCGATGCGCAGATCATCGGGTTCATTTCaa                                                                                                   >  5:24013/1‑136 (MQ=255)
                                                              catcgaaaaactgctgaacagcgaagagatggcgctgttagatagcgcggcgtcagaagtggcgtaattctcccgatttcctcAATTGaaa                                                                                                                                         <  1:153292‑M1/8‑1 (MQ=255)
                                                              catcgaaaaactgctgaacagcgaagagatggcgctgttagatagcgcggcgtcagaagtggcgtaattctcccgatttcctcAATTGaaa                                                                                                                                         >  2:153292‑M1/84‑91 (MQ=255)
                                                                tAAGGCACCTTCCTGAAAACAAGTTGATCTCGTTATCGGCAAGGAGGGGGGAAACTTTATTGCTGATGCCACCCGCCGCGAAAtt                                                                                                                                             >  2:62600/1‑85 (MQ=255)
                                                                tAAGGCACCTTCCTGAAAACAAGTTGATCTCGTTATCGGCAAGGAGGGGGGAAACTTTATTGCTGATGCCACCCGCCGCGAAAtt                                                                                                                                             <  1:62600/85‑1 (MQ=255)
                                                                    aaaactgctgaacagcgaagagatggcgctgttagatagcgcggcgtcagaagtggcgtaattctcccgatttcctcAATTGAAATAAAAAACCCGATGCGCAGATCATCGGGTTCATTTCAATTGa                                                                                               >  1:134674‑M1/78‑127 (MQ=255)
                                                                    aaaactgctgaacagcgaagagatggcgctgttagatagcgcggcgtcagaagtggcgtaattctcccgatttcctcAATTGAAATAAAAAACCCGATGCGCAGATCATCGGGTTCATTTCAATTGa                                                                                               <  2:134674‑M1/50‑1 (MQ=255)
                                                                       ccTTCCTGAAAACAAGTTGATCTCGTTATCGGCAAGGAGGGGGGAAACTTTATTGCTGATGCCACCCGCCGCGAAATTGAAATAAAAAACCCGATGCGCAGATCATCGGGTTCATTTCAATTGAGGAAATCGGGAGAATTACGCCACt                                                                       >  1:70423/1‑148 (MQ=255)
                                                                                      gTTGATCTCGTTATCGGCAAGGAGGGGGGAAACTTTATTGCTGATGCCACCCGCCGCGAAATTGAAATAAAAAACCCGATGCGCAGATCATCGGGTTCATTTCAATTGAGGAAATCGGGAGAATTACGCCACTTCTGACGCCGCGCTa                                                        >  5:105055/1‑148 (MQ=255)
                                                                                      gTTGATCTCGTTATCGGCAAGGAGGGGGGAAACTTTATTGCTGATGCCACCCGCCGCGAAATTGAAATAAAAAACCCGATGCGCAGATCATCGGGTTCATTTCAATTGAGGAAATCGGGAGAATTACGCCACTTCTGACGCCGCGCTAt                                                       <  6:118953/149‑1 (MQ=255)
                                                                                       ttGATCTCGTTATCGGCAAGGAGGGGGGAAACTTTATTGCTGATGCCACCCGCCGCGAAATTGAAATAAAAAACCCGATGCGCAGATCATCGGGTTCATTTCAATTGAGGAAATCGGGAGAATTACGCCACtt                                                                      <  6:88827/133‑1 (MQ=255)
                                                                                       ttGATCTCGTTATCGGCAAGGAGGGGGGAAACTTTATTGCTGATGCCACCCGCCGCGAAATTGAAATAAAAAACCCGATGCGCAGATCATCGGGTTCATTTCAATTGAGGAAATCGGGAGAATTACGCCACtt                                                                      >  5:88827/1‑133 (MQ=255)
                                                                                         gATCTCGTTATCGGCAAGGAGGGGGGAAACTTTATTGCTGATGCCACCCGCCGCGAAATTGAAATAAAAAACCCGATGCGCAGATCATCGGGTTCATTTCAATTGAGGAAATCGGGAGAATTACGCCACTTCTGACGCCGCCCTATCt                                                     >  1:62575/1‑148 (MQ=255)
                                                                                            ctcGTTATCGGCAAGGAGGGGGGAAACTTTATTGCTGATGCCACCCGCCGCGAAATTGAAATAAAAAACCCGATGCGCAGATCATCggg                                                                                                             <  4:19050/89‑1 (MQ=255)
                                                                                            ctcGTTATCGGCAAGGAGGGGGGAAACTTTATTGCTGATGCCACCCGCCGCGAAATTGAAATAAAAAACCCGATGCGCAGATCATCggg                                                                                                             >  3:19050/1‑89 (MQ=255)
                                                                                                ttATCGGCAAGGAGGGGGGAAACTTTATTGCTGATGCCACCCGCCGCGAAATTGAAATAAAAAACCCGATGCGCAGATCATCGGGTTCATTTCAATTGAGGAAATCGGGAGAATTACGCCACTTCTGACGCCGCGCTATCTAACAgcgc                                             >  1:91803/1‑149 (MQ=255)
                                                                                                ttATCGGCAAGGAGGGGGGAAACTTTATTGCTGATGCCACCCGCCGCGAAATTGAAATAAAAAACCCGATGCGCAGATCATCGGGTTCATTTCAATTGAGGAAATCGGGAGAATTACGCCACTTCTGACGCCGCGCTATCTAACAgcgc                                             >  6:70992/1‑149 (MQ=255)
                                                                                                               ggggAAACTTTATTGCTGATGCCACCCGCCGCGAAATTGAAATAAAAAACCCGATGCGCAGATCATCGGGTTCATTTCAATTGAGGAAATCGGGAGAATTACGCCACTTCTGACGCCGCGCTATCTAACAGCGCCATCTCTTCGCTGtt                              <  6:105055/149‑1 (MQ=255)
                                                                                                                  tcagaagtggcgtaattctcccgatttcctcAATTGAAATAAAAAACCCGATGCGCAGATCATCGGGTTCAttt                                                                                                      <  3:105082‑M1/43‑1 (MQ=255)
                                                                                                                  tcagaagtggcgtaattctcccgatttcctcAATTGAAATAAAAAACCCGATGCGCAGATCATCGGGTTCAttt                                                                                                      >  4:105082‑M1/32‑74 (MQ=255)
                                                                                                                   aaaCTTTATTGCTGATGCCACCCGCCGCGAAATTGAAATAAAAAACCCGATGCGCAGATCATCGGGTTCATTTCAATTGAGGAAATCGGga                                                                                    >  6:107133/1‑91 (MQ=255)
                                                                                                                   aaaCTTTATTGCTGATGCCACCCGCCGCGAAATTGAAATAAAAAACCCGATGCGCAGATCATCGGGTTCATTTCAATTGAGGAAATCGGga                                                                                    <  5:107133/91‑1 (MQ=255)
                                                                                                                     aCTTTATTGCTGATGCCACCCGCCGCGAAATTGAAATAAAAAACCCGATGCGCAGATCATCGGGTTCATTTCAATTGAGGAAATCGGTAGAATTACGCCACTTCTGACGCCGCGCTATCTAACAGCGCCATctct                                      <  2:599/135‑1 (MQ=255)
                                                                                                                     aCTTTATTGCTGATGCCACCCGCCGCGAAATTGAAATAAAAAACCCGATGCGCAGATCATCGGGTTCATTTCAATTGAGGAAATCGGGAGAATTACGCCACTTCTGACGCCGCGCTATCTAACAGCGCCATctct                                      >  1:599/1‑135 (MQ=255)
                                                                                                                       tttATTGCTGATGCCACCCGCCGCGAAATTGAAATAAAAAACCCGATGCGCAGATCATCGGGTTCATTTCAATTGa                                                                                               >  7:130548/1‑76 (MQ=255)
                                                                                                                       tttATTGCTGATGCCACCCGCCGCGAAATTGAAATAAAAAACCCGATGCGCAGATCATCGGGTTCATTTCAATTGa                                                                                               >  1:178183/1‑76 (MQ=255)
                                                                                                                       tttATTGCTGATGCCACCCGCCGCGAAATTGAAATAAAAAACCCGATGCGCAGATCATCGGGTTCATTTCAATTGa                                                                                               <  8:130548/76‑1 (MQ=255)
                                                                                                                       tttATTGCTGATGCCACCCGCCGCGAAATTGAAATAAAAAACCCGATGCGCAGATCATCGGGTTCATTTCAATTGa                                                                                               <  2:178183/76‑1 (MQ=255)
                                                                                                                           ttGCTGATGCCACCCGCCGCGAAATTGAAATAAAAAACCCGATGCGCAGATCATCGGGTTCATTTCAATTGAGGAAATCGGGAGAATTACGCCACTTCTGACGCCGCGCTATCTAACAGCGCCAtct                                        >  1:75809/1‑127 (MQ=255)
                                                                                                                           ttGCTGATGCCACCCGCCGCGAAATTGAAATAAAAAACCCGATGCGCAGATCATCGGGTTCATTTCAATTGAGGAAATCGGGAGAATTACGCCACTTCTGACGCCGCGCTATCTAACAGCGCCAtct                                        <  2:75809/127‑1 (MQ=255)
                                                                                                                                   gCCACCCGCCGCGAAATTGAAATAAAAAACCCGATGCGCAGATCATCGGGTTCATTTCAATTGAGGAAATCGGGAGAATTACGCCACTTCTGACGCCGCGCTATCTAAc                                                  >  5:82232/1‑109 (MQ=255)
                                                                                                                                   gCCACCCGCCGCGAAATTGAAATAAAAAACCCGATGCGCAGATCATCGGGTTCATTTCAATTGAGGAAATCGGGAGAATTACGCCACTTCTGACGCCGCGCTATCTAAc                                                  <  6:82232/109‑1 (MQ=255)
                                                                                                                                      aCCCGCCGCGAAATTGAAATAAAAAACCCGATGCGCAGATCATCGGGTTCATTTCAATTGAGGAAATCGGGAGAATTACGCCACTTCTGACGCCGCGCTATCTAACAGCGCCAtct                                        >  3:6714/1‑116 (MQ=255)
                                                                                                                                      aCCCGCCGCGAAATTGAAATAAAAAACCCGATGCGCAGATCATCGGGTTCATTTCAATTGAGGAAATCGGGAGAATTACGCCACTTCTGACGCCGCGCTATCTAACAGCGCCAtct                                        <  4:6714/116‑1 (MQ=255)
                                                                                                                                           ccgcGAAATTGAAATAAAAAACCCGATGCGCAGATCATCGGGTTCATTTCAATTGAGGAAATCGGGAGAATTACGCCACTTCTGACGCCGCGCTATCTAACAGCGCCATCTCTTCGCTGTTCAGCAGTTTTTCGATGTTCACCAGAATc  <  2:70423/149‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                              
TGTAACAGTGCGAATACCCCCAGCACCATTACCAACAGCGTGACTACGCGGATACGGTTAATCATAAGGCACCTTCCTGAAAACAAGTTGATCTCGTTATCGGCAAGGAGGGGGGAAACTTTATTGCTGATGCCACCCGCCGCGAAATTGAAATAAAAAACCCGATGCGCAGATCATCGGGTTCATTTCAATTGAGGAAATCGGGAGAATTACGCCACTTCTGACGCCGCGCTATCTAACAGCGCCATCTCTTCGCTGTTCAGCAGTTTTTCGATGTTCACCAGAATC  >  NC_000913/1972627‑1972914

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