New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913_3_bme_pgi = 454269135 (0.790)8 (0.190) 8/268 4.1 19.1% intergenic (+57/‑425) fimE/fimA tyrosine recombinase/inversion of on/off regulator of fimA/major type 1 subunit fimbrin (pilin)
?NC_000913_3_bme_pgi = 4542987 35 (0.840)intergenic (+353/‑129) fimE/fimA tyrosine recombinase/inversion of on/off regulator of fimA/major type 1 subunit fimbrin (pilin)
Rejected: Coverage evenness skew score above cutoff.
Rejected: Frequency below/above cutoff threshold.

GTTATACCGCCAGTAATGCTGCTCGTTTTGCCGGATTATGGGAAAGAAATAATCTCATAAACGAAAAATTAAAAAGAGAAGAGGTTTGATTTAACTTATTGATAATAAAGTTAAAAAAACAAATAAATACAAGACAATTGGGGCCAAACTGTCCATATCATAAATAAGTTACGTATTTTTTCTCAAGCATAAAAATATTAAAAAACGACAAAAAGCATCTAACTGTTTGATATGTAAATTATTTCTATTGTAAATTAATTTCACATCACCTCCGCTATATGTAAAGCTAA  >  NC_000913_3_bme_pgi/4542546‑4542835
                                                                                                                                                 |                                                                                                                                                
gTTAGACCGCCAGTAATGCTGCTCGTTTTGCCGGATTATGGGAAAGAAATAATCTCATAAACGAAAAATTAAAAAGAGAAGAGGTTTGATTTAACTTATTGATAATAAAGTTAAAAAAACAAATAAATACAAGACAATTGGGGCCAAAc                                                                                                                                               <  1:225644/149‑1 (MQ=255)
   ataCCGCCAGTAATGCTGCTCGTTTTGCCGGATTATGGGAAAGAAATAATCTCATAAACGAAAAATTAAAAAGAGAATAGGTTTGATTTAACTTATTGATAATAAAGTTAAAAAAACAAATAAATACAAGACAATTGGGGCCAAACTGt                                                                                                                                            <  2:345370/149‑1 (MQ=255)
   acaCCGCCAGTAATGCTGCGCGTTTTGCCGGATTATGGGAAAGAAATAATCTCATAAGCGAAAAATTAAAAAGAGAAGAGGTTTGATTTAACTTATTGATAATAAAGTTAAAAAAACAAATAAATACAAGACAATTGGGGCCAAACTGt                                                                                                                                            <  1:61128/147‑1 (MQ=255)
                  ctgctCGTCTTGCCGGATTATGGGAAAGAAATAATCTCATAAACGAAAAATTAAAAAGAGAAGAGGTTTGATTTAACTTATTGATAATAAAGTTAAAAAAACAAATAAATACAAGACAATTGGGGCCAAACTGTCCATATCataaataa                                                                                                                             <  1:684082/149‑1 (MQ=255)
                              ccGGATTATGGGAAAGAAATAATCTCATAAACGAAAAATTAAAAAGAGAAGAGGTTTGATTTAACTTATTGATAATAAAGTTAAAAAAACAAATAAATACAAGACAATTGGGGCCattttga                                                                                                                                            >  2:633663/1‑116 (MQ=255)
                              ccGGATTATGGGAAAGAAATAATCTCATAAACGAAAAATTAAAAAGAGAAGAGGTTTGATTTAACTTATTGATAATAAAGTTAAAAAAACAAATAAATACAAGACAATTGGGGCCattttga                                                                                                                                            <  1:633663/122‑7 (MQ=255)
                               cgttttATGGGAAAGAACACATGTCAAAACAGATATATTAAAGAGAAATGACGTGTTAGCTAATGTTTTTATAAGAAATGTAAAAAATCAAATATCTAGTGTACAATTGGGGCCTATCTGTCCATATCATAAATAAGTTACGTAttttt                                                                                                                <  2:729788/145‑1 (MQ=255)
                                 gATTATGGGAAAGAAATAATCTCATAAACGAAAAATTAAAAAGAGAAGAGGTTTGATTTAACTTATTGATAATAAAGTTAAAAAAACAAATAAATACAAGACAATTGGGGCCAAACTGTCCATATCATAAATAAGTTACGTATTTTTtc                                                                                                              >  1:647968/1‑149 (MQ=255)
                                 gATTATGGGAAAGAAATAATCTCATAAACGAAAAATTAAAAAGAGAAGAGGTTTGATTTAACTTATTGATAATAAAGTTAAAAAAACAAATAAATACAAGACAATTGGGGCCAAACTGTCCATATCATAAATAAGTTACGTATTTTTtc                                                                                                              >  2:363515/1‑149 (MQ=255)
                                     actggAAAGAACTAATCTCATGTGCGAAAAAATAAAAAGAGAATATGTTTTATTTAAATTATTTATAATAAAGTTAAAAAAACAAATAAATACAAGACAATTGGGGCCTATTTGTCCATATCATAAATAAGTTACGTATTTTTTCTCaa                                                                                                          <  2:200418/146‑1 (MQ=255)
                                              aaaTAATCTCATAAACGAAAAATTAAAAAGAGAAGAGGTTTGATTTAACTTATTGATAATAAAGTTAAAAAAACAAATAAATACAAGACAATTGGGGCCAAACTGTCCATATCATAAATAAGTTACGTATTTTTTCTCAAGCATaaaaa                                                                                                 <  1:363515/149‑1 (MQ=255)
                                                 taatCTCATAAACGAAAAATTAAAAAGAGAAGAGGTTTGATTTAACTTATTGATAATAAAGTTAAAAAAACAAATAAATACAAGACAATTGGGGCCAAACTGTCCATATCATAAATAAGTTACGTATTTTTTCTCAAGCATAAAAatat                                                                                              <  1:545927/149‑1 (MQ=255)
                                                   atCTCATAAACGAAAAATTAAAAAGAGAAGAGGTTTGATTTAACTTATTGATAATAAAGTTAAAAAAACAAATAAATACAAGACAATTGGGGCCAAACTGTCCATATCATAAATAAGTTACGTATTTTTTCTCAAGCATAAAAATATTa                                                                                            <  1:270579/149‑1 (MQ=255)
                                                    tctcATAAACGAAAAATTAAAAAGAGAAGAGGTTTGATTTAACTTATTGATAATAAAGTTAAAAAAACAAATAAATACAAGACAATTGGGGCCAAACTGTCCATATCATAAATAAGTTACGTATTTTTTCTCAAGCATAAAAATATTaa                                                                                           >  1:36470/1‑149 (MQ=255)
                                                              gAAAAATTAAAAAGAGAAGAGGATTGATTTAACTTATTAATAATAAAGTTAAAAAAACAAATAAATACAGGACAAATAGGGCCAAAATGTCCATATAAAAAAGAGGTTACGTATTTTTTCTCAAGCGTAAAAATATTAACAAACGaaaa                                                                                 >  2:370407/1‑149 (MQ=255)
                                                                   aTTAAAAAGAGAAGAGGTTTGATTTAACTTATTGATAATAAAGTTAAAAAAACAAATAAATACAAGACAATTGGGGCCAAACTGTCCATATCATAAATAAGTTACGTATTTTTTCTCAAGCATAAAAATATTAAAAAACGACAAAAAGc                                                                            >  1:113638/1‑149 (MQ=255)
                                                                     tAAAAAGAGAAGAGGTTTGATTTAACTTATTGATAATAAAGTTAAAAAAACAAATAAATACAAGACAATTGGGGCCAAACTGTCCATATCATAAATAAGTTACGTATTTTTTCTCAAGCATAAAAATATTAAAAAACGACAAAAAGCAt                                                                          >  1:697365/1‑149 (MQ=255)
                                                                        aaaGAGAAGAGGTTTGATTTAACTTATTGATAATAAAGTTAAAAAAACAAATAAATACAAGACAATTGGGGCCAAACTGTCCATATCATAAATAAGTTACGTATTTTTTCTCAAGCATAAAAATATTAAAAAACGACAAAAAGCATCTa                                                                       >  2:493670/1‑149 (MQ=255)
                                                                          agagaagagGTTTGATTTAACTTATTGATAATAAAGTTAAAAAAACAAATAAATACAAGGCATTTAGGGCAAAACTGTCCATATCATACATAAGTTACGTATTTTTTCTCAAGCATAAAAATATTAAAAAACGACAAACACCACCTAAc                                                                     >  1:161739/1‑149 (MQ=255)
                                                                              aagagGTTTGATTTAACTTATTGATAATAAAGTTAAAAAAACAAATAAATACAAGACAATTGGGGCCAAACTGTCCATATCATAAATAAGTTACGTATTTTTTCTCAAGCATAAAAATATTAAAAAACGACAAAAAGCATCTAACTGtt                                                                 >  1:299800/1‑149 (MQ=255)
                                                                                  ggTTTGATTTAACTTATTGATAATAAAGTTAAAAAAACAAATAAATACAAGACAATTGGGGCCAAACTGTCCATATCATAAATAAGTTACGTATTTTTTCTCAAGCATAAAAATATTAAAAAACGACAAAAAGCATCTAACTGTTTGat                                                             <  2:647968/149‑1 (MQ=255)
                                                                                   gTTTGATTTAACTTATTGATAATAAAGTTAAAAAAACAAATAAATACAAGACAATTGGGGCCAAACTGTCCATATCATAAATAAGTTACGTATTTTTTCTCAAGCATAAAAATATTAAAAAAAGACAAAAAGAATANAAATTTTTgttt                                                            >  1:87336/1‑146 (MQ=255)
                                                                                    tttGATTTAACTTATTGATAATAAAGTTAAAAAAACAAATAAATACAAGACAATTGGGGCCAAACTGTCCATATCATAAATAAGTTACGTATTTTTTCTCAAGCATAAAAATATTAAAAAACGACAAAAAGCATCTAACTGTTTGatat                                                           <  2:113638/149‑1 (MQ=255)
                                                                                         tttAACTTATTGATAATAAAGTTAAAAAAACAAATAAATACAAGACAATTGGGGCCAAACTGTCCATATCATAAATAAGTTACGTATTTTTTCTCAAGCATAAAAATATTAAAAAACGACAAAAAGCATCTAACTGTTTGATATGTaaa                                                      <  1:582719/149‑1 (MQ=255)
                                                                                                  ttGATAATAAAGTTAAAAAAACAAATAAATACAAGACAATTGGGGCCAAACTGTCCATATCATAAATAAGTTACGTATTTTTTCTCAAGCATAAAAATATTAAAAAACGACAAAAAGCATCTAACTGTTTGATATGTAAATTATTTCTa                                             <  1:726348/149‑1 (MQ=255)
                                                                                                  ttGATAATAAAGTTAAAAAAACAAATAAATACAAGACAATTGGGGCCAAACTGTCCATATCATAAATAAGTTACGTATTTTTTCTCAAGCATAAAAATATTAAAAAACGACAAAAAGCATCTAACTGTTTGATATGTAAATTATTTCTa                                             <  2:299800/149‑1 (MQ=255)
                                                                                                   tGATAATAAAGTTAAACAAACAAATAAATACAAGACAATAGGGGCCAAACTGTCCATATCATAAATAAGTTACGTATTTTTTCTCAAGCATAAAAATAATAAAAAACGACAAGAAGCATCTAACTGTTTTATATATAAACTATCTCCAt                                            >  2:311441/1‑149 (MQ=255)
                                                                                                   tGATAATAAAGTTAAAAAAACAAATAAATACAAGACAATTGGGGCCAAACTGTCCATATCATAAATAAGTTACGTATTTTTTCTCAAGCATAAAAATATTAAAAAACGACAAAAAGCATCTAACTGTTTGATATGTAAATTATTTCTAt                                            >  2:723084/1‑149 (MQ=255)
                                                                                                         taaAGTTAAAAAAACAAATAAATACAAGACAATTGGGGCCAAACTGTCCATATCATAAATAAGTTACATATTTTTTCTCAAGCATAAAAATATTAAAAAACGACAAAAAGCATCTAACTGTTTTATATGTAATTTATTTCTATTGTaaa                                      >  2:65328/1‑149 (MQ=255)
                                                                                                                  aaaaaCAAATAAATACAAGACAATTGGGGCCAAACTGTCCATATCATAAATAAGTTACGTATTTTTTCTCAAGCATAAAAATATTAAAAAACGACAAAAAGCATCTAACTGTTTGATATGTAAATTATTTCTATTGTAAATTAATTTca                             <  1:493670/149‑1 (MQ=255)
                                                                                                                               taCAAGACAATTGGGGCCAGACTGTCGATATCATAAATAAGTTACGTATTTTTTCTCAAGCATAAAAATATTAAAAAAGAACAAAATGCATCTAACTGTTTGATGTGTAAATTATTTATGTTGTAAATTAATTTTAAATCACCTCCgtt                >  2:251845/1‑147 (MQ=255)
                                                                                                                                         ttGGGGCCAAACTGTCCATATCATAAATAAGTTACGTATTTTTTCTCAAGCATAAAAATATTAAAAAACGCCAAAAAGCATCTAACTGTTTGATATGTAAATTATTTCTATTGTAAATTAATTTCACATCACCTCCGCTATATGTAAAg      >  1:501862/1‑149 (MQ=255)
                                                                                                                                             agCCAAACTGTCCATATCATAAATAAGTTACGTATTTTTTCTCAAGCATAAAAATATTAAAAAACGACAAAAAGCATCTAACTGTTTGATATGTAAATTATTTCTATTGTAAATTAATTTCACATCACCTCCGCTATATGTAAAGCTaa  <  1:178953/148‑1 (MQ=255)
                                                                                                                                              gCCAAACTGTCCATATCATAAATAAGTTACGTATTTTTTCTCAAGCATAAAAATATTaaaaaa                                                                                       <  1:205869/63‑1 (MQ=255)
                                                                                                                                              gCCAAACTGTCCATATCATAAATAAGTTACGTATTTTTTCTCAAGCATAAAAATATTaaaaaa                                                                                       >  2:205869/1‑63 (MQ=255)
                                                                                                                                               ccAAACTGTCCATATCATAAATAAGTTACGTATTTTTTCTCAAGCATAAAAATATTAAAAAACGACAAAAAGCATCTAACTGTTTGATATGTAAATTATTTCTATTGTAAATTAATTTCACATCACCTCCGCTATATGTAAAg      <  2:437766/143‑1 (MQ=255)
                                                                                                                                               ccAAACTGTCCATATCATAAATAAGTTACGTATTTTTTCTCAAGCATAAAAATATTAAAAAACGACAAAAAGCATCTAACTGTTTGATATGTAAATTATTTCTATTGTAAATTAATTTCACATCACCTCCGCTATATGTAAAg      >  1:437766/1‑143 (MQ=255)
                                                                                                                                                 |                                                                                                                                                
GTTATACCGCCAGTAATGCTGCTCGTTTTGCCGGATTATGGGAAAGAAATAATCTCATAAACGAAAAATTAAAAAGAGAAGAGGTTTGATTTAACTTATTGATAATAAAGTTAAAAAAACAAATAAATACAAGACAATTGGGGCCAAACTGTCCATATCATAAATAAGTTACGTATTTTTTCTCAAGCATAAAAATATTAAAAAACGACAAAAAGCATCTAACTGTTTGATATGTAAATTATTTCTATTGTAAATTAATTTCACATCACCTCCGCTATATGTAAAGCTAA  >  NC_000913_3_bme_pgi/4542546‑4542835

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.