New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913 4370205 =13 (0.440)4 (0.140) 4/280 2.3 18.1% coding (208/1257 nt) yjeH putative transporter
?NC_000913 4370316 = 24 (0.860)coding (97/1257 nt) yjeH putative transporter
Rejected: Frequency below cutoff threshold.

GAACCAAACGCCATACCGACGAAGTGCGCGACGCCGCCTGCGCTGGGATAGTGGCGACCCAGAATCGCAAACACAATCGCAATCGGGAACACTAAGATAATCAAAACGGGCCACGCCCACAGGCTGTTATTGCCCGCTACCAGCGCAGCTAACGCAGGAACGGCAAACACGCCAGTGCCTAATAATGACGTCGATAGCAGGCCAATGCCCTGGGCCAGCCCCAGTTCTTGTTTGAGTCCACTCATGGGTTGATGTCCGATTGCGCCCAAATTTTGGGCAACTGCGTA  >  NC_000913/4370168‑4370454
                                                                                                                                                    |                                                                                                                                          
ttattaggcactggcgtgtttgccgttcctgcgttagCTGCGCTGGGATAGTGGCGACCCAGAATCGCAAACACAATCGCAATCGGGAACACTAAGATAATCAAAACGGGCCACGCCCACAGGCTGTTATTGCCCGCTACCAGCGCAGCTa                                                                                                                                          <  1:145285‑M1/114‑1 (MQ=255)
 aaCCAAACGCCATACCGACGAAGTGCGCGACGCCGCCTGCGCTGGGATAGTGGCGACCCAGAATCGCAAACACAATCGCAATCGGGGACACTAAGATAATCAAAACGGGCCACGCCCACAGGCTGTTATTGCCCGCTACCAGCGCAGCTa                                                                                                                                          <  2:159895/150‑1 (MQ=255)
    taggcactggcgtgtttgccgttcctgcgttagctgcgctggtagcgggcaataacagcctgtgggcgtggcccgttttgattatcttagtgttcccgattgcgattgtgtttgcgattctgggtcgccactatcccagcgcagCTAACGc                                                                                                                                      <  2:159025‑M2/7‑1 (MQ=255)
     aaaCGCCATACCGACGAAGTGCGCGACGCCGCCTGCGCTGGGATAGTGGCGACCCAGAATCGCAAACACAATCGCAATCGGGAACACTAAGATAATCAAAACGGGCCACGCCCACAGGCTGTTATTGCCCGCTACCAGCGCAGCTAACGCa                                                                                                                                     <  2:385636/151‑1 (MQ=255)
        cGCCATACCGACGAAGTGCGCGACGCCGCCTGCGCTGGGATAGTGGCGACCCAGAATCGCAAACACAATCGCAATCGGGAACACTAAGATAATCAAAACGGGCCACGCCCACAGGCTGTTATTGCCCGCTACCAGCGCAGCTAACGCAGGa                                                                                                                                  >  2:402174/1‑151 (MQ=255)
            ggcgtgtttgccgttcctgcgttagCTGCGCTGGGATAGTGGCGACCCAGAATCGCAAACACAATCGCAATCGGGAACACTAAGATAATCAAAACGGGCCACGCCCACAGGCTGTTATTGCCCGCTACCAGCGCAGCTAACGCAGGAACgg                                                                                                                              <  1:116184‑M1/126‑1 (MQ=255)
                      ccgttcctgcgttagctgcgctggtagcgggcaataacagcctgtgggcgtggcccgttttgattatcttagtgttcccgattgcgattgtgtttgcgattctgggtcgccactatcccagcgcagCTAACGCAGGAACGGCAAACACGcc                                                                                                                    <  1:116184‑M2/25‑1 (MQ=255)
                       gTGCGCGACGCCGCCTGCGCTGGGATAGTGGCGACCCAGAATCGCAAACACAATCGCAATCGGGAACACTAAGATAATCAAAACGGGCCACGCCCACAGGCTGTTATTGCCCGCTACCAGCGCAGCTAACGCAGGAACGGCAAACACGCCa                                                                                                                   >  1:142890/1‑151 (MQ=255)
                              gcgttagCTGCGCTGGGATAGTGGCGACCCAGAATCGCAAACACAATCGCAATCGGGAACACTAAGATAATCAAAACGGGCCACGCCCACAGGCTGTTATTGCCCGCTACCAGCGCAGCTAACGCAGGAACGGCAAACACGCCAGTGCCta                                                                                                            <  2:159025‑M1/144‑1 (MQ=255)
                                 ccgccTGCGCTGGGATAGTGGCGACCCAGAATCGCAAACACAATCGCAATCGGGAACACTAAGATAATCAAAACGAGCCACGCCCACAGGCTGTTATTGCCCGCTACCAGCGCAGCTAACGCAGGAACGGCAAACACGCCAGTGCCtaat                                                                                                          <  1:449748/150‑1 (MQ=255)
                                  tagctgcgctggtagcgggcaataacagcctgtgggcgtggcccgttttgattatcttagtgttcccgattgcgattgtgtttgcgattctgggtcgccactatcccagcgcagCTAACGCAGGAACGGCAAACACGCCAGTGCCtaataa                                                                                                        <  1:145285‑M2/37‑1 (MQ=255)
                                  cgccTGCGCTGGGATAGTGGCGACCCAGAATCGCAAACACAATCGCAATCGGGAACACTAAGATAATCAAAACGGGCCACGCCCACAGGCTGTTATTGCCCGCTACCAGCGCAGCTAACGCAGGAACGGCAAACACGCCAGTGCCtaataa                                                                                                        <  2:41741/151‑1 (MQ=255)
                                       GCGCTGGGATAGTGGCGACCCAGAATCGCAAACACAATCGCAATCGGGAACACTAAGATAATCAAAACGGGCCACGCCCACAGGCTGTTATTGCCCGCTACCAGCGCAGCTAACGCAGGAACGGCAAACACGCCAGTGCCTAATAATGACg                                                                                                   <  2:357234‑M1/151‑1 (MQ=255)
                                             gtagcgggaaataacagcctgtgggggtggcccgttttgattatcttagtgtccccgattgcaattgtgtttgcgattctgggtcgccactatcccagcgcagCTAACGCAGGACCGGCAAACACGCCAGTGCCTAATAATGACGTCGAt                                                                                              <  1:3246‑M2/47‑1 (MQ=255)
                                               aTAGTGGCGACCCAGAATCGCAAACACAATCGCAATCGGGAACACTAAGATAATCAAAACGGGCCACGCCCACAGGCTGTTATTGCCCGCTACCAGCGCAGCTAACGCAGGAACGGCAAACACGCCAGTGCCTAATAATGACGTCGATAGc                                                                                           <  2:76484/151‑1 (MQ=255)
                                               aTAGTGGCGACCCAGAATCGCAAACACAATCGCAATCGGGAACACTAAGATAATCAAAACGGGCCACGCCCACAGGCTGTTATTGCCCGCTACCAGCGCAGCTAACGCAGGAACGGCAAACACGCCAGTGCCTAATAATGACGTCGATAGc                                                                                           <  2:366869/151‑1 (MQ=255)
                                                        aCCCAGAATCGCAAACACAATCGCAATCGGGAACACTAAGATAATCAAAACGGGCCACGCCCACAGGCTGTTATTGCCCGCTACCAGCGCAGCTAACGCAGGAACGGCAAACACGCCAGTGCCTAATAATGACGTCGATAGCAGGCCAATg                                                                                  >  2:73506/1‑151 (MQ=255)
                                                        aCCCAGAATCGCAAACACAATCGCAATCGGGAACACTAAGATAATCAAAACGGGCCACGCCCACAGGCTGTTATTGCCCGCTACCAGCGCAGCTAACGCAGGAACGGCAAACACGCCAGTGCCTAATAATGACGTCGATAGCAGGCCAATg                                                                                  >  2:338908/1‑151 (MQ=255)
                                                         cccAGAATCGCAAACACAATCGCAATCGGGAACACTAAGATAATCAAAACGGGCCACACCCACAGGCTGTTATTGCCCGCTACCAGCGCAGCTAACGCAGGAACGGCAAACACGCCAGTGCCTAATAATGACGTCGATAGCAGGCCAATg                                                                                  >  1:145368/1‑150 (MQ=255)
                                                              aaTCGCAAACACAATCGCAATCGGGAACACTAAGATAATCAAAACGGGCCACGCCCACAGGCTGTTATTGCCCGCTACCAGCGCAGCTAACGCAGGAACGGCAAACACGCCAGTGCCTAATAATGACGTCGATAGCAGGCCAATGccc                                                                               <  1:46477/148‑1 (MQ=255)
                                                                tCGCAAACACAATCGCAATCGGGAACACTAAGATAATCAAAACGGGCCACGCCCACAGGCTGTTATTGCCCGCTACCAGCGCAGCTAACGCAGGAACGGCAAACACGCCAGTGCCTAATAATGACGTCGATAGCAGGCCAATGCCCTGGGc                                                                          <  2:53813/151‑1 (MQ=255)
                                                                 cGCAAACACAATCGCAATCGGGAACACTAAGATAATCAAAACGGGCCACGCCCACAGGCTGTTATTGCCCGCTACCAGCGCAGCTAACGCAGGAACGGCAAACACGCCAGTGCCTAATAATGACGTCGATAGCAGGCCAATGCCCTGGGcc                                                                         <  2:137091/151‑1 (MQ=255)
                                                                       cacaATCGCAATCGGGAACACTAAGATAATCAAAACGGGCCACGCCCACAGGCTGTTATTGCCCGCTACCAGCGCAGCTAACGCAGGAACGGCAAACACGCCAGTGCCTAATAATGACGTCGATAGCAGGCCAATGCCCTGGGCCAGccc                                                                    >  1:461386/1‑150 (MQ=255)
                                                                           atcgcaatcgGGAACACTAAGATAATCAAAACGGGCCACGCCCACAGGCTGTTATTGCCCGCTACCAGCGCAGCTAACGCAGGAACGGCAAACACGCCAGTGCCTAATAATGACGTCGATAGCAGGCCAATGCCCTGGGCCAGCCCCAGtt                                                               >  2:240010/1‑151 (MQ=255)
                                                                                         cacTAAGATAATCAAAACGGGCCACGCCCACAGGCTGTTATTGCCCGCTACCAGCGCAGCTAACGCAGGAACGGCAAACACGCCAGTGCCTAATAATGACGTCGATAGCAGGCCAATGCCCTGGGCCAGCCCCAGTTCTTGTTTGAGTCCa                                                 >  2:172243/1‑151 (MQ=255)
                                                                                         cacTAAGATAATCAAAACGGGCCACGCCCACAGGCTGTTATTGCCCGCTACCAGCGCAGCTAACGCAGGAACGGCAAACACGCCAGTGCCTAATAATGACGTCGATAGCAGGCCAATGCCCTGGGCCAGCCCCAGTTCTTGTTTGAGTCCa                                                 >  1:241811/1‑151 (MQ=255)
                                                                                         cacTAAGATAATCAAAACGGGCCACGCCCACAGGCTGTTATTGCCCGCTACCAGCGCAGCTAACGCAGGAACGGCAAACACGCCAGTGCCTAATAATGACGTCGATAGCAGGCCAATGCCCTGGGCCAGCCCCAGTTCTTGTTTGAGTCCa                                                 >  2:377988/1‑151 (MQ=255)
                                                                                         cacTAAGATAATCAAAACGGGCCACGCCCACAGGCTGTTATTGCCCGCTACCAGCGCAGCTAACGCAGGAACGGCAAACACGCCAGTGCCTAATAATGACGTCGATAGCAGGCCAATGCCCTGGGCCAGCCCCAGTTCTTGTTTGAGTCCa                                                 >  2:42440/1‑151 (MQ=255)
                                                                                                  aaTCAAAACGGGCCACGCCCACAGGCTGTTATTGCCCGCTACCAGCGCAGCTAACGCAGGAACGGCAAACACGCCAGTGCCTAATAATGACGTCGATAGCAGGCCAATGCCCTGGGCCAGCCCCAGTTCTTGTTTGAGTCCACTCATGGGt                                        <  1:171121/151‑1 (MQ=255)
                                                                                                         aCGGGCCACGCCCACAGGCTGTTATTGCCCGCTACCAGCGCAGCTAACGCAGGAACGGCAAACACGCCAGTGCCTAATAATGACGTCGATAGCAGGCCAATGCCCTGGGCCAGCCCCAGTTCTTGTTTGAGTCCACTCATGGGTTGATGTc                                 <  1:166518/151‑1 (MQ=255)
                                                                                                                   cccACAGGCTGTTATTGCCCGCTACCAGCGCAGCTAACGCAGGAACGGCAAACACGCCAGTGCCTAATAATGACGTCGATAGCAGGCCAATGCCCTGGGCCAGCCCCAGTTCTTGTTTGAGTCCACTCATGGGTTGATGTCCGATTgcgc                        >  1:21479/1‑150 (MQ=255)
                                                                                                                           cTGTTATTGCCCGCTACCAGCGCAGCTAACGCAGGAACGGCAAACACGCCAGTGCCTAATAATGACGTCGATAGCAGGCCAATGCCCTGGGCCAGCCCCAGTTCTTGTTTGAGTCCACTCATGGGTTGATGTCCGATTGCGCCCAAAtttt               >  2:158755/1‑151 (MQ=255)
                                                                                                                                        ctataagagcAGCTAATACACGAACGGCTAACCTGCGAGAGCCTAATTGTGAAGTCGATGGAAGAACAATGCGCTGGGCCGGCCCACAATCTTGCGTGAGTCTACGCATGGGATGATGTCAGAATGCGCCCAAATTTTGGGCAACTGCGTa  >  1:116978/9‑151 (MQ=255)
                                                                                                                                                    |                                                                                                                                          
GAACCAAACGCCATACCGACGAAGTGCGCGACGCCGCCTGCGCTGGGATAGTGGCGACCCAGAATCGCAAACACAATCGCAATCGGGAACACTAAGATAATCAAAACGGGCCACGCCCACAGGCTGTTATTGCCCGCTACCAGCGCAGCTAACGCAGGAACGGCAAACACGCCAGTGCCTAATAATGACGTCGATAGCAGGCCAATGCCCTGGGCCAGCCCCAGTTCTTGTTTGAGTCCACTCATGGGTTGATGTCCGATTGCGCCCAAATTTTGGGCAACTGCGTA  >  NC_000913/4370168‑4370454

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