New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913 4542682 =30 (0.560)3 (0.060) 3/116 3.1 7.1% intergenic (+49/‑433) fimE/fimA regulator for fimA/type 1 fimbriae major subunit
?NC_000913 4542996 = 53 (1.150)intergenic (+363/‑119) fimE/fimA regulator for fimA/type 1 fimbriae major subunit
Rejected: Coverage evenness skew score above cutoff.
Rejected: Frequency below/above cutoff threshold.

AGAGAAGAGGTTTGATTTAACTTATTGATAATAAAGTT‑AAAAAAACAAATAAATACAAGACAATTGGGGCCAAACTGTCCATATCATAAATAAGTTACGTATTTTTTCTCAAGCATAAAAATATTAAAAA  >  NC_000913/4542619‑4542748
                                                                |                                                                  
agagaagagGTTTGATTTAACTTATTGATAATAAAGTT‑AAAAAAACAAATAAATACAAGACAATTGGG                                                                >  1:1219706‑M1/1‑68 (MQ=255)
agagaagagGTTTGATTTAACTTATTGATAATAAAGTT‑AAAAAAACAAATAAATACAAGACAATTGGG                                                                >  1:1469096‑M1/1‑68 (MQ=255)
agagaagagGTTTGATTTAACTTATTGATAATAAAGTT‑AAAAAAACAAATAAATACAAGACAATTGGG                                                                <  1:1521346‑M1/68‑1 (MQ=255)
  agaagagGTTTGATTTAACTTATTGATAATAAAGTT‑AAAAAAACAAATAAATACAAGACAATTGGGGC                                                              <  1:1530541‑M1/68‑1 (MQ=255)
     agagGTTTGATTTAACTTATTGATAATAAAGTT‑AAAAAAACAAATAAATACAAGACAATTGGGGCCaa                                                           <  1:2625169/68‑1 (MQ=255)
     agagGTTTGATTTAACTTATTGATAATAAAGTT‑AAAAAAACAAATAAATACAAGACAATTGGGGCCaa                                                           <  1:2248002/68‑1 (MQ=255)
     agagGTTTGATTTAACTTATTGATAATAAAGTT‑AAAAAAACAAATAAATACAAGACAATTGGGGCCAt                                                           <  1:3490393‑M1/68‑2 (MQ=255)
     agagGTTTGATTTAACTTATTGATAATAAAGTT‑AAAAAAACAAATAAATACAAGACAATTGGGGCCAt                                                           <  1:3222449‑M1/68‑2 (MQ=255)
      gagGTTTGATTTAACTTATTGATAATAAAGTTAAAAAAAACAAATAAATACAAGACAATTAGGGCCAt                                                           <  1:1440058‑M1/68‑2 (MQ=255)
        ggTTTGATTTAACTTATTGATAATAAAGTT‑AAAAAAACAAATAAATACAAGACAATTGGGGCCAtttt                                                        >  1:1086091‑M1/1‑64 (MQ=255)
        ggTTTGATTTAACTTATTGATAATAAAGTT‑AAAAAAACAAATAAATACAAGACAATTGGGGCCAAACt                                                        >  1:855774/1‑68 (MQ=255)
          tttGATTTAACTTATTGATAATAAAGTT‑AAAAAAACAAATAAATACAAGACAATTGGGGCCAAACTGt                                                      >  1:2078094/1‑68 (MQ=255)
          tttGATTTAACTTATTGATAATAAAGTT‑AAAAAAACAAATAAATACAAGACAATTGGGGCCAAACTGt                                                      >  1:1585138/1‑68 (MQ=255)
          tttGATTTAACTTATTGATAATAAAGTT‑AAAAAAACAAATAAATACAAGACAATGGGGGCCAAACTGt                                                      <  1:859520/68‑1 (MQ=255)
               tttAACTTATTGATAATAAAGTT‑AAAAAAACAAATAAATACAAGACAATTGGGGCCAAACTGTCCata                                                 >  1:903998/1‑68 (MQ=255)
               tttAACTTATTGATAATAAAGTT‑AAAAAAACAAATAAATACAAGACAATTGGGGCCAAACTGTCCata                                                 >  1:1594678/1‑68 (MQ=255)
               tttAACTTATTGATAATAAAGTT‑AAAAAAACAAATAAATACAAGACAATTGGGGCCAAACTGTCCata                                                 >  1:3307766/1‑68 (MQ=255)
               tttAACTTATTGATAATAAAGTT‑AAAAAAACAAATAAATACAAGACAATTGGGGCCAAACTGTCCata                                                 <  1:2253831/68‑1 (MQ=255)
               tttAACTTATTGATAATAAAGTT‑AAAAAAACAAATAAATACAAGACAATTGGGGCCAAACTGTCCata                                                 >  1:1416916/1‑68 (MQ=255)
                 tAACTTATTGATAATAAAGTT‑AAAAAAACAAATAAATACAAGACAATTGGGGCCAAACTGTCCATATc                                               <  1:434154/68‑1 (MQ=255)
                  aaCTTATTGATAATAAAGTT‑AAAAAAACAAATAAATACAAGACAATTGGGGCCAttttgactcataga                                              <  1:150281‑M1/68‑15 (MQ=255)
                   aCTTATTGATAATAAAGTT‑AAAAAAACAAATAAATACAAGACAATTGGGGCCAAACTGTCCATATCat                                             >  1:1348392/1‑68 (MQ=255)
                    cTTATTGCTAATAAAGTT‑AAAAAAACAAATAAATACAAGACAATTGGGGCCAAACTGTCCATATCata                                            >  1:595890/1‑68 (MQ=255)
                    cTTATTGATAATAAAGTT‑AAAAAAACAAATAAATACAAGACAATTGGGGCCAAACTGTCCATATcctg                                            >  1:3791460/1‑65 (MQ=255)
                    cTTATTGATAATAAAGTT‑AAAAAAACAAATAAATACAAGACAATTGGGGCCAAACTGTCCATATCata                                            >  1:2989017/1‑68 (MQ=255)
                        ttGATAATAAAGTT‑AAAAAAACAAATAAATACAAGACAATTGGGGCCAAACTGTCCATATCataaata                                        <  1:1095473/68‑1 (MQ=255)
                        ttGATAATAAAGTT‑AAAAAAACAAATAAAGACAAGACAATTGGGGCCAttttgactcatagaggaaag                                        <  1:3606176‑M1/68‑21 (MQ=255)
                          gATAATAAATTT‑AAAAAAACAAATAAATACAAGACAATTGGGGCCAAACTGTCCATATCATAAATAAg                                      <  1:1904747/68‑1 (MQ=255)
                          gATAATAAAGTT‑AAAAAAACAAATAAATACAAGACAATTGGGGCCAAACTGTCCATATCATAAATAAg                                      >  1:1559030/1‑68 (MQ=255)
                           ataataAAGTT‑AAAAAAACAAATAAATACAAGACAATTGGGGCCAAACTGTCCATATCATAAATAAGt                                     >  1:4046401/1‑68 (MQ=255)
                              ataaAGTT‑AAAAAAACAAATAAATACAAGACAATTGGGGCCAAACTGTCCATATCATAAATAAGTTAc                                  >  1:844473/1‑68 (MQ=255)
                                 aaGTT‑AAAAAAACAAATAAATACAAGACAATTGGGGCCAttttgactcatagaggaaagcatcgcgga                               <  1:3119207‑M1/68‑30 (MQ=255)
                                  aGTT‑AAAAAAACAAATAAATACAAGACAATTGGGGCCAAACTGTCCATATCATAAATAAGTTACGTAt                              >  1:1973530/1‑68 (MQ=255)
                                          aaaaCAAATAAATACAAGACAATTGGGGCCAttttgactcatagaggaaagcatcgcggacaaacttt                       >  1:1106251‑M1/1‑31 (MQ=255)
                                                aataaataCAAGACAATTGGGGCCAttttgactcatagaggaaagcatcgcggacaaactttttcagt                 <  1:1081972‑M1/68‑44 (MQ=255)
                                                     ataCAAGACAATTGGGGCCAAACTGTCCATATCATAAATAAGTTACGTATTTTTTCTCAAGCATaaaa            >  1:1608804/1‑68 (MQ=255)
                                                       aCAAGACAATTGGGGCCAAACTGTCCATATCATAAATAAGTTACGTATTTTTTCTCAAGCATAAAAat          >  1:74884/1‑68 (MQ=255)
                                                       aCAAGACAATTGGGGCCAAACTGTCCATATCATAAATAAGTTACGTATTTTTTCTCAAGCATAAAAat          <  1:2000851/68‑1 (MQ=255)
                                                            aCAATTGGGGCCAttttgactcatagaggaaagcatcgcggacaaactttttcagtttatttgttggc     <  1:1880005‑M1/68‑56 (MQ=255)
                                                            aCAATTGGGGCCAttttgactcatagaggaaagcatcgcggacaaactttttcagtttatttgttggc     <  1:2736274‑M1/68‑56 (MQ=255)
                                                            aCAATTGGGGCCAttttgactcatagaggaaagcatcgcggacaaactttttcagtttatttgttggc     <  1:1776621‑M1/68‑56 (MQ=255)
                                                            aCAATTGGGGCCAttttgactcatagaggaaagcatcgcggacaaactttttcagtttatttgttggc     <  1:1369200‑M1/68‑56 (MQ=255)
                                                            aCAATTGGGGCCAAACTGTCCATATCATAAATAAGTTACGTATTTTTTCTCAAGCATAAAAATATTaa     >  1:2967156/1‑68 (MQ=255)
                                                            aCAATTGGGGCCAAACTGTCCATATCATAAATAAGTTACGTATTTTTTCTCAAGCATAAAAATATTaa     >  1:195149/1‑68 (MQ=255)
                                                               aTTGGGGCCAAACTGTCCATATCATATATAAGTTACGTATTTTTTCTCAAGCATAAAAATATTaaaaa  >  1:3905730/1‑68 (MQ=255)
                                                               aTTGGGGCCAAACTGTCCATATCATAAATAAGTTACGTATTTTTTCTCAAGCATAAAAATATTaaaaa  >  1:2845174/1‑68 (MQ=255)
                                                               aTTGGGGCCAAACTGTCCATATCATAAATAAGTTACGTATTTTTTCTCAAGCATAAAAATATTaaaaa  >  1:3934951/1‑68 (MQ=255)
                                                                |                                                                  
AGAGAAGAGGTTTGATTTAACTTATTGATAATAAAGTT‑AAAAAAACAAATAAATACAAGACAATTGGGGCCAAACTGTCCATATCATAAATAAGTTACGTATTTTTTCTCAAGCATAAAAATATTAAAAA  >  NC_000913/4542619‑4542748

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