New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913 1207790 =6 (0.320)13 (0.910) 11/100 0.2 77.6% coding (290/630 nt) stfP e14 prophage; uncharacterized protein
?NC_000913 1209619 = 3 (0.210)pseudogene (1/501 nt) stfE pseudogene, e14 prophage; side tail fiber protein fragment family;Phage or Prophage Related

AAGCCGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCT  >  NC_000913/1207725‑1207857
                                                                 |                                                                   
aaGCCGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTG                                                                   <  1:159778‑M1/68‑1 (MQ=255)
 taaatccccctggtaaaagcaggaatgatgccatttgcttgttatttcttcattttcataaacgTTgg                                                                  >  1:116254‑M1/65‑68 (MQ=255)
 aGCCGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCACTGTGCTAACGGCTCTCGAAAACCTTGG                                                                  >  2:533695‑M1/1‑68 (MQ=255)
  gCCGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGT                                                                 <  2:220743‑M1/68‑1 (MQ=255)
  gCCGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGT                                                                 <  2:192026‑M1/68‑1 (MQ=255)
  gCCGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGT                                                                 >  1:469591‑M1/1‑68 (MQ=255)
    atccccctggtaaaagaaggaatgatgccatttgcttgttatttcttcattttcataaacgTTGGttt                                                               <  1:658010‑M1/7‑1 (MQ=255)
     tccccctggtaaaagcaggaatgatgccatttgcttgttatttcttcattttcataaacgTTGGTTTg                                                              <  1:585078‑M1/8‑1 (MQ=255)
       aaGAATCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGTCTCTCGATAACCTTGTTTTGGG                                                            <  1:563580‑M1/68‑1 (MQ=255)
             ccGTTTGGCGATATCAAATCGGATGGCACTGTGCAACCGGCTCTCGAAAACCTTGGTTTGGGAGAAGG                                                      >  2:372892‑M1/1‑68 (MQ=255)
             ccGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGG                                                      <  2:534382‑M1/68‑1 (MQ=255)
              cGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGa                                                     <  2:565899/68‑1 (MQ=255)
                   gcaggaatgatgccatttgcttgttatttcttcattttcataaacgTTGGTTTGGGAGAAGGAGCaaa                                                >  2:219436‑M1/47‑68 (MQ=255)
                    gCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGctctgca                                               <  2:522272‑M1/68‑8 (MQ=255)
                     cGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCAAAAc                                              <  1:207326/68‑1 (MQ=255)
                     aggaatgatgccatttgcttgttatttcttcattttcataaacgTTGGTTTGGGAGAAGGAGCAAAAc                                              >  2:538544‑M1/45‑68 (MQ=255)
                        tagcaaaTCGGAGGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGctctgcattac                                           <  2:454821‑M1/65‑12 (MQ=255)
                        aatgatgccatttgcttgttatttcttcattttcataaacgTTGGTTTGGGAGAAGGAGCAAAACTCa                                           <  1:478022‑M1/27‑1 (MQ=255)
                          tgatgccatttgcttgttatttcttcattttcataaacgTTGGTTTGGGAGAAGGAGCAAAACTCAAt                                         <  1:512389‑M1/29‑1 (MQ=255)
                          tCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCAAt                                         <  1:627555/68‑1 (MQ=255)
                           cAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGctctgcattacctg                                        >  2:228106‑M1/1‑54 (MQ=255)
                              gccatttgcttgttatttcttcattttcataaacgTTGGTTTGGGAGAAGGAGCAAAACTCAATgcag                                     <  2:8117‑M1/33‑1 (MQ=255)
                                    tgcttgttatttcttcattttcataaacgTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACgg                               >  2:96188‑M1/30‑68 (MQ=255)
                                        tgttatttcttcattttcataaacgTTGGTTTGGGAGAAGTAGCAAAACTCAATGCAGCAACGGCTAc                           >  1:185376‑M1/26‑68 (MQ=255)
                                        tgttatttcttcattttcataaacgTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTAc                           >  1:275566‑M1/26‑68 (MQ=255)
                                        aCTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGctctgcattacctgttggtgtcccggt                           >  1:191062‑M1/1‑41 (MQ=255)
                                        aCTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGctctgcattacctgttggtgtcccggt                           >  1:502133‑M1/1‑41 (MQ=255)
                                        aCTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGctctgcattaccagttggtgtcccggt                           >  2:475247‑M1/1‑41 (MQ=255)
                                         cTGCGCAAACGGCTCTCGAAAACCTTGGTTCGGGAGAAGGAGCAAAACTCAACGCAGCAACGGCTACa                          <  1:15160/68‑1 (MQ=37)
                                             gCAAACGGCTCTCGAAAACCGTGGTTTGGGAGAAGGctctgcattacctgttggtgtcccggttccgt                      >  2:646248‑M1/1‑36 (MQ=255)
                                                  cGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACAATAGGACGc                 >  2:523779/1‑68 (MQ=255)
                                                   cattttcataaacgTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCa                >  2:548642‑M1/15‑68 (MQ=255)
                                                              acgTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATa     >  2:106485‑M1/4‑68 (MQ=255)
                                                              acgTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATa     >  2:404344‑M1/4‑68 (MQ=255)
                                                              aCCTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATa     <  2:138094/68‑1 (MQ=255)
                                                               cgTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAg    >  2:637565‑M1/3‑68 (MQ=255)
                                                                gTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGc   >  2:156492‑M1/2‑68 (MQ=255)
                                                                cTTGGTTTGGGAGAAGGctctgcattacctgttggtgtcccggttccgtggccttcagccactccgcc   <  2:344479‑M1/68‑52 (MQ=255)
                                                                 TTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCt  >  2:414763‑M1/1‑68 (MQ=255)
                                                                 TTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGCTTCATAGCt  >  1:18670‑M1/1‑68 (MQ=255)
                                                                 |                                                                   
AAGCCGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCT  >  NC_000913/1207725‑1207857

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