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

CGACCAACAGTACGTCTGAAACGCTGGCGGCAACGCCAAAGGCAGTAAAATCAGCCTATGACAATGCAGAGAAACGTCTGCAGAAAGACCAGAACGGCGCTGATATACCCGATAAGGGATGCTTCCTGAACAACATTAAC  >  AP009048/1432183‑1432322
                                                                     |                                                                      
cGACCAACAGTACGTCTGAAACGCTGGCGGCAACGCCAAAGGCAGTAAAATCAGCCTATGACAATGCAGA                                                                        <  1:4508302‑M1/70‑1 (MQ=255)
  aCCAACAGTACGTCTGAAACGCTGGCGGCAACGCCAAAGGCAGTAAAATCAGCCTATGACAATGCAGAGaa                                                                     <  1:3472424/71‑1 (MQ=255)
       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:2116584/71‑1 (MQ=255)
              tCTGAAACGCTGGCGGCAACGCCAAAGGCAGTAAAATCAGCCTATGACAATGCAGAGAAACGTCTGCAGaa                                                         <  1:3133152/71‑1 (MQ=255)
                tGAAACGCTGGCGGCAACGCCAAAGGCAGTAAAATCAGCCTATGACAATGCAgaga                                                                      <  1:2183596/56‑1 (MQ=255)
                tGAAACGCTGGCGGCAACGCCAAAGGCAGTAAAATCAGCCTATGACAATGCAGA                                                                        >  1:1446425‑M1/1‑54 (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)
                           cggcAACGCCAAAGGCAGTAAAATCAGCCTATGACAATGCAGAcgtttctctgcatt                                                          >  1:4611856‑M1/1‑43 (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)
                                aCGCCAAAGGCAGTAAAATCAGCCTATGACAATGCAgaga                                                                      <  1:3678366/40‑1 (MQ=255)
                                aCGCCAAAGGCAGTAAAATCAGCCTATGACAATGCAGAGAAACGTCTGCAg                                                           >  1:1634737/1‑51 (MQ=255)
                                 cGCCAAAGGCAGTAAAATCAGCCTATGACAATGCAGAcgtttctctgcatt                                                          >  1:3522782‑M1/1‑37 (MQ=255)
                                 cGCCAAAGGCAGTAAAATCAGCCTATGACAATGCAGAGAAACGTCTGCAGAAAGACCAGAACGGCGCTGa                                       >  1:1961393/1‑70 (MQ=255)
                                 cGCCAAAGGCAGTAAAATCAGCCTATGACAATGCAGAGAAACGTCTGCAGAAAGACCAGAACGGCGCTGa                                       >  1:5247166/1‑70 (MQ=255)
                                   ccAAAGGCAGTAAAATCAGCCTATGACAATGCAGAGAAACGt                                                                 <  1:4297932‑M2/42‑2 (MQ=255)
                                     aaaGGCAGTAAAATCAGCCTATGACAATGCAGAcgtttctctgc                                                             >  1:1844267‑M1/1‑33 (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)
                                                                 gCAGAcgtttctctgcattgtcataggctgatttt                                          <  1:3789544‑M1/35‑31 (MQ=255)
                                                                 gCAGAGAAACGtctgcattgtcataggctgattttactgccttt                                 >  1:1844267‑M2/1‑11 (MQ=255)
                                                                  cAGAcgtttctctgcattgtcataggctgattttactgcctttggcgttgccgccagcgtttcagacgtac     >  1:3053952‑M1/1‑4 (MQ=255)
                                                                    gAcgtttctctgcattgtcataggctgattttactgcctttggcgttgccgccagcgtttcagacgtactg   <  1:997153‑M1/71‑70 (MQ=255)
                                                                     agaAACGTCTGCAGAAAGACCAGAACGGCGCTGATATACCCGATAAGGGATGCTTCCTGAACAACATTAAc  >  1:2613313/1‑71 (MQ=255)
                                                                     acgtttctctgcattgtcataggctgattttactgcctttggcgttgccgccagcgtttcagacgta      <  1:4934580‑M1/67‑67 (MQ=255)
                                                                     acgtttctctgcattgtcataggctgattttactgcctttgg                               <  1:4297932‑M1/42‑42 (MQ=255)
                                                                     |                                                                      
CGACCAACAGTACGTCTGAAACGCTGGCGGCAACGCCAAAGGCAGTAAAATCAGCCTATGACAATGCAGAGAAACGTCTGCAGAAAGACCAGAACGGCGCTGATATACCCGATAAGGGATGCTTCCTGAACAACATTAAC  >  AP009048/1432183‑1432322

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.