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

CCGATCAGGAGTGCCGAGAGCAGCGCCAAGGTGCCGATCAAGGCCAGGGGGATGCCGAAACCGGGAGCGTCGGTGTCGATCAGGATCAATGCGCCGACCACAAAGGCAACGATGCCGCCGAAGCCGACCACGCCGAAACTGGGCAAGAAGGCTTCGGCGATCATGAAGGTGATGCCAAGCAGGATCAGCGCCACCCCGGCAAAACTCACCGGCAACAGCTGCAGCGCATACAGCGCCAGCAGCAGGCAGATACCACCGACCACACCGCCCACCGCAGAACCTGGGTTCATGAATTCGA  >  NC_002947/250112‑250409                                                                                                                                                     |                                                                                                                                                      ccGATCAGGAGTGCCGAGAGCAGCGCCAAGGTGCCGATCAAGGCCAGGGGGATGCCGAAACCGGGAGCGTCGGTGTCGATCAGGATCAATGCGCCGACCACAAAGGCAACGATGCCGCCGAAGCCGACCACGCCGAAACTGGGCAAGAAg                                                                                                                                                      >  1:24549‑M2/1‑149 (MQ=255)     tCAGGAGTGCCGAGAGCAGCGCCAAGGTGCCGATCAAGGCCAGGGGGATGCCGAACCCGGGAGCGGCGGTGTCGATCAGCATCAATGCTCCAACCACAAGGGCAACGATGGCGCCGAAGCCGACCACGCCGAAACTGGGCAAGAAggctt                                                                                                                                                  >  1:929425‑M2/1‑145 (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                                                                                                                                   <  1:113483/150‑1 (MQ=255)                    gcagcGCCAAGGTGCCGATCAAGGCCAGGGGGATGCCGAAACCGGGAGCGTCGGTGTCGATCAGGATCAATGCGCCGACCACAAAGGCAACGATGCCGCCGAAGCCGACCACGCCGAAACTGGGCAAGAAGGCTTCGGCGATCATGAAgg                                                                                                                                   <  2:1183538/150‑1 (MQ=255)                        ggcCAAGGGGCCGAGAAAGGCCAGGGGGATGGCGAAACCGGGAGCGTCGGTGGCGATCAGGATAAATGCGCGGACCACAAAGGCAACGAGGCCGCCGAAGCGGACAACGCCGAAACGGGGCAAGAAGGCTTCGGCGATCATGAAGGTGAt                                                                                                                               <  2:209566/149‑1 (MQ=255)                           cAAGGTGCCGATCAAGGTCAGGGGGATGACGAAATCGGGAGCGTCTGTGTCTTTCAGGATTAATGCGCCGACAACAAAGTCAACGATGCCGCCGAAGCCGACCACGCCGAAACTGGGTAAGAAGTCTTCGGCGATCATGAAGGTGATGcc                                                                                                                            <  1:289898/150‑1 (MQ=255)                            aaGGTGCCGATCAAGGCCAGGGGGATGCCGAAACCGGGAGCGTCGGTGTCGATCAGGATCAATGCGCCGACCACAAAGGCAACGATGCCGCCGAAGCCGACCACGCCGAAACTGGGCAAGAAGGCTTCGGCGATCATGAAGGTGATGCCa                                                                                                                           >  2:496369/1‑150 (MQ=255)                              ggTGCCGATCAAGGCCAGGGGGATGCCGAAACCGGGAGCGTCGGTGTCGATCAGGATCAATGCGCCGACCACAAAGGCAACGATGCCGCCGAAGCCGACCACGCCGAAACTGGGCAAGAAGGCTTCGGCGATCATGAAGGTGATGCCAAg                                                                                                                         <  1:723654/150‑1 (MQ=255)                              ggTGCCGATCAAGGCCAGGGGGATGCCGAAACCGGGAGCGTCGGTGTCGATCAGGATCAATGCGCCGACCACAAAGGCAACGATGCCGCCGAAGCCGACCACGCCGAAACTGGGCAAGAAGGCTTCGGCGATCATGAAGGTGATGCCAAg                                                                                                                         <  1:496369/150‑1 (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                                                                                                           >  2:355966/1‑132 (MQ=255)                                                              cGGGAGCGTCGGTGTCGATCAGGATCAATGCGCCGACCACAAAGGCAACGATGCCGCCGAAGCCGACCACGCCGAAACTGGGCAAGAAGGCTTCGGCGATCATGAAGGTGATGCCAAGCAGGATCAGCGCCa                                                                                                           <  1:355966/132‑1 (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                                                      >  1:758795/1‑147 (MQ=255)                                                                                                    acaAAGGCAACGATGCCGCCGAAGCCGACCACGCCGAAACTGGGCAAGAAGGCTTCGGCGATCATGAAGGTGATGCCAAGCAGGATCAGCGCCACCCCGGCAAAACTCACCGGCAACAGCTGCAGCGCATACAGCGCCAGCAgcagg                                                      <  2:758795/147‑1 (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)                                                                                                                      ccgAAGCCGACCACGCCGAAACTGGGCAAGAAGGCTTCGGCGATCATGAAGGTGATGCCAAGCAGGATCAGCGCCACCCCGGCAAAACTCACCGGCAACAGCTGCAGCGCATACAGCGCCAGCAGCAGGCAGATACCACCGACCACACCg                                 >  1:355127/1‑150 (MQ=255)                                                                                                                       cgAAGCCGACCACGCCGAAACTGGGCAAGAAggcttcggcggcatcgttgcctttgtggtcggcgcattgatcctgatcgacaccgacgctcccggtttcggcatccccctggccttgatcggcaccttggcgctgctctcggcactcct                                <  1:676336‑M2/150‑120 (MQ=255)                                                                                                                           gCCGACCACGCCGAAACTGGGCAAGAAGGCTTCGGCGATCATGAAGGTGATGCCAAGCAGGATCAGCGCCACCCCGGCAAAACTCACCGGCAACAGCTGCAGCGCATACAGCGCCAGCAGCAGGCAGATACCACCGACCACACCGCCCAc                            <  1:650036/150‑1 (MQ=255)                                                                                                                                      cGAAACTGGGCAAGAAGGCTTCGGCGATCATGAAGGTGATGCCAAGCAGGATCAGCGCCACCCCGGCAAAACTCACCGGCAACAGCTGCAGCGCATACAGCGCCAGCAGCAGGCAGATACCACCGACCACACCGCCCACCGCAGAACCTg                 <  1:862470/150‑1 (MQ=255)                                                                                                                                        aaaCTGGGCAAGAAggcttcggcggcatcgttgcctttgtggtcggcgcattgatcctgatcgacaccgacgctcccggtttcggcatccccctggccttgatcggcaccttggcgctgctctcggcactcctgatcggtggcgtgctgg               >  1:646404‑M2/1‑14 (MQ=255)                                                                                                                                              ggCAAGAAGGCTTCGGCGATCATGAAGGTGATGCCAAGCAGGATCAGCGCCACCCCGGCAAAACTCACCGGCAACAGCTGCAGCGCATACAGCGCCAGCAGCAGGCAGATACCACCGACCACACCGCCCACCGCAGAACCTGGGTTCATg         <  1:1089934/150‑1 (MQ=255)                                                                                                                                                  agaagGCTTCGGCGATCATGAAGGTGATGCCAAGCAGGATCAGCGCCACCCCGGCAAAACTCACCGGCAACAGCTGCAGCGCATACAGCGCCAGCAGCAGGCAGATACCACCGACCACACCGCCCACCGCAGAACCTGGGTTCATGAAtt     <  2:849751/150‑1 (MQ=255)                                                                                                                                                   gaagGCTTCGGCGAGTATGAAGGTGATGCCAAGCAGGATCAGCGCCACCCCGGCAAAACTCACCGGGTACTGCTGCAGCGCATACAGCGCCAGCAGCAGGCAGATACCACCGACCACACCGCCCACCGCAGAACCTGGGTTCATGAATTc    <  1:380151/150‑1 (MQ=255)                                                                                                                                                     aggcttcggcggcatcgttgcctttgtggtcggcgcattgatcctgatcgacaccgacgctcccggtttcggcatccccctggccttgatcggcaccttggcgctgctctcggcactcctgatcggtggcgtgctgggcatggctctcaa  >  2:676336‑M2/1‑1 (MQ=255)                                                                                                                                                     agGCTTCGGCGATCATGAAGGTGATGCCAAGCAGGATCAGCGCCACCCCGGCAAAACTCACCGGCAACAGCTGCAGCGCATACAGCGCCAGCAGCAGGCAGATACCACCGACCACACCGCCCACCGCAGAACCTGGGTTCATGAATTCGa  <  2:1065193/150‑1 (MQ=255)                                                                                                                                                     |                                                                                                                                                      CCGATCAGGAGTGCCGAGAGCAGCGCCAAGGTGCCGATCAAGGCCAGGGGGATGCCGAAACCGGGAGCGTCGGTGTCGATCAGGATCAATGCGCCGACCACAAAGGCAACGATGCCGCCGAAGCCGACCACGCCGAAACTGGGCAAGAAGGCTTCGGCGATCATGAAGGTGATGCCAAGCAGGATCAGCGCCACCCCGGCAAAACTCACCGGCAACAGCTGCAGCGCATACAGCGCCAGCAGCAGGCAGATACCACCGACCACACCGCCCACCGCAGAACCTGGGTTCATGAATTCGA  >  NC_002947/250112‑250409

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.