New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913_3_pae_tpiA 1207790 =9 (0.250)19 (0.590) 17/250 1.1 69.1% coding (290/630 nt) stfP e14 prophage; uncharacterized protein
?NC_000913_3_pae_tpiA 1209619 = 9 (0.280)pseudogene (1/501 nt) stfE pseudogene, e14 prophage; side tail fiber protein fragment family;Phage or Prophage Related

TGGTGGAGGCATCCGGTGCCAGCCTGGAGAAGGGGCGGCACGACCAGTTACTTACCGCACTTCGCGCGCTGCTGTTAAGCCGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCGGTTATGATTGGTGGTATTGAGCAATCAGTAATCATTCAGTGGGGGTGGAATGCCGCAAAAGCATCTGCCTCTGG  >  NC_000913_3_pae_tpiA/1207649‑1207937
                                                                                                                                             |                                                                                                                                                   
tgcacagttcagcgcattacgttgtagatctgtgttctggtcatttgttgatacgcgtacatagccaataagcatgttaaatccccctggtaaaagcaggaatgatgccatttgcttgttatttcttcattttcataaacgTTGGTTTg                                                                                                                                              >  1:445474‑M1/142‑149 (MQ=255)
     agttcagcgcattacgttgtagatctgtgttctggtcatttgttgatacgcgtacatagccaataagcatgttaaatccccctggtaaaagcaggaatgatgccatttgcttgttatttcttcattttcataaacgTTGGTTTGGgaga                                                                                                                                         <  1:548458‑M1/13‑1 (MQ=255)
       ttcagcgcattacgttgtagatctgtgttctggtcatttgttgatacgcgtacatagccaataagcatgttaaatccccctggtaaaagcaggaatgatgccatttgcttgttatttcttcattttcataaacgTTGGTTTGGGAGAAg                                                                                                                                       >  1:266713‑M1/135‑149 (MQ=255)
         cagcgcattacgttgtaggtctgtgttctggtcatttgttggtacgcgtacatagccaataagcatgttaaatccccctggtaaaagcaggaatgatggcatttgcttgttatttcttctttttcataaacgTTGGTTTGGGAGAAGGa                                                                                                                                     >  1:34643‑M1/133‑149 (MQ=255)
          agcgcattacgttgtagatctgtgttctggtcatttgttgatacgcgtacatagccaataagcatgttaaatccccctggtaaaagcaggaatgatgccatttgcttgttatttcttcattttcataaacgTTGGTTTGGGAGAAGGAg                                                                                                                                    >  1:54530‑M1/132‑149 (MQ=255)
                 tacgttgtagatctgtgttctggtcatttgttgatacgcgtacatagccaataagcatgttaaatccccctggtaaaagcaggaatgatcccttttgcgtgttatttcttcattttcataaacgTTGGTTTGGGAGAAGGAGCAAAACt                                                                                                                             <  1:467242‑M1/25‑1 (MQ=255)
                   cgttgtagatctgtgttctggtcatttgttgatacgcgtacatagccaataagcatgttaaatccccctggtaaaagcaggaatgatgccatttgcttgttatttcttcattttcataaacgTTGGTTTGGGAGAAGGAGCAAAACTCa                                                                                                                           >  2:568493‑M1/123‑149 (MQ=255)
                    gttgtagatctgtgttctggtcatttgttgatacgcgtacatagccaataagcatgttaaatccccctggtaaaagcaggaatgatgccatttgcttgttatttcttcattttcataaacgTTGGTTTGGGAGAAGGAGCAAAACTCaa                                                                                                                          >  1:62002‑M1/122‑149 (MQ=255)
                    aGCCTGGAGAAGGGGCGGCACGACCAGTTACTTACCGCACTTCGCGCGCTGCTGTTAAGCCGCAAGCATCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAGGGctctgcattacc                                                                                                                          >  1:96449‑M1/1‑137 (MQ=255)
                     ttgtagatctgtgttctggtcatttgttgatacgcgtacatagccaataagcatgttaaatccccctggtaaaagcaggaatgatgccatttgcttgttatttcttcattttcataaacgTTGGTTTGGGAGAAGGAGCAAAactcgat                                                                                                                         <  1:249086‑M1/29‑4 (MQ=255)
                        tagatctgtgttctggtcatttgttgatacgcgtacatagccaataagcatgttaaatccccctggtaaaagcaggaatgatgccatttgcttgttatttcttcattttcataaacgTTGGTTTGGGAGAAGGAGCAAAACTCAATgca                                                                                                                      >  2:222078‑M1/118‑149 (MQ=255)
                        tGGCGAAGTGGCGGCACGCCAAGTTCTTAACGGCACTTCGCGCGCTGCTGTTAAGCCGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCATTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGctctgcattacctgtt                                                                                                                      <  2:152276‑M1/149‑17 (MQ=255)
                                gggCCGCACGACCAGTTATCTCTCGCACTTCGCGCGCTGCGGTTAAGCCGCAAGAATCCGTTTGGCGATATCAAAGCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTGGGGAGAAGGctctgcattacctgttggtgtccc                                                                                                              <  1:356062‑M1/149‑25 (MQ=255)
                                       aCGACCAGTTACTTACCGCACTTCGCGCGCTGCTGTTAAGCCGCAAGAATCCGTTTGGCGATCTCAAATCGGCTGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGAGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACAttc                                                                                                       >  2:465011/1‑148 (MQ=255)
                                       aCGACCAGTTACTTACCGCACTTCGCGCGCTGCTGTTAAGCCGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGctctgcattacctgttggtgtcccggttccg                                                                                                       >  1:367582‑M1/1‑118 (MQ=255)
                                       aCGACCAGTTACTTACCGCACTTCGCGCGCTGCTGTTAAGCCGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTa                                                                                                       >  2:261733/1‑149 (MQ=255)
                                         gACCAGTTACTTACCGCACTTCGCGCGCTGCTGTTAAGCCGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGctctgcattacctgttggtgtcccggttccgtg                                                                                                     <  1:136429‑M1/149‑34 (MQ=255)
                                            ttgttgatacgcgtacatagccaataagcatgttaaatccccctggtaaaagcaggaatgatgccatttgcttgttatttcttcattttcataaacgTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACg                                                                                                  >  1:362257‑M1/98‑149 (MQ=255)
                                            cAGTTACTTACCGCACTTCGCGCGCTGCTGTTAAGCCGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGctctgcattacctgttggtgtcccggttccgtggcc                                                                                                  <  2:540828‑M1/149‑37 (MQ=255)
                                             aGTTACTTACCGCACTTCGCGCGCTGCTGTTAAGCCGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGctctgcattacctgttggtgtcccggttccgtggcct                                                                                                 >  1:16769‑M1/1‑112 (MQ=255)
                                                     cgcgtacatagccaataagcatgttaaacccccctggcaaaagcaggaatgatgccatttgcttgttatttcttcattttcataaacgTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGttt                                                                                         >  1:246162‑M1/89‑149 (MQ=255)
                                                          aCTTCGCGCGCTGCTGTTAAGCCGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGctctgcattacctgttggtgtcccggttccgtggccttcagccactccgc                                                                                    >  2:142637‑M1/1‑99 (MQ=255)
                                                              cgcgcgCTGCTGTTAAGCCGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGctctgcattacctgttggtgtcccggttccgtggccttcagccactccgccaac                                                                                <  1:183330‑M1/149‑55 (MQ=255)
                                                                      agcatgttaaatccccctggtaaaagcaggaatgatgccatttgcttgttatttcttcattttcataaacgTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCGGTTa                                                                        >  2:408051‑M1/72‑149 (MQ=255)
                                                                       ctgTTAAGCCGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCGGTTAt                                                                       >  1:379609/1‑149 (MQ=255)
                                                                       ctgTTAAGCCGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCGGTTAt                                                                       <  2:630533/149‑1 (MQ=255)
                                                                          tgttaaatccccctggtaaaagcaggaatgatgccatttgcttgttatttcttcattttcataaacgTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCGGTTATGAt                                                                    >  1:349680‑M1/68‑149 (MQ=255)
                                                                               ccGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGctctgcattacctgttggtgtcccggttccgtggccttcagccactccgccaacaggctggctgaaatgca                                                               >  2:86990‑M1/1‑78 (MQ=255)
                                                                               ccGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCGGTTATGATtggtg                                                               >  2:598190/1‑149 (MQ=255)
                                                                                cGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGctctgcattacctgttggtgtcccggttccgtggccttcagccactccgccaacaggctggctgaaatgcaa                                                              >  1:247209‑M1/1‑77 (MQ=255)
                                                                                     cctggtaaaagcaggaatgatgccatttgcttgttatttcttcattttcataaacgTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCGGTTATGATTGGTGGTATTg                                                         >  2:548458‑M1/57‑149 (MQ=255)
                                                                                      ctggtaaaagcaggaatgatgccatttgcttgttatttcttcattttcataaacgTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTATGACACACCGGTTTCATAGCTATACCGGTTAtgattggtggtcttca                                                        <  2:89043‑M1/94‑6 (MQ=255)
                                                                                      ctggtaaaagcaggaatgatgccatttgcttgttatttcttcattttcataaacgTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCGGTTATGATTGGTGGTATTGa                                                        <  2:188967‑M1/94‑1 (MQ=255)
                                                                                      aaTCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGctctgcattacctgttggtgtcccggttccgtggccttcagccactccgccaacaggctggctgaaatgcaatggtgc                                                        <  2:16769‑M1/149‑79 (MQ=255)
                                                                                            tttGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGctctgcattacctgttggtgtcccggttccgtggccttcagccactccgccaacaggctggctgaaatgcaatggtgccgcttt                                                  <  1:94456‑M1/149‑85 (MQ=255)
                                                                                             ttttcgttatCAAATCTTATGTCACTGTCCAAAGGGCTTTCGACACCCTTGGTGTGGGAGAATTAGCAAAACTCAATGCAGCCACGGCTACATTAGGACGAACCGGTTTCATAGCTATACCGGTTATGATTGGTTGTATTGAGCAATCa                                                 <  2:383699/142‑1 (MQ=255)
                                                                                              agcaggaatgatgccatttgcttgttatttcttcattttcataaacgTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCGGTTATGATTGGTGGTATTGAGCAatcgg                                                <  1:106240‑M1/102‑3 (MQ=255)
                                                                                               ggCGATATCAAATCGGATGACACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGctctgcattacctgttggtgtcccggttccgtggccttcagccactccgccaacaggctggctgaaatgcaatggtgccgctttttc                                               >  2:539053‑M1/1‑62 (MQ=255)
                                                                                                caggaatgatgccatttgcttgttatttcttcattttcataaacgTTGGtt                                                                                                                                                <  2:221597‑M1/6‑1 (MQ=255)
                                                                                                caggaatgatgccatttgcttgttatttcttcattttcataaacgTTGGtt                                                                                                                                                >  1:221597‑M1/46‑51 (MQ=255)
                                                                                                  ggaatgatgccatttgcttgttatttcttcattttcataaacgTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCGGTTATGATTGGTGGTATTGAGCAATCAGTAAt                                            >  2:467242‑M1/44‑149 (MQ=255)
                                                                                                  gATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGctctgcattacctgtagatatcccgggaccgtggcctgcagccactccgccaacaggctggctgaaatgccatggtgccgctttttctgc                                            >  2:70376‑M1/1‑59 (MQ=255)
                                                                                                      tCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGctctgcattacctgttggtgtcccggttccgtggccttcagccactccgccaacaggctggctgaaatgcaatggtgccgctttttctgctgca                                        >  1:553074‑M1/1‑55 (MQ=255)
                                                                                                         aaTCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGctctgcattacctgttggtgtcccggttccgtggccttcagccactccgccaacaggctggctgaaatgcaatggtgccgctttttctgctgaagaa                                     >  1:219071‑M1/1‑52 (MQ=255)
                                                                                                                  gCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGctctgcattacctgttggtgtcccggttccgtggccttcagccactccgccaacaggctggctgaaatgcaatggtgccgctttttctgctgaagaatacccggaa                            <  1:450257‑M1/149‑107 (MQ=255)
                                                                                                                   cACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGctctgcattacctgttggtgtcccggttccgtggccttcagccactccgccaacaggctggctgaaatgcaatggtgccgctttttctgct                                           <  1:31224‑M1/133‑92 (MQ=255)
                                                                                                                   cACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGctctgcattacctgttggtgtcccggttccgtggccttcagccactccgccaacaggctggctgaaatgcaatggtgccgctttttctgct                                           >  2:31224‑M1/1‑42 (MQ=255)
                                                                                                                        atttcttcattttcataaacgTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCGGTTATGATTGGTGGTATTGAGCAATCAGTAATCATTCAGTGGGGGTGGAATGcc                      >  1:188967‑M1/22‑149 (MQ=255)
                                                                                                                          cAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGctctgcattacctgttggtgtcccggttccgtggccttcagccactccgccaacaggctggctgaaatgcaatggtgccgctttttctgctgaagaatacccggaactggcaaa                    <  1:411439‑M1/149‑115 (MQ=255)
                                                                                                                              cGGCTCTCGCAAACCTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCGGTTATGATTGGTGGTATTGAGCAATCAGTAATCATTCAGTGGGGGTGGAATGCCGCaaaa                <  2:394131/149‑1 (MQ=255)
                                                                                                                                 ctctCGAAAACCTTGGTTTGGGAGAAGGctctgcattacctgttggtgtcccggttccgtggccttcagccactccgccaacaggctggctgaaatgcaatggtgccgctttttctgctgaagaatacccggaactggcaaaggcttat             <  1:142637‑M1/149‑122 (MQ=255)
                                                                                                                                  tctcGAAAACCTTGGTTTGGGAGAAGGctctgcattacctgttggtgtcccggttccgtggccttcagccactccgccaacaggctggctgaaatg                                                                 <  2:221928‑M1/96‑70 (MQ=255)
                                                                                                                                  tctcGAAAACCTTGGTTTGGGAGAAGGctctgcattacctgttggtgtcccggttccgtggccttcagccactccgccaacaggctggctgaaatg                                                                 >  1:221928‑M1/1‑27 (MQ=255)
                                                                                                                                  tctcGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCGGTTATGATTGGTGGTATTGAGCAATCAGTAATCATTCAGTGGGGCTGGAATGCCGCAAAAGccg            >  2:267891/1‑147 (MQ=255)
                                                                                                                                      ataaacgTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGttt                                                                                         <  1:249869‑M1/61‑1 (MQ=255)
                                                                                                                                      ataaacgTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGttt                                                                                         >  2:249869‑M1/8‑68 (MQ=255)
                                                                                                                                         aaCCTTGGTTTGGGAGAAGGctctgcattacctgttggtgtcccggttccgtggccttcagccactccgccaacaggctggctgaaatgcaatggtgccgctttttctgctgaagaatacccggaactggcaaaggcttatccgacaaa     <  2:367582‑M1/149‑130 (MQ=255)
                                                                                                                                            cTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCGGTTATGATTGGTGGTATTGAGCAATCAGTAATCATTCAGTGGGGGTGGAATGCCGCAAAAGCATCTGCCTCTgg  >  2:243857/1‑149 (MQ=255)
                                                                                                                                             |                                                                                                                                                   
TGGTGGAGGCATCCGGTGCCAGCCTGGAGAAGGGGCGGCACGACCAGTTACTTACCGCACTTCGCGCGCTGCTGTTAAGCCGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCGGTTATGATTGGTGGTATTGAGCAATCAGTAATCATTCAGTGGGGGTGGAATGCCGCAAAAGCATCTGCCTCTGG  >  NC_000913_3_pae_tpiA/1207649‑1207937

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.