New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913_3_pae_pgi = 337626023 (0.780)4 (0.150) 4/258 3.4 16.7% intergenic (+24/+19) dcuD/sspB putative transporter/ClpXP protease specificity enhancing factor
?NC_000913_3_pae_pgi = 3376265 19 (0.710)intergenic (+29/+14) dcuD/sspB putative transporter/ClpXP protease specificity enhancing factor
Rejected: Coverage evenness skew score above cutoff.
Rejected: Frequency below/above cutoff threshold.

GTGGTGGTCGTTACGTCAGGAATTGCAGGCGTTTCGCCTTTTGCGGTGGTGAAGCGGACAGCGATCCCCATGGCAGTCGGTTTCGTGGTGAATATGATTGCCACAATCACGCTATTTTATTAAGTCATTAAAAAGACAAAACAGGCCGCCTGGGCCTGTTTTGTATTACTTCACAACGCGTAATGCCGGTCGACCACCGCGTGGTGGCTGCGGAGGTTCATCGTCAGGATGAGTGTCATCATCGTGATCTGGCTTGTCGCCATCAATAACCGACATAACGGTTTCGTTGTCT  >  NC_000913_3_pae_pgi/3376114‑3376405
                                                                                                                                                  |                                                                                                                                                 
gtggtggtCGTTACGTCAGGAATTGCAGGCGTTTCGCCTTTTGCGGTGGTGAAGCGGACAGCGATCCCCATGGCAGTCGGTTTCGTGGTGAATATGATTGCCACAATCACGCTATTTTATTAAGTCATTAAAAAGACAAAACAGgccgc                                                                                                                                                 >  1:310863/1‑149 (MQ=255)
    tggtCGTTACGTCAGGAATTGCAGGCGTTTCGCCTTTTGCGGTGGTGAAGCGGACAGCGATCCCCATGGCAGTCGGTTTCGTGGTGAATATGATTGCCACAATCACGCTATTTTATTAAGTCATTAAAAAGACAAAACAGGCCGCCTgg                                                                                                                                             <  1:392615/149‑1 (MQ=255)
      gtCGTTACGTCAGGAATTGCAGGCGTTTCGCCTTTTGCGGTGGTGAAGCGGACAGCGATCCCCATGGCAGTCGGTTTCGTGGTGAATATGATTGCCACAATCACGCTATTTTATTAAGTCATTAAAAAGACAAAACAGGCCGCCTGGGc                                                                                                                                           <  2:310863/149‑1 (MQ=255)
             cGTCCGGACTTGCAGGCGTGTCGGCTGTTGCGGTGGTGAAGCGGACAGCGATCCCCATGGCAGTCGGTTTCGTGGTGAAGATGATTGCCACAGTCACGCTATTTTATTAAGTCATTAAAAAGACAAAACAGGCCGCCTGGGCCTGtttt                                                                                                                                    <  1:293153/149‑1 (MQ=255)
               tCAGGGATTGGAGGCGTCTCCGCTTGTGCGGGGGTGAAGCGGACAGCTGTCCCCATGGCAGTCGGTTTCGTGGTGAATATGATTGCCACAATCACGCTATTTTATTAAGTCATTAAAAAGACAAAACAGGCCCAGGCGGCCTGTTTTGt                                                                                                                                  <  2:40301/149‑1 (MQ=255)
                ccggAATTGCAGGCGTTTCGCCTTTTGCGGTGGTGAAGCGGACAGCGATCCCCATGGCAGTCGGTTTCGTGGTGAATATGATTGCCACAATCACGCTATTTTATTAAGTCATTAAAAAGACAAAACAGGCCGCCTGGGCCTGTTTTGTa                                                                                                                                 <  1:34533/147‑1 (MQ=255)
                       tgcggtcGTTTCGCCTTCTGCGGTGGGGAAGCGGACAGCGATCCCCATGGCAGTCGGTTTCGTGGTGAATATGATTGGCTCAATCATGCTATTTTATTAAGTCATTAAAAAGACAAAACAGGCCGCCTGGGCCTGTTTTGTATTACTTc                                                                                                                          <  2:204865/143‑1 (MQ=255)
                                    ccTTTTGCGGTGGTGAAGCGGACAGCGATCCCCATGGCAGTCGGTTTCGTGGTGAATATGATTGCCACAATCACGCTATTTTATTAAGTCATTAAAAAGACAAAACAGGCCGCCTGGGCCTGTTTTGTATTACTTCACAACGCGTAATg                                                                                                             >  1:515669/1‑149 (MQ=255)
                                              tggtgAAGCGGACAGCGATCCCCATGGCAGTCGGTTTCGTGGTGAATATGATTGCCACAATCACGCTATTTTATTAAGTCATTAAAAAGACAAAACAGGCCGCCTGGGCCTGTTTTGTATTACTTCACAACGCGTAATGCCGGTCGacc                                                                                                   >  1:8392/1‑149 (MQ=255)
                                                           aGCGATCCCCATGGCAGTCGGTTTCGTGGTGAATATGATTGCCACAATCACGCTATTTTATTAAGTCATTAAAAAGACAAAACAGGCCGCCTGGGCCTGTTTTGTATTACTTCACAACGCGTAATGCCGGTCGACCACCGCGTGGTGGc                                                                                      >  1:268997/1‑149 (MQ=255)
                                                                  cccATGGCAGTCGGTTTCGTGGTGAATATGATTGCCACAATCACGCTATTTTATTAAGTCATTAAAAAGACAAAACAGGCCGCCTGGGCCTGTTTTGTATTACTTCACAACGCGTAATGCCGGTCGACCACCGCGTGGTGGCTGCGGAg                                                                               <  1:289863/149‑1 (MQ=255)
                                                                            tCGGTTTCGTGGTGAATATGATTGCCACAATCACGCTATTTTATTAAGTCATTAAAAAGACAAAACAGGCCGCCTGGGCCTGTTTTGTATTACTTCACAACGCGTAATGCCGGTCGACCACCGCGTGGTGGCTGCGGAGGTTCATCGTc                                                                     >  1:625/1‑149 (MQ=255)
                                                                                   cGTGGTGAATATGATTGCCACAATCACGCTATTTTATTAAGTCATTAAAAAGACAAAACAGGCCGCCTGGGCCTGTTTTGTATTACTTCACAACGCGTAATGCCGGTCGACCACCGCGTGGTGGCTGCGGAGGTTCATCGTCAGGATGa                                                              <  2:387504/149‑1 (MQ=255)
                                                                                   cGTGGTGAAGATGATTGCCACAATCACGCTATTTTATTAAGTCATTAAAAAGACAAAACAGGCCGCCTGGGCCTGTTTTGTATTACTTCACAACGCGTAATGCCGGTCGACCACCGCGTGGTGGCTGCGGCGGTTCTTCGTCAGGATGa                                                              <  2:322247/149‑1 (MQ=255)
                                                                                     tggtgAATATGATTGCCACAAGCACGCTATTTTATTAAGTCATTAAAAAGACAAAACAGGCCGCCTGGGCCTGTTTTGTATTACTTCGCATGGCGTAATGCCGGTCGACCACCGCGTGGTGGCTGCGGAGGTTCATCGTCAGGATGAgt                                                            <  1:334390/149‑1 (MQ=255)
                                                                                              tGATTGCCACAATCACGCTATTTTATTAAGTCATTAAAAAGACAAAACAGGCCGCCTGGGCCTGTTTTGTATTACTTCACAACGCGTAATGCCGGTCGACCACCGCGTGGTGGCTGCGGAGGTTCATCGTCAGGATGAGTGTCatcatc                                                   <  2:332755/149‑1 (MQ=255)
                                                                                                     cacaATTAAGCTATTTTATGAAGTCATTAAAAAGACAAACCAGGCCGCCTGGGCTTGGTTTGTATTACTTCCCAACGCGTACTGCCGGTCGACCCCCTCGTGGTGGCTGCGGAGGATCAGCGTCAGGATTAGTGTCATCATCGTGATCt                                            <  2:141460/149‑1 (MQ=255)
                                                                                                            aCGCTATTTTATTAAGTCATTAAAAAGACAAAACAGGCCGCCTGGGCCTGTTTTGTATTACTTCACAACGCGTAATGCCGGTCGACCACCGCGTGGTGGCTGCGGAGGTTCATCGTCAGGATGAGTGTCATCATCGTGATCTGGCTTGt                                     <  2:8392/149‑1 (MQ=255)
                                                                                                             cGCTATGTTATGAAGTCATTAGAAAGACAAGATCGGCCGCCTGGGCCTGTGTTGTATTACTTCACAACGCGTAATGCCGGGCGACCACCGCGTGGTGGCTGCGGGGGTTCATCGTCAGGATGAGTGTCATCATCGTGATCTGGCTTGTc                                    <  2:151642/149‑1 (MQ=255)
                                                                                                                            tCATTAAAAAGACAAAACAGGCCGCCTGGGCCTGTTTTGTATTACTTCACAACGCGTAATGCCGGTCGACCACCGCGTGGTGGCTGCGGAGGTTCATCGTCAGGATGAGTGTCATCATCGTGATCTGGCTTGTCGCCATCAATAACCGa                     >  2:249922/1‑149 (MQ=255)
                                                                                                                            tCATTAAAAAGAAAAAATATGCGGTCTGGGCCTGTTTTGTATTACTTAACAACGCGTAATTCCGGTCGACCACCGCGTGGTGGCTGCGGAGGTCCATCGTCAGGATGAGTGTCATCATCGTGATCTGGCTTGTCGCCATCAATAACCGa                     <  1:142571/149‑1 (MQ=255)
                                                                                                                                  aaaaGACAAAACAGGCCGCCTGGGCCTGTTTTGTATTACTTCACAACGCGTAATGCCGGTCGACCACCGCGTGGTGGCTGCGGAGGTTCATCGTCAGGATGAGTGTCATCATCGTGATCTGGCTTGTCGCCATCAATAACCGACATAAc               >  2:38272/1‑149 (MQ=255)
                                                                                                                                          aaaCAGGCCGCCTGGGCCTGTTTTGTATTACTTCACAACGCGTAATGCCGGTCGACCACCGCGTGGTGGCTGCGGAGGTTCATCGTCAGGATGAGTGTCATCATCGTGATCTGGCTTGTCGCCATCCATAACCAACCAACCGGttcgct       >  2:197881/1‑145 (MQ=255)
                                                                                                                                               ggCCGCCTGGGCCTGTTTTGTATTACTTCACAACGCGTAATGCCGGTCGACCACCGCGTGGTGGCTGCGGAGGTTCATCGTCAGGATGAGTGTCATCATCGTGATCTGGCTTGTCGCCATCAATAACCGACATAACGGTTTCGTTGTCt  <  1:299425/149‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                 
GTGGTGGTCGTTACGTCAGGAATTGCAGGCGTTTCGCCTTTTGCGGTGGTGAAGCGGACAGCGATCCCCATGGCAGTCGGTTTCGTGGTGAATATGATTGCCACAATCACGCTATTTTATTAAGTCATTAAAAAGACAAAACAGGCCGCCTGGGCCTGTTTTGTATTACTTCACAACGCGTAATGCCGGTCGACCACCGCGTGGTGGCTGCGGAGGTTCATCGTCAGGATGAGTGTCATCATCGTGATCTGGCTTGTCGCCATCAATAACCGACATAACGGTTTCGTTGTCT  >  NC_000913_3_pae_pgi/3376114‑3376405

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