New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913_3_pae_tpiA = 12078059 (0.280)22 (0.780) 20/252 0.5 81.5% coding (305/630 nt) stfP e14 prophage; uncharacterized protein
?NC_000913_3_pae_tpiA = 1209602 2 (0.070)pseudogene (18/501 nt) stfE pseudogene, e14 prophage; side tail fiber protein fragment family;Phage or Prophage Related

ATCAGGCAATTTATTTGTCGGATAAGCCTTTGCCAGTTCCGGGTATTCTTCAGCAGAAAAAGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGACCAGAACACAGATCTACAACGTAATGCGCT  >  NC_000913_3_pae_tpiA/1209461‑1209749
                                                                                                                                             |                                                                                                                                                   
aTCAGGCAATTTATTTGTCGGATAAGCCTTTGCCAGTTTCGGGTATTCTTCAGCAGAAAAAGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGccttctc                                                                                                                                              >  1:193919‑M2/1‑142 (MQ=255)
                                  aGTTCCGGGTATTCTTCAGCAGAAAAAGCGGCACCATTGCAATTCAGCCAGCCTGTTGACGGCGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGccttctcccaaaccaaggttttcgagagccgtttgcacagt                                                                                                            <  1:42170‑M2/149‑42 (MQ=255)
                                          gTATTCTTCAGCAGAAAAAGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGccttctcccaaaccaaggttttcgagagccgtttgcacagtgccatc                                                                                                      >  1:249951‑M2/1‑100 (MQ=255)
                                          gTATTCTTCAGCAGAAAAAGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGccttctcccaaaccaaggttttcgagagccgtttgcacagtgccatc                                                                                                      <  2:249951‑M2/147‑48 (MQ=255)
                                               cttcAGCAGAAAAAGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGccttctcccaaaccaaggttttcga                                                                                                                            <  2:59049‑M2/120‑26 (MQ=255)
                                               cttcAGCAGAAAAAGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGccttctcccaaaccaaggttttcga                                                                                                                            >  1:59049‑M2/1‑95 (MQ=255)
                                                ttcAACAGAAAAAGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTAGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGccttctcccaaaccaaggttttcgagagccgtttgcacagtgccatccgatttga                                                                                              <  2:552878‑M2/149‑56 (MQ=255)
                                                 gtagcagaaGAAGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGccttctcccaaaccaaggttttcgagagccgtttgcacagtgccatccgatttgat                                                                                             >  2:371510‑M2/3‑93 (MQ=255)
                                                   agcagAAAAAGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGccttctcccaaaccaaggttttcgagagccgtttgcacagtgccatccgatttgatat                                                                                           <  2:84092‑M2/149‑59 (MQ=255)
                                                                  aCCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGccttctcccaaaccaaggttttcgagagccgtttgcacagtgccatccgatttgatatcgccaaacggattct                                                                            <  2:218762‑M2/149‑74 (MQ=255)
                                                                                   agccTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGccttctcccaaaccaag                                                                                                                                    <  2:292814‑M2/76‑18 (MQ=255)
                                                                                   agccTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGccttctcccaaaccaag                                                                                                                                    >  1:292814‑M2/1‑59 (MQ=255)
                                                                                       tGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGccttctcccaaaccaaggttttcgagagccgtttgcacagtgccatccgatttgatatcggcaaacggattcttgcggcttaacagcagcgcgc                                                       <  1:48463‑M2/149‑95 (MQ=255)
                                                                                         ttGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGccttctcccaaaccaaggttttcgagagccgtttgcacagtgccatccgatttgatatcgccaaacggattcttgcggcttaacagcagcgcgcga                                                     <  2:514765‑M2/149‑97 (MQ=255)
                                                                                              ggAGTGGCTGAAGGCCAGGGAACCGGGACACCAACAGGTAATGCAGAGccttctcccaaaccaaggttttcgagagccgtttgcacagtgccatccgatttgatatcgccaaacggattcttgcggcttaacagcagcgcgcgaagtgc                                                >  1:552878‑M2/1‑48 (MQ=255)
                                                                                              ggAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGccttctcccaaaccaaggttttcgagagccgtttgcacagtgccatccgatttgatatcgccaaacggattcttgcggcttaacagcagcgcgcgaagtgc                                                >  2:414702‑M2/1‑48 (MQ=255)
                                                                                                    gCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTAt                                          >  1:352915/1‑149 (MQ=255)
                                                                                                      tGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGccttctcccaaaccaaggttttcgagagccgtttgcacagtgccatccgatttgatatcgccaaacggattcttgcggcttaacagcagcgcgcgaagtgcggtaagta                                        <  1:383953‑M2/149‑110 (MQ=255)
                                                                                                      tGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGccttctcccaaaccaaggttttcgagagccgtttgcacagtgccatccgatttgatatcgccaaacggattcttgcggcttaacagcagcgcgcgaagtgcggt                                             >  1:384525‑M2/1‑40 (MQ=255)
                                                                                                      tGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGccttctcccaaaccaaggttttcgagagccgtttgcacagtgccatccgatttgatatcgccaaacggattcttgcggcttaacagcagcgcgcgaagtgcggt                                             <  2:384525‑M2/144‑105 (MQ=255)
                                                                                                                      gggACACCAACAGGTAATGCAGAGccttctcccaaaccaagtttttcgagagccgtttgcacagtgccatccgatttgatatcgccaaacgcattcttgcggcttaacagcagcgcgcgaagtgcggtaagtaactggtcgtgccgccc                        >  1:230589‑M2/1‑24 (MQ=255)
                                                                                                                          aaccAACGGGTAATGAAGGGccttctgccacgccaaggttgtggagaggcgtgtgcacagtgccatccgatttgatatcgccaaacggattcttgcggcgtaacagcagcgcgcgaagtgcggtaagtaactggtcgtgccgccccttc                    >  2:82205‑M2/2‑20 (MQ=255)
                                                                                                                             cAACAGGTAATTCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCCTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGACCAGAACACAGATCt                 <  1:539794/149‑1 (MQ=255)
                                                                                                                                cAGGTAATGCACAGccttctcccaaaccaaggttttcgagagccgtttgcacagtgccatccgatttgatatcgccaaacggagtcttgcggctgaacagcagcgcgcgaagtgcggtaagtaactggtcgtgccgccccttctccagg              >  1:84926‑M2/1‑14 (MQ=255)
                                                                                                                                            agccttctcctaaaccaaggttttcgagagccgtttgcacagtgacaacctgtcgctgagcgccaaacggattcttgcggcttaacagcagcgcgcgaagtgcggtaagtaactggtcgtgccgccccttctccaggctggcaccggat  >  1:444645‑M2/1‑2 (MQ=255)
                                                                                                                                             |                                                                                                                                                   
ATCAGGCAATTTATTTGTCGGATAAGCCTTTGCCAGTTCCGGGTATTCTTCAGCAGAAAAAGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGACCAGAACACAGATCTACAACGTAATGCGCT  >  NC_000913_3_pae_tpiA/1209461‑1209749

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.