New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_002947	4389430 =	32 (0.530)	6 (0.100)	6/256	NT	16.0%	coding (1311/1956 nt)	PP_3865	tail protein
	?	NC_002947	4389463 =	32 (0.550)					coding (1278/1956 nt)	PP_3865	tail protein

GCCAAACACGCTTTTAATGAATTCCATGAACGGCGTGGTGTAGGCCTTGATCAGCTCCCACAGGGCGCTGAAGAACGGCCCGACCTTCTCCCAGTTGGCGACGATCAGACCGGCCGCCGCCGCG‑ATCCCCGCGGCGATCAGGCCAATCGGCGTGGCCG‑CGAAGGCCACACCGAGCACCCGCGTGGCCACTGTGGCCGCGAACACC‑GCCGCACGCAGCGCGGTGAAGGCCGCCCCGGCAATGACGATGCCGCGCACCAACTGCGGGTTGTCGCTGATCATCTGCGCCATATGGC  >  NC_002947/4389319‑4389611                                                                                                                                                  |                                                                                                                                                       gtgcgcggcatcgtcattgccggggcggccttcaccgcgctgcgtgcggcggtgttcgcggccacagtggccacgcgggtgctcggtgtggccttcgcggccacgccgattGGCCGCCGCCGCG‑ATCCCCGCGGCGATCAGGCCAATc                                                                                                                                                     >  8:210598‑M1/112‑148 (MQ=255)    aaaCACGCTTTTAATGAATTCCATGAACGGCGTGGTGTAGGCCTTGATCAGCTCCCACAGGGCGCTGAAGAACGGCCCGACCTTCTCCCAGTTGGCGACGATCAGACCGGCCGCCGCCGCGCATCCCCGCGGCGATCAGGCCAATCGGc                                                                                                                                                  <  1:17442/149‑1 (MQ=255)     aaCACGCTTTTAATGAATTCCATGAACGGCGTGGTGTAGGCCTTGATCAGCTCCCACAGGGCGCTGAAGAACGGCCCGACCTTCTCCCAGTTGGCGACGATCAGACCGGCCGCCGCCGCG‑ATCCCCGCGGCGATCAGGCCAATCGGCGt                                                                                                                                                <  7:36529/149‑1 (MQ=255)     aaCACGCTTTTAATGAATTCCATGAACGGCGTGGTGTAGGCCTTGATCAGCTCCCACAGGGCGCTGAAGAACGGCCCGACCTTCTCCCAGTTGGCGACGATCAGACCGGCCGCCGCCGCG‑ATCCCCGCGGCGATCAGGCCAATCGGCGt                                                                                                                                                <  7:302681/149‑1 (MQ=255)     aaCACGCTTTTAATGAATTCCATGAACGGCGTGGTGTAGGCCTTGATCAGCTCCCACAGGGCGCTGAAGAACGGCCCGACCTTCTCCCAGTTGGCGACGATCAGACCGGCCGCCGCCGCG‑ATCCCCGCGGCGATCAGGCCAATCGGCGt                                                                                                                                                <  7:251429/149‑1 (MQ=255)            ttttAATGAATTCCATGAACGGCGTGGTGTAGGCCTTGATCAGCTCCCACAGGGCGCTGAAGAACGGCCCGACCTTCTCCCAGTTGGCGACGATCAGACCGGCCGCCGCCGCG‑ATCCCCGCGGCGATCAGGCCAATCGGCGTGGCcg‑c                                                                                                                                         >  1:124147/1‑148 (MQ=255)               tAATGAATTCCATGAACGGCGTGGTGTAGGCCTTGATCAGCTCCCACAGGGCGCTGAAGAACGGCCCGACCTTCTCCCAGTTGGCGACGATCAGACCGGCCGCCGCCGCG‑ATCCCCGCGGCGATCAGGCCAATCGGCGTGGCCG‑CGaa                                                                                                                                      <  2:124147/148‑1 (MQ=255)                aaTGAATTCCATGAACGGCGTGGTGTAGGCCTTGATCAGCTCCCACAGGGCGCTGAAGAACGGCCCGACCTTCTCCCAGTTGGCGACGATCAGACCGGCCGCCGCCGCG‑CACCCCCCGGCGCTCAGGCCCATCGGCGTGGCCG‑CGAAgg                                                                                                                                    >  7:153889/1‑149 (MQ=255)                 aTGAATTCCATGAACGGCGTGGTGTAGGCCTTGATCAGCTCCCACAGGGCGCTGAAGAACGGCCCGACCTTCTCCCAGTTGGCGACGATCAGACCGGCCGCCGCCGCG‑ATCCCCGCGGGGATCAGGGCAATACGCGTGGCCG‑CGCAGGc                                                                                                                                   >  1:295197/1‑149 (MQ=255)                        ccATGAACGGCGTGGTGTAGGCCTTGATCAGCTCCCACAGGGCGCTGAAGAACGGCCCGACCTTCTCCCAGTTGGCGACGATCAGACCGGCCGCCGCCGCG‑ATCCCCGCGGCGATCAGGCCAATCGGCGTGGCCG‑CGAAGGCCACAcc                                                                                                                             <  2:274557/148‑1 (MQ=255)                          aTGAACGGCGTGGTGTAGGCCTTGATCAGCTCCCACAGGGCGCTGAAGAACGGCCCGACCTTCTCCCAGTTGGCGACGATCAGACCGGCCGCCGCCGCG‑ATCCCCGCGGCGATCAGGCCAATCGGCGTGGCCG‑CGAAGGCCACACCGAg                                                                                                                          >  4:309141/1‑149 (MQ=255)                                            gCCTTGATCAGCTCCCACAGGGCGCTGAAGAACGGCCCGACCTTCTCCCAGTTGGCGACGATCAGACCGGCCGCCGCCGCG‑ATCCCCGCGGCGATCAGGCCAATCGGCGTGGCCG‑CGAAGGCCACACCGAGCACCCGCGTGGCCACtgt                                                                                                        >  5:277704/1‑149 (MQ=255)                                                  aTCAGCTCCCACAGGGCGCTGAAGAACGGCCCGACCTTCTCCCAGTTGGCGACGATCAGACCGGCCGCCGCCGCG‑ATCCCCGCGGCGATCAGGCCAATCGGCGTGGCCG‑CGAAGGCCACACCGAGCACCCGCGTGGCCACTGTGGCcgc                                                                                                  <  1:176348/149‑1 (MQ=255)                                                    cAGCTCCCACAGGGCGCTGAAGAACGGCCCGACCTTCTCCCAGTTGGCGACGATCAGACCGGCCGCCGCCGCG‑ATCCCCGCGGCGATCAGGCCAATCGGCGTGGCCG‑CGAAGGCCACACCGAGCACCCGCGTGGCCACTGTGGCcgcg                                                                                                 >  7:28566/1‑148 (MQ=255)                                                     aGCTCCCACAGGGCGCTGAAGAACGGCCCGACCTTCTCCCAGTTGGCGACGATCAGACCGGCCGCCGCCGCG‑ATCCCCGCGGCGATCAGGCCAATCGGCGTGGCCG‑CGAAGGCCACACCGAGCACCCGCGTGGCCACTGTGGCCGCGaa                                                                                               <  5:297615/149‑1 (MQ=255)                                                       cTCCCACAGGGCGCTGAAGAACGGCCCGACCTTCTCCCAGTTGGCGACGATCAGACCGGCCGCCGCCGCG‑ATCCCCGCGGCGATCAGGCCAATCGGCGTGGCCG‑CGAAGGCCACACCGAGCACCCGCGTGGCCACTGTGGCCGCGAaca                                                                                             >  5:116478/1‑149 (MQ=255)                                                            acaGGGCGCTGAAGAACGGCCCGACCTTCTCCCAGTTGGCGACGATCAGACCGGCCGCCGCCGCG‑ATCCCCGCGGCGATCAGGCCAATCGGCGTGGCCG‑CGAAGGCCACACCGAGCACCCGCGTGGCCACTGTGGCCGCGAACAcc‑gcc                                                                                       <  3:137417/149‑1 (MQ=255)                                                                tgcgcTGAATAACGGCCCGACCTTCTCCCAGTTGGCGACGATCAGACCGGCAGCCGCCGCG‑ATCCCCGCGGCGATCAGGCCAATCGGCGTGGCCG‑CGAAGGCCACACCGAGCACCCGCGTGGCCACTGTGGCCGCGAACACC‑GCcgcac                                                                                   <  8:276896/148‑1 (MQ=255)                                                                       aagaaCGGCCCGACCTTCTCCCAGTTGGCGACGATCAGACCGGCCGCCGCCGCG‑ATCCCCGCGGCGATCAGGCCAATCGGCGTGGCCG‑CGAAGGCCACACCGAGCACCCGCGTGGCCACTGTGGCCGCGAACACC‑GCCGCACGCAgcgc                                                                            <  6:277704/149‑1 (MQ=255)                                                                                   aCCTTCTCCCAGTTGGCGACGATCAGACCGGCCGCCGCCGCG‑ATCCCCGCGGCGCTCAGGCCAATAGGCGTGGCCG‑CGGAGGCCCCACCGAGCCCCCGCGGGGGCCCAGGGGGCGCCCACACC‑CCCGCCCGCCGCGCGGGGAAGgcccc                                                                >  2:296795/1‑147 (MQ=255)                                                                                   aCCTTCTCCCAGTTGGCGACGATCAGACCGGCCGCCGCCGCG‑ATACCCGCGGCGATCAGGCCAATCGGCGTGGCCG‑CGAAGGCCACACCGAGCACCCGCGTGGCCACTGTGGCCGCGAACACA‑GCCGCACGCAGCGCGGTGAAGgccgc                                                                >  5:375222/1‑149 (MQ=255)                                                                                        ctcCCAGTTGGCGACGATCAGACCGGCCGCCGCCGCG‑ATCCCCGCGGCGATCAGGCCAATCGGCGTGGCCG‑CGAAGGCCACACCGAGCACCCGCGTGGCCACTGTGGCCGCGAACACCAGCCGCACGCAGCGCGGTGAAGGCCGccc                                                              <  1:214249/147‑1 (MQ=255)                                                                                           ccAGTTGGCGACGATCAGACCGGCCGCCGCCGCG‑ATCCCCGCGGCGATCAGGCCAATCGGCGTGGCCG‑CGAAGGCCACACCGAGCACCCGCGTGGCCACTGTGGCCGCGAACACC‑GCCGCACGCAGCGCGGTGAAGGCCGCCCCgg                                                           <  6:116478/146‑1 (MQ=255)                                                                                                   cgacgaTCAGACCGGCCGCCGCCGCG‑ATCCCCGCGGCGATCAGGCCAATCGGCGTGGCCG‑CGAAGGCCACACCGAGCACCCGCGTGGCCACTGTGGCCGCGAACACC‑GCCGCACGCAGCGCGGTGAAGGCCGCCCCGGCAATGACGATg                                                >  3:80726/1‑149 (MQ=255)                                                                                                    gacgaTCAGACCGGCCGCCGCCGCA‑AACCCCGCCGCGATCAAGCCAAAACCCCTGGACGAGGACGACCAACCCCCCCACCC‑CGAGGCCACATTGGCCGCGAACACC‑CACGCAAGCAGCGCGCTGAAGCCCCACCCGGCAAAGCCGAttc                                               >  6:147084/1‑147 (MQ=255)                                                                                                        aTCAGACCGGCCGCCGCCGCG‑ATCCCCGCGGCGATCAGGCCAATCGGCGTGGCCG‑CGAAGGCCACACCGAGCACCCGCGTGGCCACTGTGGCCGCGAACACC‑GCCGCACGCAGCGCGGTGAAGGCCGCCCCGGCAATGACGATGCcgcg                                           <  8:153889/149‑1 (MQ=255)                                                                                                             gattggcctgatcgccgcggggatcgcggcggcggccAATCGGCGTGGCCG‑CGAAGGCCACACCGAGCACCCGCGTGGCCACTGTGGCCGCGAACACC‑GCCGCACGCAGCGCGGTGAAGGCCGCCCCGGCAATGACGATGCCGCGCAc                                        >  8:210598‑M2/38‑148 (MQ=255)                                                                                                              ccGGCCGCCGCCGCG‑ATCCCCGCGGCGATCAGGCCAATCGGCGTGGCCG‑CGAAGGCCACACCGAGCACCCGCGTGGCCACTGTGGCCGCGAACACC‑GCCGCACGCAGCGCGGTGAAGGCCGCCCCg                                                            >  2:68983/1‑126 (MQ=255)                                                                                                              ccGGCCGCCGCCGCG‑ATCCCCGCGGCGATCAGGCCAATCGGCGTGGCCG‑CGAAGGCCACACCGAGCACCCGCGTGGCCACTGTGGCCGCGAACACC‑GCCGCACGCAGCGCGGTGAAGGCCGCCCCg                                                            <  1:68983/126‑1 (MQ=255)                                                                                                                gcctgatcgccgcggggatcgcggcggtcaggccAATCGGCGTGGCCG‑CGAAGGCCACACCGAGCACCCGCGTGGCCACTGTGGCCGCGAACACC‑GCCGCACGCAGCGCGGTGAAGGCCGCCCCGGCAATGACGATGCCGCGCACCAAc                                    >  3:299135‑M2/35‑149 (MQ=255)                                                                                                                  tgatcgccgcG‑ATCCCCGCGGCGATCAGGCCAATCGGCGTGGCCG‑CGAAGGCCACACCGAGCACCCGCGTGGCCACTGTGGCCGCGAACACC‑GCCGCACGCAGCGCGGTGAAGGCCGCCCCGGCAATGACGATGCCGCGCACCAACTGc                                 >  6:47594/5‑149 (MQ=255)                                                                                                                   tgatcgccgcggggatcgcggcggtcaggccAATCGGCGTGGCCG‑CGAAGGCCACACCGAGCACCCGCGTGGCCACTGTGGCCGCGAACACC‑GCCGCACGCAGCGCGGTGAAGGCCGCCCCGGCAATGACGATGCCGCGCACCAACTGc                                 >  8:82391‑M2/32‑149 (MQ=255)                                                                                                                   cctgatcgccgcggggatcgcggcggcggccAATCGGCGTGGCCG‑CGAAGGCCACACCGAGCACCCGCGTGGCCACTGTGGCCGCGAACACC‑GCCGCACGCAGCGCGGTGa                                                                       >  4:222867‑M2/32‑111 (MQ=255)                                                                                                                   cctgatcgccgcggggatcgcggcggcggccAATCGGCGTGGCCG‑CGAAGGCCACACCGAGCACCCGCGTGGCCACTGTGGCCGCGAACACC‑GCCGCACGCAGCGCGGTGa                                                                       <  3:222867‑M2/80‑1 (MQ=255)                                                                                                                    gccgccgcG‑ATCCCCGCGGCGATCAGGCCAATCGGCGTGGCCG‑CGAAGGCCACACCGAGCACCCGCGTGGCCACTGTGGCCGCGAACACC‑GCCGCACGCAGCGCGGTGAAGGCCGCCCCGGCAATGACGATGCCGCGCACCAACTGCgg                               >  2:224258/1‑149 (MQ=255)                                                                                                                           cg‑ATCCCCGCGGCGATCAGGCCAATCGGCGTGGCCG‑CGAAGGCCACACCGAGCACCCGCGTGGCCACTGTGGCCGCGAACACC‑GCCGCACGCAGCGCGGTGAAGGCCGCCCCGGCAATGACGATGCCGCGCACCAACTGCGGGTTGTCg                        <  8:28566/149‑1 (MQ=255)                                                                                                                                gggatcgcggcggcggccAATCGGCGTGGCCG‑CGAAGGCCACACCGAGCACCCGCGTGGCCACTGTGGCCGCGAACACC‑GCCGCACGCAGCGCGGTGAAGGCCGCCCCGGCAATGACGATGCCGCGCACCAACTGCGGGTTGTCGCTGa                    <  7:210598‑M2/131‑1 (MQ=255)                                                                                                                                    gcggcgATCAGGCCAATCGGCGTGGCCG‑CGAAGGCCACACCGAGCACCCGCGTGGCCACTGTGGCCGCGAACACC‑GCCGCACGCAGCGCGGTGAAGGCCGCCCCGGCaa                                                        <  2:19594/109‑1 (MQ=255)                                                                                                                                    gcggcgATCAGGCCAATCGGCGTGGCCG‑CGAAGGCCACACCGAGCACCCGCGTGGCCACTGTGGCCGCGAACACC‑GCCGCACGCAGCGCGGTGAAGGCCGCCCCGGCaa                                                        >  1:19594/1‑109 (MQ=255)                                                                                                                                          aTCAGGCCAATCGGCGTGGCCG‑CGAAGGCCACACCGAGCACCCGCGTGGCCACTGTGGCCGCGAACACC‑GCCGCACGCAGCGCGGTGAAGGCCGCCCCGGCAATGACGATGCCGCGCACCAACTGCGGGTTGTCGCTGATCATCTgcgc          >  8:144483/1‑149 (MQ=255)                                                                                                                                               gccAATCGGCGTGGCCG‑CGAAGGCCACACCGAGCACCCGCGTGGCCACTGTGGCCGCGAACACC‑GCCGCACGCAGCGCGGTGAAGGCCGCCCCGGCAATGACGATGCCGCGCACCAACTGCGGGTTGTCGCTGATCATCTGCGCCatat     >  2:155266‑M2/4‑149 (MQ=255)                                                                                                                                                 cAATCGGCGTGGCCG‑CGAAGGCCACACCGAGCACCCGCGTGGCCACTGTGGCCGCGAACACC‑GCCGCACGCAGCGCGGTGAAGGCCGCCCCGGCAATGACGATGCCGCGCACCAACTGCGGGTTGTCGCTGATCATCTGCGCCATATgg   <  4:80726‑M2/148‑1 (MQ=255)                                                                                                                                                  aaTCGGCGTGGCCG‑CGAAGGCCACACCGAGCACCCGCGTGGCCACTGTGGCCGCGAACACC‑GCCGCACGCAGCGCGGTGAAGGCCGCCCCGGCAATGACGATGCCGCGCACCAACTGCGGGTTGTCGCTGATCATCTGCGCCATATGGc  >  8:208147/1‑149 (MQ=255)                                                                                                                                                  aaTCGGCGTGGCCG‑CGAAGGCCACACCGAGCACCCGCGTCGCCACTGTGGCCGCGAACACC‑GCCGCACGCAGCGCGGTGAAGGCCGCCCCGGCAATGACGATGCCGCGCACCAACTGCGGGTTGTCGCTGATCATCTGCGCCATATGGc  >  2:140862/1‑149 (MQ=255)                                                                                                                                                  |                                                                                                                                                       GCCAAACACGCTTTTAATGAATTCCATGAACGGCGTGGTGTAGGCCTTGATCAGCTCCCACAGGGCGCTGAAGAACGGCCCGACCTTCTCCCAGTTGGCGACGATCAGACCGGCCGCCGCCGCG‑ATCCCCGCGGCGATCAGGCCAATCGGCGTGGCCG‑CGAAGGCCACACCGAGCACCCGCGTGGCCACTGTGGCCGCGAACACC‑GCCGCACGCAGCGCGGTGAAGGCCGCCCCGGCAATGACGATGCCGCGCACCAACTGCGGGTTGTCGCTGATCATCTGCGCCATATGGC  >  NC_002947/4389319‑4389611

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