New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913_3_pae_tpiA 4542598 =27 (0.920)8 (0.290) 8/264 2.2 24.9% intergenic (+49/‑433) fimE/fimA tyrosine recombinase/inversion of on/off regulator of fimA/major type 1 subunit fimbrin (pilin)
?NC_000913_3_pae_tpiA 4542912 = 23 (0.840)intergenic (+363/‑119) fimE/fimA tyrosine recombinase/inversion of on/off regulator of fimA/major type 1 subunit fimbrin (pilin)

TTCCTCATTCTTCTCTCCAAAAACCACCTCATGCAATATAAACATCTATAAATAAAGATAACAATAGAATATTAAGCCAACAAATAAACTGAAAAAGTTTGTCCGCGATGCTTTCCTCTATGAGTCAAAATGGCCCCAAATGTTTCATCTTTTGGGGGAAAACTGTGCAGTGTTGGCAGTCAAACTCGTTTACAAAACAAAGTGTACAGAACGACTGCCCATGTCGATTTAGAAATAGTTTTTTGAAAGGAAAGCAGCATGAAAATTAAAACTCTGG  >  NC_000913_3_pae_tpiA/4542773‑4543049
                                                                                                                                           |                                                                                                                                         
ttCCTCATTCTTCTCTCCAAAAACCACCTCATGCAATATAAACATCTATAAATAAAGATAACAATAGAATATTAAGCCAACAAATAAACTGAAAAAGTTTGTCCGCGATGCTTTCCTCTATGAGTCAAAATGGCCCCAAATGTTTCATc                                                                                                                                  >  1:123991/1‑149 (MQ=255)
agcggaggtgatgtgaaattaatttacaatagaaataatttacatatcaaacagttagatgctttttgtcgttttttaatataaaaagaatgaagaaggaatacgtaacttatttatgatatggacagtttggccccaaATGTTTCATc                                                                                                                                  >  1:505687‑M2/140‑149 (MQ=255)
 gcggaggtgatgtgaaattaatttacaatagaaataatttacatatcaaacagttagatgctttttttttttttttaatatttttatgctggataaaaaatacgtaacttatttatgatatggacagtttggccccaaATGTTTCATCt                                                                                                                                 >  2:250203‑M2/139‑149 (MQ=255)
    tCATTCTTCTCTCCAAAAACCACCTAATGCAATATAAACATCTATAAATAAAGATAACAATAGAATATTAAGCCAACAAATAAACTGAAAAAGTTTGTCCGCGATGCTTTCCTCTATGAGTCAAAATGGCCCCAAATATTTCATCtttt                                                                                                                              >  2:59646/1‑149 (MQ=255)
     cATTCTTCTCTCCAAAAACCACCTCATGCAATATAAACATCGATAAATAAAGATACCACTCGCCTATCCGTCCAACAAATAAACTGAAAAAGTTTGTCCGCGATGCTTTCCTCTATGAGTCAAAATGGCCCCAAATGTTTCATCTTTTg                                                                                                                             >  1:228872/1‑149 (MQ=255)
      aTTCTTCTCTCCAAAAACCACCTCATGCAATATAAACCTCTATAAATAAAGATAACAATAGAATATTAAGCCAACAAATAAACTGAAAAAGTTTGTCCGCGATGCTTTCCTCTATGAGTCAAAATGGCCCCAAATGTTTCATCTTTTgg                                                                                                                            >  2:436667/1‑149 (MQ=255)
       gtgatgtgaaattaatttacaatagaaataatttacatatcaaacagttagatgctttttgtcgttttttaatatttttatgcttgagaaaaaatacgtaacttatttatgatatggacagtttggccccaaATGTTTCATCTTTTggg                                                                                                                           <  1:475553‑M2/17‑1 (MQ=255)
             tctgCAAAAACCACCCCATGCAATATAAACATCTAGAAATAAAGATAACAAGAGAAGATTAAGCCAACAAATAAACTGAAAAAGTTTGTCCGCGATGCTTTCCTCTATGAGTCAAAATGGCCCCAAATGTTTCATCTTTTGGGGGaaaa                                                                                                                     <  2:126740/149‑1 (MQ=255)
               tcCAAAAACCACCTCATGCAATATAAACATCTATAAATAAAGATAACAATAGCATATTAAGCCAACAAATAAACTGAAAAAGTTTGTCCGCGATGCTTTCCTCTATGAGTCAAAATGGCCCCAAAGGTTTAATCTTTTGGGGGAAAACt                                                                                                                   >  1:222027/1‑149 (MQ=255)
                 attaatttacaatagaaataatttacatatcaaacagttagatgctttttgtcgttttttaatatttttatgcttgagaaaaaatacgtaacttatttatgatatggacagtttggccccaaATGTTTCATCTTTTTGGGGAAAACtgt                                                                                                                 >  2:146216‑M2/123‑149 (MQ=255)
                  aaaaaCCACCTCATGCAATATAAACATCTATAAATAAAGATAACAATAGAATATTAAGCCAACAAATAAACTGAAAAAGTTTGTCCGCGATGCTTTCCTCTATGAGTCAAAATGGCCCCAAATGTTTCATCTTTTGGGGGAAAACtgtg                                                                                                                <  1:3129/149‑1 (MQ=255)
                         accTCATGCAATATAAACATCTATAAATAAAGATAACAATAGAATATTAAGCCAACAAATAAACTGAAAAAGTTTGTCCGCGATGCTTTCCTCTATGAGTCAAAATGGCCCCAAATGTTTCATCTTTTGGGGGAAAACTGTGCAGTGtt                                                                                                         <  2:270235/149‑1 (MQ=255)
                           cgcATGGAATATAAACAGCTAAAAAAAGAGATGACGATATAATATTAAGCAAACAAATAAACTTAAAAAGTTTGTCCGCGATGCTTTCCTCTATGAGTGAAAATGGCCCGAATGGTTTTATCTGTTGGGGGAAAACTGTGCAGTGTTgg                                                                                                       <  2:237948/147‑1 (MQ=255)
                           cTCATGCAATATAAACTTCTATAAATAAAGATAACAATAGAATATTAAGCCAACAAATAAACTGAAAAAGTTTGTCCGCGATGCTTTCCTCTATGAGTCAAAATGGCCCCAAATGTTTCATCTTTTGGGGGAAAACTGTGCAGTGTTgg                                                                                                       <  1:350877/149‑1 (MQ=255)
                                                   ataaaGATAACAATAGAATATTAAGCCAACAAATAAACGGAAAAAGTTTGTCGGCGATGCTTTCCTCCATGAGTCAAAATGGCACCAAATGTTTCATCTTTTGGGTGAAAACTGTGCAGTGTAGGCAGTCAAACTCGATGaaaaaacaa                                                                               >  2:67510/1‑149 (MQ=255)
                                                             cAATAGAATATTAAGCCAACAAATAAACTGAAAAAGTTTGTCCGCGATGCTTTCCTCTATGAGTCAAAATGGCCCCAAATGTTTCATCTTTTGGGGGAAAACTGTGCAGTGTTGGCAGTCAAACTCGTTTACAAAACAAAGTGTACAGa                                                                     >  2:32165/1‑149 (MQ=255)
                                                              aaTAGAATATTAAGCCAACAAATAAACTGAAAAAGTTTGTCCGCGATGCTTTCCTCTATGAGTCAAAATGGCCCCAAATGTTTCATCTTTTGGGGGAAAACTGTGCAGTGTTGGCAGTCAAACTCGTTTACAAAACAAAGTGTACAGaa                                                                    >  1:282119/1‑149 (MQ=255)
                                                              aaTAGAATATTAAGCCAACAAATAAACTGAAAAAGTTTGTCCGCGATGCTTTCCTCTATGAGTCAAAATGGCCCCAAATGTTTCATCTTTTGGGGGAAAACTGTGCAGTGTTGGCAGTCAAACTCGTTTACAAAACAAAGTGTACAGaa                                                                    >  2:271219/1‑149 (MQ=255)
                                                                    tcgttttttaatatttttatgcttgagaaaaaatacgtaacttatttatgatatggacagtttggccccaaATGTTTCATCTTTTGGGGGAAAACTGTGCAGTGTTGGCAGTCAAACTCGTTTACAAAACAAAGTGTACAGAACGACTg                                                              >  2:181232‑M2/72‑149 (MQ=255)
                                                                                caaATAAACTGAAAAAGTTTGTCCGCGATGCTTTCCTCTATGAGTCAAAATGGCCCCAAATGTTTCATCTTTTGGGGGAAAACTGTGCAGTGTTGGCAGTCAAACTCGTTTACAAAACAAAGTGTACAGAACGACTGCCCATGTCGAtt                                                  <  2:185862/149‑1 (MQ=255)
                                                                                      aaCTGAAAAAGTTTGTCCGCGATGCTTTCCTCTATGAGTCAAAATGGCCCCAAATGTTTCATCTTTTGGGGGAAAACTGTGCAGTGTTGGCAGTCAAACTCGTTTACAAAACAAAGTGTACAGAACGACTGCCCATGTCGATTTAGaaa                                            <  1:32165/149‑1 (MQ=255)
                                                                                             agaaaaaatacgtaacttatttatgatatggacagtttggccccaaATGTTTCATCTTTTGGGGGAAAACTGTGCAGTGTTGGCAGTCAAACTCGTTTACAAAACAAAGTGTACAGAACGACTGCCCATGTCGATTTAGAAATAGtttt                                     <  1:181232‑M2/103‑1 (MQ=255)
                                                                                                            tGCTTTCCTCTATGAGTCAAAATGGCCCCAAATGTTTCATCTTTTGGGGGAAAACTGTGCAGTGTTGGCAGTCAAACTCGTTTACAAAACAAAGTGTACAGAACGACTGCCCATGTCGATTTAGAAATAGTTTTTTGAAAGGAAagcag                      <  2:282119/149‑1 (MQ=255)
                                                                                                            cttatttatgatatggacagtttggccccaaATGTTTCATCTTTTGGGGGCACCCTGAGCAACGTTGGCAGTCAAACTCGTTTACAAAACAAAGTGTACAGAACGACTGCCCATGTCGATTTAGAAATAGTTTTTTGAAAGGAAagcag                      <  1:250203‑M2/118‑1 (MQ=255)
                                                                                                             gCTTTCCTCTATGAGTCAAAATGGCCCCAAATGTTTCATCTTTTGGGGGAAAACTGTGCAGTGTTGGCAGTCAAAc                                                                                              <  2:136236/76‑1 (MQ=255)
                                                                                                             gCTTTCCTCTATGAGTCAAAATGGCCCCAAATGTTTCATCTTTTGGGGGAAAACTGTGCAGTGTTGGCAGTCAAAc                                                                                              >  1:136236/1‑76 (MQ=255)
                                                                                                              cTTTCCTCTATGAGTCAAAATGGCCCCAAATGTTTCATCTTTTGGGGGAAAACTGTGCAGTGTTGGCAGTCAAACTCGTTTACAAAACAAAGTGTACAGAACGACTGCCCATGTCGATTTAGAAATAGTTTTTTGAAAGGAAAgcagca                    >  1:285042/1‑149 (MQ=255)
                                                                                                               tttCCTCTATGAGTCAAAATGGCCCCAAATGTTTCATCTTTTGGGGGAAAACTGTGCAGTGTTGGCAGTCAAACTCGTTTACAAAACAAAGTGTACAGAACGACTGCCCATGTCGATTTAGAAATAGTTTTTTGAAAGGAAAGCAGCAt                   <  2:363524/149‑1 (MQ=255)
                                                                                                                          aGTCAAAATGGCCCCAAATGTTTCATCTTTTGGGGGAAAACTGTGCAGTGTTGGCAGTCAAACTCGTTTACAAAACAAAGTGTACAGAACGACTGCCCATGTCGATTTAGAAATAGTTTTTTGAAAGGAAAGCAGCATGAAAATTaaaa        <  1:280085/149‑1 (MQ=255)
                                                                                                                          aGTCAAAATGGCCCCAAATGTTTCATCTTTTGGGGGAAAACTGTGCAGTGTTGGCAGTCAAACTCGTTTACAAAACAAAGTGTACAGAACGACTGCCCATGTCGATTTAGAAATAGTTTTTTGAAAGGAAAGCAGCATGAAAATTaaaa        <  1:271219/149‑1 (MQ=255)
                                                                                                                                tttggcaccaaATGTTTCATCTTTTGGGGGAAAACTGTGCAGTGTTGGCAGTCAAACTCGTTTACAAAACAAAGTGTACAGAACGACTGCCCATGTCGATTTAGAAATAGTTTTTTGAAAGGAAAGCAGCATGAAAATTAAAACTCTgg  <  1:270361‑M2/138‑1 (MQ=255)
                                                                                                                                           |                                                                                                                                         
TTCCTCATTCTTCTCTCCAAAAACCACCTCATGCAATATAAACATCTATAAATAAAGATAACAATAGAATATTAAGCCAACAAATAAACTGAAAAAGTTTGTCCGCGATGCTTTCCTCTATGAGTCAAAATGGCCCCAAATGTTTCATCTTTTGGGGGAAAACTGTGCAGTGTTGGCAGTCAAACTCGTTTACAAAACAAAGTGTACAGAACGACTGCCCATGTCGATTTAGAAATAGTTTTTTGAAAGGAAAGCAGCATGAAAATTAAAACTCTGG  >  NC_000913_3_pae_tpiA/4542773‑4543049

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