New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913_3_bme_pgi = 12078057 (0.210)16 (0.540) 14/248 1.2 69.3% coding (305/630 nt) stfP e14 prophage; uncharacterized protein
?NC_000913_3_bme_pgi = 1209602 8 (0.270)pseudogene (18/501 nt) stfE pseudogene, e14 prophage; side tail fiber protein fragment family;Phage or Prophage Related

ATCCGGTGCCAGCCTGGAGAAGGGGCGGCACGACCAGTTACTTACCGCACTTCGCGCGCTGCTGTTAAGCCGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCGGTTATGATTGGTGGTATTGAGCAATCAGTAATCATTCAGTGGGGGTGGAATGCCGCAAAAGCATCTGCCTCTGGGGGGGATGGAAATACA  >  NC_000913_3_bme_pgi/1207659‑1207953
                                                                                                                                                  |                                                                                                                                                    
aTCCGGTGCCAGCCTGGAGAAGGGGCGCAACGTTCTTTTACTTACCGCACTTCGCGCGCTGCTGTTAAGCCGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGGGAGGGct                                                                                                                                                    >  2:430032‑M1/1‑147 (MQ=255)
 tCCGGTGCCAGCCTGGAGAAGGGGCGGCACGACCAGTTACTTACCGCACTTCGCGCGCTGCTGTTAAGCCGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAATCCTTGGTTTGGGAGAAGGctc                                                                                                                                                   <  2:526583‑M1/149‑4 (MQ=255)
   cGGTGCCAGCCTGGAGAAGGGGCGGCACGACCAGTTACTTACCGCACTTCGCGCGCTGCTGTTAAGCCGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGctctg                                                                                                                                                 <  1:28804‑M1/149‑6 (MQ=255)
     gTGCCAGCCTGGAGAAGGGGCGGCACGACCAGTTACTTACCGCACTTCGCGCGCTGCTGTTAAGCCGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGctctgca                                                                                                                                               >  2:153964‑M1/1‑142 (MQ=255)
      tGCCAGCCTGGAGAAGGGGCGGCATGACCAGTTACTTACCGCACTTCGCGCGCTGCTGTTAAGCCGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCAAAAc                                                                                                                                              <  2:598240/149‑1 (MQ=255)
       tacgttgtagacctgtgttctggtcatttgttgatacgcgtacatagccaataagcatgttaaatccccctggtaaaagcaggaatgatgccatttgcttgttatttcttcattttcataaacgTTGGTTTGGGAGAAGGAGCAAAACt                                                                                                                                             <  2:269801‑M1/25‑1 (MQ=255)
              tagatctgtgttctggtcatttgttgatacgcgtacatagccaataagcatgttaaatccccctggtaaaagcaggaatgatgccatttgcttgttatttcttcattttcataaacgTTGGTTTGGGAGAAGGAGCAAAACTCAATgca                                                                                                                                      >  1:471200‑M1/118‑149 (MQ=255)
              tGGAGAAGGGGCGGCACGACCAGTTACTTACCGCACTTCGCGCGCTGCTGTTAAGCCGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCAATgca                                                                                                                                      >  1:304998/1‑149 (MQ=255)
              tGGAGAAGGGGCGGCACGACCAGTTACGTACCGCACTTCGCGCGCTGCTGTTAAGCCGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGctctgcattacctgtt                                                                                                                                      <  1:252243‑M1/149‑17 (MQ=255)
                   ttttgtttttgtcatttgtggatatgcttacatagccaataagcatgttaactccccctgttaaaagcaggaatgatgccatttgcttgttatttcttcattttcataaacgTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAAc                                                                                                                                 >  2:28210‑M1/113‑149 (MQ=255)
                           tggtcatttgttgatacgcgtacatagccaataagcatgttaaatccccctggtaaaagcaggaatgatgccatttgcttgttatttcttcattttcataaacgTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACAt                                                                                                                         <  1:429143‑M1/45‑1 (MQ=255)
                                   tgtggatacgcgtacatagccaataagcatgttaaatccccctggtaaaagcaggaatgatgccatttgcttgttatttcttcattttcataaacgTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGc                                                                                                                 >  1:269801‑M1/97‑149 (MQ=255)
                                         ttacCGCACTTCGCGCGCTGCTGTTAAGCCGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGATGCACCGGt                                                                                                           >  1:281842/1‑149 (MQ=255)
                                           acCGCACGTCGCGCGCTGCTGTTAAGCCGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGTGAGAAGGctctgcattacctgttggtgtcccggttccgtggccttcagccac                                                                                                         <  2:562403‑M1/149‑46 (MQ=255)
                                              gCACTTCGCGCGCTGCTGTTAAGCCGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGctctgcattacctgttggtgtcccggttccgtggccttcagccactcc                                                                                                      >  1:188727‑M1/1‑101 (MQ=255)
                                                  ttCGCGCGCTGCTGTTAAGCCGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGctctgcattacctgttggtgtcccggttccgtggccttcagccactccgcca                                                                                                  >  1:146599‑M1/1‑97 (MQ=255)
                                                   ttgtgtgctgttgttaatCCGCAAGAATCCGTTTGGCGAGATCAAATCGGATCGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGctctgcattacctgttggtgtcccggttccgtggccttcagccactccgccaa                                                                                                 <  1:466787‑M1/143‑54 (MQ=255)
                                                                tgttaaatccccctggtaaaagcaggaatgatgccatttgcttgttatttcttcattttcataaacgTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCGGTTATGAt                                                                                    <  2:345057‑M1/82‑1 (MQ=255)
                                                                 gttaaatccccctggtaaaagcaggaatgatgccatttgcttgttatttcttcattttcataaacgTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCGGTCATGAtt                                                                                   <  2:556145‑M1/83‑1 (MQ=255)
                                                                   aGCCGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGctctgcattacctgttggtgtcccggttccgtggccttcagccactccgccaacaggctggctgaaatg                                                                                 >  1:75712‑M1/1‑80 (MQ=255)
                                                                    aaatcctgcaggtaaaagcaggaaggtggcgatttgattgttatttcttcatttctataaacgTTGGTTTGGGAGAAGGCGCAAAACTCAATGCAGCAACGGCTACATTAGGACTCACCGGTTTCATAGCTATACCGGTTATGATtggt                                                                                >  2:562958‑M1/64‑149 (MQ=255)
                                                                                 aaaagcaggaatgatgccatttgctttttatttcttcatttttattaacgTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCGGTTATGATTGGTGGTATTGAGCAAt                                                                   >  1:13948‑M1/51‑149 (MQ=255)
                                                                                    tGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGATGCACCGGTTTCATAGCTATACCGGTTATGATTGGTGGTATTGAGCAATCAg                                                                <  2:281842/149‑1 (MQ=255)
                                                                                     gcaggaatgatgccatttgcttgttatttcttcattttcataaacgTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCGGTTATGATTGGTGGTATTGAGCAATCAGt                                                               <  2:100887‑M1/103‑1 (MQ=255)
                                                                                            tCAAATCGGATGGCACTGTGCAAACGTCTCTCGAAAACCTTGGTTTGGGAGAAGGctctgcattacctgttggtgtcccggttccgtggccttcagccactccgccaacaggctggctgaaatgcaatggtgccgctttttctgctgaa                                                        <  1:317185‑M1/149‑95 (MQ=255)
                                                                                            tCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGctctgcattacctgttggtgtcccggttccgtggccttcagccactccgccaacaggctggctgaaatgcaatggtgccgctttttctgctgaa                                                        <  1:279311‑M1/149‑95 (MQ=255)
                                                                                            gtaaattGGATGGCACAGTGGAAACGGCTCTCGAAAACCTTGGTTGGGGATAAGGctctgcattacctgttggtgtcccggttccgtggccttcagccactccgccaacaggctggctgaaatgcaatggtgccgctttttctgctgaa                                                        <  2:437065‑M1/147‑95 (MQ=255)
                                                                                                   atttgcttgttagttcttgagtttcataaacgTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCGGTTATGATTGGTGGTATTGAGCAATCAGTAATCATTCAGTggg                                                 >  2:582416‑M1/33‑149 (MQ=255)
                                                                                                    tttgcttgttatttcttcattttcataaacgTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCGGTTATGATTGGTGGTATTGAGCAATCAGTAATCATTCAGTgggg                                                >  1:345057‑M1/32‑149 (MQ=255)
                                                                                                              tgCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGctctgcattacctgttggtgtcccggttccgtggccttcagccactccgccaacaggctggctgaaatgcaatggtgccgctttttctgctgaagaatacccggaactggca                                      >  1:410375‑M1/1‑37 (MQ=255)
                                                                                                                 aaaCGGCTCTCGAAAACCTTGGTTTGGGAGAAGGctctgcattacctgttggtgtcccggttccgtggccttcagccactccgccaacaggctggctgaaatgcaatggtgccgctttttctgctgaagaatacccggaactggcaaag                                   >  1:399767‑M1/1‑34 (MQ=255)
                                                                                                                        tctcGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCGGTTATGATTGGTGGTATTGAGCAATCAGTAATCATTCAGTGGGGGTGGAATGCCGCAAAAGCAt                            <  2:304998/149‑1 (MQ=255)
                                                                                                                         ctcGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCGGTTATGATTGTTGGTATTGAGCAATCAGTAATCATTCAGTGGGGGTGGAATGCCGCAAAAGCATc                           <  2:595592/149‑1 (MQ=255)
                                                                                                                          tcGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCGGTTATGATTGGTGGTATTGAGCAATCAGTAATCATTCAGTGGGGGTGGAATTCCGCAAAAGCATCt                          >  1:98854/1‑149 (MQ=255)
                                                                                                                            ataaacgTTGGTTTGGGAGAAGGGGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCGGTTATGATTGGTGGTATTGAGCAATCAGTAATCATTCCGTGTGGGTGGAATGCCGAAAAcgcatatcc                        <  2:539016‑M1/142‑5 (MQ=255)
                                                                                                                            ataaacgTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCGGTTATGATTGGTGGTATTGAGCAATCAGTAATCATTCAGTGGGGGTGGAATGCCGCAAAAGCATCTGc                        <  2:277054‑M1/142‑1 (MQ=255)
                                                                                                                                  cTTGGTTTGGGAGAAGGctctgcattacctgttggtgtcccggttccgtggccttcagccactccgccaacaggctggctgaaatgcaatggtgccgctttttctgctgaagcataccgnggactggcaaagtctgntacgacagataa                  >  1:165647‑M1/1‑17 (MQ=255)
                                                                                                                                   ttGGTTTGGGAGAAGGctttgcattacctgttggtgtcccggttccgtggccttcagccactccgccaacaggctggctgaaatgcaatggtgccgctttttctgctgaagaatacccggaactggcaaaggcttatccgacaaataaa                 <  1:153964‑M1/149‑134 (MQ=255)
                                                                                                                                     GGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTTGGACGCACCGGTTTCATAGCTATACCGGTTATGATTGGTGGTATTGAGCAATCAGTAATCATTCAGTGGGGGTGGAATGCCGCAAAAGCATCTGCCTCTggggg               >  2:474873‑M1/1‑149 (MQ=255)
                                                                                                                                     GGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCGGTTATGATTGGTGGTATTGAGCAATCAGTAATCATTCAGTGGGGGTGGAATGCCGCAAAAGCATCTGCCTCTggggg               >  2:523385‑M1/1‑149 (MQ=255)
                                                                                                                                       tttGGGAGAAGGctctgcattacctgttggtgtccgggttccgtggccttcagccactccgccaacaggctggctgaaatgcaatggtgccgctttttctgctgaagaatacccggaactggcaaaggcttatccgacaaataaattgc             <  2:266199‑M1/149‑138 (MQ=255)
                                                                                                                                        ttGGGAGAAGGcgctgaattacctgtgggtgtcccggttccgtggccttcagccactccgccaacaggctggctgaaatgcaatggtgccgctttttctgctgaagaatacccggaactggcaaaggcttatccgacaaataaattgcc            <  1:230429‑M1/149‑139 (MQ=255)
                                                                                                                                          gggAGAAGGctctgcattacctgttggtgtcccggttccgtggccttcagccactccgccaacaggctggctgaaatgcaatggtgccgctttttctgctgaagaatacccggaactggcaaaggcttatccgacaaataaattgcctg          >  1:110993‑M1/1‑9 (MQ=255)
                                                                                                                                                  gctctgcattacctgttggtgtcccggttccgtggccttcagccactccgccaacaggctggctgaaatgcaatggtgccgctttttctgctgaagaatacccggaactggcaaaggcttatccgacaaataaattgcctgatttacgt  <  2:146599‑M1/149‑149 (MQ=255)
                                                                                                                                                  GAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAACTATACCGGTTATGATTGGTGGTATTGAGCAATCAGTAATCATTCAGTGGGGGTGGAATGCCGCAAACGCATCTGCGTCTGGGGGGGATGGaaatcca  <  2:97668‑M1/149‑4 (MQ=255)
                                                                                                                                                  |                                                                                                                                                    
ATCCGGTGCCAGCCTGGAGAAGGGGCGGCACGACCAGTTACTTACCGCACTTCGCGCGCTGCTGTTAAGCCGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCGGTTATGATTGGTGGTATTGAGCAATCAGTAATCATTCAGTGGGGGTGGAATGCCGCAAAAGCATCTGCCTCTGGGGGGGATGGAAATACA  >  NC_000913_3_bme_pgi/1207659‑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.