New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_002947 4383136 =46 (0.810)5 (0.090) 3/244 NT 9.9% coding (719/804 nt) PP_3858 tail fiber protein
?NC_002947 4383173 = 48 (0.900)coding (682/804 nt) PP_3858 tail fiber protein

ATCGCACGGAACATCGGAATGGAAGAACGCCAGTCGCGACGGGCTGAAGGAGTACCCCATGATCAGTTCCCTTTGGCATGCCAGGTGTAGGCGGTTGAGGTGTTACCACCCGAGTTGAGGATGGTCACCGAAGAGTTCGTCCACCGCAGGATGCGCACAGTATTGCCCTCGGCCGCCGCATCCTCCAGGTTGGTCAGCCCCACATACCTGCATGCGGTAGGAAAGCTCACCGGGAAGGGGATAGTCTGCTCATTGGGGCCTGCCGCCATGTTGAACACACCGTACTGTTGGACT  >  NC_002947/4382989‑4383282
                                                                                                                                                   |                                                                                                                                                  
aTCGCACGGAACATCGGAATGGAAGAACGCCAGTCGCGACGGGCTGAAGGAGTACCCCATGATCAGTTCCCTTTGGCATGCCAGGTGTAGGCGGTTGAGGTGTTACCACCCGAGTTGAGGATGGTCACCGAAGAGTTCGTCCACCGCa                                                                                                                                                    <  2:40215/148‑1 (MQ=255)
     gcgtcgtgggagcgcgatggctcaaccggcgcagtccaacagtacggtgtgttcaacatggcggcaggccccaatgagcagactatccccttcccggtgagctttcctaccgcatgcaggtatgtggggctgaccaacctggAGGATgc                                                                                                                                              >  5:91759‑M1/143‑149 (MQ=255)
               ggAATGGAAGAACGCCAGTCGCGACGGGCTGAAGGAGTACCCCATGATCAGTTCCCTTTGGCATGCCAGGTGTAGGCGGTTGAGGTGTTACCACCCGAGTTGAGGATGGTCACCGAAGAGTTCGTCCACCGCAGGATGCGCACAGTAtt                                                                                                                                    <  1:247591/149‑1 (MQ=255)
                gAATGGAAGAACGCCAGTCGCGACGGGCTGAAGGAGTACCCCATGATCAGTTCCCTTTGGCATGCCAGGTGTAGGCGGTTGAGGTGTTACCACCCGAGTTGAGGATGGTCACCGAAGAGTTCGTCCACCGCAGGATGCGCACAGTATTg                                                                                                                                   <  1:302878/149‑1 (MQ=255)
                      aagaaCGCCAGTCGCGACGGGCTGAAGGAGTACCCCATGATCAGTTCCCTTTGGCATGCCAGGTGTAGGCGGTTGAGGTGTTACCACCCGAGTTGAGGATGGTCACCGAAGAGTTCGTCCACCGCAGGATGCGca                                                                                                                                           >  4:87387/1‑135 (MQ=255)
                                      gtccaacagtacggtgtgttcaacatggcggcaggccccaatgagcagactatccccttcccggtgagctttcctaccgcatgcaggtatgtggggctgaccaacctggAGGATGCGCACAGTATTGCCCTCGGTCGCCGCATCCTCCa                                                                                                             >  2:338570‑M1/110‑149 (MQ=255)
                                      gtccaacagtacggtgtgttcaacatggcggcaggccccaatgagcagactatccccttcccggtgagctttcctaccgcatgcaggtatgtggggctgaccaacctggAGGATGCGCACAGTATTGCCCTCGGCCGCCGCATCCTCCa                                                                                                             >  5:180030‑M1/110‑149 (MQ=255)
                                      gtacaacagtacggtgtgttcaacatggcggcaggccccaatgagcagactatccccttcccggtgagctttcctaccgcatgcaggtatgtggggctgaccaacctggAGGATGCGCACAGTATTGCCCTCGGTCGCCGCATCCTCCa                                                                                                             >  2:333369‑M1/110‑149 (MQ=255)
                                       cgtgCTGAACGAGTACCCCATGATCAGTTCCCTTTGGCATGCCAGGGGTAGGCGGTTGAGGTGTTACCACCCGAGTTGAGCATGGTGACCGAAGAGTTCGTCCACCGCAGGATGCGCACAGTATTGCCCTCGGCCGCCGCATCCTCCAg                                                                                                            <  8:75660/146‑1 (MQ=255)
                                       cGGGCTGAAGGAGTACCCCATGATCAGTTCCCTTTGGCATGCCAGGTGTAGGCGGTTGAGGTGTTACCACCCGAGTTGAGGATGGTCACCGAAGAGTTCGTCCACCGCAGGATGCGCACAGTATTGCCCTCGGCCGCCGCATCCTCCAg                                                                                                            <  8:28390/149‑1 (MQ=255)
                                         ggCTGAAGGAGTACCCCATGATCAGTTCCCTTTGGCATGCCAGGTGTAGGCGGTTGAGGTGTTACCACCCGAGTTGAGGATGGTCACCGAAGAGTTCGTCCACCGCAGGATGCGCACAGTATTGCCCTCGGCCGCCGCATCCTCCAGGt                                                                                                          >  6:319978/1‑149 (MQ=255)
                                           cTGAAGGAGTACCCCATGATCAGTTCCCTTTGGCATGCCAGGTGTAGGCGGTTGAGGTGTTACCACCCGAGTTGAGGATGGTCACCGAAGAGTTCGTCCACCGCAGGATGCGCACAGTATTGCCCTCGGCCGCCGCATCCTCCAGGTTg                                                                                                        <  1:204450/149‑1 (MQ=255)
                                           cTGAAGGAGTACCCCATGATCAGTTCCCTTTGGCATGCCAGGTGTAGGCGGTTGAGGTGTTACCACCCGAGTTGAGGATGGTCACCGAAGAGTTCGTCCACCGCAGGATGCGCACAGTATTGCCCTCGGCCGCCGCATCCTCCAGGTTg                                                                                                        <  1:125270/149‑1 (MQ=255)
                                            tGAAGGAGTACCCCATGATCAGTTCCCTTTGGCATGCCAGGTGTAGGCGGTTGAGGTGTTACCACCCGAGTTGAGGATGGTCACCGAAGAGTTCGTCCACCGCAGGATGCGCACAGTATTGCCCTCGGCCGCCGCATCCTCCAGGTTgg                                                                                                       >  7:188087/1‑149 (MQ=255)
                                            tGAAGGAGTACCCCATGATCAGTTCCCTTTGGCATGCCAGGTGTAGGCGGTTGAGGTGTTACCACCCGAGTTGAGGATGGTCACCGAAGAGTTCGTCCACCGCAGGATGCGCACAGTATTGCCCTCGGCCGCCGCATCCTCCAGGTTgg                                                                                                       >  4:108812/1‑149 (MQ=255)
                                              aaGGAGTACCCCATGATCAGTTCCCTTTGGCATGCCAGGTGTAGGCGGTTGAGGTGTTACCACCCGAGTTGAGGATGGTCACCGAAGAGTTCGTCCACCGCAGGATGCGCACAGTATTGCCCTCGGCCGCCGCATCCTCCAGGTTGGTc                                                                                                     <  4:326273/149‑1 (MQ=255)
                                               aGGAGTACCCCATGATCAGTTCCCTTTGGCATGCCAGGTGTAGGCGGTTGAGGTGTTACCACCCGAGTTGAGGATGGTCACCGAAGAGTTCGTCCACCGCAGGATGCGCACAGTATTGCCCTCGGCCGCCGCATCCTCCAGGTTGGt                                                                                                      >  8:291418/1‑147 (MQ=255)
                                                ggAGTACCCCATGATCAGTTCCCTTTGGCATGCCAGGTGTAGGCGGTTGAGGTGTTACCACCCGAGTTGAGGATGGTCACCGAAGAGTTCGTCCACCGCAGGATGCGCACAGTATTGCCCTCGGCCGCCGCATCCTCCAGGTTGGTCAg                                                                                                   >  2:362799/1‑149 (MQ=255)
                                                        ccATGATCAGTTCCCTTTGGCATGCCAGGTGTAGGCGGTTGAGGTGTTACCACCCGAGTTGAGGATGGTCACCGAAGAGTTCGTCCACCGCAGGATGCGCACAGTATTGCCCTCGGCCGCCGCATCCTCCAGGTTGGTCAGCCCcaca                                                                                            >  2:76998/1‑148 (MQ=255)
                                                        ccATGATCAGTTCCCTTTGGCATGCCAGGTGTAGGCGGTTGAGGTGTTACCACCCGAGTTGAGGATGGTCACCGAAGAGTTCGTCCACCGCAGGATGCGCACAGTATTGCCCTCGGCCGCCGCATCCTCCAGGTTGGTCAGCCCcaca                                                                                            >  5:85312/1‑148 (MQ=255)
                                                        ccATGATCAGTTCCCTTTGGCATGCCAGGTGTAGGCGGTTGAGGTGTTACCACCCGAGTTGAGGATGGTCACCGAAGAGTTCGTCCACCGCAGGATGCGCACAGTATTGCCCTCGGCCGCCGCATCCTCCAGGTTGGTCAGCCCCACAt                                                                                           <  1:362799/149‑1 (MQ=255)
                                                            acatggcggcaggccccaatgagcagactatccccttcccggtgagctttcctaccgcatgcaggtatgtggggctgaccaacctggAGGATGCGCACAGTATTGCCCTCGGCCGCCGCAtcctc                                                                                                               >  7:233276‑M1/88‑125 (MQ=255)
                                                            acatggcggcaggccccaatgagcagactatccccttcccggtgagctttcctaccgcatgcaggtatgtggggctgaccaacctggAGGATGCGCACAGTATTGCCCTCGGCCGCCGCAtcctc                                                                                                               <  8:233276‑M1/38‑1 (MQ=255)
                                                                  ttCCCTTTGGCATGCCAGGTGTAGGCGGTTGAGGTGTTACCACCCGAGTTGAGGATGGTCACCGAAGAGTTCGTCCACCGCAGGATGCGCACAGTATTGCCCTCGGCCGCCGCATCCTCCAGGTTGGTCAGCCCCACATACCTGCATGc                                                                                 >  4:355019/1‑149 (MQ=255)
                                                                     ccTTTGGCATGCCAGGTGTAGGCGGTTGAGGTGTTACCACCCGAGTTGAGGATGGTCACCGAAGAGTTCGTCCACCGCAGGATGCGCACAGTATTGCCCTCGGCCGCCGCATCCTCCAGGTTGGTCAGCCCCACATACCTGCATGCgg                                                                               >  4:151491/1‑148 (MQ=255)
                                                                     ccTTTGGCATGCCAGGTGTAGGCGGTTGAGGTGTTACCACCCGAGTTGAGGATGGTCACCGAAGAGTTCGTCCACCGCAGGATGCGCACAGTATTGCCCTCGGCCGCCGCATCCTCCAGGTTGGTCAGCCCCACATACCTGCATGCgg                                                                               <  4:190937/148‑1 (MQ=255)
                                                                      cTTTGGCATGCCAGGTGTAGGCGGTTGAGGTGTTACCACCCGAGTTGAGGATGGTCACCGAAGAGTTCGTCCACCGCAGGATGCGCACAGTATTGCCCTCGGCCGCCGCATCCTCGAGGTTGGTCAGCCCCACATACCTGCATGCGGTa                                                                             <  3:151491/149‑1 (MQ=255)
                                                                      cTTTGGCATGCCAGGTGTAGGCGGTTGAGGTGTTACCACCCGAGTTGAGGATGGTCACCGAAGAGTTCGTCCACCGCAGGATGCGCACAGTATTGCCCTCGGCCGCCGCATCCTCCAGGTTGGTCAGCCCCACATACCTGCATGCGGTa                                                                             <  1:76998/149‑1 (MQ=255)
                                                                         tGGCATGCCAGGTGTAGGCGGTTGAGGTGTTACCACCCGAGTTGAGGATGGTCACCGAAGAGTTCGTCCACCGCAGGATGCGCACAGTATTGCCCTCGGCCGCCGCATCCTCCAGGTTGGTCAGCCCCACATACCTGCATGCGGTAg                                                                            <  6:85312/147‑1 (MQ=255)
                                                                           gCATGCCAGGTGTAGGCGGTTGAGGTGTTACCACCCGAGTTGAGGATGGTCACCGAAGAGTTCGTCCACCGCAGGATGCGCACAGTATTGCCCTCGGCCGCCGCATCCTCCAGGTTGGTCAGCCCCACATACCTGCATGCGGTAGGaaa                                                                        <  8:188087/149‑1 (MQ=255)
                                                                            cATGCCAGGTGTAGGCGGTTGAGGTGTTACCACCCGAGTTGAGGATGGTCACCGAAGAGTTCGTCCACCGCAGGATGCGCACAGTATTGCCCTCGGCCGCCGCATCCTCCAGGTTGGTCAGCCCGACATACCTGCATGCGGTAGGAAAg                                                                       <  3:108812/149‑1 (MQ=255)
                                                                                ccAGGTGTAGGCGGTTGAGGTGTTACCACCCGAGTTGAGGATGGTCACCGAAGAGTTCGTCCACCGCAGGATGCGCACAGTATTGCCCTCGGCCGCCGCATCCTCCAGGTTGGTCAGCCCCACATACCTGCATGCGGTAGGAAAGCTCa                                                                   >  8:70913/1‑149 (MQ=255)
                                                                                   ggTGTAGGCGGTTGAGGTGTTACCACCCGAGTTGAGGATGGTCACCGAAGAGTTCGTCCACCGCAGGATGCGCACAGTATTGCCCTCGGCCGCCGCATCCTCCAGGTTGGTCAGCCCCACATACCTGCATGCGGTAGGAAAGCTCACCg                                                                <  7:291418/149‑1 (MQ=255)
                                                                                            ggTTGAGGTGTTACCACCCGAGTTGAGGATGGTCACCGAAGAGTTCGTCCACCGCAGGATGCGCACAGTATTGCCCTCGGCCGCCGCATCCTCCAGGTTGGTCAGCCCCACATACCTGCATGCGGTAGGAAAGCTCACCGGGAAGGGGa                                                       <  5:319978/149‑1 (MQ=255)
                                                                                             gTTGAGGTGTTACCACCCGAGTTGAGGATGGTCACCGAAGAGTTCGTCCACCGCAGGATGCGCACAGTATTGCCCTCGGCCGCCGCATCCTCCAGGt                                                                                                          <  8:155860/97‑1 (MQ=255)
                                                                                                       tACCACCCGAGTTGAGGATGGTCACCGAAGAGTTCGTCCACCGCAGGATGCGCACAGTATTGCCCTCGGCCGCCGCATCCTCCAGGTTGGTCAGCCCCACATACCTGCATGCGGTAGGAAAGCTCACCGGGAAGGGGATAGTCTGCTCa                                            >  5:119439/1‑149 (MQ=255)
                                                                                                             ccGAGTTGAGGATGGTCACCGAAGAGTTCGTCCACCGCAGGATGCGCACAGTATTGCCCTCGGCCGCCGCATCCTCCAGGTTGGTCAGCCCCACATACCTGCATGCGGTAGGAAAGCTCACCGGGAAGGGGATAGTCTGCTCATTg                                         >  8:328312/1‑146 (MQ=255)
                                                                                                             ccGAGTTGAGGATGGTCACCGAAGAGTTCGTCCACCGCAGGATGCGCACAGTATTGCCCTCGGCCGCCGCATCCTCCAGGTTGGTCAGCCCCACATACCTGCATGCGGTAGGAAAGCTCACCGGGAAGGGGATAGTCTGCTCATTg                                         <  7:328312/146‑1 (MQ=255)
                                                                                                                     aGGATGGTCACCGAAGAGTTCGTCCACCGCAGGATGCGCACAGTATTGCCCTCGGCCGCCGCATCCTCCAGGTTGGTCAGCCCCACATACCTGCATGCGGTAGGAAAGCTCACCGGGAAGGGGATAGTCTGCTCATTGGGGCCTgccgc                              >  1:53856/1‑149 (MQ=255)
                                                                                                                      ggATGGTCACCGAAGAGTTCGTCCACCGCAGGATGCGCACAGTATTGCCCTCGGCCGCCGCATCCTCCAGGTTGGTCAGCCCCACATACCTGCATGCGGTAGGAAAGCTCACCGGGAAGGGGATAGTCTGCTCATTGGGGCCTgccgc                              >  7:33427/1‑148 (MQ=255)
                                                                                                                       gATGGTCACCGAAGAGTTCGTCCACCGCAGGATGCGCACAGTATTGCCCTCGGCCGCCGCATCCTCCAGGTTGGTCTGCCCCACATACCTGCATGCGGTAGGAAAGCTCACCGGGAAGGGGATAGTCTGCTCATTGGGGCCTGCCGCCa                            >  3:239207/1‑149 (MQ=255)
                                                                                                                         tGGTCACCGAAGAGGTCGTCCACCGCAGGATGCGCACAGTATTGCCCGCGGCCGCCGCATCCTCCAGGTTGGTCAGCCCCACATACCTGCATGCGGTAGGAAAGCTCACCGGGAAGGGGATAGTCTGCTCATTGGGGCCTGCCGCCATg                          <  3:272809/149‑1 (MQ=255)
                                                                                                                            tCACCGAAGAGTTCGTCCACCGCAGGATGCGCACAGTATTGCCCTCGGCCGCCGCATCCTCCAGGTTGGTCAGCCCCACATACCTGCATGCGGTAGGAAAGCTCACCGGGAAGGGGATAGTCTGCTCATTGGGGCCTGCCGCCATGtt                        >  4:348622/1‑148 (MQ=255)
                                                                                                                               ccGAAGAGTTCGTCCACCGCAGGATGCGCACAGTATTGCCCTCGGCCGCCGCATCCTCCAGGTTGGTCAGCCCCACATACCTGCATGCGGTAGGAAAGCTCACCGGGAAGGGGATAGTCTGCTCATTgg                                        <  6:51246/129‑1 (MQ=255)
                                                                                                                               ccGAAGAGTTCGTCCACCGCAGGATGCGCACAGTATTGCCCTCGGCCGCCGCATCCTCCAGGTTGGTCAGCCCCACATACCTGCATGCGGTAGGAAAGCTCACCGGGAAGGGGATAGTCTGCTCATTgg                                        >  5:51246/1‑129 (MQ=255)
                                                                                                                                   agagTTCGTCCACCGCAGGATGCGCACAGTATTGCCCTCGGCCGCCGCATCCTCCAGGTTGGTCAGCCCCACATACCTGCATGCGGTAGGAAAGCTCACCGGGAAGGGGATAGTCTGCTCATTGGGGCCTGCCGCCATGTTGAacacac                <  3:348622/149‑1 (MQ=255)
                                                                                                                                   agagTTCGTCCACCGCAGGATGCGCACAGTATTGCCCTCGGCCGCCGCATCCTCCAGGTTGGTCAGCCCCACATACCTGCATGCGGTAGGAAAGCTCACCGGGAAGGGGATAGTCTGCTCATTGGGGCCTGCCGCCATGTTGAacacac                >  8:61014/1‑149 (MQ=255)
                                                                                                                                       ttCGTCCACCGCAGGATGCGCACAGTATTGCCCTCGGCCGCCGCATCCTCCAGGTTGGTCAGCCCCACATACCTGCATGCGGTAGGAAAGCTCACCGGGAAGGGGATAGTCTGCTCATTGGGGCCTGCCGCCATGTTGAACACACCGTa            >  1:54393/1‑149 (MQ=255)
                                                                                                                                          gTCCACCGCAGGATGCGCACAGTATTGCCCTCGGCCGCCGCATCCTCCAGGTTGGTCAGCCCCACATACCTGCATGCGGTAGGAAAGCTCACCGGGAAGGGGATAGTCTGCTCATTGGGGCCTGCCGCCATGTTGAACACACCGTACTg         >  3:162466/1‑149 (MQ=255)
                                                                                                                                           tCCACCGCAGGATGCGCACAGTATTGCCCTCGGCCGCCGCATCCTCCAGGTTGGTCAGCCCCACATACCTGCATGCGGTAGGAAGGCTCACCGGGAAGGGGATAGTCTGCTCATTGGGGCCTGCCGCCATGTTGAACACACCGTACTGt        >  6:302201/1‑149 (MQ=255)
                                                                                                                                            ccACCGCAGGATGCGCACAGTATTGCCCTCGGCCGCCGCATCCTCCAGGTTGGTCAGCCCCACATACCTGCATGCGGTAGGAAAGCTCACCGGGAAGGGGATAGTCTg                                                <  2:336164/108‑1 (MQ=255)
                                                                                                                                            ccACCGCAGGATGCGCACAGTATTGCCCTCGGCCGCCGCATCCTCCAGGTTGGTCAGCCCCACATACCTGCATGCGGTAGGAAAGCTCACCGGGAAGGGGATAGTCTg                                                >  1:336164/1‑108 (MQ=255)
                                                                                                                                             cACCGCAGGATGCGCACAGTATTGCCCTCGGCCGCCGCATCCTCCAGGTTGGTCAGCCCCACATACCTGCATGCGGTAGGAAAGCTCACCGGGAAGGGGATAGTCTGCTCATTGGGGCCTGCCGCCATGTTGAACACACCGTACTGtt       <  7:70913/148‑1 (MQ=255)
                                                                                                                                               ccGCAGGATGCGCACAGTATTGCCCTCGGCCGCCGCATCCTCCAGGTTGGTCAGCCCCACATACCTGCATGCGGTAGGAAGGCTCACCGGGAAGGGGATAGTCTGCTCATTGGGGCCTGCCGCCATGTTGAACACACCGTACTGTTGGa    <  5:302201/149‑1 (MQ=255)
                                                                                                                                                 gCAGGATGCGCACAGTATTGCCCTCGGCCGCCGCATCCTCCAGGTTGGTCAGCCCCACATACCTGCATGCGGTAGGAAAGCTCACCGGGAAGGGGATAGTCTGCTCATTGGGGCCTGCCGCCATGTTGAACACACCGTACTGTTGGACt  <  8:33427/149‑1 (MQ=255)
                                                                                                                                                  cAGGATGCGCACAGTATTGCCCTCGGCCGCCGCATCCTCCAGGTTGGTCAGCCCCACATACCTGCATGCGGTAGGAAAGCTCACCGGGAAGGGGATAGTCTGCTCATTGGGGCCTGCCGCCATGTTGAACACACCGTACTGTTg      <  7:78397/144‑1 (MQ=255)
                                                                                                                                                  cAGGATGCGCACAGTATTGCCCTCGGCCGCCGCATCCTCCAGGTTGGTCAGCCCCACATACCTGCATGCGGTAGGAAAGCTCACCGGGAAGGGGATAGTCTGCTCATTGGGGCCTGCCGCCATGTTGAACACACCGTACTGTTg      >  8:78397/1‑144 (MQ=255)
                                                                                                                                                   |                                                                                                                                                  
ATCGCACGGAACATCGGAATGGAAGAACGCCAGTCGCGACGGGCTGAAGGAGTACCCCATGATCAGTTCCCTTTGGCATGCCAGGTGTAGGCGGTTGAGGTGTTACCACCCGAGTTGAGGATGGTCACCGAAGAGTTCGTCCACCGCAGGATGCGCACAGTATTGCCCTCGGCCGCCGCATCCTCCAGGTTGGTCAGCCCCACATACCTGCATGCGGTAGGAAAGCTCACCGGGAAGGGGATAGTCTGCTCATTGGGGCCTGCCGCCATGTTGAACACACCGTACTGTTGGACT  >  NC_002947/4382989‑4383282

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.