New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_002947	3435481 =	37 (0.820)	3 (0.070)	3/234	NT	8.5%	coding (350/981 nt)	gpJ	baseplate/tail fiber protein
	?	NC_002947	3435521 =	30 (0.700)					coding (390/981 nt)	gpJ	baseplate/tail fiber protein

CTGGTGAACCGCGCCGCTGATTACGGTGTGGAGAAGCTGGTTATCCGCCCGGCTGATCCTGATGCGGTGCCGCCGGTAGAGGCGGTAATGGAGGGCGACGAAGCGCTGCGCTACCGCACGCGGCTGTCGCTTGAGGCCTTGTCGGTAGCCGGCAGCAGCGGCGCCTATGAGTATCACGGGCTGAGCGCGTCGGCCGAGCTGACCAACGTGTCGGTCGATTCGCCCCGGTTTGTTGAGGCGGAGTTGACCCCCGAGGTACGGGCCCAGCTGCCGGCCGGGGCCATTGTGGTGGTCTGCGA  >  NC_002947/3435372‑3435670                                                                                                                                                      |                                                                                                                                                      tacctcgggggtcaactccgcctcaacaaaccggggcgaatcgaccgacacgttggtcagctcggccgacgcgctcagcccgtgatactcataggcgccgctgctgccgGCTACCGCACGCGGCTGTCGCTTGAGGCCTTGTCGGTAGcc                                                                                                                                                       >  2:657161‑M1/110‑150 (MQ=255)        aCCGCGCCGCTGATTACGGTGTGGAGAAGCTGGTTATCCGCCCGGCTGATCCTGATGCGGTGCCGCCGGTAGAGGCGGTAATGGAGGGCGACGAAGCGCTGCGCTACCGCACGCGGCTGTCGCTTGAGGCCTTGTCGGTAGCCGgcagca                                                                                                                                                >  1:525960/1‑150 (MQ=255)             gccgcTGATTACGGTGTGGAGAAGCTGGTTATCCGCCCGGCTGATCCTGATGCGGTGCCGCCGGTAGAGGCGGTAATGGAGGGCGACGAAGCGCTGCGCTACCGCACGCGGCTGTCGCTTGAGGCATTGTCGGTATCCGGCAGCAgcggc                                                                                                                                           <  2:412076/150‑1 (MQ=255)               cgcTGATTACGGTGTGGAGAAGCTGGTTATCCGCCCGGCTGATCCTGATGCGGTGTCGCCGGTAGAGGCGGTAATGGAGGGCGACGAAGCGCTGCGCTACCGCACGCGGCTGTCGCTTGAGGCCTTGTCGGTAGCCGGCAGCAGCGgcgc                                                                                                                                         >  1:330427/1‑150 (MQ=255)               cgcTGATTACGGTGTGGAGAAGCTGGTTATCCGCCCGGCTGATCCTGATGCGGTGCCGCCGGTAGAGGCGGTAATGGAGGGCGACGAAGCGCTGCGCTACCGCACGCGGCTGTCGCTTGAGGCCTTGTCGGTAGCCGGCAGCAGCGgcgc                                                                                                                                         <  1:736859/150‑1 (MQ=255)                   gATTACGGTGTGGAGAAGCTGGTTATCCGCCCGGCTGCTGTTGGTCCAGGTCCGCCTGTAGAGGCGGTAATGGAGTGCGACGAAGCGCTGCGCTACCGCACGCGGCTGTCGCTTGAGGCCTTGTCGGTAGCCGGCAGCAGCGGCGCCTAt                                                                                                                                     <  1:1173207/150‑1 (MQ=255)                    aTTACGGTGTGGAGAAGCTGGTTATCCGCCCGGCTGATCCTGATGCGGTGCCGCCGGTAGAGGCGGTAATGGAGGGCGACGAAGCGCTGCGCTACCGCACGCGGCTGTCGCTTGAGGCCTTGTCGGTAGCCGGCAGCAGCGGCGCCTATg                                                                                                                                    >  1:717981/1‑150 (MQ=255)                     ttACGGTGTGGAGAAGCTGGTTATCCGCCCGGCTGATCCTGATGCGGTGCCGCCGGTAGAGGCGGTAATGGAGGGCGACGAAGCGCTGCGCTACCGCACGCGGCTGTCGCTTGAGGCCTTGTCGGTAGCCGGCAGCAGCGGCGCCTATGa                                                                                                                                   <  2:525960/150‑1 (MQ=255)                         ggTGTGGAGAAGCTGGTTATCCGCCCGGCTGATCCTGATGCGGTGCCGCCGGTAGAGGCGGTAATGGAGGGCGACGAAGCGCTGCGCTACCGCACGCGGCTGTCGCTTGAGGCCTTGTCGGTAGCCGGCAGCAGCGGCGCCTATGAGTAt                                                                                                                               <  2:1060425/150‑1 (MQ=255)                              ggAGAAGCTGGTTATCCGCCCGGCTGATCCTGATGCGGTGCCGCCGGTAGAGGCGGTAATGGAGGGCGACGAAGCGCTGCGCTACCGCACGCGGCTGTCGCTTGAGGCCTTGTCGGTAGCCGGCAGCAGCGGCGCCTATGAGTATCACgg                                                                                                                          <  1:261887/150‑1 (MQ=255)                                 gaAGCTGGTTATCCGCCCGGCTGATCCTGATGCGGTGCCGCCGGTAGAGGCGGTAATGGAGGGCGACGAAGCGCTGCGCTACCGCACGCGGCTGTCGCTTGAGGCCTTGTCGGTAGCCGGCAGCAGCGGCGCCTa                                                                                                                                      <  1:551358/135‑1 (MQ=255)                                 gaAGCTGGTTATCCGCCCGGCTGATCCTGATGCGGTGCCGCCGGTAGAGGCGGTAATGGAGGGCGACGAAGCGCTGCGCTACCGCACGCGGCTGTCGCTTGAGGCCTTGTCGGTAGCCGGCAGCAGCGGCGCCTa                                                                                                                                      >  2:551358/1‑135 (MQ=255)                                     cTGGTTATCCGCCCGGCTGATCCTGATGCGGTGCCGCCGGTAGAGGCGGTAATGGAGGGCGACGAAGCGCTGCGCTACCGCACGCGGCTGTCGCTTGAGGCCTTGTCGGTAGCCGGCAGCAGCGGCGCCTATGAGTATCACGGGCTGAgc                                                                                                                   >  2:209802/1‑150 (MQ=255)                                         ttATCCGCCCGGCTGATCCTGATGCGGTGCCGCCGGTAGAGGCGGTAATGGAGGGCGACGAAGCGCTGCGCTACCGCACGCGGCTGTCGCTTGAGGCCTTGTCGGTAGCCGGCAGCAGCGGCGCCTATGAGTATCACGGGCTGAGCGCGt                                                                                                               <  1:543497/150‑1 (MQ=255)                                            tCCGCCCGGCTGATCCTGATGCGGTGCCGCCGGTAGAGGCGGTAATGGAGGGCGACGAAGCGCTGCGCTACCGCACGCGGCTGTCGCTTGAGGCCTTGTCGGTAGCCGGCAGCAGCAGCCCCTATGACTATCACGTGCTGAGCGCGTCgg                                                                                                            >  2:120124/1‑150 (MQ=255)                                             ccGCCCGGCTGATCCTTATGCGGTGCCGCCGGTAGAGGCGGTAATGGAGGGCGACGAAGCGCTGCGCTACCGCACGCGGCTGTCGCTTGAGGCCTTGTCGGTAGCCGGCAGCAGCGGCGCCTATGAGTATCACGGGCTGAGCGCGTCGGc                                                                                                           <  2:530623/150‑1 (MQ=255)                                                cccGGCTGAGCCTGATGCGGTGCCGCCGGTAGAGGCGGTAATGGAGGGCGACGAAGCGCTGCGCTCCCGCACGCGGCTGTCGCTTGAGGCCTTGTCGGTAGCCGGCAGCAGCGGCGCCTATGAGTATCACGGGCTGAGCGCGTCGGCCGa                                                                                                        <  2:463913/150‑1 (MQ=255)                                                 ccGGCTGATCCTGATGCGGTGCCGCCGGTAGAGGCGGTAATGGAGGGCGACGAAGCGCTGCGCTACCGCACGCGGCTGTCGCTTGAGGCCTTGTCGGTAGCCGGCAGCAGCGGCGCCTATGAGTATCACGGGCTGAGCGCGTCGGCCGAg                                                                                                       <  1:922095/150‑1 (MQ=255)                                                        aTCCTGATGCGGTGCCGCCGGTAGAGGCGGTAATGGAGGGCGACGAAGCGCTGCGCTACCGCACGCGGCTGTCGCTTGAGGCCTTGTCGGTAGCCGGCAGCAGCGGCGCCTATGAGTATCACGGGCTGAGCGCGTCGGCCGa                                                                                                        >  2:371433/1‑142 (MQ=255)                                                        aTCCTGAGGCGGTGCCGCCGGTAGAGGCGGTAATGGAGGGCGACGAAGCGCTGCGCTACCGCACGCGGCTGTCGCTTGAGGCCTTGTCGGTAGCCGGCAGCAGCGGCGCCTATGAGTATCACGGGCTGAGCGCGTCGGCCGa                                                                                                        <  1:371433/142‑1 (MQ=255)                                                          cagctcggccgacgcgctcagcccgtgatactcataggcgccgctgctgccgGCTACCGCACGCGGCTGTCGCTTGAGGCCTTGTCGGTAGCCGgca                                                                                                                                                   >  1:778397‑M1/53‑97 (MQ=255)                                                          cagctcggccgacgcgctcagcccgtgatactcataggcgccgctgctgccgGCTACCGCACGCGGCTGTCGCTTGAGGCCTTGTCGGTAGCCGgca                                                                                                                                                   <  2:778397‑M1/45‑1 (MQ=255)                                                                           ggTAGAGGCGGTAATGGAGGGCGACGAAGCGCTGCGCTACCGCACGCGGCTGTCGCTTGAGGCCTTGTCGGTAGCCGGCAGCAGCGGCGCCTATGAGTATCACGGGCTGAGCGCGTCGGCCGAGCTGACCAACGTGTCGGTCGATTCGcc                                                                             >  2:446026/1‑150 (MQ=255)                                                                            gTAGAGGCGGTAATGGAGGGCGACGAAGCGCTGCGCTACCGCACGCGGCTGTCGCTTGAGGCCTTGTCGGTAGCCGGCAGCAGCGGCGCCTATGAGTATCACGGGCTGAGCGCGTCGGCCGAGCTGa                                                                                                   >  2:1181813/1‑127 (MQ=255)                                                                            gTAGAGGCGGTAATGGAGGGCGACGAAGCGCTGCGCTACCGCACGCGGCTGTCGCTTGAGGCCTTGTCGGTAGCCGGCAGCAGCGGCGCCTATGAGTATCACGGGCTGAGCGCGTCGGCCGAGCTGa                                                                                                   <  1:1181813/127‑1 (MQ=255)                                                                                  gCGGTAATGGAGGGCGACGAAGCGCTGCGCTACCGCACGCGGCTGTCGCTTGAGGCCTTGTCGGTAGCCGGCAGCAGCGGCGCCTATGAGTATCACGGGCTGAGCGCGTCGGCCGAGCTGACCAACGTGTCGGTCGATTCGCCCCGGttt                                                                      >  2:1018977/1‑150 (MQ=255)                                                                                    ggTAATGGAGGGCGACGAAGCGCTGCGCTACCGCACGCGGCTGTCGCTTGAGGCCTTGTCGGTAGCCGGCAGCAGCGGCGCCTATGAGTATCACGGGCTGAGCGCGTCGGCCGAGCTGACCAACGTGTCGGTCGATTCGCCCCGGCttgt                                                                    <  2:777183/150‑1 (MQ=255)                                                                                     gTAATGGAGGGCGACGAAGCGCTGCGCTACCGCACGCGGCTGTCGCTTGAGGCCTTGTCGGTAGCCGGCAGCAGCGGCGCCTATGAGTATCACGGGCTGAGCGCGTCGGCCGAGCTGACCAACGTGTCGGTCGATTCGCCCCGGTttgtt                                                                   <  2:717981/150‑1 (MQ=255)                                                                                                         tgccggctaccgacaaggcctcaagcgacagccgcgtgcggtagcCGGCAGCAGCGGCGCCTATGAGTATCACGGGCTGAGCGCGTCGGCCGAGCTg                                                                                                    <  2:778397‑M2/52‑1 (MQ=255)                                                                                                         tgccggctaccgacaaggcctcaagcgacagccgcgtgcggtagcCGGCAGCAGCGGCGCCTATGAGTATCACGGGCTGAGCGCGTCGGCCGAGCTg                                                                                                    >  1:778397‑M2/46‑97 (MQ=255)                                                                                                             ggctaccgacaaggcctcaagcgacagccgcgtgcggtagcCGGCAGCAGCGGCGCCTATGAGTATCACGGGCTGAGCGCGTCGGCCGAGCTGACCAACGTGTCGGTCGATTCGCCCCGGTTTGTTGAGGCGGAGTTGACCCCCGAGGTa                                           >  2:657161‑M2/42‑150 (MQ=255)                                                                                                                            cTGTCGCTTGAGGCCTTGTCGGTAGCCGGCAGCAGCGGCGCCTATGAGTATCACGGGCTGAGCGCGTCGGCCGAGCTGACCAACGTGTCGGTCGATTCGCCCCGGTTTGTTGAGGCGGAGTTGACCCCCGAGGTACGGGCCCAGCTgccg                            >  2:604672/1‑150 (MQ=255)                                                                                                                            cTGTCGCTTGAGGCCTTGTCGGTAGCCGGCAGCAGCGGCGCCTATGAGTATCACGGGCTGAGCGCGTCGGCCGAGCTGACCAACGTGTCAGGCGATTCGCCCCGGGTGGGTGATTCGGAGTTGACCCCCGAGGTACGGGCCCAGCTgccg                            >  2:1040724/1‑150 (MQ=255)                                                                                                                                cGCTTGAGGCCTTGTCGGTAGCCGGCAGCAGCGGCGCCTATGAGTATCACGGGCTGAGCGCGTCGGCCGAGCTGACCAACGTGTCGGTCGATTCGCCCCGGTTTGTTGAGGCGGAGTTGACCCCCGAGGTACGGGCCCAGCTgccggccg                        <  2:448082/150‑1 (MQ=255)                                                                                                                                cGCTTGAGGCCTTGTCGGTAGCCGGCAGCAGCGGCGCCTATGAGTATCACGGGCTGAGCGCGTCGGCCGAGCTGACCAACGTGTCGGTCGATTCGCCCCGGTTTGTTGAGGCGGAGTTGACCCCCGAGGTACGGGCCCAGCTgccggccg                        <  1:446026/150‑1 (MQ=255)                                                                                                                                        gCCTCGTCGGTAGCCGGCAGCAGCGGCGCCTATGAGTATCACGGGCTGAGCGCGTCGGCCGAGCTGACCAACGTGTCGGTCGATTCGCCCCGGTTTGTTGAGGCGGAGTTGACCCCCGAGGTACGGGCCCAGCTGCCGGCCGGGGCCAtt                <  1:180392/150‑1 (MQ=255)                                                                                                                                           ttGTCGGTAGCCGGCAGCAGCGGCGCCTATGAGTATCACGGGCTGAGCGCGTCGGCCGAGCTGACCAACGTGTCGGTCGATTCGCCCCGGTTTGTTGAGGCGGAGTTGACCCCCGAGGTACGGGCCCAGCTGCCGGCCGGGGCCATtgtg             >  1:731595/1‑150 (MQ=255)                                                                                                                                            tGTCGGTAGCCGGCAGCAGCGGCGCCTATGAGTATCACGGGCTGAGCGCGTCGGCCGAGCTGACCAACGTGTCGGTCGATTCGCCCCGGTTTGTTGAGGCGGAGTTGACCCCCGAGGTACGGGCCCAGCTGCCGGCCGGGGCCATTgtgg            <  2:330427/150‑1 (MQ=255)                                                                                                                                              tCGGTAGCCGGCAGCAGCGGCGCCTATGAGTATCACGGGCTGAGCGCGTCGGCCGAGCTGACCAACGTGTCGGTCGAGTCGCCCCGGTTTGTTGAGGCGTAGTTGACCCCCGAGGTACGGGCCCAGCTGCCGGCCGGGGCCATTgtggtg          <  1:120124/150‑1 (MQ=255)                                                                                                                                               cggtagcCGGCAGCAGCGGCGCCTATGAGTATCACGGGCTGAGCGCGTCGGCCGAGCTGACCAACGTGTCGGTCGATTCGCCCCGGTTTGTTGAGGCGGAGTTGACCCCCGAGGTACGGGCCCAGCTGCCGGCCGGGGCCATTGtggtgg         <  1:906227‑M2/143‑1 (MQ=255)                                                                                                                                                      cGGCAGCAGCGGCGCCTATGAGTATCACGGGCTGAGCGCGTCGGCCGAGCTGACCAACGTGTCGGTCGATTCGCCCCGGTTTGTTGAGGCGGAGTTGACCCCCGAGGTACGGGCCCAGCTGCCGGCCGGGGCCATTGTGGTGGTCTGCGa  <  1:458401/150‑1 (MQ=255)                                                                                                                                                      |                                                                                                                                                      CTGGTGAACCGCGCCGCTGATTACGGTGTGGAGAAGCTGGTTATCCGCCCGGCTGATCCTGATGCGGTGCCGCCGGTAGAGGCGGTAATGGAGGGCGACGAAGCGCTGCGCTACCGCACGCGGCTGTCGCTTGAGGCCTTGTCGGTAGCCGGCAGCAGCGGCGCCTATGAGTATCACGGGCTGAGCGCGTCGGCCGAGCTGACCAACGTGTCGGTCGATTCGCCCCGGTTTGTTGAGGCGGAGTTGACCCCCGAGGTACGGGCCCAGCTGCCGGCCGGGGCCATTGTGGTGGTCTGCGA  >  NC_002947/3435372‑3435670

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