New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913 608007 =NA (NA)7 (0.270) 5/224 2.3 15.6% noncoding (1/1345 nt) IS186 repeat region
?NC_000913 609352 = 38 (1.480)intergenic (+175/+107) insL‑2/entD IS186/IS421 transposase/phosphopantetheinyl transferase EntD
Rejected: Frequency below/above cutoff threshold.

GTGTTTAACGGTGCTGGGATTTACGCTTCTGGTCGGAGACTCGCTGTGTGAGTTTACGGTGAAGGAACGTAATATTGAGTTTAAGGCTGTTCTCGCTTACGAACCGAAGAAGTAGCCGTTGTGCGGGGAGTAATCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTACTCTCACGATAACTGGTCAGCAATTCTGGCCCATATTGGTAAGCCCGAAGAACTGGATACTTC  >  NC_000913/607874‑608135
                                                                                                                                     |                                                                                                                                
gtgtTTAACGGTGCTGGGATTTACGCTTCTGGTCGGAGACTCGCTGTGTGAGTTTACGGTGAAGGAACGTAATATTGAGTTTAAGGCTGTTCTCGCTTACGAACCGAAGAAGTAGCCGTTGTGCGGGGAGTAATCCCATa                                                                                                                            <  1:155198/140‑1 (MQ=255)
    ttAACGGTGCTGGGATTTACGCTTCTGGTCGGAGACTCGCTGTGTGAGTTTACGGTGAAGGAACGTAATATTGAGTTTAAGGCTGTTCTCGCTTACGAACCGAAGAAGTAGCCGTTGTGCGGGGATTACTCCCATAAgcg                                                                                                                        <  1:189080/140‑1 (MQ=255)
     tAACGGTGCTGGGATTTACGCTTCTGGTCGGAGACTCGCTGTGTGAGTTTACGGTGAAGGAACGTAATATTGAGTTTAAGGCTGTTCTCGCTTACGAACCGAAGAAGTAGCCGTTGTGCGGGGAGTAATCCCATAAgcgc                                                                                                                       <  1:8833/140‑1 (MQ=255)
      aaCGGTGCTGGGATTTACGCTTCTGGTCGGAGACTCGCTGTGTGAGTTTACGGTGAAGGAACGTAATATTGAGTTTAAGGCTGTTCTCGCTTACGAACCGAAGAAGTAGCCGTTGTGCGGGGAGTAATCCCATAAGCGCt                                                                                                                      <  1:98915/140‑1 (MQ=255)
           tGCTGGGATTTACGCTTCTGGTCGGAGACTCGCTGTGTGAGTTTACGGTGAAGGAACGTAATATTGAGTTTAAGGCTGTTCTCGCTTACGAACCGAAGAAGTAGCCGTTGTGCGGGGAGTAATCCCATAAGCGCTAACtt                                                                                                                 >  2:10643/1‑140 (MQ=255)
           tGCTGGGATTTACGCTTCTGGTCGGAGACTCGCTGTGTGAGTTTACGGTGAAGGAACGTAATATTGAGTTTAAGGCTGTTCTCGCTTACGAACCGAAGAAGTAGCCGTTGTGCGGGGAGTAATCCCATAAGCGCTAACtt                                                                                                                 >  2:50186/1‑140 (MQ=255)
           tGCTGGGATTTACGCTTCTGGTCGGAGACTCGCTGTGTGAGTTTACGGTGAAGGAACGTAATATTGAGTTTAAGGCTGTTCTCGCTTACGAACCGAAGAAGTAGCCGTTGTGCGGGGAGTAATCCCATAAGCGCTAACtt                                                                                                                 >  2:384342/1‑140 (MQ=255)
           tGCTGGGATTTACGCTTCTGGTCGGAGACTCGCTGTGTGAGTTTACGGTGAAGGAACGTAATATTGAGTTTAAGGCTGTTCTCGCTTACGAACCGAAGAAGTAGCCGTTGTGCGGGGAGTAATCCCATAAGCGCTAACtt                                                                                                                 <  1:50186/140‑1 (MQ=255)
             cTGGGATTTACGCTTCTGGTCGGAGACTCGCTGTGTGAGTTTACGGTGAAGGAACGTAATATTGAGTTTAAGGCTGTTCTCGCTTACGAACCGAAGAAGTAGCCGTTGTGCGGGGAGTAATCCCATAAGCGCTAACTTaa                                                                                                               >  1:162689/1‑140 (MQ=255)
                ggATTTACGCTTCTGGTCGGAGACTCGCTGTGTGAGTTTACGGTGAAGGAACGTAATATTGAGTTTAAGGCTGTTCTCGCTTACGAACCGAAGAAGTAGCCGTTGTGCGGGGAGTAATCCCATAAGCGCTAACTTAAggg                                                                                                            <  1:10643/140‑1 (MQ=255)
                   tttACGCTTCTGGTCGGAGACTCGCTGTGTGAGTTTACGGTGAAGGAACGTAATATTGAGTTTAAGGCTGTTCTCGCTTACGAACCGAAGAAGTAGCCGTTGTGCGGGGAGTAATCCCATAAGCGCTAACTTAAGGGTtg                                                                                                         <  1:394981/140‑1 (MQ=255)
                   tttACGCTTCTGGTCGGAGACTCGCTGTGTGAGTTTACGGTGAAGGAACGTAATATTGAGTTTAAGGCTGTTCTCGCTTACGAACCGAAGAAGTAGCCGTTGTGCGGGGAGTAATCCCATAAGCGCTAACTTAAGGGTtg                                                                                                         <  2:289232/140‑1 (MQ=255)
                           tCTGGTCGGAGACTCGCTGTGTGAGTTTACGGTGAAGGAACGTAATATTGAGTTTAAGGCTGTTCTCGCTTACGAACCGAAGAAGTAGCCGTTGTGCGGGGAGTAATCCCATAAGCGCTAACTTAAGGGTTGTGGTATTa                                                                                                 >  2:237960/1‑140 (MQ=255)
                               gTCGGAGACTCGCTGTGTGAGTTTACGGTGAAGGAACGTAATATTGAGTTTAAGGCTGTTCTCGCTTACGAACCGAAGAAGTAGCCGTTGTGCGGGGAGTAATCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGcc                                                                                             <  1:371387/140‑1 (MQ=255)
                               gTCGGAGACTCGCTGTGTGAGTTTACGGTGAAGGAACGTAATATTGAGTTTAAGGCTGTTCTCGCTTACGAACCGAAGAAGTAGCCGTTGTGCGGGGAGTAATCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGcc                                                                                             >  2:28637/1‑140 (MQ=255)
                                                             aaGGAACGTAATATTGAGTTTAAGGCTGTTCTCGCTTACGAACCGAAGAAGTAGCCGTTGTGCGGGGAGTAATCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTAct                                                               <  2:111743/140‑1 (MQ=255)
                                                                             aGTTTAAGGCTGTTCTCGCTTACGAACCGAAGAAGTAGCCGTTGTGCGGGGAGTAATCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTACTCTCACGATAACTGGTc                                               >  2:412812/1‑140 (MQ=255)
                                                                               tttAAGGCTGTTCTCGCTTACGAACCGAAGAAGTAGCCGTTGTGCGGGGAGTAATCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTACTCTCACGATAACTGGTCAg                                             >  2:187051/1‑140 (MQ=255)
                                                                               tttAAGGCTGTTCTCGCTTACGAACCGAAGAAGTAGCCGTTGTGCGGGGAGTAATCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTACTCTCACGATAACTGGTCAg                                             >  2:304929/1‑140 (MQ=255)
                                                                                  aaGGCTGTTCTCGCTTACGAACCGAAGAAGTAGCCGTTGTGCGGGGAGTAATCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTACTCTCACGATAACTGGTCAg                                             <  2:186490/137‑1 (MQ=255)
                                                                                  aaGGCTGTTCTCGCTTACGAACCGAAGAAGTAGCCGTTGTGCGGGGAGTAATCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTACTCTCACGATAACTGGTCAg                                             >  1:186490/1‑137 (MQ=255)
                                                                                  aaGGCTGTTCTCGCTTACGAACCGAAGAAGTAGCCGTTGTGCGGGGAGTAATCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTACTCTCACGATAACTGGTCAGCaa                                          >  2:168550/1‑140 (MQ=255)
                                                                                         ttCTCGCTTACGAACCGAAGAAGTAGCCGTTGTGCGGGGAGTAATCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTACTCTCACGATAACTGGTCAGCAATTCTGGc                                   <  1:384342/140‑1 (MQ=255)
                                                                                         ttCTCGCTTACGAACCGAAGAAGTAGCCGTTGTGCGGGGAGTAATCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTACTCTCACGATAACTGGTCAGCAATTCTGGc                                   <  1:28637/140‑1 (MQ=255)
                                                                                           ctcGCTTACGAACCGAAGAAGTAGCCGTTGTGCGGGGAGTAATCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTACTCTCACGATAACTGGTCAGCAATTCTGGccc                                 <  2:78639/140‑1 (MQ=255)
                                                                                              gCTTACGAACCGAAGAAGTAGCCGTTGTGCGGGGAGTAATCCCATAAGCGCTATCTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTACTCTCACGATAACTGGTCAGCAATTCTGGCCCata                              >  1:487747/1‑140 (MQ=255)
                                                                                                   cGAACCGAAGAAGTAGCCGTTGTGCGGGGAGTAATCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTACTCTCACGATAACTGGTCAGCAATTCTGGc                                   >  1:85967/1‑130 (MQ=255)
                                                                                                   cGAACCGAAGAAGTAGCCGTTGTGCGGGGAGTAATCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTACTCTCACGATAACTGGTCAGCAATTCTGGc                                   <  2:85967/130‑1 (MQ=255)
                                                                                                               gTAGCCGTTGTGCGGGGAGTAATCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTACTCTCACGATAACTGGTCAGCAATTCTGGCCCATATTGGTAAGCCCgaagaa             <  1:166461/140‑1 (MQ=255)
                                                                                                                          gCGGGGAGTAATCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTACTCTCACGATAACTGGTCAGCAATTCTGGCCCATATTGGTAAGCCCGAAGAACTGGATACTTc  >  2:408908/1‑140 (MQ=37)
                                                                                                                                     |                                                                                                                                
GTGTTTAACGGTGCTGGGATTTACGCTTCTGGTCGGAGACTCGCTGTGTGAGTTTACGGTGAAGGAACGTAATATTGAGTTTAAGGCTGTTCTCGCTTACGAACCGAAGAAGTAGCCGTTGTGCGGGGAGTAATCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTACTCTCACGATAACTGGTCAGCAATTCTGGCCCATATTGGTAAGCCCGAAGAACTGGATACTTC  >  NC_000913/607874‑608135

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