New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913 = 454269052 (1.020)4 (0.090) 4/110 2.3 8.0% intergenic (+57/‑425) fimE/fimA regulator for fimA/type 1 fimbriae major subunit
?NC_000913 = 4542986 47 (1.080)intergenic (+353/‑129) fimE/fimA regulator for fimA/type 1 fimbriae major subunit
Rejected: Frequency below/above cutoff threshold.

AAGAGAAGAGGTTTGATTTAACTTATTGATAATAAAGTTAAAAAAACAAATAAATACAAGACAATTGGGGCCA‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  >  NC_000913/4542618‑4542690
‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑ttggggccaTTTTGACTCATAGAGGAAAGCATCGCGGACAAACTTTTTCAGTTTATTTGTTGGCTTAATAT  <  NC_000913/4542986‑4542925
                                                                                                                                       
AAGAGAAGAGGTTTGATTTAACTTATTGATAATAAAGTTAAAAAAACAAATAAATACAAGACAAT                                                                        <  1:707754/65‑1
 AGAGAAGAGGTTTGATTTAACTTATTGATAATAAAGTTAAAAAAACAAATAAATACAAGACAATT                                                                       >  1:3643072/1‑65
    GAAGAGGTTTGATTTAACTTATTGATAATAAAGTTAAAAAAACAAATAAATACAAGACAATTGGG                                                                    <  1:1953480/65‑1
      AGAGGTTTGATTTAACTTATTGATAATAAAGTTAAAAAAACAAATAAATACAAGACAATTGGGGC                                                                  >  1:1636969/1‑65
      AGAGGTTTGATTTAACTTATTGATAATAAAGTTAAAAAAACAAATAAATACAAGACAATTGGGGC                                                                  >  1:3215608/1‑65
           TTTGATTTAACTTATTGATAATAAAGTTAAAAAAACAAATAAATACAAGACAATTGGGGCCATTT                                                             >  1:2541938/1‑65
                        ATTGATAATAAAGTTAAAAAAACAAATAAATACAAGACAATTGGGGCCATTTTGACTCATAGAGG                                                <  1:2250630/65‑1
                          TGATAATAAAGTTAAAAAAACAAATAAATACAAGACAATTGGGGCCATTTTGACTCATAGAGGAA                                              >  1:1016866/1‑65
                                           AAACAAATAAATACAAGACAATTGGGGCCATTTTGACTCATAGAGGAAAGCATCGCGGACAAACT                             <  1:2657675/65‑1
                                                                  gtGGCCATTTTGACTCATAGAGGAAAGCATCGCGGACAAACTTTTTCAGTTTATTTGTTGGCTTA      <  1:3351030/63‑1
                                                                    GGCCATTTTGACTCATAGAGGAAAGCATCGCGGACAAACTTTTTCAGTTTATTTGTTGGCTTAAT    <  1:2282263/65‑1
                                                                    GGCCATTTTGACTCATAGAGGAAAGCATCGCGGACAAACTTTTTCAGTTTATTTGTTGGCTTAAT    <  1:3483101/65‑1
                                                                    GGCCATTTTGACTCATAGAGGAAAGCATCGCGGACAAACTTTTTCAGTTTATTTGTTGGCTTAAT    <  1:3210042/65‑1
                                                                    GGCCATTTTGACTCATAGAGGAAAGCATCGCGGACAAACTTTTTCAGTTTATTTGTTGGCTTAAT    <  1:3660579/65‑1
                                                                      CCATTTTGACTCATAGAGGAAAGCATCGCGGACAAACTTTTTCAGTTTATTTGTTGGCTTAATAT  <  1:3117478/65‑1
                                                                      CCATTTTGACTCATAGAGGAAAGCATCGCGGACAAACTTTTTCAGTTTATTTGTTGGCTTAATAT  <  1:3707667/65‑1
                                                                      CCATTTTGACTCATAGAGGAAAGCATCGCGGACAAACTTTTTCAGTTTATTTGTTGGCTTAATAT  <  1:3579907/65‑1
                                                                      CCATTTTGACTCATAGAGGAAAGCATCGCGGACAAACTTTTTCAGTTTATTTGTTGGCTTAATAT  <  1:658507/65‑1
                                                                      CCATTTTGACTCATAGAGGAAAGCATCGCGGACAAACTTTTTCAGTTTATTTGTTGGCTTAATAT  >  1:726930/1‑65
                                                                      CCATTTTGACTCATAGAGGAAAGCATCGCGGACAAACTTTTTCAGTTTATTTGTTGGCTTAATAT  <  1:3393916/65‑1
                                                                      CCATTTTGACTCATAGAGGAAAGCATCGCGGACAAACTTTTTCAGTTTATTTGTTGGCTTAATAT  <  1:3220743/65‑1
                                                                      CCATTTTGACTCATAGAGGAAAGCATCGCGGACAAACTTTTTCAGTTTATTTGTTGGCTTAATAT  <  1:2434233/65‑1
                                                                      CCATTTTGACTCATAGAGGAAAGCATCGCGGACAAACTTTTTCAGTTTATTTGTTGGCTTAATAT  <  1:1466873/65‑1
                                                                      CCATTTTGACTCATAGAGGAAAGCATCGCGGACAAACTTTTTCAGTTTATTTGTTGGCTTAATAT  <  1:138185/65‑1
                                                                      CCATTTTGACTCATAGAGGAAAGCATCGCGGACAAACTTTTTCAGTTTATTTGTTGGCTTAATAT  <  1:1294578/65‑1
                                                                      CCATTTTGACTCATAGAGGAAAGCATCGCGGACAAACTTTTTCAGTTTATTTGTCGGCTTAATAT  <  1:2141121/65‑1
                                                                      CCATTTTGACGCATAGAGGAAAGCATCGCGGACAAACTTTTTCAGTTTATTTGTTGGCTTAATAT  <  1:209133/65‑1
                                                                                                                                       
AAGAGAAGAGGTTTGATTTAACTTATTGATAATAAAGTTAAAAAAACAAATAAATACAAGACAATTGGGGCCA‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  >  NC_000913/4542618‑4542690
‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑ttggggccaTTTTGACTCATAGAGGAAAGCATCGCGGACAAACTTTTTCAGTTTATTTGTTGGCTTAATAT  <  NC_000913/4542986‑4542925

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