New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NZ_CP009273 2442890 =27 (0.900)3 (0.100) 3/162 2.1 10.2% coding (639/822 nt) yfcO fimbrial adhesin YfcO
?NZ_CP009273 2442931 = 27 (0.930)coding (598/822 nt) yfcO fimbrial adhesin YfcO
Rejected: Frequency below/above cutoff threshold.

TGTCCGTTATAGGGTGTCGATTTGGCGTGACCGAGACTTCCGCCTTGTCACATCCCGACGCGGTAAAAGTTATCATCTGTGCAACCCAGCGGCAGTAACCACCTGTTTCATTCATCTGCGTACTGTCGCACAACTTCCCATCTTTCATCAAATACGTTGAAAGCGTTTTCTCAACCAC  >  NZ_CP009273/2442802‑2442979
                                                                                        |                                                                                         
tGTCCGTTATAGGGTGTCGATTTGGCGTGACCGAGACTTCCGCCTTGTCACATCCCGACGCGGTAAAAGTTATCATCTGTGCAACCCAGc                                                                                          >  2:160739/1‑90 (MQ=255)
tGTCCGTTATAGGGTGTCGATTTGGCGTGACCGAGACTTCCGCCTTGTCACATCCCGACGCGGTAAAAGTTATCATCTGTGCAACCCAGc                                                                                          <  1:823694/90‑1 (MQ=255)
agtaagtcaaaactgctcaccctgcagtcgacgggcggtgtggttgagaaaacgctttcaacgtatttgatgaaagatgggaagttgtGc                                                                                          <  2:696041‑M1/2‑1 (MQ=255)
agtaagtcaaaactgctcaccctgcagtcgacgggcggtgtggttgagaaaacgctttcaacgtatttgatgaaagatgggaagttgtGc                                                                                          <  1:764763‑M1/2‑1 (MQ=255)
agtaagtcaaaactgctcaccctgcagtcgacgggcggtgtggttgagaaaacgctttcaacgtatttgatgaaagatgggaagttgtGc                                                                                          >  2:385537‑M1/89‑90 (MQ=255)
agtaagtcaaaactgctcaccctgcagtcgacgggcggtgtggttgagaaaacgctttcaacgtatttgatgaaagatgggaagttgtGc                                                                                          <  1:55333‑M1/2‑1 (MQ=255)
agtaagtcaaaactgctcaccctgcagtcgacgggcggtgtggttgagaaaacgctttcaacgtatttgatgaaagatgggaagttgtGc                                                                                          <  1:60211‑M1/2‑1 (MQ=255)
     gTTATAGGGTGTCGATTTGGCGTGACCGAGACTTCCGCCTTGTCACATCCCGACGCGGTAAAAGTTATCATCTGTGCAACCCAGCGGCAg                                                                                     >  1:839883/1‑90 (MQ=255)
         aaactgctcaccctgcagtcgacgggcggtgtggttgagaaaacgctttcaacgtatttgatgaaagatgggaagttgtGCGGCAgtacg                                                                                 <  2:410725‑M1/11‑3 (MQ=255)
              tgtCGATTTGGCGTGACCGAGACTTCCGCCTTGTCACATCCCGACGCGGTAAAAGTTATCATCTGTGCAACCCAGCGGCAGTAACCACCt                                                                            >  2:473532/1‑90 (MQ=255)
              tgtCGATTTGGCGTGACCGAGACTTCCGCCTTGTCACATCCCGACGCGGTAAAAGTTATCATCTGTGCAACCCAGCGGCAGTAACCACCt                                                                            >  2:140621/1‑90 (MQ=255)
              tgtCGATTTGGCGTGACCGAGACTTCCGCCTTGTCACATCCCGACGCGGTAAAAGTTATCATCTGTGCAACCCAGCGGCAGTAACCACCt                                                                            >  2:748930/1‑90 (MQ=255)
               gtCGATTTGGCGTGACCGAGACTTCCGCCTTGTCACATCCCGACGCGGTAAAAGTTATCATCTGTGCAACCCAGCGGCAGTAACCACCTg                                                                           >  1:592488/1‑90 (MQ=255)
               gtCGATTTGGCGTGACCGAGACTTCCGCCTTGTCACATCCCGACGCGGTAAAAGTTATCATCTGTGCAACCCAGCGGCAGTAACCACCTg                                                                           >  1:800859/1‑90 (MQ=255)
               ctcaccctgcagtcgacgggcggtgtggttgagaaaacgctttcaacgtatttgatgaaagatgggaagttgtGc                                                                                          <  2:246277‑M1/2‑1 (MQ=255)
               ctcaccctgcagtcgacgggcggtgtggttgagaaaacgctttcaacgtatttgatgaaagatgggaagttgtGc                                                                                          >  1:246277‑M1/74‑75 (MQ=255)
                tCGATTTGGCGTGACCGAGACTTCCGCCTTGTCACATCCCGACGCGGTAAAAGTTATCATCTGTGCAACCCAGCGGCAGTAACCACCTGt                                                                          <  2:592709/90‑1 (MQ=255)
                    tttGGCGTGACCGAGACTTCCGCCTTGTCACATCCCGACGCGGTAAAAGTTATCATCTGTGCAACCCAGCGGCAGTAACCACCTGTttca                                                                      >  1:408002/1‑90 (MQ=255)
                       ggCGTGACCGAGACTTCCGCCTTGTCACATCCCGACGCGGTAAAAGTTATCATCTGTGCAACCCAGCGGCAGTAACCACCTGTttcattc                                                                   <  1:140621/90‑1 (MQ=255)
                            gACCGAGACTTCCGCCTTGTCACATCCCGACGCGGTAAAAGTTATCATCTGTGCAACCCAGCGGCAGTAACCACCTGTTTCATTCATCTg                                                              >  2:613248/1‑90 (MQ=255)
                            gACCGAGACTTCCGCCTTGTCACATCCCGACGCGGTAAAAGTTATCATCTGTGCAACCCAGCGGCAGTAACCACCTGTTTCATTCATCTg                                                              >  2:731335/1‑90 (MQ=255)
                                       ccGCCTTGTCACATCCCGACGCGGTAAAAGTTATCATCTGTGCAACCCAGCGGCAGTAACCACCTGTTTCATTCATCTGCGTACTGTCGc                                                   >  2:858240/1‑90 (MQ=255)
                                       ccGCCTTGTCACATCCCGACGCGGTAAAAGTTATCATCTGTGCAACCCAGCGGCAGTAACCACCTGTTTCATTCATCTGCGTACTGTCGc                                                   <  1:829921/90‑1 (MQ=255)
                                           cTTGTCACATCCCGACGCGGTAAAAGTTATCATCTGTGCAACCCAGCGGCAGTAACCACCTGTTTCATTCATCTGCGTACTGTCGCACaa                                               <  2:839883/90‑1 (MQ=255)
                                               gaaaacgctttcaacgtatttgatgaaagatgggaagttgtGc                                                                                          <  2:594158‑M1/2‑1 (MQ=255)
                                               gaaaacgctttcaacgtatttgatgaaagatgggaagttgtGc                                                                                          >  1:594158‑M1/42‑43 (MQ=255)
                                                   acgctttcaacgtatttgatgaaagatgggaagttgtGc                                                                                          <  1:71133‑M1/2‑1 (MQ=255)
                                                   acgctttcaacgtatttgatgaaagatgggaagttgtGc                                                                                          >  2:71133‑M1/38‑39 (MQ=255)
                                                        gACGCGGTAAAAGTTATCATCTGTGCAACCCAGCGGCAGTAACCACCTGTTTCATTCATCTGCGTACTGTCGCACAACTTCCCATCTTtc                                  >  1:22374/1‑90 (MQ=255)
                                                            cgGTAAAAGTTATCATCTGTGCAACCCAGCGGCAGTAACCACCTGTTTCATTCATCTGCGTACTGTCGCACAACTTCCCATCTTtcatca                              >  2:699742/1‑90 (MQ=255)
                                                                atttgatgaaagatgggaagttgtGCGGCAGTAACCACCTGTttcattc                                                                   <  1:345912‑M1/25‑1 (MQ=255)
                                                                atttgatgaaagatgggaagttgtGCGGCAGTAACCACCTGTttcattc                                                                   >  2:345912‑M1/25‑49 (MQ=255)
                                                                      tATCATCTGTGCAACCCAGCGGCAGTAACCACCTGTTTCATTCATc                                                                <  2:509722/46‑1 (MQ=255)
                                                                      tATCATCTGTGCAACCCAGCGGCAGTAACCACCTGTTTCATTCATc                                                                >  1:509722/1‑46 (MQ=255)
                                                                      tATCATCTGTGCAACCCAGCGGCAGTAACCACCTGTTTCATTCATCTGCGTACTGTCGCACAACTTCCCATCTTTCATCAAATACGTTGa                    >  1:98988/1‑90 (MQ=255)
                                                                             tgtgCAACCCAGCGGCAGTAACCACCTGTTTCATTCATCTGCGTACTGTCGCACAACTTCCCATCTTTCATCAAATACGTTGAAAGCGtt             >  1:358180/1‑90 (MQ=255)
                                                                              gtgCAACCCAGCGGCAGTAACCACCTGTTTCATTCATCTGCGTACTGTCGCACAACTTCCCATCTTTCATCAAATACGTTGAAAGCGttt            >  2:618495/1‑90 (MQ=255)
                                                                              gtgCAACCCAGCGGCAGTAACCACCTGTTTCATTCATCTGCGTACTGTCGCACAACTTCCCATCTTTCATCAAATACGTTGAAAGCGttt            >  1:619256/1‑90 (MQ=255)
                                                                                   aCCCAGCGGCAGTAACCACCTGTTTCATTCATCTGCGTACTGTCGCACAACTTCCCATCttt                                   <  2:413050/62‑1 (MQ=255)
                                                                                   aCCCAGCGGCAGTAACCACCTGTTTCATTCATCTGCGTACTGTCGCACAACTTCCCATCttt                                   >  1:413050/1‑62 (MQ=255)
                                                                                      cAGCGGCAGTAACCACCTGTTTCATTCATCTGCGTACTGTCGCACAACTTCCCATCTTTCATCAAATACGTTGAAAGCGTTTTCTCAAcc    >  2:589906/1‑90 (MQ=255)
                                                                                        GCGGCAGTAACCACCTGTTTCATTCATCTGCGTACTGTCGCACAACTTCCCATCTTTCATCAAATACGTTGAAAGCGTTTTCTCAACcac  <  2:418950‑M1/90‑1 (MQ=255)
                                                                                        GCGGCAGTAACCACCTGTTTCATTCATCTGCGTACTGTCGCACAACTTCCCATCTTTCATCAAATACGTTGAAAGCGTTTTCTCAACcac  <  1:246976‑M1/90‑1 (MQ=255)
                                                                                        |                                                                                         
TGTCCGTTATAGGGTGTCGATTTGGCGTGACCGAGACTTCCGCCTTGTCACATCCCGACGCGGTAAAAGTTATCATCTGTGCAACCCAGCGGCAGTAACCACCTGTTTCATTCATCTGCGTACTGTCGCACAACTTCCCATCTTTCATCAAATACGTTGAAAGCGTTTTCTCAACCAC  >  NZ_CP009273/2442802‑2442979

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