New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913_3_pae_tpiA = 449816219 (0.490)4 (0.110) 4/274 5.2 19.4% intergenic (+306/‑26) intB/insC1 pseudogene, integrase homology;IS, phage, Tn; Phage‑related functions and prophages; KpLE2 phage‑like element; P4‑like integrase/IS2 repressor TnpA
?NC_000913_3_pae_tpiA = 4498168 15 (0.400)intergenic (+312/‑20) intB/insC1 pseudogene, integrase homology;IS, phage, Tn; Phage‑related functions and prophages; KpLE2 phage‑like element; P4‑like integrase/IS2 repressor TnpA
Rejected: Coverage evenness skew score above cutoff.
Rejected: Frequency below/above cutoff threshold.

AGGGTGATGCGGGATAACTGACGTAAGGAGTGCTGCAGGGATTTGATGGTCCCTTGTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGTTTTTAGACTAGTCATTGGAGAACAGATGATTGATGTCTTAGGGCCGGAGAAACGCAGACGGCGTACCACACAGGAAAAGA  >  NC_000913_3_pae_tpiA/4498041‑4498242
                                                                                                                               |                                                                          
aGGGTGATGCGGGATAACTGACGTAAGGAGTGCTGCAGGGATTTGATGGTCCCTTGTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGTTTTTAGACTAGTCATTGGAGAACAgat                                                       <  2:137933/149‑1 (MQ=255)
 gggTGATGCGGGATAACTGACGTAAGGAGTGCTGCAGGGATTTGATGGTCCCTTGTAGGCTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAATCTTAATACTAGTCTTTAGACTAGTCATTGGAAAACAgatg                                                      >  1:39925/1‑149 (MQ=255)
           ggATAACTGACGTAAGGAGTGCTGCAGGGATTTGATGGTCCCTTGTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGTTTTTAGACTAGTCATTGGAGAACAGATGATTGATGTCt                                            >  2:294194/1‑149 (MQ=255)
           ggATAACTGACGTAAGGAGTGCTGCAGGGATTTGATGGTCCCTTGTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGTTTTTAGACTAGTCATTGGAGAACAGATGATTGATGTCt                                            >  2:294423/1‑149 (MQ=255)
            gATAACTGACGTAAGGAGTGCTGCAGGGATTTGATGGTCCCTTGTAGACTGGTCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGTTTTTAGACTAGTCATTGGAGAACAGATGATTGATGTCtt                                           <  2:550105/149‑1 (MQ=255)
            gATAACTGACGTAAGGAGTGCTGCAGGGATTTGATGGTCCCTTGTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGTTTTTAGACTAGTCATTGGAGAACAGATGATTGATGTCtt                                           <  2:21259/149‑1 (MQ=255)
               aaCTGACGTAAGGAGTGCTGCAGGGATTTGATGGTCCCTTGTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGTTTTTAGACTAGTCATTGGAGAACAGATGATTGATGTCTTAgg                                        <  2:45464/149‑1 (MQ=255)
                           gAGTGCTGCAGGGATTTGATGGTCCCTTGTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTGAATAAGTGATAGTCTTAATACTAGTTTTTAGACTAGTCATTGGAGAACAGATGATTGATGTCTTAGGGCCGGAGGAACg                            >  1:561305/1‑149 (MQ=255)
                                    aGGGATTTGATGGTCCCTTGTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGTTTTTAGACTAGTCATTGGAGAACAGATGATTGATGTCTTAGGGCCGGAGAAACGCAGACGGCg                   <  1:681220/149‑1 (MQ=255)
                                        aTTTGATGGTCCCTTGTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGTTTTTAGACTAGTCATTGGAGAACAGATGATTGATGTCTTAGGGCCGGAGAAACGCAGACGGCGTAcc               >  1:447620/1‑149 (MQ=255)
                                            gATGGTCCCTTGTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGTTTTTAGAc                                                                          >  1:582565/1‑86 (MQ=255)
                                            gATGGTCCCTTGTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGTTTTTAGAc                                                                          <  2:582565/86‑1 (MQ=255)
                                              tGGTCCCTTGTAGACTGGCCCCCTTAATCTTCAGACAATCAATATCACTTAAATAAGTGATAGTCTTAATACTAGTTTTTAGACTAGTCATTGGAGAACAGATGATTGATGTCTTAGGGCCGGAGAAACGCAGACGGCGTACCACACAg         <  1:258124/149‑1 (MQ=255)
                                              tGGTCCCTTGTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGTTTTTAGACTAGTCATTGGAGAACAGATGATTGATGTCTTAGGGCCGGAGAAACGCAGACGGCGTACCACACAg         >  2:327167/1‑149 (MQ=38)
                                                   ccTTGTAGACTGGCCCCCTGAATCTCCAGACAACCAATAGCACTTAAATAAGTGATAGTCTTAATACTAGTTTTTAGACTAGTCATTGGAGAACAGATGATTGATGTCTTAGGGCCGGAGAAACGCAGACGGCGTACCACACAGGaaaa    <  2:599910/149‑1 (MQ=18)
                                                     ttGTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGTTTTTAGACTAGTCATTGGAGAACAGATGATTGATGTCTTAGGGCCGGAGAAACGCAGACGGCGTACCACACAGGAAAAGa  >  1:358297/1‑149 (MQ=32)
                                                     ttGTAGACTGGCCCCCTGAATCTCCAGACAACCAATAACAATTAAATAAGTGATAGTCTTAATACTAGTTTTTAGACTAGTAAATGGAGAACAGATTATTGATGTCTAAAGTCAGGAGACACGCAGACGACGTACGACACAgacacagc  >  2:14453/1‑142 (MQ=11)
                                                                                                                               |                                                                          
AGGGTGATGCGGGATAACTGACGTAAGGAGTGCTGCAGGGATTTGATGGTCCCTTGTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGTTTTTAGACTAGTCATTGGAGAACAGATGATTGATGTCTTAGGGCCGGAGAAACGCAGACGGCGTACCACACAGGAAAAGA  >  NC_000913_3_pae_tpiA/4498041‑4498242

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