New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913_3_hsa_tpiA 1207790 =10 (0.210)39 (0.910) 34/250 0.3 84.0% coding (290/630 nt) stfP e14 prophage; uncharacterized protein
?NC_000913_3_hsa_tpiA 1209619 = 6 (0.140)pseudogene (1/501 nt) stfE pseudogene, e14 prophage; side tail fiber protein fragment family;Phage or Prophage Related

ATAAGCCTTTGCCAGTTCCGGGTATTCTTCAGCAGAAAAAGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGACCAGAACACAGATCTACAACGTAATGCGCTGAACTGTGCAGGATG  >  NC_000913_3_hsa_tpiA/1209482‑1209764
                                                                                                                                         |                                                                                                                                                 
aTAAGCCTTTGCCAGTTCCGGGTATTCTTCAGCAGAAAAAGCGGCGCAATTGCTTTTCACGAAGCCTGTTGGAAGTGTGCCCGCCTCCTTCAGACCCGCCACATCAAGAGGCAAAAGAGAACCTTCTCCCAAACCAAGGttttcgagag                                                                                                                                        >  2:201689/1‑142 (MQ=255)
         tGCCAGTTCCGGGTATTCTTCAGCAGAAAAAGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAAt                                                                                                                               <  1:50050/149‑1 (MQ=255)
         ggcattccacccccactgaatgattactgattgctcaataccacaaatcctatcctttatagctaggaaagctgtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaagcaaCGTTGATGAAAAGTAAGAAAt                                                                                                                               >  1:382002‑M2/129‑149 (MQ=255)
          gCCAGTTCCGGGTATTCTTCAGCAGAAAAAGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACg                                                                                                                                                  >  2:323259/1‑129 (MQ=255)
          gCCAGTTCCGGGTATTCTTCAGCAGAAAAAGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACg                                                                                                                                                  <  1:323259/129‑1 (MQ=255)
           cattccacccccactgaatgattactgattgctcaataccaccaatcataaccggtatagctatgaaaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATaa                                                                                                                             <  2:239252‑M2/23‑1 (MQ=255)
                   ccccactgaatgattactgattgctcaataccaccaatcataaccggtatagctatgaaaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCaaa                                                                                                                     >  1:618882‑M2/119‑149 (MQ=255)
                          gaatgattactgattgctcaataccaccaatcataaccggtatagctatgaaaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGcatc                                                                                                              >  1:159659‑M2/112‑149 (MQ=255)
                           aatgattactgattgctcaataccaccaatcataaccggtatagctatgaaaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATAAAGAAATAACAAGCAAATGGcatca                                                                                                             >  2:357927‑M2/111‑149 (MQ=255)
                            atgattactgattgctcaataccaccaatcataaccggtatagctatgaaaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGcatcat                                                                                                            <  2:618882‑M2/40‑1 (MQ=255)
                              gattactgattgctcaataccaccaatcataaccggtatagctatgaaaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATgaag                                                                                                                                   >  2:142308‑M2/108‑124 (MQ=255)
                              gattactgattgctcaataccaccaatcataaccggtatagctatgaaaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATgaag                                                                                                                                   <  1:142308‑M2/17‑1 (MQ=255)
                              gattactgattgctcaataccaccaatcataaccggtatagctatgaaaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTc                                                                                                          <  1:783248‑M2/42‑1 (MQ=255)
                               attactgattgctcaataccaccaatcataaccggtatagctatgaaaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaatccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTcc                                                                                                         <  1:568278‑M2/43‑1 (MQ=255)
                                           tcaataccaccaatcataaccggtatagctatgaaaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAg                                                                                             >  1:345116‑M2/95‑149 (MQ=255)
                                               cATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAggggg                                                                                         >  1:227765/1‑149 (MQ=255)
                                                   accaatcataaccggtatagctatgaaaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGAttt                                                                                     >  2:535147‑M2/87‑149 (MQ=255)
                                                   accaatcataaccggtatagctatgaaaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGAttt                                                                                     >  1:609790‑M2/87‑149 (MQ=255)
                                                     caatcataaccggtatagctatgaaaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTaa                                                                                   <  1:39698‑M2/65‑1 (MQ=255)
                                                       atcataaccggtatagctatgaaaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACa                                                                                 >  1:33259‑M2/83‑149 (MQ=255)
                                                              ccggtatagctatgaaaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTAt                                                                          <  1:659784‑M2/74‑1 (MQ=255)
                                                                ggtatagctatgaaaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccgaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTg                                                                        <  2:518653‑M2/76‑1 (MQ=255)
                                                                ggtatagctatgaaaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTg                                                                        <  2:512064‑M2/76‑1 (MQ=255)
                                                                     agctatgaaaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaatcaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTAt                                                                   <  2:679355‑M2/81‑1 (MQ=255)
                                                                     agctatgaaaccggtgcgtcctaatgtacccgttgctgcattgagttttgctccttctcccaaatcaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTAt                                                                   <  2:679240‑M2/81‑1 (MQ=255)
                                                                      gctatgaaaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGcatcat                                                                                                            <  1:256333‑M2/40‑1 (MQ=255)
                                                                      gctatgaaaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGcatcat                                                                                                            >  2:256333‑M2/68‑107 (MQ=255)
                                                                       gCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGt                                                                 <  2:87207/149‑1 (MQ=255)
                                                                       ctatgaaaccggtgcgtcctaatgtagcccttgctgtattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGt                                                                 <  1:638322‑M2/83‑1 (MQ=255)
                                                                        tatgaaaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGc                                                                      <  1:310433‑M2/78‑1 (MQ=255)
                                                                        tatgaaaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGc                                                                      >  2:310433‑M2/66‑143 (MQ=255)
                                                                        tatgaaaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTa                                                                >  1:315712‑M2/66‑149 (MQ=255)
                                                                         atgaaaccggtgcgtcctaatgtagccgttgctggagtgagttttgctccttctcacaaaccagCGTTTAGGAAAATTAGGAAATAACAAGCAAATGGCATCATTCCTGCTTGTACCAGGGGGATTTAACATGCATGTTGGCCATCTAc                                                               >  2:435410‑M2/65‑149 (MQ=255)
                                                                            aaaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTAcgc                                                             >  1:581038‑M2/62‑148 (MQ=255)
                                                                            aaaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTAcgc                                                             <  2:581038‑M2/87‑1 (MQ=255)
                                                                             aaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTgg                                                                       <  1:333577‑M2/77‑1 (MQ=255)
                                                                             aaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTgg                                                                       >  2:333577‑M2/61‑137 (MQ=255)
                                                                                      gtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACaaa                                                  <  2:345116‑M2/98‑1 (MQ=255)
                                                                                           aatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGAcc                                             <  2:159659‑M2/103‑1 (MQ=255)
                                                                                               tagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGACCAGaa                                         <  2:550768‑M2/107‑1 (MQ=255)
                                                                                                         tgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGACCAGAACACAGATCTa                               >  2:334947‑M2/33‑149 (MQ=255)
                                                                                                            attgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGACCAGAACACAGATCTACaa                            <  1:535147‑M2/120‑1 (MQ=255)
                                                                                                             ttgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGACCAGAACACAGATCTACAAc                           >  2:554651‑M2/29‑149 (MQ=255)
                                                                                                             ttgagttttgccccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGACCAGAACACAGATCTACAAc                           <  1:282975‑M2/121‑1 (MQ=255)
                                                                                                                       ctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGACCAAAACACAGATCTACAACGTAATGCACt                 >  1:378916‑M2/19‑149 (MQ=255)
                                                                                                                           ttCTCCAAAACCAACCTGTATGTAAATGTAGAAGTAACAAGCATATGGCATCATTCCTGCTTTTACCAGTGGGATTTAACATGCTTATTGGCTGTGTACGCGTATCAACAAATAACCAGAACACAGATGTACAACGTAATGCTCTCAAc             >  2:693821/1‑149 (MQ=255)
                                                                                                                                ccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGACCAGAACACAGATCTACAACGTAATGCGCTGAACTGTGc        >  1:769897‑M2/10‑149 (MQ=255)
                                                                                                                                 caaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAAt                                                 <  2:98745‑M2/99‑1 (MQ=255)
                                                                                                                                 caaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAAt                                                 >  1:98745‑M2/9‑107 (MQ=255)
                                                                                                                                 caaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGACCAGAACACAGATCTACAACGTAATGGGCTGAACTGTGCa       >  1:397577‑M2/9‑149 (MQ=255)
                                                                                                                                 cAAACCAACGTTTATGAAAATGAAGAAATACCAAGCAAATGGCATCATTCCTGCTTTTACGAGGGGACTTTAACATGCTTATTGGCTATGTACGGGTATCAGCAAATGACCAGATCACAGATCTACAAGGTAGTGCGCGGAAGAGTGcc       >  2:353806/1‑148 (MQ=255)
                                                                                                                                    accaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGACCAGAACACAGATCTACAACGTAATGCGCTGAACTGTGCAGGa    <  2:22912‑M2/144‑1 (MQ=255)
                                                                                                                                      caaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATGAACAAATGACCAGAACACAGATCTGCAACGTAATGCGCTGAACTGTGCAGGATg  >  1:727535‑M2/4‑149 (MQ=255)
                                                                                                                                         |                                                                                                                                                 
ATAAGCCTTTGCCAGTTCCGGGTATTCTTCAGCAGAAAAAGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGACCAGAACACAGATCTACAACGTAATGCGCTGAACTGTGCAGGATG  >  NC_000913_3_hsa_tpiA/1209482‑1209764

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.