New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913_3_pae_tpiA 1207790 =7 (0.170)26 (0.710) 23/254 0.7 84.9% coding (290/630 nt) stfP e14 prophage; uncharacterized protein
?NC_000913_3_pae_tpiA 1209619 = 3 (0.080)pseudogene (1/501 nt) stfE pseudogene, e14 prophage; side tail fiber protein fragment family;Phage or Prophage Related

GTCGGATAAGCCTTTGCCAGTTCCGGGTATTCTTCAGCAGAAAAAGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGACCAGAACACAGATCTACAACGTAATGCGCTGAACTGTGCAGGATGCGA  >  NC_000913_3_pae_tpiA/1209477‑1209767
                                                                                                                                              |                                                                                                                                                    
gTCGGATAAGCCTTTGCCAGTTCCGGGTATTCTTCAGCAGAAAAAGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTTGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTAt                                                                                                                                                <  2:313093/149‑1 (MQ=255)
    atgcttttgcggcattccacccccactgaatgattactgattgctcaataccaccaatcataaccggtatagctatgaaaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGaaa                                                                                                                                            <  1:429438‑M2/11‑1 (MQ=255)
    atgcttttgcggcattccacccccactgaatgattactgattgctcaataccaccaatcataaccggtatagctatgaaaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGGTTATGaaa                                                                                                                                            >  1:580703‑M2/139‑149 (MQ=255)
     tgcggtggcggcattccacccccactgaatgattactgattgctcaataccaccaatcataaccggtatagctatgaaaccggtgcgtcctaatgtagccgtggctgcattgagttttgctccttctcccaaaccaaCGTTTATGaaaa                                                                                                                                           <  2:305763‑M2/12‑1 (MQ=255)
              ggcattccacccccactgaatgattactgattgctcaataccaccaatcataaccggtatagctatgaaaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGttt                                                                                                                                                  >  2:98786‑M2/129‑133 (MQ=255)
              ggcattccacccccactgaatgattactgattgctcaataccaccaatcataaccggtatagctatgaaaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGttt                                                                                                                                                  <  1:98786‑M2/5‑1 (MQ=255)
                 attccacccccactgaatgattactgattgcgcaataccaccaatcataaccggtatagctatgaaaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAAc                                                                                                                               <  2:2024‑M2/24‑1 (MQ=255)
                 actccatcgccactgaatgatcactgatggcgcacttccaccaatcataaccggtatagctctgaaaccggtgcgtcctaatgtagccgttgctgcattgagttttccttcttctcccaaaccaaCGTTTATGAAAATGAAGAAatacc                                                                                                                               <  2:164360‑M2/24‑3 (MQ=255)
                     caccccccctgaaggattactgattgctcaatagcaccaatcataaccggtatagctatgaaaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccatCGTTTATGAAAATGAAGAAATAACAAGc                                                                                                                           <  2:282268‑M2/28‑1 (MQ=255)
                                          gattgctcaataccacgaatcataaccggtatagctatgaaaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCttt                                                                                                      >  1:629382‑M2/101‑149 (MQ=255)
                                          aaaGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCttt                                                                                                      <  2:462410/149‑1 (MQ=255)
                                                tcaataccaccaatcataaccggtatagctatgaaaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAg                                                                                                >  1:123914‑M2/95‑149 (MQ=255)
                                                  aataccaccaatcataaccggtatagctatgaaaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAggg                                                                                              >  1:252837‑M2/93‑149 (MQ=255)
                                                   ataccaccaatcataaccggtatagctatgaaaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAgggg                                                                                             <  1:311244‑M2/58‑1 (MQ=255)
                                                        accaatcataaccggtatagctatgaaaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGAttt                                                                                        >  1:337594‑M2/87‑149 (MQ=255)
                                                          caatcataaccggtatagttatgaaaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaatcaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTaa                                                                                      <  2:385580‑M2/65‑1 (MQ=255)
                                                          caatcataaccggtatagctatgaaaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctccgaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTaa                                                                                      <  1:410411‑M2/65‑1 (MQ=255)
                                                               ataaccggtatagctatgaaaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGc                                                                                 <  2:123914‑M2/70‑1 (MQ=255)
                                                               agaaccggtagagctatgaaaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGc                                                                                 <  1:317861‑M2/70‑1 (MQ=255)
                                                                        atagctatgaaaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACGAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCt                                                                        >  2:548095‑M2/71‑149 (MQ=255)
                                                                            ctatgaaaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTCATGAAAAGGAAGAAATAACACGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTAATGGCTATGt                                                                    >  1:524151‑M2/67‑149 (MQ=255)
                                                                               tgaaaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGGTTTTACGAGGGGGATTTAACATGCTTATTGGCTATGTAcg                                                                 >  1:18568‑M2/64‑149 (MQ=255)
                                                                                    ccggtgcgtcctaatgtagccgttgccgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTAt                                                            <  1:43691‑M2/91‑1 (MQ=255)
                                                                                        aaGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAAAAATAACAAGCAAATGGCATCACCCGGTCTCTCTCCAGTGGGATTTAACATGCTTGTTGGCTATGTCCGCGTATcacc                                                        >  1:20903/1‑147 (MQ=255)
                                                                                            tcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAAt                                                    >  2:602894‑M2/51‑149 (MQ=255)
                                                                                                 atgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTg                                                                                                          >  1:48623‑M2/46‑90 (MQ=255)
                                                                                                 atgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTg                                                                                                          <  2:48623‑M2/45‑1 (MQ=255)
                                                                                                    tagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGACCAGaa                                            <  2:584299‑M2/107‑1 (MQ=255)
                                                                                                                  ttgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGACCAGAACACAGATCTAAAAc                              >  1:109856‑M2/29‑149 (MQ=255)
                                                                                                                                   tcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGACCAGAACACAGATCTACAACGTAATGCGCTGAACtgt             <  2:331107‑M2/138‑1 (MQ=255)
                                                                                                                                     ccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTc                                                                                                             >  1:608492‑M2/10‑51 (MQ=255)
                                                                                                                                     ccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTc                                                                                                             <  2:608492‑M2/42‑1 (MQ=255)
                                                                                                                                     ccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGACCAGAACACAGATCTACAACGTAATGCGCTGAAGTGTGc           >  1:237499‑M2/10‑149 (MQ=255)
                                                                                                                                     ccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGACCAGAACACAGATCTACAACGTAATGCGCTGAACTGTGc           >  2:11908‑M2/10‑149 (MQ=255)
                                                                                                                                       aaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGACCAGAACACAGATCTACAACGTAATGCGCTGAACTGTGCAg         <  1:602894‑M2/142‑1 (MQ=255)
                                                                                                                                          ccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGACCAGAACACAGATCTACAACGTAATGCGCTGAACTGTGCAGGAt      >  1:705008‑M2/5‑149 (MQ=255)
                                                                                                                                              cGTTTATGAAAATGAAGAAATAACAATCAAATCGCATCACTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGACCAGAACTCAGATCTACAACGTAATGCGCTGAACTGTGCAGGATGCGa  <  1:597440/149‑1 (MQ=255)
                                                                                                                                              |                                                                                                                                                    
GTCGGATAAGCCTTTGCCAGTTCCGGGTATTCTTCAGCAGAAAAAGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGACCAGAACACAGATCTACAACGTAATGCGCTGAACTGTGCAGGATGCGA  >  NC_000913_3_pae_tpiA/1209477‑1209767

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