New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_002947	= 3437057	41 (0.760)	4 (0.080)	4/240	NT	9.4%	coding (183/879 nt)	PP_3056	tail fiber protein
	?	NC_002947	= 3437081	39 (0.780)					coding (207/879 nt)	PP_3056	tail fiber protein

CCTGGCTGGAGTCTGGAATGCCACCAGTACAGGGCTTTCTGCCGAGATCACCCACATCGCGCTCGGTACGTCGGGTTACACCCCGACCAACGAGCAAACCAGCTTGCGCGCCCAGGTGGCCAAGTATCCCATTGCTGGCGGTGAGCGGCTGAGTGACAGCCTGATTCACCTCACCGCCATCGCTGATGGCCCGGCCGCGTTCTGGGTGCGTGAAATCGGCTTTCTGTTAGCCGATGGCACGCTGTTGGCGGTTTGGTCGCATGCGACTGACGCTCTGACCTACAAGCCG  >  NC_002947/3436913‑3437201                                                                                                                                                 |                                                                                                                                                 ccTGGCTGGAGTCTGGAATGCCACCAGTACAGGGCTTTCTGCCGAGATCACCCACATCGCGCTCGGTACGTCGGGTTACACCCCGACCAACGAGCAAACCAGCTTGCGCGCCCAGGTGGCCAAGTATCCCATTGCTGGCGGTGAGCGGc                                                                                                                                              <  4:75889/149‑1 (MQ=255)   tggctggAGTCTGGAATGCCACCAGTACAGGGCTTTCTGCCGAGATCACCCACATCGCGCTCGGTACGTCGGGTTACACCCCGACCAACGAGCAAACCAGCTTGCGCGCCCAGGTGGCCAAGTATCCCATTGCTGGCGGTGAGCGGCt                                                                                                                                             >  1:77685/1‑148 (MQ=255)     gctggAGTCTGGAATGCCACCAGTACAGGGCTTTCTGCCGAGATCACCCACATCGCGCTCGGTACGTCGGGTTACACCCCGACCAACGAGCAAACCAGCTTGCGCGCCCAGGTGGCCAAGTATCCCATTGCTGGCGGTGAGCGGCTGa                                                                                                                                           <  6:340562/148‑1 (MQ=255)      ctggAGTCTGGAATGCCACCAGTACAGGGCTTTCTGCCGAGATCACCCACATCGCGCTCGGTACGTCGGGTTACACCCCGACCAACGAGCAAACCAGCTTGCGCGCCCAGGTGGCCAAGTATCCCATTGCTGGCGGTGAGCGGCTGAGt                                                                                                                                         <  2:240324/149‑1 (MQ=255)          aGTCTGGAATGCCACCAGTACAGGGCTTTCTGCCGAGATCACCCACATCGCGCTCGGTACGTCGGGTTACACCCCGACCAACGAGCAAACCAGCTTGCGCGCCCAGGTGGCCAAGTATCCCATTGCTGGCGGTGAGCGGCTGAGTGACa                                                                                                                                     <  6:220108/149‑1 (MQ=255)          aGTCTGGAATGCCACCAGTACAGGGCTTTCTGCCGAGATCACCCACATCGCGCTCGGTACGTCGGGTTACACCCCGACCAACGAGCAAACCAGCTTGCGCGCCCAGGTGGCCAAGTATCCCATTGCTGGCGGTGAGCGGCTGAGTGACa                                                                                                                                     <  3:8801/149‑1 (MQ=255)              tGGAATGCCACCAGTACAGGGCTTTCTGCCGAGATCACCCACATCGCGCTCGGTACGTCGGGTTACACCCCGACCAACGAGCAAACCAGCTTGCGCGCCCAGGTGGCCAAGTATCCCATTGCTGGCGGTGAGCGGCTGAGTGACAGCCt                                                                                                                                 >  3:164585/1‑149 (MQ=255)              tGGAATGCCACCAGTACAGGGCTTTCTGCCGAGATCACCCACATCGCGCTCGGTACGTCGGGTTACACCCCGACCAACGAGCAAACCAGCTTGCGCGCCCAGGTGGCCAAGTATCCCATTGCTGGCGGTGAGCGGCTGAGTGACAGCCt                                                                                                                                 >  8:60956/1‑149 (MQ=255)                     ccaccaGTACAGGGCTTTCTGCCGAGATCACCCACATCGCGCTCGGTACGTCGGGTTACACCCCGACCAACGAGCAAACCAGCTTGCGCGCCCAGGTGGCCAAGTATCCCATTGCTGGCGGTGAGgtgaatcaggctgtcactcagccg                                                                                                                          >  3:115683‑M1/1‑125 (MQ=255)                     ccaccaGTACAGGGCTTTCTGCCGAGATCACCCACATCGCGCTCGGTACGTCGGGTTACACCCCGACCAACGAGCAAACCAGCTTGCGCGCCCAGGTGGCCAAGTATCCCATTGCTGGCGGTGAGCGGCTGAGTGAc                                                                                                                                      <  5:208592/137‑1 (MQ=255)                     ccaccaGTACAGGGCTTTCTGCCGAGATCACCCACATCGCGCTCGGTACGTCGGGTTACACCCCGACCAACGAGCAAACCAGCTTGCGCGCCCAGGTGGCCAAGTATCCCATTGCTGGCGGTGAGCGGCTGAGTGAc                                                                                                                                      >  6:208592/1‑137 (MQ=255)                          aGTACAGGGCTTTCTGCCGAGATCACCCACATCGCGCTCGGTACGTCGGGTTACACCCCGACCAACGAGCAAACCAGCTTGCGCGCCCAGGTGGCCAAGTATCCCATTGCTGGCGGTGAGCGGCTGAGTGACAGCCTGATTCACctcac                                                                                                                     >  4:323725‑M2/1‑144 (MQ=255)                              cAGGGCTTTCTGCCGAGATCACCCACATCGCGCTCGGTACGTCGGGTTACACCCCGACCAACGAGCAAACCAGCTTGCGCGCCCAGGTGGCCAAGTATCCCATTGCTGGCGGTGAGgtgaatcaggctgtcactcagccgctcaccgcc                                                                                                                 <  4:115683‑M1/149‑34 (MQ=255)                                          ccGAGATCACCCACATCGCGCTCGGTACGTCGGGTTACACCCCGACCAACGAGCAAACCAGCTTGCGCGCCCAGGTGGCCAAGTATCCCATTGCTGGCGGTGAGCGGCTGAGTGACAGCCTGATTCACCTCACCGCCATCGCTGATGGc                                                                                                     >  3:123475/1‑149 (MQ=255)                                             agaTCACCCACATCGCGCTCGGTACGTCGGGTTACA‑CCCGACCAACGAGCAAACCAGCTTGCGCGCCCAGGTGGCCAAGTATCCCATTGCTGGCGGTGAGCGGCTGAGTGACAGCCTGATTCACCTCACCGCCATCGCTGATGGCCCgg                                                                                                 <  2:194713/149‑1 (MQ=255)                                                    ccACATCGCGCTCGGTACGTCGGGTTACACCCCGACCAACGAGCAAACCAGCTTGCGCGCCCAGGTGGCCAAGTATCCCATTGCTGGCGGTGAGCGGCTGAGTGACAGCCTGATTCACCTCACCGCCATCGCTGATGGCCCGGCCGCGt                                                                                           >  6:48843/1‑149 (MQ=255)                                                         tCGCGCTCGGTACGTCGGGTTACACCCCGACCAACGAGCAAACCAGCTTGCGCGCCCAGGTGGCCAAGTATCCCATTGCTGGCGGTGAGCGGCTGAGTGACAGCCTGATTCACCTCACCGCCATCGCTGATGGCCCGGCCTCGTTGTgg                                                                                      <  4:164585/149‑1 (MQ=255)                                                         tCGCGCTCGGTACGTCGGGTTACACCCCGACCAACGAGCAAACCAGCTTGCGCGCCCAGGTGGCCAAGTATCCCATTGCTGGCGGTGAGCGGCTGAGTGACAGCCTGATTCACCTCACCGCCATCGCTGATGGCCCGGCCGTGTTCTgg                                                                                      <  2:77685/149‑1 (MQ=255)                                                          cgcgCTCGGTACGTCGGGTTACACCCCGACCAACGAGCAAACCAGCTTGCGCGCCCAGGTGGCCAAGTATCCCATTGCTGGCGGTGAGCGGCTGAGTGACAGCCTGATTCACCTCACCGCCATCGCTGATGGCCCGGCCGCGTTCTggg                                                                                     >  3:5230/1‑149 (MQ=255)                                                           gcgcTCGGTACGTCGGGTTACACCCCGACCAACGAGCAAACCAGCTTGCGCGCCCAGGTGGCCAAGTATCCCATTGCTGGCGGTGAGCGGCTGAGTGACAGCCTGATTCACCTCACCGCCATCGCTGATGGCCCGGCCGCGTTCTGGGt                                                                                    <  4:53786/149‑1 (MQ=255)                                                            cgcTCGGTACGTCGGGTTACACCCCGACCAACGAGCAAACCAGCTTGCGCGCCCAGGTGGCCAAGTATCCCATTGCTGGCGGTGAGgtgaatcaggctgtcactcagccgctcaccgccag                                                                                                               <  2:340829‑M1/121‑36 (MQ=255)                                                            cgcTCGGTACGTCGGGTTACACCCCGACCAACGAGCAAACCAGCTTGCGCGCCCAGGTGGCCAAGTATCCCATTGCTGGCGGTGAGgtgaatcaggctgtcactcagccgctcaccgccag                                                                                                               >  1:340829‑M1/1‑86 (MQ=255)                                                                 ggTACGTCGGGTTACACCCCGACCAACGAGCAAACCAGCTTGCGCGCCCAGGTGGCCAAGTATCCCATTGCTGGCGGTGAGCGGCTGAGTGACAGCCTGATTCACCTCACCGCCATCGCTGATGGCCCGGCCGCGTTCTGGGTGCGTg                                                                               >  6:227465/1‑148 (MQ=255)                                                                      gTCGGGTTACACCCCGACCAACGAGCAAACCAGCTTGCGCGCCCAGGTGGCCAAGTATCCCATTGCTGGCGGTGAGCGGCTGAGTGACAGCCTGATTCACCTCACCGCCATCGCTGATGGCCCGGCCGCGTTCTGGGt                                                                                    <  2:16320/138‑1 (MQ=255)                                                                      gTCGGGTTACACCCCGACCAACGAGCAAACCAGCTTGCGCGCCCAGGTGGCCAAGTATCCCATTGCTGGCGGTGAGCGGCTGAGTGACAGCCTGATTCACCTCACCGCCATCGCTGATGGCCCGGCCGCGTTCTGGGt                                                                                    >  1:16320/1‑138 (MQ=255)                                                                        cGGGTTACACCCCGACCAACGAGCAAACCAGCTTGCGCGCCCAGGTGGCCAAGTATCCCATTGCTGGCGGTGAGCGGCTGAGTGACAGCCTGATtcac                                                                                                                          <  2:334520/98‑1 (MQ=255)                                                                        cGGGTTACACCCCGACCAACGAGCAAACCAGCTTGCGCGCCCAGGTGGCCAAGTATCCCATTGCTGGCGGTGAGCGGCTGAGTGACAGCCTGATtcac                                                                                                                          >  1:334520/1‑98 (MQ=255)                                                                           gTTACACCCCGACCAACGAGCAAACCAGCTTGCGCGCCCAGGTGGCCAAGTATCCCATTGCTGGCGGTGAGCGGCTGAGTGACAGCCTGATTCACCTCACCGCCATCGCTGATGGCCCGGCCGCGTTCTGGGTGCGTGAAATCGGCttt                                                                    <  5:48843/149‑1 (MQ=255)                                                                                  cccGACCAACGAGCAAACCAGCTTGCGCGCCCAGGTGGCCAAGTATCCCATTGCTGGCGGTGAGCGGCTGAGTGACAGCCTGATTCACCTCACCGCCATCGCTGATGGCACGGCCGCGTTCTGGGTGCGTGAAATCGGCTTTCTGTTAg                                                             >  4:208621/1‑149 (MQ=255)                                                                                        cAACGAGCAAACCAGCTTGCGCGCCCAGGTGGCCAAGTATCCCATTGCTGGCGGTGAGCGGCTGAGTGACAGCCTGATTCACCTCACCGCCATCGCTGATGGCCCGGCCGCGTTCTGGGTGCGTGAAATCGGCTTTCTGTTAGCCGATg                                                       <  3:208621/149‑1 (MQ=255)                                                                                            gAGCAAACCAGCTTGCGCGCCCAGGTGGCCAAGTATCCCATTGCTGGCGGTGAGCGGCTGAGTGACAGCCTGATTCACCTCACCGCCATCGCTGATGGCCCGGCCGCGTTCTGGGTGCGTGAAATCGGCTTTCTGTTAGCCGATGGCAc                                                   >  5:62174/1‑149 (MQ=255)                                                                                            gAGCAAACCAGCTTGCGCGCCCAGGTGGCCAAGTATCCCATTGCTGGCGGTGAGCGGCTGAGTGACAGCCTGATTCACCTCACCGCCATCGCTGATGGCCCGGCCGCGTTCTGGGTGCGTGAAATCGGCTTTCTGTTAGCCGATGGCAc                                                   >  3:62441/1‑149 (MQ=255)                                                                                              gCAAACCAGCTTGCGCGCCCAGGTGGCCAAGTATCCCATTGCTGGCGGTGAGCGGCTGAGTGACAGCCTGATTCACCTCACCGCCATCGCTGATGGCCCGGCCGCGTTCTGGGTGCGTGAAATCGGCTTTCTGTTAGCCGATGGCACGc                                                 >  5:63585/1‑149 (MQ=255)                                                                                                aaaCCAGCTTGCGCGCCCAGGTGGCCAAGTATCCCATTGCTGGCGGTGAGCGGCTGAGTGACAGCCTGATTCACCTCACCGCCATCGCTGATGGCCCGGCCGCGTTCTGGGTGCGTGAAATCGGCTTTCTGTTAGCCGATGGCACGCTg                                               >  7:250548/1‑149 (MQ=255)                                                                                                        ttGCGCGCCCAGGTGGCCAAGTATCCCATTGCTGGCGGTGAGCGGCTGAGTGACAGCCTGATTCACCTCACCGCCATCGCTGATGGCCCGGCCGCGTTCTGGGTGCGTGAAATCGGCTTTCTGTTAGCCGATGGCACGCTGTTGGCGGt                                       >  6:278720/1‑149 (MQ=255)                                                                                                        ttGCGCGCCCAGGTGGCCAAGTATCCCATTGCTGGCGGTGAGCGGCTGAGTGACAGCCTGATTCACCTCACCGCCATCGCTGATGGCCCGGCCGCGTTCTGGGTGCGTGAAATCGGCTTTCTGTTAGCCGATGGCACGCTGTTGGCGGt                                       >  1:314978/1‑149 (MQ=255)                                                                                                            gcgcCCAGGTGGCCAAGTATCCCATTGCTGGCGGTGAGCGGCTGAGTGACAGCCTGATTCACCTCACCGCCATCGCTGATGGCCCGGCCGCGTTCTGGGTGCGTGAAATCGGCTTTCTGTTAGCCGATGGCACGCTGTTGGCGGTTTgg                                   <  2:314978/149‑1 (MQ=255)                                                                                                             cgcCCAGGTGGCCAAGTATCCCATTGCTGGCGGTGAGCGGCTGAGTGACAGCCTGATTCACCTCACCGCCATCGCTGATGGCCCGGCCGCGTTCTGGGTGCGTGAAATCGGCTTTCTGTTAGCCGATGGCACGCTGTTGGCGGTTTGGt                                  <  5:227465/149‑1 (MQ=255)                                                                                                                ccAGGTGGCCAAGTATCCCATTGCTGGCGGTGAGCGGCTGAGTGACAGCCTGATTCACCTCACCGCCATCGCTGATGGCCCGGCCGCGTTCTGGGTGCGTGAAATCGGCTTTCTGTTAGCCGATGGCACGCTGTTGGCGGTTTGGTCg                                >  2:155024/1‑148 (MQ=255)                                                                                                                   ggTGGCCAAGTATCCCATTGCTGGCGGTGAGCGGCTGAGTGACAGCCTGATTCACCTCACCGCCATCGCTGATGGCCCGGCCGCGTTCTGGGTGCGTGAAATCGGCTTTCTGTTAGCCGATGGCACGCTGTTGGCGGTTTGGTCGCATg                            <  8:250548/149‑1 (MQ=255)                                                                                                                               tCCCATTGCTGGCGGTGAGCGGCTGAGTGACAGCCTGATTCACCTCACCGCCATCGCTGATGGCCCGGCAGCGTTCTGGGTGCGTGAAATCGGCTTTCTGTTAGCCGATGGCACGCTGTTGGCGGTTTGGTCGCATGCGACTGACGctc                >  4:196783/1‑149 (MQ=255)                                                                                                                                   aTTGCTGGCGGTGAGCGGCTGAGTGACAGCCTGATTCACCTCACCGCCATCGCTGATGGCCCGGCCGCGTTCTGGGTGCGTGAAATCGGCTTTCTGTTAGCCGATGGCACGCTGTTGGCGGTTTGGTCGCATGCGACTGACGCTCTGAc            <  5:278720/149‑1 (MQ=255)                                                                                                                                      gCTGGCGGTGAGCGGCTGAGTGAGAGCCTGATTCACCTCACCGCCATCGCTGATGGCCCGGCCGCGTTCTGGGAGCGTGAAATCGGCTTTCTGTTAGCCGATGGCACGCTGTTGGCGGTTTGGTCGCATGCGACTGACGCTCTGACCt          <  7:60956/148‑1 (MQ=255)                                                                                                                                       cTGGCGGTGAGCGGCTGAGTGACAGCCTGATTCACctcaccgccagcaatgggatacttggccacctgggcgcgcaagctggtttgctcgttggtcggggtgtaacccgacgtaccgagcg                                    <  2:340829‑M2/121‑87 (MQ=255)                                                                                                                                       cTGGCGGTGAGCGGCTGAGTGACAGCCTGATTCACctcaccgccagcaatgggatacttggccacctgggcgcgcaagctggtttgctcgttggtcggggtgtaacccgacgtaccgagcg                                    >  1:340829‑M2/1‑35 (MQ=255)                                                                                                                                       cTGGCGGTGAGCGGCTGAGTGACAGCCTGATTCACCTCACCGCCATCGCTGATGGCCCGGCCGCGTTCTGGGTGCGTGAAATCGGCTTTCTGTTAGCCGATGGCACGCTGTTGGCGGTTTGGTCGCATGCGACTGACGCTCTGAcc           >  5:277347/1‑146 (MQ=255)                                                                                                                                       cTGGCGGTGAGCGGCTGAGTGACAGCCTGATTCACCTCACCGCCATCGCTGATGGCCCGGCCGCGTTCTGGGTGCGTGAAATCGGCTTTCTGTTAGCCGATGGCACGCTGTTGGCGGTTTGGTCGCATGCGACTGACGCTCTGAcc           <  6:277347/146‑1 (MQ=255)                                                                                                                                       cTGGCGGTGAGCGGCTGAGTGACAGCCTGATTCACCTCACCGCCATCGCTGATGGCCCGGCCGCGTTCTGGGTGCGTGAAATCGGCTTTCTGTTAGCCGATGGCACGCTGTTGGCGGTTTGGTCGCATGCGACTGACGCTCTGACCTa         <  4:123475/148‑1 (MQ=255)                                                                                                                                        tGGCGGTGAGCGGCTGAGTGACAGCCTGATTCACCTCACCGCCATCGCTGATGGCCCGGCCGCGTTCTGGGTGCGTGAAATCGGCTTTCTGTTAGCCGATGGCACGCTGTTGGCGGTTTGGTCGCATGCGACTGACGCTCTGACCTAc        >  3:258132/1‑148 (MQ=255)                                                                                                                                         ggCGGTGAGCGGCTGAGTGACAGCCTGATTCACctcaccgccagcaatgggatacttggccacctgggcgcgcaagctggtttgctcgttggtcggggtgtaacccgacgtaccgagcgcgatgtgggtgatctcggcagaaagccctg      <  4:115683‑M2/149‑117 (MQ=255)                                                                                                                                          gCGGTGAGCGGCTGAGTGACAGCCTGATTCACCTCACCGCCATCGCTGATGGCCCGGCCGCGTTCTGGGTGCGTGAAATCGGCTTTCTGTTAGCCGATGGCACGCTGTTGGCGGTTTGGTCGCATGCGACTGACGCTCTGACCTACAAg     >  3:201056/1‑149 (MQ=255)                                                                                                                                             gtgaggtgaatcaggctgtcactcagccgctcaccgccagcaatgggatacttggccacctgggcgcgcaagctggtttgctcgttggtcggggtgtaacccgacgtaccgagcgcgatgtgggtgatctcggcagaaagccctgtact  >  4:323725‑M1/1‑5 (MQ=255)                                                                                                                                                 |                                                                                                                                                 CCTGGCTGGAGTCTGGAATGCCACCAGTACAGGGCTTTCTGCCGAGATCACCCACATCGCGCTCGGTACGTCGGGTTACACCCCGACCAACGAGCAAACCAGCTTGCGCGCCCAGGTGGCCAAGTATCCCATTGCTGGCGGTGAGCGGCTGAGTGACAGCCTGATTCACCTCACCGCCATCGCTGATGGCCCGGCCGCGTTCTGGGTGCGTGAAATCGGCTTTCTGTTAGCCGATGGCACGCTGTTGGCGGTTTGGTCGCATGCGACTGACGCTCTGACCTACAAGCCG  >  NC_002947/3436913‑3437201

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