New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913	= 1432648	NA (NA)	41 (0.790)	16/280	1.3	12.8%	coding (238/576 nt)	tfaR	Rac prophage; putative tail fiber assembly protein
	?	NC_000913	1634576 =	294 (5.380)					coding (310/576 nt)	tfaQ	Qin prophage; putative tail fibre assembly protein
Rejected: Frequency below cutoff threshold.

GCCTCATACAGGTCTGCCAGCAAACAGTACCGATATTGCACCGCCAGATATTCCGGCTGGCTTCGTGGCTGTTTTCAACAGTGATGAGTCATCGTGGCATCTCGTTGAAGATCATCGGGGTAAAACGGTTTATGACGTGGCTTCCGGCAACGCGTTATTTATTTCTGAACTCGGTCCGTTACCGGAAAATGTTACCTGGTTATCGCCGGAAGGGGAGTTTCAGAAGTGGAACGGCACAGCCTGGGTGAAGGATACGGAAGCAGAAAAACTGTTCCGGATCCGGGAGGCGG  >  NC_000913/1432509‑1432798                                                                                                                                            |                                                                                                                                                       gCCTCATACAGGTCTGCCAGCAAACAGTACCGATATTGCACCGCCAGATATTCCGGCTGGCTTCGTGGCTGTTTTCAACAGTGATGAGTCATCGTGGCATCTCGTTGAAGATCATCGGGGTAAAACGGTTTATGACGTGGCTTCCGGCaa                                                                                                                                              >  1:198191/1‑150 (MQ=31) gCCTCATACAGGTCTGCCAGCAAACAGTACCGATATTGCACCGCCAGATATTCCGGCTGGCTTCGTGGCTGTTTTCAACAGTGATGAGTCATCGTGGCATCTCGTTGAAGATCATCGGGGTAAAACGGTTTATGACGTGGCTTCCGGCaa                                                                                                                                              >  2:501507/1‑150 (MQ=31)  ccTCATACAGGTCTGCCAGCAAACAGTACCGATATTGCACCGCCAGATATTCCGGCTGGCTTCGTGGCTGTTTTCAACAGTGATGAGTCATCGTGGCATCTCGTTGAAGATCATCGGGGTAAAACGGTTTATGACGTGGCTTCCGGCAAcg                                                                                                                                            <  1:254098/151‑1 (MQ=31)    tCATACAGGTCTGCCAGCAAACAGTACCGATATTGCACCGCCAGATATTCCGGCTGGCTTCGTGGCTGTTTTCAACAGTGATGAGTCATCGTGGCATCTCGTTGAAGATCATCGGGGTAAAACGGTTTATGACGTGGCTTCCGGCAAcgcg                                                                                                                                          >  2:850208/1‑151 (MQ=31)    tCATACAGGTCTGCCAGCAAACAGTACCGATATTGCACCGCCAGATATTCCGGCTGGCTTCGTGGCTGTTTTCAACAGTGATGAGTCATCGTGGCATCTCGTTGAAGATCATCGGGGTAAAACGGTTTATGACGTGGCTTCCGGCAAcgcg                                                                                                                                          >  2:447067/1‑151 (MQ=31)      aTACAGGTCTGCCAGCAAACAGTACCGATATTGCACCGCCAGATATTCCGGCTGGCTTCGTGGCTGTTTTCAACAGTGATGAGTCATCGTGGCATCTCGTTGAAGATCATCGGGGTAAAACGGTTTATGACGTGGCTTCCGGCAACGCGtt                                                                                                                                        >  2:458314/1‑151 (MQ=31)            gTCTGCCAGCAAACAGTACCGATATTGCACCGCCAGATATTCCGGCTGGCTTCGTGGCTGTTTTCAACAGTGATGAGTCATCGTGGCATCTCGTTGAAGATCATCGGGGTAAAACGGTTTATGACGTGGCTTCCGGCAACGCGttattt                                                                                                                                    <  1:225549/149‑1 (MQ=31)              cTGCCAGCAAACAGTACCGATATTGCACCGCCAGATATTCCGGCTGGCTTCGTGGCTGTTTTCAACAGTGATGAGTCATCGTGGCATCTCGTTGAAGATCATCGGGGTAAAACGGTTTATGACGTGGCTTCCGGCAACGCGTTatttattt                                                                                                                                >  2:563271/1‑151 (MQ=31)               tGCCAGCAAACAGTACCGATATTGCACCGCCAGATATTCCGGCTGGCTTCGTGGCTGTTTTCAACAGTGATGAGTCATCGTGGCATCTCGTTGAAGATCATCGGGGTAAAACGGTTTATGACGTGGCTTCCGGCAACGCGTTATTTATTTc                                                                                                                               <  2:460257/151‑1 (MQ=31)                     cAAACAGTACCGATATTGCACCGCCAGATATTCCGGCTGGCTTCGTGGCTGTTTTCAACAGTGATGAGTCATCGTGGCATCTCGTTGAAGATCATCGGGGTAAAACGGTTTATGACGTGGCTTCCGGCAACGCGTTATTTATTTCTGAAc                                                                                                                          <  1:800522/150‑1 (MQ=31)                           gTACCGGTCTTGCACCGCCAGATATTCCGGCTGGTTTCGTGGCTGTTTTAAAAAGTGATGCTGCATCATGCCATCTCGCTGACGATAATCCTGGCAGAACCGGTTATGACGTGGCTTGCGGCAACGCGTCATTTCTTTCTGAACTCGGTc                                                                                                                    >  1:655654/1‑150 (MQ=2)                           gTACCGATATTGCACCGCCAGATATTCCGGCTGGCTTCGTGGCTGTTTTCAACAGTGATGAGTCATCGTGGCATCTCGTTGAAGATCATCGGGGTAAAACGGTTTATGACGTGGCTTCCGGCAACGCGTTATTTATTTCTGAACTCGGTc                                                                                                                    >  1:517637/1‑150 (MQ=31)                                 atatTGCACCGCCAGATATTCCGGCTGGCTTCGTGGCTGTTTTCAACAGTGATGAGTCATCGTGGCATCTCGTTGAAGATCATCGGGGTAAAACGGTTTATGACGTGGCTTCCGGCAACGCGTTATTTATTTCTGAACTCGGTCCGTTAc                                                                                                              >  2:518426/1‑150 (MQ=31)                                 atatTGCACCGCCAGATATTCCGGCTGGCTTCGTGGCTGTTTTCAACAGTGATGAGTCATCGTGGCATCTCGTTGAAGATCATCGGGGTAAAACGGTTTATGACGTGGCTTCCGGCAACGCGTTATTTATTTCTGAACTCGGTCCGTTAc                                                                                                              >  1:671117/1‑150 (MQ=31)                                 atatTGCACCGCCAGATATTCCGGCTGGCTTCGTGGCTGTTTTCAACAGTGATGAGTCATCGTGGCATCTCGTTGAAGATCATCGGGGTAAAACGGTTTATGACGTGGCTTCCGGCAACGCGTTATTTATTTCTGAACTCGGTCCGTTAc                                                                                                              >  1:419542/1‑150 (MQ=31)                                           gCCAGATATTCCGGCTGGCTTCGTGGCTGTTTTCAACAGTGATGAGTCATCGTGGCATCTCGTTGAAGATCATCGGGGTAAAACGGTTTATGACGTGGCTTCCGGCAACGCGTTATTTATTTCTGAACTCGGTCCGTTACCGGAAAATGtt                                                                                                   >  2:67420/1‑151 (MQ=31)                                             cAGATATTCCGGCTGGCTTCGTGGCTGTTTTCAACAGTGATGAGTCATCGTGGCATCTCGTTGAAGATCATCGGGGTAAAACGGTTTATGACGTGGCTTCCGGCAACGCGTTATTTATTTCTGAACTCGGTCCGTTACCGGAAAATGTTa                                                                                                  >  2:713889/1‑150 (MQ=31)                                                atatTCCGGCTGGCTTCGTGGCTGTTTTCAACAGTGATGAGTCATCGTGGCATCTCGTTGAAGATCATCGGGGTAAAACGGTTTATGACGTGGCTTCCGGCAACGCGTTATTTATTTCTGAACTCGGTCCGTTACCGGAAAATGTTACCt                                                                                               <  1:710657/150‑1 (MQ=31)                                                       ggctggctTCGTGGCTGTTTTCAACAGTGATGAGTCATCGTGGCATCTCGTTGAAGATCATCGGGGTAAAACGGTTTATGACGTGGCTTCCGGCAACGCGTTATTTATTTCTGAACTCGGTCCGTTACCGGAAAATGTTACCTGGTTATCg                                                                                       >  1:148063/1‑151 (MQ=31)                                                             ctTCGTGGCTGTTTTCAACAGTGATGAGTCATCGTGGCATCTCGTTGAAGATCATCGGGGTAAAACGGTTTATGACGTGGCTTCCGGCAACGCGTTATTTATTTCTGAACTCGGTCCGTTACCGGAAAATGTTACCTGGTTATCGCCGGaa                                                                                 >  2:264637/1‑151 (MQ=31)                                                             ctTCGTGGCTGTTTTCAACAGTGATGAGTCATCGTGGCATCTCGTTGAAGATCATCGGGGTAAAACGGTTTATGACGTGGCTTCCGGCAACGCGTTATTTATTTCTGAACTCGGTCCGTTACCGGAAAATGTTACCTGGTTATCGCCGGaa                                                                                 >  2:867750/1‑151 (MQ=31)                                                                  tGGCTGTTTTCAACAGTGATGAGTCATCGTGGCATCTCGTTGAAGATCATCGGGGTAAAACGGTTTATGACGTGGCTTCCGGCAACGCGTTATTTATTTCTGAACTCGGTCCGTTACCGGAAAATGTTACCTGGTTATCGCCGGAAgggg                                                                             <  2:497461/150‑1 (MQ=31)                                                                    gCTGTTTTCAACAGTGATGAGTCATCGTGGCATCTCGTTGAAGATCATCGGGGTAAAACGGTTTATGACGTGGCTTCCGGCAACGCGTTATTTATTTCTGAACTCGGTCCGTTACCGGAAAATGTTACCTGGTTATCGCCGGAAGGGGAg                                                                           <  2:91691/150‑1 (MQ=31)                                                                        ttttCAACAGTGATGAGTCATCGTGGCATCTCGTTGAAGATCATCGGGGTAAAACGGTTTATGACGTGGCTTCCGGCAACGCGTTATTTATTTCTGAACTCGGTCCGTTACCGGAAAATGTTACCTGGTTATCGCCGGAAGGGGAGTTTc                                                                       <  2:486818/150‑1 (MQ=31)                                                                         tttCAACAGTGATGAGTCATCGTGGCATCTCGTTGAAGATCATCGGGGTAAAACGGTTTATGACGTGGCTTCCGGCAACGCGTTATTTATTTCTGAACTCGGTCCGTTACCGGAAAATGTTACCTGGTTATCGCCGGAAGGGGAGTTTCAg                                                                     <  2:697893/151‑1 (MQ=31)                                                                            cAACAGTGATGAGTCATCGTGGCATCTCGTTGAAGATCATCGGGGTAAAACGGTTTATGACGTGGCTTCCGGCAACGCGTTATTTATTTCTGAACTCGGTCCGTTACCGGAAAATGTTACCTGGTTATCGCCGGAAGGGGAGTTTCAGAAg                                                                  <  1:213409/151‑1 (MQ=31)                                                                                 gTGATGAGTCATCGTGGCATCTCGTTGAAGATCATCGGGGTAAAACGGTTTATGACGTGGCTTCCGGCAACGCGTTATTTATTTCTGAACTCGGTCCGTTACCGGAAAATGTTACCTGGTTATCGCCGGAAGGGGAGTTTCAGAAGTGGaa                                                             >  2:409454/1‑151 (MQ=31)                                                                                   gctctttCATCGTGGCATCTCGTTGAAGCTCATCGGGGCAAAACGGGTTATGACGTGGCTTACCGCAACGCGTTATTTATTTCTGAACTCGGTCCGTTACCGGAAAATGTTACCTGGTTATCGCCGGCAGGGGAGTTTCAGAAGGGGAACg                                                           >  1:312863/7‑151 (MQ=11)                                                                                    atgaGTCATCGTGGCATCTCGTTGAAGATCATCGGGGTAAAACGGTTTATGACGTGGCTTCCGGCAACGCGTTATTTATTTCTGAACTCGGTCCGTTACCGGAAAATGTTACCTGGTTATCGCCGGAAGGGGAGTTTCAGAAGTGGAACgg                                                          <  2:14704/151‑1 (MQ=31)                                                                                      gaGTCATCGTGGCATCTCGTTGAAGATCATCGGGGTAAAACGGTTTATGACGTGGCTTCCGGCAACGCGTTATTTATTTCTGAACTCGGTCCGTTACCGGAAAATGTTACCTGGTTATCGCCGGAAGGGGAGTTTCAGAAGTGGAACGGca                                                        >  2:541428/1‑151 (MQ=31)                                                                                        gTCATCGTGGCATCTCGTTGAAGATCATCGGGGTAAAACGGTTTATGACGTGGCTTCCGGCAACGCGTTATTTATTTCTGAACTCGGTCCGTTACCGGAAAATGTTACCTGGTTATCGCCGGAAGGGGAGTTTCAGAAGTGGAACGGcaca                                                      >  1:120731/1‑151 (MQ=31)                                                                                         tCATCGTGGCATCTCGTTGAAGATCATCGGGGTAAAACGGTTTATGACGTGGCTTCCGGCAACGCGTTATTTATTTCTGAACTCGGTCCGTTACCGGAAAATGTTACCTGGTTATCGCCGGAAGGGGAGTTTCAGAAGTGGAACGGcac                                                       <  2:611140/149‑1 (MQ=31)                                                                                         tCATCGTGGCATCTCGTTGAAGATCATCGGGGTAAAACGGTTTATGACGTGGCTTCCGGCAACGCGTTATTTATTTCTGAACTCGGTCCGTTACCGGAAAATGTTACCTGGTTATCGCCGGAAGGGGAGTTTCAGAAGTGGAACGGCACAg                                                     >  2:634491/1‑151 (MQ=31)                                                                                           aTCGTGGCATCTCGTTGAAGATCATCGGGGTAAAACGGTTTATGACGTGGCTTCCGGCAACGCGTTATTTATTTCTGAACTCGGTCCGTTACCGGAAAATGTTACCTGGTTATCGCCGGAAGGGGAGTTTCAGAAGTGGAACGGCACAGcc                                                   <  2:356660/151‑1 (MQ=31)                                                                                           aTCGTGGCATCTCGTTGAAGATCATCGGGGTAAAACGGTTTATGACGTGGCTTCCGGCAACGCGTTATTTATTTCTGAACTCGGTCCGTTACCGGAAAATGTTACCTGGTTATCGCCGGAAGGGGAGTTTCAGAAGTGGAACGGCACAGcc                                                   <  2:520855/151‑1 (MQ=31)                                                                                            tCGTGGCATCTCGTTGAAGATCATCGGGGTAAAACGGTTTATGACGTGGCTTCCGGCAACGCGTTATTTATTTCTGAACTCGGTCCGTTACCGGAAAATGTTACCTGGTTATCGCCGGAAGGGGAGTTTCAGAAGTGGAACGGCACAGc                                                    >  1:432644/1‑149 (MQ=31)                                                                                            tCGTGGCATCTCGTTGAAGATCATCGGGGTAAAACGGTTTATGACGTGGCTTCCGGCAACGCGCTATTTATTTCTGAACTCGGTCCGTTACCGGAAAATGTTACCTGGTTATCGCCGGAAGGGGAGTTTCAGAAGTGGAACGGCACAGc                                                    >  2:297183/1‑149 (MQ=14)                                                                                               tGGCATCTCGTTGAAGATCATCGGGGTAAAACGGTTTATGACGTGGCTTCCGGCAACGCGTTATTTATTTCTGAACTCGGTCCGTTACCGGAAAATGTTACCTGGTTATCGCCGGAAGGGGAGTTTCAGAAGTGGAACGGCACAGCCTggg                                               <  2:308638/151‑1 (MQ=31)                                                                                               tGGCATCTCGTTGAAGATCATCGGGGTAAAACGGTTTATGACGTGGCTTCCGGCAACGCGTTATTTATTTCTGAACTCGGTCCGTTACCGGAAAATGTTACCTGGTTATCGCCGGAAGGGGAGTTTCAGAAGTGGAACGGCACAGCCTggg                                               <  2:120731/151‑1 (MQ=31)                                                                                                  cATCTCGTTGAAGATCATCGGGGTAAAACGGTTTATGACGTGGCTTCCGGCAACGCGTTATTTATTTCTGAACTCGGTCCGTTACCGGAAAATGTTACCTGGTTATCGCCGGAAGGGGAGTTTCAGAAGTGGAACGGCACAGCCTGGGTg                                             <  2:452379/150‑1 (MQ=31)                                                                                                  cATCTCGTTGAAGATCATCGGGGTAAAACGGTTTATGACGTGGCTTCCGGCAACGCGTTATTTATTTCTGAACTCGGTCCGTTACCGGAAAATGTTACCTGGTTATCGCCGGAAGGGGAGTTTCAGAAGTGGAACGGCACAGCCTGGGTg                                             <  1:262872/150‑1 (MQ=31)                                                                                                  cATCTCGTTGAAGATCATCGGGGTAAAACGGTTTATGACGTGGCTTCCGGCAACGCGTTATTTATTTCTGAACTCGGTCCGTTACCAGAAAATGTTACCTGGTTATCGCCGGAAGGGGAGTTTCAGAAGTGGAACGGCACAGCCTGGGTg                                             <  2:311782/150‑1 (MQ=14)                                                                                                     ctcGTTGAAGATCATCGGGGTAAAACGGTTTATGACGTGGCTTCCGGCAACGCGTTATTTATTTCTGAACTCGGTCCGTTACCGGAAAATGTTACCTGGTTATCGCCGGAAGGGGAGTTTCAGAAGTGGAACGGCACAGCCTGGGTGAAgg                                         <  1:77006/151‑1 (MQ=31)                                                                                                      tcGTTGAAGATCATCGGGGTAAAACGGTTTATGACGTGGCTTCCGGCAACGCGTTATTTATTTCTGAACTCGGTCCGTTACCGGAAAATGTTACCTGGTTATCGCCGGAAGGGGAGTTTCAGAAGTGGAACGGCACAGCCTGGGTGAAGGa                                        <  1:42469/151‑1 (MQ=31)                                                                                                       cGTTGAAGATCATCGGGGTAAAACGGTTTATGACGTGGCTTCCGGCAACGCGTTATTTATTTCTGAACTCGGTCCGTTACCGGAAAATGTTACCTGGTTATCGCCGGAAGGGGAGTTTCAGAAGTGGAACGGCACAGCCTGGGTGAAGGAt                                       <  1:793920/151‑1 (MQ=31)                                                                                                            aaGATCATCGGGGTAAAACGGTTTATGACGTGGCTTCCGGCAACGCGTTATTTATTTCTGAACTCGGTCCGTTACCGGAAAATGTTACCTGGTTATCGCCGGAAGGGGAGTTTCAGAAGTGGAACGGCACAGCCTGGGTGAAGGATACGGa                                  <  1:265126/151‑1 (MQ=31)                                                                                                            aaGATCATCGGGGTAAAACGGTTTATGACGTGGCTTCCGGCAACGCGTTATTTATTTCTGAACTCGGTCCGTTACCGGAAAATGTTACCTGGTTATCGCCGGAAGGGGAGTTTCAGAAGTGGAACGGCACAGCCTGGGGGAAGGATACGGa                                  <  2:340619/151‑1 (MQ=14)                                                                                                                 catcGGGGTAAAACGGTTTATGACGTGGCTTCCGGCAACGCGTTATTTATTTCTGAACTCGGTCCGTTACCGGAAAATGTTACCTGGTTATCGCCGGAAGGGGAGTTTCAGAAGTGGAACGGCACAGCCTGGGTGAAGGATACGGAAGCAg                             >  2:543816/1‑151 (MQ=31)                                                                                                                      gggTAAAACGGTTTATGACGTGGCTTCCGGCAACGCGTTATTTATTTCTGAACTCGGTCCGTTACCGGAAAATGTTACCTGGTTATCGCCGGAAGGGGAGTTTCAGAAGTGGAACGGCACAGCCTGGGTGAAGGATACGGAAGCAGaaaaa                        >  1:817917/1‑151 (MQ=31)                                                                                                                      gggTAAAACGGTTTATGACGTGGCTTCCGGCAACGCGTTATTTATTTCTGAACTCGGTCCGTTACCGGAAAATGTTACCTGGTTATCGCCGGAAGGGGAGTTTCAGAAGTGGAACGGCACAGCCTGGGTGAAGGATACGGAAGCAGaaaaa                        >  1:819991/1‑151 (MQ=31)                                                                                                                                gTTTATGACGTGGCTTCCGGCAACGCGTTATTTATTTCTGAACTCGGTCCGTTACCGGAAAATGTTACCTGGTTATCGCCGGAAGGGGAGTTTCAGAAGTGGAACGGCACAGCCTGGGTGAAGGATACGGAAGCAGAAAAACTGTTCCGGa              >  1:648338/1‑151 (MQ=31)                                                                                                                                 tttATGACGTGGCTTCCGGCAACGCGTTATTTATTTCTGAACTCGGTCCGTTACCGGAAAATGTTACCTGGTTATCGCCGGAAGGGGAGTTTCAGAAGTGGAACGGCACAGCCTGGGTGAAGGATACGGAAGCAGAAAAACTGTTCCGGAt             <  2:720218/151‑1 (MQ=31)                                                                                                                                 tttATGACGTGGCTTCCGGCAACGCGTTATTTATTTCTGAACTCGGTCCGTTACCGGAAAATGTTACCTGGTTATCGCCGGAAGGGGAGTTTCAGAAGTGGAACGGCACAGCCTGGGTGAAGGATACGGAAGCAGAAAAACTGTTCCGGAt             <  1:195069/151‑1 (MQ=31)                                                                                                                                            gCTTCCGGCAACGCGTTATTTATTTCTGAACTCGGTCCGTTACCGGAAAATGTTACCTGGTTATCGCCGGAAGGGGAGTTTCAGAAGTGGAACGGCACAGCCTGGGTGAAGGATACGGAAGCAGAAAAACTGTTCCGGATCCGGGAGGCgg  >  1:356981/1‑151 (MQ=31)                                                                                                                                            |                                                                                                                                                       GCCTCATACAGGTCTGCCAGCAAACAGTACCGATATTGCACCGCCAGATATTCCGGCTGGCTTCGTGGCTGTTTTCAACAGTGATGAGTCATCGTGGCATCTCGTTGAAGATCATCGGGGTAAAACGGTTTATGACGTGGCTTCCGGCAACGCGTTATTTATTTCTGAACTCGGTCCGTTACCGGAAAATGTTACCTGGTTATCGCCGGAAGGGGAGTTTCAGAAGTGGAACGGCACAGCCTGGGTGAAGGATACGGAAGCAGAAAAACTGTTCCGGATCCGGGAGGCGG  >  NC_000913/1432509‑1432798

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