New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? AP009048 = 143225244 (0.750)3 (0.060) 3/90 NT 7.4% coding (1490/3363 nt) stfR predicted tail fiber protein
?AP009048 = 1432258 36 (0.700)coding (1496/3363 nt) stfR predicted tail fiber protein

CAGTACGTCTGAAACGCTGGCGGCAACGCCAAAGGCAGTAAAATCAGCCTATGACAATGCAGAGAAACGTCTGCAGAAAGACCAGAACGGCGCTGATATACCCGATAAGGGATGCTTCCTGAACAACATTAACGCGGT  >  AP009048/1432190‑1432327
                                                                    |                                                                     
cAGTACGTCTGAAACGCTGGCGGCAACGCCAAAGGCAGTAAAATCAGCCTATGACAATGCAGAGAAACGtc                                                                     <  1:997153‑M2/71‑3 (MQ=255)
  gTACGTCTGAAACGCTGGCGGCAACGCCAAAGGCAGTAAAATCAGCCTATGACAATGCAGAGAAACGtctg                                                                   >  1:3053952‑M2/1‑67 (MQ=255)
   tACGTCTGAAACGCTGGCGGCAACGCCAAAGGCAGTAAAATCAGCCTATGACAATGCAGAGAAACGt                                                                      <  1:4934580‑M2/67‑2 (MQ=255)
      gTCTGAAACGCTGGCGGCAACGCCAAAGGCAGTAAAATCAGCCTATGACAATGCAGAGAAACGTCTGCAGa                                                               <  1:491775/71‑1 (MQ=255)
      gTCTGAAACGCTGGCGGCAACGCCAAAGGCAGTAAAATCAGCCTATGACAATGCAGAGAAACGTCTGCAGa                                                               <  1:3451987/71‑1 (MQ=255)
       tCTGAAACGCTGGCGGCAACGCCAAAGGCAGTAAAATCAGCCTATGACAATGCAGAGAAACGTCTGCAGaa                                                              <  1:3133152/71‑1 (MQ=255)
       tCTGAAACGCTGGCGGCAACGCCAAAGGCAGTAAAATCAGCCTATGACAATGCAGAGAAACGTCTGCAGaa                                                              <  1:2116584/71‑1 (MQ=255)
           aaaCGCTGGCGGCAACGCCAAAGGCAGTAAAATCAGCCTATGACAATGCAGAGAAACga                                                                      >  1:2063060/1‑58 (MQ=255)
           aaaCGCTGGCGGCAACGCCAAAGGCAGTAAAATCAGCCTATGACAATGCAGAGAAACGTCTGCAGAAAg                                                            <  1:768518/69‑1 (MQ=255)
           aaaCGCTGGCGGCAACGCCAAAGGCAGTAAAATCAGCCTATGACAATGCAGAGAAACGTCTGCAGAAAg                                                            <  1:422733/69‑1 (MQ=255)
            aaCGCTGGCGGCAACGCCAAAGGCAGTAAAATCAGCCTATGACAATGCAGAGAAACGTCTGCAGAAAGa                                                           >  1:3876677/1‑69 (MQ=255)
             aCGCTGGCGGCAACGCCAAAGGCAGTAAAATCAGCCTATGACAATGCAGAGAAACGTCTGCAGAAAGAc                                                          >  1:854352/1‑69 (MQ=255)
              cGCTGGCGGCAACGCCAAAGGCAGTAAAATCAGCCTATGACAATGCAGAGAAACGTCTGCAGAAAGAc                                                          >  1:718918/1‑68 (MQ=255)
               gCTGGCGGCAACGCCAAAGGCAGTAAAATCAGCCTATGACAATGCAGAGAAACGTCTGCAGAAAGACCAg                                                       >  1:1756116/1‑70 (MQ=255)
                cTGGCGGCAACGCCAAAGGCAGTAAAATCAGCCTATGACAATGCAGAGAAACGTCTGCAGAAAGACCAGaa                                                     <  1:1197486/71‑1 (MQ=255)
                 tGGCGGCAACGCCAAAGGAAGTAAAATCAGCCTATGACAATGCAGAGAAACGTCTGCAGAAAGACCAGaa                                                     <  1:3783373/70‑1 (MQ=255)
                   gcggcAACGCCAAAGGCAGTAAAATCAGCCTATGACAATGCAGAGAAACGTCTGCAGAAAGACCAGAACgg                                                  >  1:278715/1‑71 (MQ=255)
                   gcggcAACGCCAAAGGCAGTAAAATCAGCCTATGACAATGCAGAGAAACGTCTGCAGAAAGACCAGAACgg                                                  >  1:4126340/1‑71 (MQ=255)
                    cggcAACGCCAAAGGCAGTAAAATCAGCCTATGACAATGCAGAGAAACGTCTGCAGAAAGACCAGAACGgc                                                 <  1:1369425/71‑1 (MQ=255)
                      gcAACGCCAAAGGCAGTAAAATCAGCCTATGACAATGCAGAGAAACGTCTGCAGAAAGACCAGAACGgcgc                                               <  1:4997133/71‑1 (MQ=255)
                       cAACGCCAAAGGCAGTAAAATCAGCCTATGACAATGCAGAGAAACGTCTGCAGAAAGa                                                           >  1:1458761/1‑58 (MQ=255)
                       cAACGCCAAAGGCAGTAAAATCAGCCTATGACAATGCAGAGAAACGTCTGCAGAAAGACCAGAACGGCGCt                                              >  1:3773156/1‑71 (MQ=255)
                        aaCGCCAAAGGCAGTAAAATCAGCCTATGACAATGCAGAGAAACGTCTGCAGAAAg                                                            <  1:4938359/56‑1 (MQ=255)
                        aaCGCCAAAGGCAGTAAAATCAGCCTATGACAATGCAGAGAAACGTCTGCAGAAAGACCAGaa                                                     <  1:2199438/63‑1 (MQ=255)
                         aCGCCAAAGGCAGTAAAATCAGCCTATGACAATGCAGAGAAACGTCTGCAg                                                                >  1:1634737/1‑51 (MQ=255)
                          cGCCAAAGGCAGTAAAATCAGCCTATGACAATGCAGAGAAACGTCTGCAGAAAGACCAGAACGGCGCTGa                                            >  1:5247166/1‑70 (MQ=255)
                          cGCCAAAGGCAGTAAAATCAGCCTATGACAATGCAGAGAAACGTCTGCAGAAAGACCAGAACGGCGCTGa                                            >  1:1961393/1‑70 (MQ=255)
                            ccAAAGGCAGTAAAATCAGCCTATGACAATGCAGAGAAACGt                                                                      <  1:4297932‑M2/42‑2 (MQ=255)
                               aaGGCAGTAAAATCAGCCTATGACAATGCAGAGAAACGTCTGCAGAAAGACCAGAACGGCGCTGATATAcc                                      >  1:1043800/1‑71 (MQ=255)
                                 ggCAGTAAAATCAGCCTATGACAATGCAGAGAAACGTCTGCAGa                                                               <  1:1646337/44‑1 (MQ=255)
                                   cAGTAAAATCAGCCTATGACAATGCAGAGAAACGTCTGCAGAAAGACCAGAACGGCGCTGATATACCCGAt                                  <  1:2110134/71‑1 (MQ=255)
                                   cAGTAAAATCAGCCTATGACAATGCAGAGAAACGTCTGCAGAAAGACCAGAACGGCGCTGATATACCCGAt                                  <  1:4250907/71‑1 (MQ=255)
                                     gTAAAATCAGCCTATGACAATGCAGAGAAACGTCTGCAGAAAGACCAGAACGGCGCTGATATACCCGATaa                                <  1:2134304/71‑1 (MQ=255)
                                      tAAAATCAGCCTATGACAATGCAGAGAAACGTCTGCAGAAAGACCAGAACGGCGCTGATATACCCGATAAg                               <  1:3231179/71‑1 (MQ=255)
                                       aaaaTCAGCCTATGACAATGCAGAGAAACGtctgc                                                                  <  1:3789544‑M2/35‑6 (MQ=255)
                                       aaaaTCAGCCTATGACAATGCAGAGAAACGTCTGCAGAAAGACCAGAACGGCGCTGATATACCCGATAAgg                              >  1:2651565/1‑71 (MQ=255)
                                        aaaTCAGCCTATGACAATGCAGAGAAACGTCTGCAGAAAGACCAGAACGGCGCTGATATACCCGATAAgg                              >  1:3750443/1‑70 (MQ=255)
                                         aaTCAGCCTATGACAATGCAGAGAAACGTCTGCAGAAAGa                                                           <  1:2830695/40‑1 (MQ=255)
                                             aGCCTATGACAATGCAGAGAAACGTCTGCAGAAAGACCAGAACGGCGCTGATATACCCGATAAGGGATGc                         <  1:5229577/70‑1 (MQ=255)
                                             aGCCTATGACAATGCAGAGAAACGTCTGCAGAAAGACCAGAACGGCGCTGATATACCCGATAAGGGATGCt                        >  1:669063/1‑71 (MQ=255)
                                             aGCCTATGACAATGCAGAGAAACGTCTGCAGAAAGACCAGAACGGCGCTGATATACCCGATAAGGGATGCt                        <  1:2913592/71‑1 (MQ=255)
                                              gCCTATGACAATGCAGAGAAACGTCTGCAGAAAGACCAGAACGGCGCTGa                                            <  1:3487961/50‑1 (MQ=255)
                                               ccTATGACAATGCAGAGAAACGTCTGCAGAAAGACCAGAACGGCGCTGATATACCCGATAAGGGATGCTTc                      >  1:1308609/1‑71 (MQ=255)
                                                   tGACAATGCAGAGAAACGTCTGCAGAAAGACCAGAACGGCGCTGATATACCCGATAAGGGATGCTTCCTg                   >  1:175629/1‑70 (MQ=255)
                                                     aCAATGCAGAGAAACGTCTGCAGAAAGACCAGAACGGCGCTGATATACCCGATAAGGGATGCTTCCTGaa                 >  1:2274107/1‑70 (MQ=255)
                                                      cAATGCAGAGAAACGTCTGCAGAAAGACCAGAACGGCGCTGa                                            <  1:4032451/42‑1 (MQ=255)
                                                       aaTGCAGAGAAACGtctgcattgtcataggctgattttactgcctttggcgttgccg                            >  1:4611856‑M2/1‑14 (MQ=255)
                                                       aaTGCAGAGAAACGtctgcattgtcataggctgattttactgcctttggcg                                  >  1:3522782‑M2/1‑14 (MQ=255)
                                                       aaTGCAGAGAAACGTCTGCAGAAAGACCAGAACGGCGCTGATATACCCGATAAgg                              <  1:434015/55‑1 (MQ=255)
                                                         tGCAGAGAAACGTCTGCAGAAAGACCAGAACGGCGCTGATATACCCGATAAGGGATGCt                        >  1:2453062/1‑59 (MQ=255)
                                                         tGCAGAGAAACGTCTGCAGAAAGACCAGAACGGCGCTGATATACCCGATAAGGGATGCTTCCTGAacaaca            >  1:4503029/1‑71 (MQ=255)
                                                          gCAGAGAAACGtctgcattgtcataggctgattttactgccttt                                      >  1:1844267‑M2/1‑11 (MQ=255)
                                                              agaAACGTCTGCAGAAAGACCAGAACGGCGCTGATATACCCGATAAGGGATGCTTCCTGAACAACATTAAc       >  1:2613313/1‑71 (MQ=255)
                                                               gaAACGTCTGCAGAAAGACCAGAACGGCGCTGATATACCCGATAAGGGATGCTTCCTGAACAACATTAAcg      <  1:351073/71‑1 (MQ=255)
                                                                   cGTCTGCAGAAAGACCAGAACGGCGCTGATATACCCGATAAGGGATGCTTCCTGAACAACATTAACGCGGt  <  1:4842217/71‑1 (MQ=255)
                                                                    |                                                                     
CAGTACGTCTGAAACGCTGGCGGCAACGCCAAAGGCAGTAAAATCAGCCTATGACAATGCAGAGAAACGTCTGCAGAAAGACCAGAACGGCGCTGATATACCCGATAAGGGATGCTTCCTGAACAACATTAACGCGGT  >  AP009048/1432190‑1432327

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