New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_002947	= 250236	36 (0.800)	5 (0.120)	4/228	NT	13.3%	coding (967/1347 nt)	nfeD	membrane bound peptidase NefD
	?	NC_002947	= 250260	32 (0.770)					coding (943/1347 nt)	nfeD	membrane bound peptidase NefD

TTGAGAGCCATGCCCAGCACGCCACCGATCAGGAGTGCCGAGAGCAGCGCCAAGGTGCCGATCAAGGCCAGGGGGATGCCGAAACCGGGAGCGTCGGTGTCGATCAGGATCAATGCGCCGACCACAAAGGCAACGATGCCGCCGAAGCCGACCACGCCGAAACTGGGCAAGAAGGCTTCGGCGATCATGAAGGTGATGCCAAGCAGGATCAGCGCCACCCCGGCAAAACTCACCGGCAACAGCTGCAGCGCATACAGCGCCAGCAGCAGGCAGATACCACCGACCACACCGCCCACCG  >  NC_002947/250088‑250385                                                                                                                                                     |                                                                                                                                                      ttGAGAGCCATGCCCAGCACGCCACCGATCAGGAGTGCCGAGAGCAGCGCCAAGGTGCCGATCAAGGCCAGGGGGATGCCGAAACCGGGAGCGTCGGTGTCGATCAGGATCAATGCGCCGACCACAAAGGCAACGATGCCGCCGAAGCCt                                                                                                                                                      >  2:676336‑M1/1‑149 (MQ=255)        ccATGCCCAGCACGCCACCGATCAGGAGTGCCGAGAGCAGCGCCAAGGTGCCGATCAAGGCCAGGGGGATGCCGAAACCGGGAGCGTCGGTGTCGATCAGGATCAATGCGCAGACCACAAAGGCAACGATGCCGCCGAAGCCGACCACGc                                                                                                                                               >  1:313229/1‑150 (MQ=255)             cccAGCACGCCACCGATCAGGAGTGCCGAGAGCAGCGCCAAGGTGCCGATCAAGGCCAGGGGGATGCCGAAACCGGGAGCGTCGGTGTCGATCAGGATCAATGCGCCGACCACAAAGGCAACGATGCCGCCGAAGCCGACCACGCCGaaa                                                                                                                                          >  2:1072929/1‑150 (MQ=255)              ccAGCACGCCACCGATCAGGAGTGCCGAGAGCAGCGCCAAGGTGCCGATCAAGGCCAGGGGGATGCCGAAACCGGGAGCGTCGGTGTCGATCAGGATCAATGCGCCGACCACAAAGGCAACGATGCCGCCGAAGCCttcttgcccagttt                                                                                                                                         >  1:646404‑M1/1‑136 (MQ=255)                         ccGATCAGGAGTGCCGAGAGCAGCGCCAAGGTGCCGATCAAGGCCAGGGGGATGCCGAAACCGGGAGCGTCGGTGTCGATCAGGATCAATGCGCCGACCACAAAGGCAACGATGCCGCCGAAGCCGACCACGCCGAAACTGGGCAAGAAg                                                                                                                              >  1:24549‑M2/1‑149 (MQ=255)                             tCAGGAGTGCCGAGAGCAGCGCCAAGGTGCCGATCAAGGCCAGGGGGATGCCGAACCCGGGAGCGGCGGTGTCGATCAGCATCAATGCTCCAACCACAAGGGCAACGATGGCGCCGAAGCCGACCACGCCGAAACTGGGCAAGAAggctt                                                                                                                          >  1:929425‑M2/1‑145 (MQ=255)                               aGGAGTGCCGAGAGCAGCGCCAAGGTGCCGATCAAGGCCAGGGGGATGCCGAAACCGGGAGCGTCGGTGTCGATCAGGATCAATGCGCCGACCACAAAGGCAACGATGCCGCCGAAGCCttcttgcccagtttcggcgtggtcggcttcg                                                                                                                        <  1:676336‑M1/150‑32 (MQ=255)                                ggAGTGCCGAGAGCAGCGCCAAGGTGCCGATCAAGGCCAGGGGGATGCCGAAACCGGGAGCGTCGGTGTCGATCAGGATCAATGCGCCGACCACAAAGGCAACGATGCCGCCGAAGCCGACCACGCCGAAACTGGGTAAGAAggattcgg                                                                                                                       <  1:849751‑M2/150‑9 (MQ=255)                                ggAGTGCCGAGAGCAGCGCCAAGGTGCCGATCAAGGCCAGGGGGATGCCGAAACCGGGAGCGTCGGTGTCGATCAGGATCAATGCGCCGACCACAAAGGCAACGATGCCGCCGAAGCCGACCACGCCGAAACTGGGCAAGAAggcttcgg                                                                                                                       >  2:455144‑M2/1‑142 (MQ=255)                                   gTGCCGAGAGCAGCGCCAAGGTGCCGATCAAGGCCAGGGGGATGCCGAAACCGGGAGCGTCGGTGTCGATCAGGATCAATGCGCCGACCACAAAGGCAACGATGCCGCCGAAGCCGACCACGCCGAAACTGGGCAAGAAGGCTTCGGCGa                                                                                                                    <  1:825797/150‑1 (MQ=255)                                     gCCGAGAGCAGCGCCAAGGTGCCGATTAAGGCCAGGGGGATGCCGAAACCGGGAGCGTCGGTGTCGATCAGGATCAATGCGCCGACCACAAAGGCAACGATGCCGCCGAAGCCGACCACGCCGAAACTGGGCAAGAAGGCTTCGGCGATc                                                                                                                  <  1:305486/150‑1 (MQ=255)                                     gCCGAGAGCAGCGCCAAGGTGCCGATCAAGGCCAGGGGGATGCCGAAACCGGGAGCGTCGGTGTCGATCAGGATCAATGCGCCGACCACAAAGGCAACGATGCCGCCGAAGCCGACCACGCCGAAACTGGGCAAGAAGGCTTCGGCGATc                                                                                                                  <  2:949074/150‑1 (MQ=255)                                            gcagcGCCAAGGTGCCGATCAAGGCCAGGGGGATGCCGAAACCGGGAGCGTCGGTGTCGATCAGGATCAATGCGCCGACCACAAAGGCAACGATGCCGCCGAAGCCGACCACGCCGAAACTGGGCAAGAAGGCTTCGGCGATCATGAAgg                                                                                                           <  2:1183538/150‑1 (MQ=255)                                            gcagcGCCAAGGTGCCGATCAAGGCCAGGGGGATGCCGAAACCGGGAGCGTCGGTGTCGATCAGGATCAATGCGCCGACCACAAAGGCAACGATGCCGCCGAAGCCGACCACGCCGAAACTGGGCAAGAAGGCTTCGGCGATCATGAAgg                                                                                                           <  1:113483/150‑1 (MQ=255)                                                ggcCAAGGGGCCGAGAAAGGCCAGGGGGATGGCGAAACCGGGAGCGTCGGTGGCGATCAGGATAAATGCGCGGACCACAAAGGCAACGAGGCCGCCGAAGCGGACAACGCCGAAACGGGGCAAGAAGGCTTCGGCGATCATGAAGGTGAt                                                                                                       <  2:209566/149‑1 (MQ=255)                                                 gcCAAGGTGCCGATCAAGGCCAGGGGGATGCCGAAACCGGGAGCGTCGGTGTCGATCAGGATCAATGCGCCGACCACAAAGGCAACGATGCCGCCGAAGCCttcttgcccagtttcggcgtggtcggcttcggcggcatcgttgcctttg                                                                                                      >  2:146530‑M1/1‑101 (MQ=255)                                                   cAAGGTGCCGATCAAGGTCAGGGGGATGACGAAATCGGGAGCGTCTGTGTCTTTCAGGATTAATGCGCCGACAACAAAGTCAACGATGCCGCCGAAGCCGACCACGCCGAAACTGGGTAAGAAGTCTTCGGCGATCATGAAGGTGATGcc                                                                                                    <  1:289898/150‑1 (MQ=255)                                                   cAAGGTGCCGATCAAGGCCAGGGGGATGCCGAAACCGGGAGCGTCGGTGTCGATCAGGATCAATGCGCCGACCACAAAGGCAACGATGCCGCCGAAGCCttcttgcccagtttcggcgtggtcggcttcggcggcatcgttgcctttgtg                                                                                                    <  2:646404‑M1/150‑52 (MQ=255)                                                   cAAGGTGCCGATCAAGGCCAGGGGGATGCCGAAACCGGGAGCGTCGGTGTCGATCAGGATCAATGCGCCGACCACAAAGGCAACGATGCCGCCGAAGCCttcttgcccagtttcggcgtggtcggcttcggcggcatcgttgcctttgtg                                                                                                    <  1:146530‑M1/150‑52 (MQ=255)                                                    aaGGTGCCGATCAAGGCCAGGGGGATGCCGAAACCGGGAGCGTCGGTGTCGATCAGGATCAATGCGCCGACCACAAAGGCAACGATGCCGCCGAAGCCGACCACGCCGAAACTGGGCAAGAAGGCTTCGGCGATCATGAAGGTGATGCCa                                                                                                   >  2:496369/1‑150 (MQ=255)                                                      ggTGCCGATCAAGGCCAGGGGGATGCCGAAACCGGGAGCGTCGGTGTCGATCAGGATCAATGCGCCGACCACAAAGGCAACGATGCCGCCGAAGCCGACCACGCCGAAACTGGGCAAGAAGGCTTCGGCGATCATGAAGGTGATGCCAAg                                                                                                 <  1:496369/150‑1 (MQ=255)                                                      ggTGCCGATCAAGGCCAGGGGGATGCCGAAACCGGGAGCGTCGGTGTCGATCAGGATCAATGCGCCGACCACAAAGGCAACGATGCCGCCGAAGCCGACCACGCCGAAACTGGGCAAGAAGGCTTCGGCGATCATGAAGGTGATGCCAAg                                                                                                 <  1:723654/150‑1 (MQ=255)                                                          ccGATCAAGGCCAGGGGGATGCCGAAACCGGGAGCGTCGGTGTCGATCAGGATCAATGCGCCGACCACAAAGGCAACGATGCCGCCGAAGCCGa                                                                                                                                                     <  2:1130314/94‑1 (MQ=255)                                                          ccGATCAAGGCCAGGGGGATGCCGAAACCGGGAGCGTCGGTGTCGATCAGGATCAATGCGCCGACCACAAAGGCAACGATGCCGCCGAAGCCGa                                                                                                                                                     >  1:1130314/1‑94 (MQ=255)                                                            gATCAAGGCCAGGGGGATGCCGAAACCGGGAGCGTCGGTGTCGATCAGGATCAATGCGCCGACCACAAAGGCAACGATGCCGCCGAAGCCGACCACGCCGAAACTGGGCAAGAAGGCTTCGGCGATCATGAAGGTGATGCCAAGCAGGAt                                                                                           >  2:1089934/1‑150 (MQ=255)                                                                aaGGCCAGGGGGATGCCGAAACCGGGAGCGTCGGTGTCGATCAGGATCAATGCGCCGACCACAAAGGCAACGATGCCGCCGAAGCCGACCACGCCGAAACTGGGCAAGAAGGCTTCGGCGATCATGAAGGTGATGCCAAGCAGGATCAgc                                                                                       <  2:992832/150‑1 (MQ=255)                                                                      aGGGGGATGCCGAAACCGGGAGCGTCGGTGTCGATCAGGATCAATGCGCCGACCACAAAGGCAACGATGCCGCCGAAGCCGACCACGCCGAAACTGGGCAAGAAGGCTTCGGCGATCATGAAGGTGATGCCAAGCAGGATCAGCGCCAcc                                                                                 >  2:274184/1‑150 (MQ=255)                                                                         gggATGCCGAAACCGGGAGCGTCGGTGTCGATCAGGATCAATGCGCCGACCACAAAGGCAACGATGCCGCCGAAGCCGACCACGCCGAAACTGGGCAAGAAGGCTTCGGCGATCATGAAGGTGATGCCAAGCAGGATCAGCGCCACCCCg                                                                              <  1:237438/150‑1 (MQ=255)                                                                            aTGCCGAAACCGGGAGCGTCGGTGTCGATCAGGATCAATGCGCCGACCACAAAGGCAACGATGCCGCCGAAGCCGACCACGCCGAAACTGGGCAAGAAGGCTTCGGCGATCATGAAGGTGATGCCAAGCAGGATCAGCGCCACCCCGGCa                                                                           <  1:1072929/150‑1 (MQ=255)                                                                               ccGAAACCGGGAGCGTCGGTGTCGATCAGGATCAATGCGCCGACCACAAAGGCAACGATGCCGCCGAAGCCGACCACGCCGAAACTGGGCAAGAAGGCTTCGGCGATCATGAAGGTGATGCCAAGCAGGATCAGCGCCACCCCGGCaaaa                                                                        <  2:504935/150‑1 (MQ=255)                                                                                    aCCGGGAGAGTCGGTGTCGATCAGGATCAATGCGCCGACCACAAAGGCAACGATGCCGCCGAAGCCGACCACGCCGAAACTGGGCAAGAAGGCTTCGGCGATCATGAAGGTGATGCCAAGCAGGATCAGCGCCACCCCGGCAAAACTCAc                                                                   <  2:40303/150‑1 (MQ=255)                                                                                      cGGGAGCGTCGGTGTCGATCAGGATCAATGCGCCGACCACAAAGGCAACGATGCCGCCGAAGCCGACCACGCCGAAACTGGGCAAGAAGGCTTCGGCGATCAt                                                                                                                <  1:196821/103‑1 (MQ=255)                                                                                      cGGGAGCGTCGGTGTCGATCAGGATCAATGCGCCGACCACAAAGGCAACGATGCCGCCGAAGCCGACCACGCCGAAACTGGGCAAGAAGGCTTCGGCGATCAt                                                                                                                >  2:196821/1‑103 (MQ=255)                                                                                      cGGGAGCGTCGGTGTCGATCAGGATCAATGCGCCGACCACAAAGGCAACGATGCCGCCGAAGCCGACCACGCCGAAACTGGGCAAGAAGGCTTCGGCGATCATGAAGGTGATGCCAAGCAGGATCAGCGCCa                                                                                   <  1:355966/132‑1 (MQ=255)                                                                                      cGGGAGCGTCGGTGTCGATCAGGATCAATGCGCCGACCACAAAGGCAACGATGCCGCCGAAGCCGACCACGCCGAAACTGGGCAAGAAGGCTTCGGCGATCATGAAGGTGATGCCAAGCAGGATCAGCGCCa                                                                                   >  2:355966/1‑132 (MQ=255)                                                                                      cGGGAGCGTCGGTGTCGATCAGGATCAATGCGCCGACCACAAAGGCAACGATGCCGCCGAAGCCGACCACGCCGAAACTGGGCAAGAAGGCTTCGGCGATCATGAAGGTGATGCCAAGCAGGATCAGCGCCACCCCGGCAAAACTCACCg                                                                 <  1:37926/150‑1 (MQ=255)                                                                                            cGTCGCTGTCTATCCGCACCAAGGCGCGGACCACAAAGGCAACGATGCCGCCGAAGCCGACCACGCCGAAACTGCGCAAGAAGCCTTCAGCGATGATGAAGGTGATCCCAAGCAGGATCAGCGCCACCCCGGCAAAACTCACCGGCAACa                                                           <  2:313229/150‑1 (MQ=255)                                                                                                ggTGTCGATCAGGATCAATGCGCCGACCACAAAGGCAACGATGCCGCCGAAGCCGACCACGCCGAAACTGGGCAAGAAGGCTTCGGCGATCATGAAGGTGATGCCAAGCAGGATCAGCGCCACCCCGGCAAAACTCACCGGCAACAGCTg                                                       <  1:274184/150‑1 (MQ=255)                                                                                                      gATCAGGATCAATGCGCCGACCACAAAGGCAACGATGCCGCCGAAGCCGACCACGCCGAAACTGGGCAAGAAGGCTTCGGCGATCATGAAGGTGATGCCAAGCAGGATCAGCGCCACCCCGGCAAAACTCACCGGCAACAGCTGCAgcgc                                                 >  1:463988/1‑150 (MQ=255)                                                                                                                aaTGCGCCGACCACAAAGGCAACGATGCCGCCGAAGCCGACCACGCCGAAACTGGGCAAGAAGGCTTCGGCGATCATGAAGGTGATGCCAAGCAGGATCAGCGCCACCCCGGCAAAACTCACCGGCAACAGCTGCAGCGCATACAGCGcc                                       >  1:695366/1‑150 (MQ=255)                                                                                                                           cacaAAGGCAACGATGCCGCCGAAGCCGACCACGCCGAAACTGGGCAAGAAggcttcggcggcatcgttgcctttgtggtcggcgcattgatcctgatcgacaccgacgctcccggtttcggcatccccctggccttgatcggcaccttg                            <  2:646404‑M2/150‑100 (MQ=255)                                                                                                                           cacaAAGGCAACGATGCCGCCGAAGCCGACCACGCCGAAACTGGGCAAGAAggcttcggcggcatcgttgcctttgtggtcggcgcattgatcctgatcgacaccgacgctcccggtttcggcatccccctggccttgatcggcaccttg                            <  1:146530‑M2/150‑100 (MQ=255)                                                                                                                            acaAAGGCAACGATGCCGCCGAAGCCGACCACGCCGAAACTGGGCAAGAAGGCTTCGGCGATCATGAAGGTGATGCCAAGCAGGATCAGCGCCACCCCGGCAAAACTCACCGGCAACAGCTGCAGCGCATACAGCGCCAGCAgcaggcag                           >  1:1065193/1‑150 (MQ=255)                                                                                                                            acaAAGGCAACGATGCCGCCGAAGCCGACCACGCCGAAACTGGGCAAGAAGGCTTCGGCGATCATGAAGGTGATGCCAAGCAGGATCAGCGCCACCCCGGCAAAACTCACCGGCAACAGCTGCAGCGCATACAGCGCCAGCAgcaggcag                           >  2:380151/1‑150 (MQ=255)                                                                                                                            acaAAGGCAACGATGCCGCCGAAGCCGACCACGCCGAAACTGGGCAAGAAGGCTTCGGCGATCATGAAGGTGATGCCAAGCAGGATCAGCGCCACCCCGGCAAAACTCACCGGCAACAGCTGCAGCGCATACAGCGCCAGCAgcagg                              <  2:758795/147‑1 (MQ=255)                                                                                                                            acaAAGGCAACGATGCCGCCGAAGCCGACCACGCCGAAACTGGGCAAGAAGGCTTCGGCGATCATGAAGGTGATGCCAAGCAGGATCAGCGCCACCCCGGCAAAACTCACCGGCAACAGCTGCAGCGCATACAGCGCCAGCAgcagg                              >  1:758795/1‑147 (MQ=255)                                                                                                                             caAAGGCAACGATGCCGCCGAAGCCGACCACGCCGAAACTGGGCAAGAAggcttcggcggcatcgttgcctttgtggtcggcgcattgatcctgatcgacaccgacgctcccggtttcggcatccccctggccttgatcggcaccttggc                          >  2:146530‑M2/1‑49 (MQ=255)                                                                                                                                aGGCAACGATGCCGCCGAAGCCGACCACGCCGAAACTGGGCAAGAAGGCTTCGGCGATCATGAAGGTGATGCCAAGCAGGATCAGCGCCACCCCGGCAAAACTCACCGGCAACAGCTGCAGCGCATACAGCGCCAGCAGCAGGCAGATac                       <  2:463988/150‑1 (MQ=255)                                                                                                                                    aaCGATGCCGCCGAAGCCGACCACGCCGAAACTGGGCAAGAAGGCTTCGGCGATCATGAAGGTGATGCCAAGCAGGATCAGCGCCACCCCGGCAAAACTCACCGGCAACAGCTGCAGCGCATACAGCGCCAGCAGCAGGCAGATaccacc                   >  2:668841/1‑150 (MQ=255)                                                                                                                                    aaCGATGCCGCCGAAGCCGACCACGCCGAAACTGGGCAAGAAGGCTTCGGCGATCATGAAGGTGATGCCAAGCAGGATCAGCGCCACCCCGGCAAAACTCACCGGCAACAGCTGCAGCGCATACAGCGCCAGCAGCAGGCAGATaccacc                   >  2:955786/1‑150 (MQ=255)                                                                                                                                              ccgAATCCttcttacccagtttcggcgtggtcggcttcggcggcatcgttgcctttgtggtcggcgcattgatcctgatcgacaccgacgctcccggtttcggcatccccctggccttgatcggcaccttggcgctgctctcggcactcc         <  1:849751‑M1/150‑143 (MQ=255)                                                                                                                                              ccgAAGCCttcttgcccagtttcggcgtggtcggcttcggcggcatcgttgcctttgtggtcggcgcattgatcctgatcgacaccgacgctcccggtttcggcatccccctggccttgatcggcaccttggcgctgctctcggcactcc         >  2:455144‑M1/1‑8 (MQ=255)                                                                                                                                              ccgAAGCCGACCACGCCGAAACTGGGCAAGAAGGCTTCGGCGATCATGAAGGTGATGCCAAGCAGGATCAGCGCCACCCCGGCAAAACTCACCGGCAACAGCTGCAGCGCATACAGCGCCAGCAGCAGGCAGATACCACCGACCACACCg         >  1:355127/1‑150 (MQ=255)                                                                                                                                               cgAAGCCGACCACGCCGAAACTGGGCAAGAAggcttcggcggcatcgttgcctttgtggtcggcgcattgatcctgatcgacaccgacgctcccggtttcggcatccccctggccttgatcggcaccttggcgctgctctcggcactcct        <  1:676336‑M2/150‑120 (MQ=255)                                                                                                                                                 aaGCCttcttgcccagtttcggcgtggtcggcttcggcgccatcgttgcccttgtggttggagcattgatgctgatcgacaccgccgctcccgggttcggcatccccctggccttgatcggcaccttggcgctgctctcggcactcctga      >  1:929425‑M1/1‑5 (MQ=255)                                                                                                                                                   gCCGACCACGCCGAAACTGGGCAAGAAGGCTTCGGCGATCATGAAGGTGATGCCAAGCAGGATCAGCGCCACCCCGGCAAAACTCACCGGCAACAGCTGCAGCGCATACAGCGCCAGCAGCAGGCAGATACCACCGACCACACCGCCCAc    <  1:650036/150‑1 (MQ=255)                                                                                                                                                     cttcttgcccagtttcggcgtggtcggcttcggcggcatcgttgcctttgtggtcggcgcattgatcctgatcgacaccgacgctcccggtttcggcatccccctggccttgatcggcaccttggcgctgctctcggcactcctgatcgg  >  1:24549‑M1/1‑1 (MQ=255)                                                                                                                                                     |                                                                                                                                                      TTGAGAGCCATGCCCAGCACGCCACCGATCAGGAGTGCCGAGAGCAGCGCCAAGGTGCCGATCAAGGCCAGGGGGATGCCGAAACCGGGAGCGTCGGTGTCGATCAGGATCAATGCGCCGACCACAAAGGCAACGATGCCGCCGAAGCCGACCACGCCGAAACTGGGCAAGAAGGCTTCGGCGATCATGAAGGTGATGCCAAGCAGGATCAGCGCCACCCCGGCAAAACTCACCGGCAACAGCTGCAGCGCATACAGCGCCAGCAGCAGGCAGATACCACCGACCACACCGCCCACCG  >  NC_002947/250088‑250385

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