New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913_3_pae_tpiA = 120780526 (0.460)14 (0.280) 13/254 3.8 32.5% coding (305/630 nt) stfP e14 prophage; uncharacterized protein
?NC_000913_3_pae_tpiA = 1209602 35 (0.690)pseudogene (18/501 nt) stfE pseudogene, e14 prophage; side tail fiber protein fragment family;Phage or Prophage Related
Rejected: Coverage evenness skew score above cutoff.

CATCCGGTGCCAGCCTGGAGAAGGGGCGGCACGACCAGTTACTTACCGCACTTCGCGCGCTGCTGTTAAGCCGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCGGTTATGATTGGTGGTATTGAGCAATCAGTAATCATTCAGTGGGGGTGGAATGCCGCAAAAGCATCTGCCTCTGGGGGGGATGGAAATACA  >  NC_000913_3_pae_tpiA/1207658‑1207953
                                                                                                                                                   |                                                                                                                                                    
cagcgcattacgttagagatctgtgttctggtcatttgttgatacgcgtacatagcgaataagcatgttaaatccccctggtaaaagcaggaatgatgccatttgcttgttatttcttcattttcataaacgTTGGTTTGGGAGAAGGa                                                                                                                                                     >  2:65458‑M1/133‑149 (MQ=255)
  tCCGGTGCCAGCCTGGAGAAGGGGCGGCACGACCAGTTACTTACCGCACTTCGCGCGCTGCTGTTAAGCCGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGc                                                                                                                                                   >  2:221722/1‑149 (MQ=255)
   cgcattacgttgtagatcttgtttctggtcatttgttgatacgcgtacatagccaataagcatgttaaatcccccttgtaaaagcaggaatgatgccatttgcttgttatttcttcattttcataaccgTTGGTTTGGGAGAaggatct                                                                                                                                                  >  2:145204‑M1/130‑146 (MQ=255)
     ggTGCCAGCCTGGAGAAGGGGCGGCACGACCAGTTACTTACCGCACTTCGCGCGCTGCTGTTAAGCCGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAggctctgc                                                                                                                                                <  2:748844/149‑7 (MQ=255)
           gttgtagatctgtgttctggtcatttgttgatacgcgtacataggcaataaacatgtcatattgtcctggtaaaagcaggaatgatgacatttgcttgttatttcttcattttcataaacgTTGGTTTGGGTGATGGAGTAAAActcga                                                                                                                                          <  2:622605‑M1/28‑3 (MQ=255)
             ccTGGAGAAGGGGCGGCACGACCAGTTACTTACCGCACTTCGCGCGCTGCTGTTAAGCCGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCAATg                                                                                                                                        <  1:651372/149‑1 (MQ=255)
              cTGGAGAAGGGGCGGCACGACCAGTTACTTACCGCACTTCGCGCGCTGCTGTTAAGCCGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCAATgc                                                                                                                                       >  1:930864/1‑149 (MQ=255)
                ggAGAAGGGGCGGCACGACCAGTTACTTACCGCACTTCGCGCGCTGCTGTTAAGCCGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCAATgcag                                                                                                                                     <  2:930864/149‑1 (MQ=255)
                  atctgtgttctggtcatttgttgatacgcgtacatagccaataagcatgttaaatccccctggtaaaagcaggaatgatgccatttgcttgttatttcttcattttcataaacgTTGGTTTGGGAGAAGGAGCAAAACTCAATgcagca                                                                                                                                   <  2:983341‑M1/35‑1 (MQ=255)
                   tctgtgttctggtcatttgttgatacgcgtacatagccaataagcatgttaaatccccctggtaaaagcaggaatgatgccatttgcttgttatttcttcattttcataaacgTTGGTTTGGGAGAAGGAGCAAAACTCCATGCAGCaa                                                                                                                                  <  2:443464‑M1/36‑1 (MQ=255)
                    aaGGGGCGGCACGACCAGTTACTTACCGCACTTCGCGCGCTGCTGTTAAGCCGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAAc                                                                                                                                 <  1:221722/149‑1 (MQ=255)
                              aCGACCAGTTACTTACCGCACTTCGCGCGCTGCTGTTAAGCCGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTa                                                                                                                       >  1:217280/1‑149 (MQ=255)
                              aCGACCAGTTACTTACCGCACTTCGCGCGCTGCTGTTAAGCCGCAAGAATCCGTTTGGCGATAACAAACGGGTTGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTGTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTa                                                                                                                       >  1:447847/1‑149 (MQ=25)
                                     gttgagacgcgtacatagccaaaacgcatgttaaacccccctgggaaaagcaggaatgaggccatttgcttgttatttcttcattttcataaacgTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCa                                                                                                                >  1:1010495‑M1/96‑149 (MQ=255)
                                     gTTACTTACCGCACTTCGCGCGCTGCTGTTAAGCCGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCa                                                                                                                >  2:827242/1‑149 (MQ=255)
                                                 aCTTCGCGCGCTGCTGTTAAGCCGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAg                                                                                                    <  2:44496/149‑1 (MQ=255)
                                                   atagccaataagcatgttaaatccccctggtaaaagcaggaatgcttcgaatttcttattgactcgttattttcataaacgTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCt                                                                                                  <  1:747023‑M1/68‑1 (MQ=255)
                                                   atagccaataagcatgttaaatccccctggtaaaagcaggaatgatgccatttgcttgttatttcttcattttcataaacgTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCt                                                                                                  <  1:176460‑M1/68‑1 (MQ=255)
                                                       cgcgcTGCTGTTAAGCCGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATAc                                                                                              <  1:908326/149‑1 (MQ=255)
                                                        caataagcatgttaaatccccctggtaaaagcaggaatgatgccatttgcttgtcatttcttcattttcataaacgATAGCAAGGAACAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATAcc                                                                                             >  2:565763‑M1/77‑149 (MQ=255)
                                                          gctgctGTTAAGCCGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCgg                                                                                           >  1:572229/1‑149 (MQ=255)
                                                          gctgctGTTAAGCCGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCgg                                                                                           <  1:482629/149‑1 (MQ=255)
                                                                     gacGCAAGGAGCCGTGTGCGGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCGGTTATGATtggt                                                                                <  2:245668/147‑1 (MQ=38)
                                                                       cGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCGGTTATGGTtggtgg                                                                              >  2:799318/1‑149 (MQ=255)
                                                                          cccctggtaaaagcaggaatgatgccatttgcttgttatttcttcattttcataaacgTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCGGTTATGATTGGTGGTAt                                                                           >  1:656966‑M1/59‑149 (MQ=255)
                                                                           ccctggtaaacgcaggaatgatgccatttgctgcttatttcttcattttcataaacgTTGGTTTGTGAGAAGAAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCGGTTATGATTGGTGGTAtt                                                                          >  1:596467‑M1/58‑149 (MQ=255)
                                                                             aaTCCGGTGGGCGATAGCAAAGCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCGGTTATGATTGGTGGTATTGa                                                                        <  1:220717/149‑1 (MQ=255)
                                                                              aTCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCGGTTATGATTGGTGGTATTGAg                                                                       >  1:356868/1‑149 (MQ=255)
                                                                                gtaaaagcaggaatgatgccatttgcttgttatttcttcattttcataaacgTTGGGTGGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCGGTTATGATTGGTGGTATTGCGCa                                                                     <  2:512965‑M1/97‑1 (MQ=255)
                                                                                 taaaagcaggaatgatgccatttggttgttatttcttcattttcataaacgTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGTACGCACCGGTTTCATAGCTATACCGGTTATGATTGGTGGTATTGAGCaa                                                                    >  2:245562‑M1/52‑149 (MQ=255)
                                                                                 cGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCaaaa                                                                                                                                               >  1:151620/1‑74 (MQ=255)
                                                                                 cGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCaaaa                                                                                                                                               <  2:151620/74‑1 (MQ=255)
                                                                                  gTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCGGTTATGATTGGTGGTATTGAGCAAt                                                                   >  1:711265/1‑149 (MQ=255)
                                                                                  aaaagcaggaatgatgccatttgcttgttatttcttcattttcataaacgTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCGGTTATGATTGGTGGTATTGAGCAAt                                                                   >  2:988447‑M1/51‑149 (MQ=255)
                                                                                   tttGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCAATGCACCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCGGTTATGATTGGTGGTATTGAGCAATc                                                                  >  1:490417/1‑149 (MQ=255)
                                                                                    ttGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAACCCTTGGTTGGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGAACGCACCCGTTTCATAGCTATACCGGTTATGATTGGTGGTATTGAGCAATCa                                                                 <  2:423921/149‑1 (MQ=255)
                                                                                    ttGGCAATATGAAATTCGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTGTCATAGCTATCCCGGTTATGATTGGTGGTATTGAGCAATCa                                                                 <  2:423910/149‑1 (MQ=38)
                                                                                    aagcaggaatgatgccatttgcttgttatttcttcattttcataaacgTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCGGTTATGATTGGTGGTATTGAGCAATCa                                                                 >  2:700364‑M1/49‑149 (MQ=255)
                                                                                     ttgCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCGGTTATGATTGGTGGTATTGAGCAATCAg                                                                <  2:322619/147‑1 (MQ=255)
                                                                                     agcaggaatgatgccatttgcttgttatttcttcattttcataaacgTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAAGACGCACCGGTTTCATAGCTATACTGGTTATGACTGGTGGTATTGAGCAATCAg                                                                <  1:870720‑M1/102‑1 (MQ=255)
                                                                                       gCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCg                                                                                            >  1:882087/1‑119 (MQ=255)
                                                                                       gCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCg                                                                                            <  2:882087/119‑1 (MQ=255)
                                                                                       caggaatgatgccatttgcttgttatttcttcattttcataaacgTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCGGTTATGATTGGTGGTATTGAGCAATCAGTa                                                              >  1:431994‑M1/46‑149 (MQ=255)
                                                                                                    atttgcttgttatttcttcattttcataaacgTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCGGTTATGATTGGTGGTATTGAGCAATCAGTAATCATTCAGTggg                                                 >  1:983341‑M1/33‑149 (MQ=255)
                                                                                                     tttgcttgttatttcttcattttcataaacgTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCGGTTATGATTGGTGGTATTGAGCAATCAGTAATCATTCAGTgggg                                                >  1:252544‑M1/32‑149 (MQ=255)
                                                                                                         cttgttatttctgcattttcataaacgTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCGGTTATGATTGGTGGTATTGAGCAATCAGTAATCATTCAGTGGGGGTgg                                            >  1:324541‑M1/28‑149 (MQ=255)
                                                                                                               atttcttcattttcataaacgTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCGGTTATGATTAGTACCGACGTTCTATCAGTAATCATTCAGTGGGGGTGGAATGcc                                      >  1:497160‑M1/22‑149 (MQ=255)
                                                                                                                    ttcattttcataaacgTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCGGTTATGATTGGTGGTATTGAGCAATCAGTAATCATTCAGTGGGGGTGGAATGCCGCaaa                                 <  1:11309‑M1/133‑1 (MQ=255)
                                                                                                                                  ccTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCGGTTATGATTGGTGGTATTGAGCAATCAGTAATCATTCAGTGGGGGTGGAATGCCGCAAAAGCATCTGCCTCTg                   <  1:827242/149‑1 (MQ=255)
                                                                                                                                   gTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCGGTTATGATTGGTGGTATTGAGCAATCAGTAATCATTCAGTGGGGGTGGAATGCCGCAAAAGCATCTGCCTCTgg                  <  1:418647‑M1/148‑1 (MQ=255)
                                                                                                                                     TGGTTGGGGAGAAGGAGCAAAACGCAATGCAGCAACGGCTACATTAGGAAACCCTCGCCTCATAGCTATACCGGTTATGATTGGTGGTAGTGAGCAATCAGTAATCATTCAGTGGGGGTGGAATGCCGCAAAAGCATCTGCCTCTgggg                >  1:233062‑M1/1‑149 (MQ=255)
                                                                                                                                           GGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTa                                                                                                                            <  2:627228‑M1/35‑1 (MQ=255)
                                                                                                                                           GGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTa                                                                                                                            >  1:627228‑M1/1‑35 (MQ=255)
                                                                                                                                            GGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACAGGTTATGATTGGTGGTATTGAGAAATAAGTAATCATTCAATAGGGGTGGAATGCCGCAAAAGACTCTGACTCTGGGGGGGATgg         <  2:698299‑M1/149‑1 (MQ=255)
                                                                                                                                             GAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCGGTTATGATTGGTGGTATTGAGCAATCAGTAATCATTCAGTGGGGGTGGAATGCCGCAAAAGCATCTGCCTCTGGGGGGGATGGa        <  2:328559‑M1/149‑1 (MQ=255)
                                                                                                                                             GAGAAGGAGCAAAACGCAATGCAGCAACGGCTACATTAGGACGCACCGGGTTCATAGCTATACCAGTTATGATTGGTTGTTTTGAGCAATCAGTAATCAGTCAGTGGGGGTGGAATGCCGCAAAAGCATCTGCCTCTGGGGGGGATGGa        >  2:349935‑M1/1‑149 (MQ=255)
                                                                                                                                               GAAGGGGCAAAACTGAATGCGGGAACGGCTACATTAGGACGCCCCGGTTTCATAGCTATACCGGTTATGATTGGTGGTATTGAGCAATCAGTAATCATTCAGGGGGGGTGGAATGCCGCAAAAGCATCTGCCTCTGGGGGGGATGGaaa      >  2:572229‑M1/1‑149 (MQ=255)
                                                                                                                                               GAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCGGTTATGATTGGTGGTATTGAGCAATCAGTAATCATTCAGTGGGGGTGGAATGCCGCAAAAGCATCTGCCTCTGGGGGGGATGGaaa      >  2:356868‑M1/1‑149 (MQ=255)
                                                                                                                                               GAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCGGTTATGATTGGTGGTATTGAGCAATCAGTAATCATTCAGTGGGGGTGGAATGCCGCAAAAGCATCTGCCTCTGGGGGGGATGGaaa      >  1:753721‑M1/1‑149 (MQ=255)
                                                                                                                                                   GAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCGGTTATGATTGGTGGTATTGAGCAATCAGTAATCATTCAGTGGGGGTGGAATGCGGCAAAAGCATCTGCCTCTGGGAGGCATGGAAATACa  <  2:536676‑M1/149‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                    
CATCCGGTGCCAGCCTGGAGAAGGGGCGGCACGACCAGTTACTTACCGCACTTCGCGCGCTGCTGTTAAGCCGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCGGTTATGATTGGTGGTATTGAGCAATCAGTAATCATTCAGTGGGGGTGGAATGCCGCAAAAGCATCTGCCTCTGGGGGGGATGGAAATACA  >  NC_000913_3_pae_tpiA/1207658‑1207953

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