New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913	1207790 =	11 (0.280)	10 (0.290)	6/244	NT	46.8%	coding (290/630 nt)	stfP	e14 prophage; uncharacterized protein
	?	NC_000913	1209619 =	13 (0.380)					pseudogene (1/501 nt)	stfE	pseudogene, e14 prophage; side tail fiber protein fragment family;Phage or Prophage Related

CAGCGTGGTGGAGGCATCCGGTGCCAGCCTGGAGAAGGGGCGGCACGACCAGTTACTTACCGCACTTCGCGCGCTGCTGTTAAGCCGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCGGTTATGATTGGTGGTATTGAGCAATCAGTAATCATTCAGTGGGGGTGGAATGCCGCAAAAGCATCTGCCTC  >  NC_000913/1207644‑1207934                                                                                                                                                   |                                                                                                                                                 caggctgcaacgttcatctctgtacgttgtagatctgtgttctggtcatttgttgatacgcgtacatagccaataagcatgttaaatccccctggtaaaagcaggaatgatgccatttgcttgttatttcttcattttcataaacgTTg                                                                                                                                                >  2:384855‑M1/147‑149 (MQ=255)   gCGTGGTGGAGGCATCCGGTGCCAGCCTGGAGAAGGGGCGGCACGACCAGTTACTTACCGCACTTCGCGCGCTGCTGTTAAGCCGCAAGAATCCGTTTGGCGATATCAAATCGGTTGGCCCTGTGCAAACGGCTCTGGAAAACCTTGGg                                                                                                                                              >  2:62030‑M1/1‑148 (MQ=255)       gcacagttcaccgcattacgttgtagatctgtgttccggtcatttgttgatacgcgtacatagccaataagcatgttaaatccccctggtaaaagcaggaatgatgccatttgcttgttatttcttcattttcataaacgTTGGTTTgg                                                                                                                                          >  1:528105‑M1/141‑149 (MQ=255)         acagttcagcgcattacgttgtagatctgtgttctggtcatttgttgatacgcgtacatagccaataagcatgttaaatccccctggtaaaagcaggaatgatgccatttgcttgttatttcttcattttcataaacgTTGGTTtgggg                                                                                                                                        <  1:118462‑M1/11‑2 (MQ=255)         acagttcagcgcattacgttgtagatctgtgttctggtcatttgttgatacgcgtacatagccaataagcatgttaaatccccctggtaaaagcaggaatgatgccatttgcttgttatttcttcattttcataaacgTTGGTTTGGga                                                                                                                                        >  1:595941‑M1/139‑149 (MQ=255)            aGGCACCCGGTGCCAGCCTGGAGAAGGGGCCGCACCACCAGTTACTTACCGCACTTCGCGCGCTGCTGTTAAGCCGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAA                                                                                                                                     <  2:48572‑M1/149‑1 (MQ=255)              tcagcgcattacgttgtagatctgtgttctggtcatttgttgatacgcgtacatagccaataagcatgttaaatccccctggtaaaagcaggaatgatgccatttgcttgttatttcttcattttcataaacgTTGGTTTGGGAGAAgg                                                                                                                                   >  2:127574‑M1/134‑149 (MQ=255)              gCATCCGGTGCCAGCCTGGAGAAGGGGCGGCACGACCAGTTACTTACCGCACTTCGCGCGCTGCTGTTAAGCCGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGG                                                                                                                                   >  1:75500‑M1/1‑149 (MQ=255)              gCATCCGGTGCCAGCCTGGAGAAGGGGCGGCACGACCAGTTACTTACCGCACTTCGCGCGCTGCTGTTAAGCCGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGG                                                                                                                                   >  1:267956‑M1/1‑149 (MQ=255)              gCATCCGGTGCCAGCCTGGAGAAGGGGCGGCACGACCAGTTACTTACCGCACTTCGCGCGCTGCTGTTAAGCCGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGG                                                                                                                                   >  1:492234‑M1/1‑149 (MQ=255)              gCATCCGGTGCCAGCCTGGAGAAGGGGCGGCACGAACCGATACTTACCGCACTTCGCGCGCTGCTGTTAAGCCGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGG                                                                                                                                   >  1:338092‑M1/1‑149 (MQ=255)               cATCCGGTGCCAGCCTGGAGAAGGGGCGGCACGACCAGTTACTTACCGCACTTCGCGCGCTGCTGTTAAGCCGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCCCGCGAAAACCTTGGTTTGGGAGAAGGa                                                                                                                                  <  2:151498/149‑1 (MQ=255)                  ccGGTGCCAGCCTGGAGAAGGGGCGGCACGACCAGTTACTTACCGCACTTCGCGCGCTGCTGTTAAGCCGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGctct                                                                                                                               >  2:223072‑M1/1‑145 (MQ=255)                         cAGCCCGGAGAAGGGGCGGCACGACCAGTTACTTACCGCACTTCGCGCGCTGCTGTTAAGCCGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGctctgcattac                                                                                                                        >  2:485187‑M1/1‑138 (MQ=255)                          gttgtagatctgtgttctggtcatttgttgatacgcgtacatagcccataagcatgttaaatccccctggtaaaagcaggaatgatgccatttgcttgttatttcttcattttcataaacgTTGGTTTGGGAGAAGGAGCAAAACTCaa                                                                                                                       >  1:471151‑M1/122‑149 (MQ=255)                          gttgtagatctgtgttctggtcatttgttgatacgcgtacatagccaataagcatgttaaattcccctggtaaaagcaggaatgatgccatttgcttgttatttcttcattttcataaacgTTGGTTTGGGAGAAGGAGCAAAACTCaa                                                                                                                       >  1:468878‑M1/122‑149 (MQ=255)                              tGGAGAAGGGGCGGCACGACCAGTTACTTACCGCACTTCGCGCGCTGCTGTTAAGCCGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCAATgca                                                                                                                   >  1:606923/1‑149 (MQ=255)                                        gcggcACGACCAGTTACTTACCGCACTTCGCGCGCTGCTGTTAAGCCGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGctctgcattacctgttggtgtcccgg                                                                                                         <  2:574072‑M1/149‑27 (MQ=255)                                            cACGACCAGTTACTTACCGCACTTCGCGCGCTGCTGTTAAGCCGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGctctgcattacctgttggtgtcccggttcc                                                                                                     >  1:80586‑M1/1‑119 (MQ=255)                                                   aGTTACTTACCGCACTTCGCGCGCTGCTGTTAAGCCGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGctctgcattacctgttggtgtcccggttccgtggcct                                                                                              >  1:628348‑M1/1‑112 (MQ=255)                                                      cacttgccGCACCTCGCGAGCTGCTGTGACGCCGCAAGAATCCGTTGGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGctctgcattacctgttggtgtcccggttccgtggccttca                                                                                           <  2:14669‑M1/148‑41 (MQ=255)                                                         ttacCCCACTTCGCGCGCTGCTGTTAAGCCGCAAGAATCCGTTTGGCGATATCCAATCGGAAGGCACCGTGCAAACGGCTCTCGAAAACCGTGGTTTGGGAGAAGGctctccattacctgttcgtgcgccggttccgcagccgtcaacc                                                                                        >  2:353761‑M1/1‑106 (MQ=255)                                                                acatagccaataagcatgttaaatccccctggtaaaagcaggagtgatgccatttgcttgttatttcttcattttcataaacgTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAg                                                                                 >  1:10183‑M1/84‑149 (MQ=255)                                                                acatagccaataagcatgttaaatccccctggtaaaagcaggaatgatgccatttgcttgttatttcttcattttcataaacgTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAg                                                                                 >  2:581710‑M1/84‑149 (MQ=255)                                                                      cgcgcTGCTGTTAAGCCGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGctctgcattacctgttggtgtcccggttccgtggccttcagccactccgccaacag                                                                           >  1:433858‑M1/1‑93 (MQ=255)                                                                         gctgctGTTAAGCCGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGAAGCACCGGGTTCATAGCTATACCgg                                                                        >  2:513306/1‑149 (MQ=255)                                                                               gTTAAGCCGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCGGTTATGa                                                                  <  2:75226/149‑1 (MQ=255)                                                                                    gcggcaaGAATCCGGTTCGCGATATCAAATCGGATGGCACTGCGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGctctgcattacctgttggtgtcccggttccgtggccttcagccactccgccaacaggctggctgaaatgc                                                             <  2:339947‑M1/146‑71 (MQ=255)                                                                                     aatccccctggtaaaagcaggaatgatgccatttgcttgttatttcttcattttcataaacgTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCGGTTATGATtggtg                                                            >  1:237202‑M1/63‑149 (MQ=255)                                                                                     aagcccccgtgtaaatggaggaatgatgccatttgcttgttatttcttcattttcataaacgTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCGGTTATGATtggtg                                                            >  2:386675‑M1/63‑149 (MQ=255)                                                                                            aacccGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCGGTTATGATTGGTGGTATTGa                                                     <  2:338092/146‑1 (MQ=255)                                                                                            aaTCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCGGTTATGATTGGTGGTATTGa                                                     <  2:267956/149‑1 (MQ=255)                                                                                             aTCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGctctgcattacctgttggtgtcccggttccgtggccttcagccactccgccaacaggctggctgaaatgcaatggtgcc                                                    >  2:641882‑M1/1‑70 (MQ=255)                                                                                                  ttgggcgaTATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGctctgcattacctgttggtgtcccggttccgtggccttcagccactccgccaacaggctggctgaaatgcaatggtgccgcttt                                               <  1:282219‑M1/146‑85 (MQ=255)                                                                                                   ttGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGctctgcattacctgttggtgtcccggttccgtggccttcagccactccgccaacaggctggctgaaatgcaatggtgccgctttt                                              >  1:401863‑M1/1‑64 (MQ=255)                                                                                                     ggCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGctctgcattacctgttggtgtcccggttccgtggccttcagccactccgccaacaggctggctgaaatgcaatggtgccgctttttc                                            <  2:425177‑M1/149‑88 (MQ=255)                                                                                                       cgctatCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCGGTTATGATTGGTGGTATTGAGCAATCAGTaa                                          <  2:487966/146‑1 (MQ=255)                                                                                                         atatCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCGGTTATGATTGGTGGTATTGAGCAATCAGTAATc                                        >  1:667687/1‑149 (MQ=255)                                                                                                             cAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGctctgcattacctgttggtgtcccggttccgtggccttcagccactccgccaacaggctggctgaaatgcaatggtgccgctttttctgctgaag                                    <  1:223072‑M1/149‑96 (MQ=255)                                                                                                               aaTCGGATGGCACTGTGCAAACGGCTCTAGAAAACCTTGGTTTGGGAGAAGGctctgcattacgtgttggtgtcccggttccgtggccttcagccactccgccaacaggctagctgaaatgcaatggtgccgctttttctgctgaagaa                                  >  1:451186‑M1/1‑52 (MQ=255)                                                                                                                  catttgcttgttatttcttccttttcataaccgTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTACACCGGTTATGATTGGTGGTATTGAGCAATCAGTAATCATTCAGTgg                               >  1:166038‑M1/34‑149 (MQ=255)                                                                                                                      tgcttgttatttcttcattttcataaacgTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCGGTTATGATTGGTGGTATTGAGCAATCAGTAATCATTCAGTGGGGGt                           >  1:562724‑M1/30‑149 (MQ=255)                                                                                                                       ggCACTGTGCAAACGGCTCTCGAAAACATTGGTTTGGGAGGAGGctctgcattgcctgttggtgtcccggttccgtggccttcagccacaccgacaacaggctggctgaaatgaaatggtgccaatttttctgctcaagaataccagga                          >  2:675199‑M1/1‑44 (MQ=255)                                                                                                                        gCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGctctgcattacctgttggtgtcccggttccgtggccttcagccactccgccaacaggctggctgaaatgcaatggtgccgctttttctgctgaagaatacccggaa                         <  2:80586‑M1/149‑107 (MQ=255)                                                                                                                        gCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCGGTTATGATTGGTGGTATTGAGCAATCAGTAATCATTCAGTGGGGGTgg                         <  2:606923/149‑1 (MQ=255)                                                                                                                         cACTGCGCAAACGGCTCTCGAAACCCTTGGTTCGGGAGAAGGAGCAAAACGCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCGGTTATGATTGGTGGTATTGAGCAATCAGTAATCATTCAGTGGGGGTGGa                        <  1:62030/149‑1 (MQ=255)                                                                                                                            tgtgCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCGGTTATGATTGGTGGTATTGAGCAATCAGTAATCATTCAGTGGGGGTGGAATg                     >  1:610485/1‑149 (MQ=255)                                                                                                                             gtgCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGctctgcattacctgttggtgtcccggttccgtggccttcagccactccgccaacaggctggctgaaatgcaatggtgccgctttttctgctgaagaatacccggaactggc                    >  2:183915‑M1/1‑38 (MQ=255)                                                                                                                                    cGGCTCTCGAAAACCTTGGTTTGGGAGAAGGctctgcattacctgttggtgtcccggttccgtggccttcagccactccgccaacaggctggctgaaatgcaatggtgccgctttttctgctgaagaatacccggaactggcaaaggct             <  1:183915‑M1/149‑119 (MQ=255)                                                                                                                                      attttcataaacgTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCGGTTATGATTGGTGGTATTGAGCAATCAGTAATCATTCAGTGGGGGTGGAATGCCGCAAAAGc           <  1:520861‑M1/136‑1 (MQ=255)                                                                                                                                         tttataaccgTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCGGTTATGATTGGTGGTATTGAGCAATCAGTAATCATTCAGTGGGGGTGGAATGCCGCAAAAgcacc        >  2:237656‑M1/11‑147 (MQ=255)                                                                                                                                            gAAAACCTTGGTTTGGGAGAAGGctctgcattacctgttggtgtcccggttccgtggccttcagccactccgccaacaggctggctgaaatgcaatggtgccgctttttctgctgaagaatacccggaactggcaaaggcttatccgac     <  1:485187‑M1/149‑127 (MQ=255)                                                                                                                                             aaaaCCTTGGTTTGGGAGAAGGctctgcattacctgttggtgtcccggttccgtggccttcagccactccgccaacagtctggctgaaatgcaatggtgccgctttttctgctgaagaatacccggaactggcaagggcttatccgaca    >  2:180194‑M1/1‑22 (MQ=255)                                                                                                                                               aaCCTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCGGTTATGATTGGTGGTATTGAGCAATCAGTAATCATTCAGTGGGGGTGGAATGCCGCAAAAGCATCCGcccc  >  2:312118/1‑147 (MQ=255)                                                                                                                                                   |                                                                                                                                                 CAGCGTGGTGGAGGCATCCGGTGCCAGCCTGGAGAAGGGGCGGCACGACCAGTTACTTACCGCACTTCGCGCGCTGCTGTTAAGCCGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCGGTTATGATTGGTGGTATTGAGCAATCAGTAATCATTCAGTGGGGGTGGAATGCCGCAAAAGCATCTGCCTC  >  NC_000913/1207644‑1207934

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