New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913_3_pae_tpiA 1207790 =4 (0.140)21 (0.820) 20/254 0.4 86.5% coding (290/630 nt) stfP e14 prophage; uncharacterized protein
?NC_000913_3_pae_tpiA 1209619 = 3 (0.120)pseudogene (1/501 nt) stfE pseudogene, e14 prophage; side tail fiber protein fragment family;Phage or Prophage Related

TTATTTGTCGGATAAGCCTTTGCCAGTTCCGGGTATTCTTCAGCAGAAAAAGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGACCAGAACACAGATCTACAACGTAATGCGCTGAACTGTGCAGGATGCGA  >  NC_000913_3_pae_tpiA/1209471‑1209767
                                                                                                                                                    |                                                                                                                                                     
ttatttGTCGGATAAGCCTTTGCCAGTTCCGGGTATTCTTCAGCAGAAAAAGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAAc                                                                                                                                                       <  2:247433/149‑1 (MQ=255)
      gcagatgcttttgcggcatcccacccccactgaatgattactgattgctcaataccaccaatcataaccggtatagctatgaaaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTAt                                                                                                                                                 <  1:170410‑M2/7‑1 (MQ=255)
                               cccactgaatgattactgattgctcaataccaccaatcataaccggtatagctatgaaaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAAt                                                                                                                        <  2:268978‑M2/32‑1 (MQ=255)
                                 atctgaaggactacggctggctcaatatcagcacacatatccgttatagctatgaaaccggtgcttcctaatgtagccgtggctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATgg                                                                                                                      <  2:155970‑M2/34‑1 (MQ=255)
                                       atgattactgattgctcaataccaccaatcataaccggtatagctatgaaaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGcatcat                                                                                                                >  1:213380‑M2/110‑149 (MQ=255)
                                             gAAAAAGCGGCACCAGTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCTACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTTTCCCAAACCATCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGc                                                                                                          <  2:353007/149‑1 (MQ=255)
                                                     ctcaataccaccaatcataaccggtatagctatgaaaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAACCAAATGGCATCATTCCTGCTTTTACCa                                                                                                  >  2:335562‑M2/96‑149 (MQ=255)
                                                      tcaataccaccaatcataaccggtatagctatgaaaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAg                                                                                                 >  1:83012‑M2/95‑149 (MQ=255)
                                                      tcaataccaccaatcataaccgataccacgcccaaaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttcccccaaaccaaCGCTTATAAAAATGAAGAAATAACAAGCAAATAGCAACATTCCTCCTTTTACCAg                                                                                                 >  2:488382‑M2/95‑149 (MQ=255)
                                                                caatccgaaccggtatagctatgaaaccggtgctttctaatgtagccgttgctgcattgagttttgctccttctcccaacccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGGTTTaa                                                                                       <  2:416165‑M2/65‑1 (MQ=255)
                                                                caatcataaccggtatagctatgaaaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTaa                                                                                       >  2:443077‑M2/85‑149 (MQ=255)
                                                                 aatcataaccggtatagctatgaaaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAAc                                                                                      <  2:83012‑M2/66‑1 (MQ=255)
                                                                  atcataaccggtatagctatgaaaccggtgcgtcctcatgtagccgttgctgcattgcgttttgctccttctcccaaaccaaCGTTTATGAAAATGCAGAAATCACAATCAAATGACATCATTCCTGCTTTTACCAGGGGGATTTAACa                                                                                     >  2:85930‑M2/83‑149 (MQ=255)
                                                                  atcataaacggtatagctatgaaacgtgtgcgtccgcatgtaggcgttgctgtattgagttttgctgcttctcccaagccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACa                                                                                     <  2:417840‑M2/67‑1 (MQ=255)
                                                                   tcataaccggtatagctatgaaaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACAt                                                                                    <  1:75130‑M2/68‑1 (MQ=255)
                                                                     ataaccggtatagctatgaaaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAAAGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGAGGGATTTAACCTGc                                                                                  >  1:478553‑M2/80‑149 (MQ=255)
                                                                        acctgtatagctatgaaaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTa                                                                               >  2:341488‑M2/77‑149 (MQ=255)
                                                                               tagctatgaaaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAATAAATAACAAGCAAATGGCATCATTCCTTCTTTTACCAGGGGGATTTAACATGCTTATTGGCTa                                                                       <  1:142444‑M2/80‑1 (MQ=255)
                                                                                   tatgaaaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTTTTGGCTATGTa                                                                   >  1:298026‑M2/66‑149 (MQ=255)
                                                                                                     taatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGAc                                                 >  2:118932‑M2/48‑149 (MQ=255)
                                                                                                         gtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGACCAGa                                             >  1:339489‑M2/44‑149 (MQ=255)
                                                                                                            gccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGACCAGaact                                          >  2:342005‑M2/41‑148 (MQ=255)
                                                                                                                 aacAACAGGGAATGCAGATCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGACCAGAACACAGAt                                     <  2:366264/147‑1 (MQ=255)
                                                                                                                   ctgcattgagttttgctccttctcccaaaccaaCGATGATGAAAATGAAGAAATAACAAGCAAATGACATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTAGCTATGTACGCGTATCAACAAATGACCAAACCACAGATCt                                   >  2:43033‑M2/34‑149 (MQ=255)
                                                                                                                                         tccgaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGACCAGAACACAGATCTACAACGTAATGCGCTGAACtgt             <  2:97103‑M2/138‑1 (MQ=255)
                                                                                                                                           ccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGACCAGAACCCAGATCTACAACGTAATGCGCTAAACTCTGc           >  2:9003‑M2/10‑149 (MQ=255)
                                                                                                                                           ccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGACCAGAACACAGAGCTACAACGTAATCCGCAGAActgtga           >  1:163032‑M2/10‑148 (MQ=255)
                                                                                                                                                 caaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGACCAGAACACAGATCTACACCGTAATGCGCTGAACTGTGCAggagg     >  1:235404‑M2/4‑147 (MQ=255)
                                                                                                                                                   aCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGACCAGAACACAGATCTACAACGTAATGCGCTGAACTGTGCAGGATGCg   >  1:3305‑M2/2‑149 (MQ=255)
                                                                                                                                                    cGTTTATGACAATGAAGAAATAACAGGCAATAAGCATCCTTAGCAACTTTACCAGGGTGATTTACCCGGTCATTGAGCGATATTCAGAAAAGAGCAAAGGTCCAGAACACAGATCTACAACGTAATGCGCTGAACTGTGCAGGATGCGa  <  1:497743/149‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
TTATTTGTCGGATAAGCCTTTGCCAGTTCCGGGTATTCTTCAGCAGAAAAAGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGACCAGAACACAGATCTACAACGTAATGCGCTGAACTGTGCAGGATGCGA  >  NC_000913_3_pae_tpiA/1209471‑1209767

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