New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_002947	88477 =	42 (0.990)	3 (0.080)	3/240	NT	6.1%	coding (841/1518 nt)	betC	choline‑sulfatase
	?	NC_002947	88524 =	53 (1.350)					coding (794/1518 nt)	betC	choline‑sulfatase

AGGCAGTGACCCGCCCCGCCGCGAAGCGCTTCGGCGCGTGGATCAGCAGCGGCACCCGCGCCGACATTTCGAACCAGTGCATTTTGTACCAGAGGCCACGCTCGCCAAGCATGTCGCCGTGGTCGCCGGAAAACACGATCAGTGTGTCATCGGCCAGGTTGCATTCCTCCAGGGTCTGCAGCAATTGGCCGATGTTGTCATCGATATAGCTGCACGCGCCGAAGTAGGCGCGGCGGGCATCGCGAACCTTATCCACAGGCAGTGGCTTGTTCCAAAGGTCGTAGACCTTCAGCAG  >  NC_002947/88330‑88624                                                                                                                                                    |                                                                                                                                                    aGGCAGTGACCCGCCCCGCCGCGAAGCGCTTCGGCGCGTGGATCAGCAGCGGCACCCGCGCCGACATTTCGAACCAGTGCATTTTGTACCAGAGGCCACGCTCGCCAAGCATGTCGCCGTGGTCGCCGGAAAACACGATCAGTGTGTc                                                                                                                                                     >  2:249156/1‑148 (MQ=255) aGGCAGTGACCCGCCCCGCCGCGAAGCGCTTCGGCGCGTGGATCAGCAGCGGCACCCGCGCCGACATTTCGAACCAGTGCATTTTGTACCAGAGGCCACGCTCGCCAAGCATGTCGCCGTGGTCGCCGGAAAACACGATCAGTGTGTc                                                                                                                                                     >  2:63251/1‑148 (MQ=255)    tgccgagttcagccaggcagaactcgacccgcattcacagcgcctgctgaaggtctacgacctttggaacaagccactgcctgtggataaggttcgcgatgcccgccgcgcctacttcggcgcgtgcagctatatcgatgacaaCATCg                                                                                                                                                 >  1:147935‑M1/145‑149 (MQ=255)    tgccgagttcagccaggcagaactcgacccgcattcacagcgcctgctgaaggtctacgacctttggaacaagccactgcctgtggataaggttcgcgatgcccgccgcgcctacttcggcgcgtgcagctatatcgatgacaaCATCg                                                                                                                                                 >  5:99085‑M1/145‑149 (MQ=255)    ggtcgagttctgcctggctgaactcgacccgcattcacagcgcctgctgaaggtctacgacctttggaacaagccactgcctgtggataaggttcgcgatgcccgccgcgcctacttcggcgcgtgcagctatatcgatgacaaCATCg                                                                                                                                                 <  1:173989‑M1/5‑1 (MQ=255)    cAGTGACCCGCCCCGCCGCGAAGCGCTTCGGCGCGTGGATCAGCAGCGGCACCCGCGCCGACATTTCGAACCAGTGCATTTTGTACCAGAGGCCACGCTCGCCAAGCATGTCGCCGTGGTCGCCGGAAAACACGATCAGTGTGTCATCg                                                                                                                                                 <  2:33674/149‑1 (MQ=255)        gagttcagccaggcagaactcgacccgcattcacagcgcctgctgaaggtctacgacctttggaacaagccactgcctgtggataaggttcgcgatgcccgccgcgcctacttcggcgcgtgcagctatatcgatgacaaCATCGGcc                                                                                                                                              <  8:174667‑M1/8‑1 (MQ=255)          ttcagccaggcagaactcgacccgcattcacagcgcctgctgaaggtcgaacgacctttggaacaagccactgcctgtggataaggttcgcgatgcccgccgcgcctacttcggcgcgtgcagctatatcgatgacaaCAtcggcaaat                                                                                                                                           <  1:63251‑M1/11‑5 (MQ=255)                ccgccgCGAAGCGCTTCGGTGCGTGGATCAGCAGCGGCACCCGCGCCGACATTTCGAACCAGTGCATTTTGTACCAGAGGCCACGCTCGCCAAGCATGTCGCCGTGGTCGCCGGAAAACACGATCAGTGTGTCATCGGCCAGGTTGCAt                                                                                                                                     >  4:107831/1‑149 (MQ=255)                  gccgcGAAGCGCTTCGGCGCGTGGATCAGCAGCGGCACCCGCGCCGACATTTCGAACCAGTGCATTTTGTACCAGAGGCCACGCTCGCCAAGCATGTCGCCGTGGTCGCCGGAAAACACGATCAGTGTGTCATCGGCCAGGTTGCATtc                                                                                                                                   <  1:148680/149‑1 (MQ=255)                            gcTTCGGCGCGTGGAACAGCAGCGGCACCCGCGCCGACATTTCGAACCAGTGCATTTTGTACCAGAGGCCACGCTCGCCAAGCATGTCGCCGTGGTCGCCGGAAAACACGATCAGTGCGTCATCGGCCAGGTTGCATTCCTCCAgg                                                                                                                            >  8:84015/1‑146 (MQ=255)                            gcTTCGGCGCGTGGAACAGCAGCGGCACCCGCGCCGACATTTCGAACCAGTGCATTTTGTACCAGAGGCCACGCTCGCCAAGCATGTCGCCGTGGTCGCCGGAAAACACGATCAGTGCGTCATCGGCCAGGTTGCATTCCTCCAgg                                                                                                                            <  7:84015/146‑1 (MQ=255)                                     cgTGGATCAGCAGCGGCACCCGCGCCGACATTTCGAACCAGTGCATTTTGTACCAGAGGCCACGCTCGCCAAGCATGTCGCCGTGGTCGCCGGAAAACACGATCAGTGTGTCATCGGCCAGGTTGCATTCCTCCAg                                                                                                                             >  4:247542/1‑136 (MQ=255)                                     cgTGGATCAGCAGCGGCACCCGCGCCGACATTTCGAACCAGTGCATTTTGTACCAGAGGCCACGCTCGCCAAGCATGTCGCCGTGGTCGCCGGAAAACACGATCAGTGTGTCATCGGCCAGGTTGCATTCCTCCAg                                                                                                                             <  3:247542/136‑1 (MQ=255)                                       tGGATCAGCAGCGGCACCCGCTCCGACATTTCGAACTAGTGCATTTTGTACCAGAGGCCACGCTCGCCAAGCATGTCGCGGTGGTCGGCGTAAAACACGATCAGTGTGTCATCGGCCAGGTTGCATTCCTCCAGGGGCTGCAGCAATTg                                                                                                              <  7:159273/149‑1 (MQ=255)                                         gATCAGCAGCGGCACCCGCGCCGACATTTCGAACCAGTGCATTTTGTACCAGAGGCCACGCTCGCCAAGCATGTCGCCGTGGTCGCCGGAAAACACGATCAGTGTGTCATCGGCCAGGTTGCATTCCTCCAGGGTCTGCAGCAATTGGc                                                                                                            >  8:15126/1‑149 (MQ=255)                                          aTCAGCAGCGGCACCCGCGCCGACATTTCGAACCAGTGCATTTTGTACCAGAGGCCACGCTCTCCAAGCATGTCGCCGTGGTCGCCGGAAAACACGATCAGTGTGTCATCTGCCAGGTTGCATTCCTCCAGGGTCTgcagc                                                                                                                   <  7:8714/141‑1 (MQ=255)                                          aTCAGCAGCGGCACCCGCGCCGACATTTCGAACCAGTGCATTTTGTACCAGAGGCCACGCTCGCCAAGCATGTCGCCGTGGTCGCCGGAAAACACGATCAGTGTGTCATCGGCCAGGTTGCATTCCTCCAGGGTCTgcagc                                                                                                                   >  8:8714/1‑141 (MQ=255)                                                     cACCCGCGCCGACATTTCGAACCAGTGCATTTTGTACCAGAGGCCACGCTCGCCAAGCATGTCGCCGTGGTCGCCGGAAAACACGATCAGTGTGTCATCGGCCAGGTTGCATTCCTCCAGGGTCTGCAGCAATTGGCCGATGTTGTCAt                                                                                                <  6:37296/149‑1 (MQ=255)                                                     aaggtctacgacctttggaacaagccactgcctgtggataaggttcgcgatgcccgccgcgcctacttcggcgcgtgcagctatatcgatgacaaCATCGGCCAGGTTGCATTCCTCCAGGGTCTGCAGCAATTGGCCGATGTTGTCa                                                                                                 >  6:113531‑M1/96‑148 (MQ=255)                                                          gcgcCGACATTTCGAACCAGTGCATTTTGTACCAGAGGCCACGCTCGCCAAGCATGTCGCCGTGGTCGCCGGAAAACACGATCAGTGTGTCATCGGCCAGGTTGCATTCCTCCAGGGTCTGCAGCAATTGGCCGATGTTGTCATCGata                                                                                           <  7:15126/149‑1 (MQ=255)                                                           cgactaCATTTCGAACCAGTGCATTTTGTACCAGATGCCACGCTCGCCAATCATGTCGCCGTTTTCGCCGGAAAACACGATCATTGTGTCATCGTCCAGGTTTCATTCCTCCAGGGTATGCAGCAATTGGCCTATGTTGTCATCGatat                                                                                          <  4:142582/144‑1 (MQ=255)                                                             ccGACATTTCGAACCAGTGCATTTTGTACCAGAGGCCACGCTCGCCAAGCATGTCGCCGTGGTCGCCGGAAAACACGATCAGTGTGTCATCGGCCAGGTTGCATTCCTCCAGGGTCTGCAGCAATTGGCCGATGTTGTCATCGATATAg                                                                                        >  4:39634/1‑149 (MQ=255)                                                              cGACATTTCGAACCAGTGCATTTTGTACCAGAGGCCACGCTCGCCAAGCATGTCGCCGTGGTCGCCGGAAAACACGATCAGTGTGTCATCGGCCAGGTTGCATTCCTCCAGGGTCTGCAGCAATTGGCCGATGTTGTCATCGATATAGc                                                                                       >  8:35352/1‑149 (MQ=255)                                                              cGACATTTCGAACCAGTGCATTTTGTACCAGAGGCCACGCTCGCCAAGCATGTCGCCGTGGTCGCCGGAAAACACGATCAGTGTGTCATCGGCCAGGTTGCATTCCTCCAGGGTCTGCAGCAATTGGCCGATGTTGTCATCGATATAGc                                                                                       >  4:132597/1‑149 (MQ=255)                                                               gACATTTCGAACCAGTGCATTTTGTACCAGAGGCCACGCTCGCCAAGCATGTCGCCGTGGTCGCCGGAAAACACGATCAGTGTGTCATCGGCCAGGTTGCATTCCTCCAGGGTCTGCAGCAATTGGCCGATGTTGTCATCGATATAGCt                                                                                      >  7:75175/1‑149 (MQ=255)                                                                aCATTTCGAACCAGTGCATTTTGTACCAGAGGCCACGCTCGCCAAGCATGTCGCCGTGGTCGCCGGAAAACACGATCAGTGTGTCATCGGCCAGGTTGCATTCCTCCAGGGTCTGCAGCAATTGGCCGa                                                                                                         <  8:29650/129‑1 (MQ=255)                                                                aCATTTCGAACCAGTGCATTTTGTACCAGAGGCCACGCTCGCCAAGCATGTCGCCGTGGTCGCCGGAAAACACGATCAGTGTGTCATCGGCCAGGTTGCATTCCTCCAGGGTCTGCAGCAATTGGCCGa                                                                                                         >  7:29650/1‑129 (MQ=255)                                                                    ttCGAACCAGTGCATTTTGTACCAGAGGCCACGCTCGCCAAGCATGTCGCCGTGGTCGCCGGAAAACACGATCAGTGTGTCATCGGCCAGGTTGCATTCCTCCAGGGTCTGCAGCAATTGGCCGATGTTGTCATCGATATAGCTGCAcg                                                                                 >  3:13460/1‑149 (MQ=255)                                                                           cAGTGCATTTTGTACCAGAGGCCACGCTCGCCAAGCATGTCGCCGTGGTCGCCGGAAAACACGATCAGTGTGTCATCGGCCAGGTTGCATTCCTCCAGGGTCTGCAGCAATTGGCCGATGTTGTCATCGATATAGCTGCACGCGCCGaa                                                                          >  4:265379/1‑149 (MQ=255)                                                                              tGCATTTTGTACCAGAGGCCACGCTCGCCAAGCATGTCGCCGTGGTCGCCGGAAAACACGATCAGTGTGTCATCGGCCAGGTTGCATTCCTCCAGGGTCTGCAGCAATTGGCCGATGTTGTCATCGATATAGCTGCACGCGCCGAAGTa                                                                       <  5:35755/149‑1 (MQ=255)                                                                                    ttGTACCAGAGGCCACGCTCGCCAAGCATGTCGCCGTGGTCGCCGGAAAACACGATCAGTGTGTCATCGGCCAGGTTGCATTCCTCCAGGGTCTGCAGCAATTGGCCGATGTTGTCATCGATATAGCTGCACGCGCCGAAGTAGgcgc                                                                  <  6:268746/148‑1 (MQ=255)                                                                                         ccAGAGGCCACGCTCGCCAAGCATGTCGCCGTGGTCGCCGGAAAACACGATCAGTGTGTCATCGGCCAGGTTGCATTCCTCCAGGGTCTGCAGCAATTGGCCGATGTTGTCATCGATATAGCTGCACGCGCCGAAGTAGGCGTGGCggg                                                            >  4:127393/1‑149 (MQ=255)                                                                                          cAGAGGCCACGCTCGCCAAGCATGTCGCCGTGGTCGCCGGAAAACACGATCAGTGTGTCATCGGCCAGGTTGCATTCCTCCAGGGTCTGCAGCAATTGGCCGATGTTGTCATCGATATAGCTGCACGCGCCGAAGTAGGCGCGGCggg                                                            >  2:58050/1‑148 (MQ=255)                                                                                            gagGCCACGCTCGCCAAGCATGTCGCCGTGGTCGCCGGAAAACACGATCAGTGTGTCATCGGCCAGGTTGCATTCCTCCAGGGTCTGCAGCAATTGGCCGATGTTGTCATCGATATAGCTGCACGCGCCGAAGTAGGCGCGGCGGGCAt                                                         >  5:124353/1‑149 (MQ=255)                                                                                              ggCCACGCTCGCCAAGCATGTCGCCGTGGTCGCCGGAAAACACGATCAGTGTGTCATCGGCCAGGTTGCATTCCTCCAGGGTCTGCAGCAATTGGCCGATGTTGTCATCGATATAGCTGCACGCGCCGAAGTAGGCGCGGCGGGCATcg                                                       <  4:13460/149‑1 (MQ=255)                                                                                               gCCACGCTCGCCAAGCATGTCGCCGTGGTCGCCGGAAAACACGATCAGTGTGTCATCGGCCAGGTTGCATTCCTCCAGGGTCTGCAGCAATTGGCCGATGTTGTCATCGATATAGCTGCACGCGCCGAAGTAGGCGCGGCGGGCATcgc                                                      >  2:173989/1‑149 (MQ=255)                                                                                                 tcgcgatgcccgccgcgcctacttcggcgcgtgcagctatatcgatgacaaCATCGGCCAGGTTGCATTCCTCCAGGGTCTGCAGCAATTGGCCg                                                                                                          <  1:175912‑M1/44‑1 (MQ=255)                                                                                                 tcgcgatgcccgccgcgcctacttcggcgcgtgcagctatatcgatgacaaCATCGGCCAGGTTGCATTCCTCCAGGGTCTGCAGCAATTGGCCg                                                                                                          >  2:175912‑M1/52‑95 (MQ=255)                                                                                                   cGCTCGCCAAGCATGTCGCCGTGGTCGCCGGAAAACACGATCAGTGTGTCATCGGCCAGGTTGCATTCCTCCAGGGTCTGCAGCAATTGGCCGATGTTGTCATCGATATAGCTGCACGCGCCGAAGTAGGCGCGGCGGGCATCGCGAAc                                                  <  3:132597/149‑1 (MQ=255)                                                                                                      tCGCCAAGCATGTCGCCGTGGTCGCCGGAAAACACGATCAGTGTGTCATCGGCCAGGTTGCATTCCTCCAGGGTCTGCAGCAATTGGCCGATGTTGTCATCGATATAGCTGCACGCGCCGAAGTAGGCGCGGCGGGCATCGCGAACCt                                                >  3:16728/1‑148 (MQ=255)                                                                                                      tCGCCAAGCATGTCGCCGTGGTCGCCGGAAAACACGATAAGTGTGTCATCGGCCAGGTTGCATTCCTCCAGGGTCAGCAGCAATTGGCCGATGTTGTCATCGATATAGCTGCACGCGCCGAAGAAGGCGCGGCGGGCATCGCGAACCt                                                >  6:156131/1‑148 (MQ=255)                                                                                                        gCCAAGCATGTCGCCGTGGTCGCCGGAAAACACGATCAGTGTGTCATCGGCCAGGTTGCATTCCTCCAGGGTCTGCAGCAATTGGCCGATGTTGTCATCGATATAGCTGCACGCGCCGAAGTAGGCGCGGCGGGCATCGCGAACCTTAt                                             >  1:145858/1‑149 (MQ=255)                                                                                                         ccAAGCATGTCGCCGTGGTCGCCGGAAAACACGATCAGTGTGTCATCGGCCAGGTTGCATTCCTCCAGGGTCTGCAGCAATTGGCCGATGTTGTCATCGATATAGCTGCACGCGCCGAAGTAGGCGCGGCGGGCATCGCGAACCTTATc                                            <  3:36925/149‑1 (MQ=255)                                                                                                          cAAGCATGTCGCCGTGGTCGCCGGAAAACACGATCAGTGTGTCATCGGCCAGGTTGCATTCCTCCAgg                                                                                                                            >  7:240049/1‑68 (MQ=255)                                                                                                          cAAGCATGTCGCCGTGGTCGCCGGAAAACACGATCAGTGTGTCATCGGCCAGGTTGCATTCCTCCAgg                                                                                                                            <  8:240049/68‑1 (MQ=255)                                                                                                                    gCCGTGGTCGCCGGAAAACACGATCAGTGTGTCATCGGCCAGGTTGCATTCCTCCAGGGTCTGCAGCAATTGGCCGATGTTGTCATCGATATAGCTGCACGCGCCGAAGt                                                                        <  3:86636/110‑1 (MQ=255)                                                                                                                    gCCGTGGTCGCCGGAAAACACGATCAGTGTGTCATCGGCCAGGTTGCATTCCTCCAGGGTCTGCAGCAATTGGCCGATGTTGTCATCGATATAGCTGCACGCGCCGAAGt                                                                        >  4:86636/1‑110 (MQ=255)                                                                                                                        tGGTCGCCGGAAAACACGATCAGTGTGTCATCGGCCAGGTTGCATTCCTCCAGGGTCTGCAGCAATTGGCCGATGTTGTCATCGATATAGCTGCACGCGCCGAAGTAGGCGCGGCGGGCATCGCGAACCTTATCCACAGGCAGTGGCtt                             <  7:81720/149‑1 (MQ=255)                                                                                                                          gTCGCCGGAAAACACGATCAGTTTGTCATCGGCCAGGTTGCATTCCTCCAGGGTCTGAAGCAATTGGCCGATGTTGTCATCGATATAGCTGCACGCGCCGAAGTAGGCGCGGCGGGCATCGCGAACCTTATCCACAGGCAGTGGCTTGt                           >  3:243299/1‑149 (MQ=255)                                                                                                                            cGCCGGAAAACACGATCAGTGTGTCATCGGCCAGGTTGCATTCCTCCAGGGTCTGCAGCAATTGGCCGATGTTGTCATCGATATAGCTGCACGCGCCGAAGTAGGCGCGGCGGGCATCGCGAACCTTATCCACAGGCAGTGGCTTGTTc                         <  4:16728/149‑1 (MQ=255)                                                                                                                                  aaaaCACGATCAGTGTGTCATCGGCCAGGTTGCATTCCTCCAGGGTCTGCAGCAATTGGCCGATGTTGTCATCGATATAGCTGCACGCGCCGAAGTAGGCGCGGCGGGCATCGCGAACCTTATCCACAGGCAGTGGCTTGTTCCAAAgg                   >  4:78782/1‑149 (MQ=255)                                                                                                                                         gATCAGTGTGTCATCGGCCAGGTTGCATTCCTCCAGGGTCTGCAGCAATTGGCCGATGTTGTCATCGATATAGCTGCACGCGCCGAAGTAGGCGCGGCGGGCATCGCGAACCTTATCCACAGGCAGTGGCTTGTTCCAAAGGTCGTAGa            >  2:65791/1‑149 (MQ=255)                                                                                                                                                 tgtCATCGGCCAGGTTGCATTCCTCCAGGGTCTGCAGCAATTGGCCGATGTTGTCATCGATATAGCTGCACGCGCCGAAGTAGGCGCGGCGGGCATCGCGAACCTTATCCACAGGCAGTGGCTTGTTCCAAAGGTCGTAGACCTTcag     >  7:93936/1‑148 (MQ=255)                                                                                                                                                   tCATCGGCCAGGTTGCATTCCTCCAGGGTCTGCAGCAATTGGCCGATGTTGTCATCGATATAGCTGCACGCGCCGAAGTAGGCGCGGCGGGCATCGCGAACCTTATCCACAGGCAGTGGCTTGTTCCAAAGGTCGTAGACCTTcagcag  <  3:107831/149‑1 (MQ=255)                                                                                                                                                    |                                                                                                                                                    AGGCAGTGACCCGCCCCGCCGCGAAGCGCTTCGGCGCGTGGATCAGCAGCGGCACCCGCGCCGACATTTCGAACCAGTGCATTTTGTACCAGAGGCCACGCTCGCCAAGCATGTCGCCGTGGTCGCCGGAAAACACGATCAGTGTGTCATCGGCCAGGTTGCATTCCTCCAGGGTCTGCAGCAATTGGCCGATGTTGTCATCGATATAGCTGCACGCGCCGAAGTAGGCGCGGCGGGCATCGCGAACCTTATCCACAGGCAGTGGCTTGTTCCAAAGGTCGTAGACCTTCAGCAG  >  NC_002947/88330‑88624

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