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)

CGTTATACCGCCAGTAATGCTGCTCGTTTTGCCGGATTATGGGAAAGAAATAATCTCATAAACGAAAAATTAAAAAGAGAAGAGGTTTGATTTAACTTATTGATAATAAAGTTAAAAAAACAAATAAATACAAGACAATTGGGGCCAAACTGTCCATATCATAAATAAGTTACGTATTTTTTCTCAAGCATAAAAATATTAAAAAACGACAAAAAGCATCTAACTGTTTGATATGTAAATTATTTCTATTGTAAATTAATTTCACATCACCTCCGCTATATGTAAAG  >  NC_000913_3_pae_tpiA/4542460‑4542746
                                                                                                                                          |                                                                                                                                                    
ccagagttttaattttcatgctgctttcctttcaaaaaactatttctaaatcgacatgggcagtcgttctgtacactttgttttgtaaacgagtttgactgccaacactgcacagttttcccccaaaagatgaaacatTTGGTGCCaaa                                                                                                                                            <  1:270361‑M1/11‑1 (MQ=255)
     taCCGCCAGTAAAGCGTCTCGCTATGACAGATTATGGGAAAGAAATAATCTTATAAAAAAAAAAATAAAAAGAGAAGAGGTTTGATTTAACTTATTGATAATAAAGTTAAAAAAACAAATAAATACAAGACAATTGGGGCCAAACTGTc                                                                                                                                       <  2:209429/149‑1 (MQ=255)
      aCCGCCAGTAATGCTGCTCGTTTTGCCGGATTATGGGAAAGAAATAATCTCATAAACGAAAAATTAAAAAGAGAAGAGGTTTGATTTAACTTATTGATAATAAAGTTAAAAAAACAAATAAATACAAGACAATTGGGGCCAAACTGTcc                                                                                                                                      <  1:71188/149‑1 (MQ=255)
         gCCAGTAATGCTGCTCGTTTTGCCGGATTATGGGAAAGAAATAATCTCATAAACGAAAAATTAAAAAGAGAAGAGGTTTGATTTAACTTATTGATAATAAAGTTAAAAAAACAAATAAATACAAGACAATTAGGGCCAAACTATCCata                                                                                                                                   >  1:207274/1‑149 (MQ=255)
             gTAATGCTGCTCGTTTTGCCGGATTATGGGAAAGAAATAATCTCATAAACGAAAAATTAAAAAGAGAAGAGGTTTGATTTAACTTATTGATAATAAAGTTAAAAAAACAAATAAATACAAGACAATTGGGGCCAAACTGTCCATATCat                                                                                                                               >  1:126740/1‑149 (MQ=255)
                 tgctgcTCGTTTTGCCGGATTATGGGAAAGAAATAATCTCATAAACGAAAAATTAAAAAGAGAAGAGGTTTGATTTAACTTATTGATAATAAAGTTAAAAAAACAAATAAATACAAGACAATTGGGGCCAAACTGTCCATATCataaat                                                                                                                           >  1:330822/1‑149 (MQ=255)
                   ctgctCGTTTTGCCGGATTATGGGAAAGAAATAATCTCATAAACGAAAAATTAAAAAGAGAAGAGGTTTGATTTAACTTATTGATAATAAAGTTAAAAAAACAAATAAATACAAGACAATTGGGGCCAAACTGTCCATATCataaataa                                                                                                                         <  2:249190/149‑1 (MQ=255)
                   ctgctCGTTTTGCCGGATTATGGGAAAGAAATAATCTCATAAACGAAAAATTAAAAAGAGAAGAGGTTTGATTTAACTTATTGATAATAAAGTTAAAAAAACAAATAAATACAAGACAATTGGGGCCAAACTGTCCATATCataaataa                                                                                                                         <  2:249215/149‑1 (MQ=255)
                    ctgctttcctttcaaaaaactatttctaaatcgacatgggcagtcgttctgtacactttgttttgtaaacgagtttgactgccaacgttgctcagggtgcccccaaaagatgaaacatTTGGGGCCAAACTGTCCATATCATAAATAAg                                                                                                                        <  1:250203‑M1/31‑1 (MQ=255)
                      ctCGTTTTGCCGGATTATGGGAAAGAAATAATCTCATAAACGAAAAATGAAAAAGCGAAGAGGTTTGATTTAACTTATTGATAATAAAGTTAAAAAAACAAATAAATACAAGACAATTGGGGCCAAACTGTCCATATCATAAATAAGtt                                                                                                                      >  2:3129/1‑149 (MQ=255)
                                   aaaactatttctaaatcgacatgggcagtcgttctgtacactttgttttgtaaacgagtttgactgccaacactgcacagttttcccccaaaagatgaaacatTTGGGGCCAAACTGTCCATATCATAAATAAGTTACGTATTTTTtct                                                                                                         <  1:181232‑M1/46‑1 (MQ=255)
                                      aTGGGAAAGAAATAATATCATAAACGAAAAATTAAAAAGAGAAGAGGTTTGATTTAACTTATTGATAATAAAGTTAAAAAAACAAATAAATACAAGACAATTGGGGCCGAACTGTCCATATCATAAATAAGTTACGTATTTTTTCTCaa                                                                                                      <  2:97796/149‑1 (MQ=255)
                                                       tcATAAACGAAAAATTAAAAAGAGAAGAGGTTTGATTTAACTTATTGATAATAAAGTTAAAAAAACAAATAAATACAAGACAATTGGGGCCAAACTGTCCATATCATAAATAAGTTACGTATTTTTTCTCAACCATAAAAATATTaaaa                                                                                     >  2:350797/1‑149 (MQ=255)
                                                         aTAAACGAAAAATTAAAAAGAGAAGAGGTTTGATTTAACTTATTGATAATAAAGTTAAAAAAACAAATAAATACAAGACAATTGGGGCCAAACTGTCCATATCATAAATAAGTTACGTATTTTTTCTCAAGCATAAAAATATTaaaaaa                                                                                   <  2:330822/149‑1 (MQ=255)
                                                            cagtcgttctgtacactttgttttgtaaacgagtttgactgccaacactgcacagttttcccccaaaagatgaaacatTTGGGGCCAAACTGTCCATATCATAAATAAGTTACGTATTTTTTCTCAAGCATAAAAATATTAAAAAACGa                                                                                >  2:181232‑M1/79‑149 (MQ=255)
                                                             aCGAAAAATTAAAAAGAGAAGAGGTTTGATTTAACTTATTGATAATAAAGTTAAAAAAACAAATAAATACAAGACAATTGGGGCCAAACTGTCCATATCATAAATAAGTTACGTATTTTTTCTCAAGCATAAAAATATTAAAAAACGAc                                                                               >  1:217163/1‑149 (MQ=255)
                                                                                agagGTAAGATATAACTTATTGATAAAAAAGTTAAAAAAACAAATAAATACAAGACAATTGGGGCCAAACTGCCCATAACATAAATAAGTTACGTATTTTTTCTCAAGCATAAAAATATTAAAAAACGACAAAAAGCATCTAACTGttt                                                            <  2:93121/149‑1 (MQ=255)
                                                                                 gagGTTTGATTTAACTTATTGATAATAAAGTTAAAAAAACAAATAAATACAAGACAATTGGGGCCAAACTGTCCATATCATAAATAAGTTACGTATTTTTTCTCAAGCATAAAAATATTAAAAAACGACAAAAAGCATTTAACTGtttt                                                           >  2:56272/1‑148 (MQ=255)
                                                                                   cgTTTGATTTAACATATTAATAATAAAGTTAAAAAAACAAAAAAATACCAGACAATTGGGGCCAAACTGTACATATCATAAATAAGTTACGTATTTTTTCTCAAGCATAAAAATATTAAAAAACGACAAAAAGCATCTAACTGTTTGat                                                         <  2:78330/148‑1 (MQ=255)
                                                                                     tttGATTTAACTTATTGATAATAAAGTTAAAAAAACAAATAAATACAAGACAATTGGGGCCAAACTGTCCATATCATAAATAAGTTACGTATTTTTTCTCAAGCATAAAAATATTAAAAAACGACAAAAAGAATCTAACTGTTTGatat                                                       >  2:280085/1‑149 (MQ=255)
                                                                                        gATTTAACTTATTGATAATAAAGTTAAAAAAACAAATAAATACAAGACAATTGGGGCCAAACTGTCCATATCATAAATAAGTTACGTATTTTTTCTCAAGCATAAAAATATTAAAAAACGACAAAAAGCATCTAACTGTTTGATATGTa                                                    >  1:270235/1‑149 (MQ=255)
                                                                                           ttAACTTATTGATAATAAAGTTAAAAAAACAAATAAATACAAGACAATTGGGGCCAAACTGTCCATATCATAAATAAGTTACGTATTTTTTCTCAAGCATAAAAATATTGTAACTCTACAAAAAGCATCTAACTGTTTGATATGTAAat                                                 <  1:510099/149‑1 (MQ=255)
                                                                                                    tGATAAAAAATTTAAAAAAATAAATAAATACAAGACAATTGGGGCCAAACTGTCCATATCATAAATAAGTTACGTCTTTTTTCTCAAGCATAAAAATATTAAAAAACGACAAAAAGCATCTAACTGTTTGATATGTAAATTATTTCTAt                                        <  2:207274/149‑1 (MQ=255)
                                                                                                        aataaAGTTAAAAAAACAAATAAATACAAGACAATTGGGGCCAAACTGTCCATATCATAAATAAGTTACGTATTTTTTCTCAAGCATAAAAATATTAAAACACGACAAAAAGCATCTAACTGTTTGata                                                        >  2:108452/1‑129 (MQ=255)
                                                                                                        aataaAGTTAAAAAAACAAATAAATACAAGACAATTGGGGCCAAACTGTCCATATCATAAATAAGTTACGTATTTTTTCTCAAGCATAAAAATATTAAAAAACGACAAAAAGCATCTAACTGTTTGata                                                        <  1:108452/129‑1 (MQ=255)
                                                                                                         ataaAGTTAAAAAAACAAATAAATACAAGACAATTGGGGCCAAACTGTCCATATCATAAATAAGTTACGTATTTTTTCTCAAGCATAAAAATATTAAAAAACGACAAAAAGCATCTAACTGTTTGATATGTAAATTATTTCTATTGTaa                                   >  2:471048/1‑149 (MQ=255)
                                                                                                          aagaaatAAAAAGACAAATAACAAGAAGTCAATTGGGGCAAAACTTTGCATATAAAAAATAAGTTACGTATTTTTTGGCAAGCATAAAAATATTAAAAAACGACAAAAAGCATCTAACTGTTTGATATGTAAATTATTTCTATTGTaaa                                  <  2:353463/143‑1 (MQ=255)
                                                                                                               acagttttccccaaaaagatgaaacatTTGGGGCCAAACTGTCCATATCATAAATAAGTTACGTATTTTTTCTCAAGCATAAAAATATTAAAAAACGACAAAAAGCATCTAACTGTTTGATATGTAAATTATTTCTATTGTAaattaat                             >  2:146216‑M1/28‑149 (MQ=255)
                                                                                                                  aaaaaaCAAATAAATACAAGACAATTGGGGCCAAACTGTCCATATCataaataa                                                                                                                         <  1:74166/54‑1 (MQ=255)
                                                                                                                  aaaaaaCAAATAAATACAAGACAATTGGGGCCAAACTGTCCATATCataaataa                                                                                                                         >  2:74166/1‑54 (MQ=255)
                                                                                                                    aaaaCAAATAAATACAAGACCATTGGGGCCAAACTGTCCAAATCATAAATAAGTTACGTATTTTTTCTCAAGCATAAAAATATTAAAAAACGACAAAAATCATCTAACTGTTTGATATGTAAATTATTTCTATTGTAAATTAATTTcac                        <  1:105173/149‑1 (MQ=255)
                                                                                                                       aCAAATAAATACAAGACAATTGGGGCCAAACTGTCCATATCATAAATAAGTTACGTATTTTTTCTCAAGCATAAAAATATTAAAAAACGACAAAAAGCATCTAACTGTTTGATATGTAAATTATTTCTATTGTAAATTAATTTCACATc                     >  1:363524/1‑149 (MQ=255)
                                                                                                                         cccaaaagatgaaacatTTGGGGCCAAACTGTCCATATCATAAATAAGTTACGTATTTTTTCTCAAGCATAAAAATATTAAAAAACGACAAAAAGCATCTAACTGTTTGATATGTAAATTATTTCTATTGTAAATTAATTTCACATCAc                   <  1:475553‑M1/132‑1 (MQ=255)
                                                                                                                               agatgaaacatTTGGGGCCAAACTGTCCATATCATAAATAAGTTACGTATTTTTTATCCAGCATAAAAATATTAAAAAAAAAAAAAAAGCATCTAACTGTTTGATATGTAAATTATTTCTATTGTAAATTAATTTCACATCACCTCCGc             >  2:250203‑M1/12‑149 (MQ=255)
                                                                                                                                gatgaaacatTTGGGGCCAAACTGTCCATATCATAAATAAGTTACGTATTCCTTCTTCATTCTTTTTATATTAAAAAACGACAAAAAGCATCTAACTGTTTGATATGTAAATTATTTCTATTGTAAATTAATTTCACATCACCTCCGCt            >  1:505687‑M1/11‑149 (MQ=255)
                                                                                                                                          TTGGGGCCAAACTGTCCATATCATAAATAAGTTACGTATTTTTTCTCAAGCATAAAAATATTAAAAAACGACAAAAAGCATCTAACTGTTTGATATGTAAATTATTTCTATTGTAAATTAATTTCACATCACCTCCGCTATATGTAAAg  >  1:350797‑M1/1‑149 (MQ=255)
                                                                                                                                          |                                                                                                                                                    
CGTTATACCGCCAGTAATGCTGCTCGTTTTGCCGGATTATGGGAAAGAAATAATCTCATAAACGAAAAATTAAAAAGAGAAGAGGTTTGATTTAACTTATTGATAATAAAGTTAAAAAAACAAATAAATACAAGACAATTGGGGCCAAACTGTCCATATCATAAATAAGTTACGTATTTTTTCTCAAGCATAAAAATATTAAAAAACGACAAAAAGCATCTAACTGTTTGATATGTAAATTATTTCTATTGTAAATTAATTTCACATCACCTCCGCTATATGTAAAG  >  NC_000913_3_pae_tpiA/4542460‑4542746

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.