New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_002947	3193191 =	32 (0.560)	3 (0.060)	3/250	NT	7.9%	coding (270/1425 nt)	prr	gamma‑aminobutyraldehyde dehydrogenase
	?	NC_002947	3193230 =	39 (0.720)					coding (231/1425 nt)	prr	gamma‑aminobutyraldehyde dehydrogenase

GTACTCGCCCGCCGCCGAGCCTTGCAGGCAGCGAACGGCGCCGGCGAAGAAGCGGAACACGTCCGACACGGCGGGCAGTTCATCCTGCAGCACGGCCGTGAACGGCTTGCCGCAGTTGTCCGCCTCCAGGCGGGCGAACACCGCGGCCTGGGCGTCGATGCGGTCGGCGAGGGTGAGCAGCAGTGCCGCACGGTCCTTGGGCGGTGTTTGCGACCAGGCTTCGAAGGCCGCATCCGCGCTTTGCACTGCGGCATGGACCTGTGCATGACTGGCCTCGTGGGTCTGTACC  >  NC_002947/3193083‑3193371                                                                                                                                                    |                                                                                                                                              catgccgcagtgcaaagcgcggatgcggccttcgaagcctggtcgcaaacaccgcccaaggaccgtgcggcactgctgctcaccctcgccgaccgcatcgacgcccagGCCGCAGTTGTCCGCCTCCAGGCGGGCGAACACCGCGGCCt                                                                                                                                              >  6:69987‑M1/109‑149 (MQ=255)  tACTCGCCCGCCGCCGAGCCTTGCAGGCAGCGAACGGCGCCGGCGAAGAAGCGGAACACGTCCGACACGGCGGGCAGTTCATCCTGCAGCACGGCCGTGAACGGCTTGCCGCAGTTGTCCGCCTCCAGGCGGGCGAACACCGCGGCCTg                                                                                                                                             <  5:92723/149‑1 (MQ=255)       gcccgccGCCGAGCCTTGCAGGCAGCGAACGGCGCCGGCGAAGAAGCGGAACACGTCCGACACGGCGGGCAGTTCATCCTGCAGCACGGCCGTGAACGGCTTGCCGCAGTTGTCCGCCTCCAGGCGGGCGAACACCGCGGCCTGGGc                                                                                                                                          <  6:175432/147‑1 (MQ=255)        cccgccGCCGAGCCTTGCAGGCAGCGAACGGCGCCGGCGAAGAAGCGGAACACGTCCGACACGGCGGGCAGTTCATCCTGCAGCACGGCCGTGAACGGCTTGCCGCAGTTGTCCGCCTCCAGGCGGGCGAACACCGCGGCCTGGGCGTc                                                                                                                                       <  4:211230/149‑1 (MQ=255)              gccgAGCCTTGCAGGCAGCGAACGGCGCCGGCGAAGAAGCGGAACACGTCCGACACGGCGGGCAGTTCATCCTGCAGCACGGCCGTGAACGGCTTGCCGCAGTTGTCCGCCTCCAGGCGGGCGAACACCGCGGCCTGGGCGTCGATGCg                                                                                                                                 <  6:138297/149‑1 (MQ=255)                       tGCAGGCAGCGAACGGCGCCGGCGAAGAAGCGGAACACGTCCGACACGGCGGGCAGTTCATCCTGCAGCACGGCCGTGAACGGCTTGCCGCAGTTGTCCGCCTCCAGGCGGGCGAACACCGCGGCCTGGGCGTCGATGCGGTCGGCGAg                                                                                                                        >  2:78644/1‑149 (MQ=255)                        gcaggcagCGAACGGCGCCGGCGAAGAAGCGGAACACGTCCGACACGGCGGGCAGTTCATCCTGCAGCACGGCCGTGAACGGCTTGCCGCAGTTGTCCGCCTCCAGGCGGGCGAACACCGCGGCCTGGGCGTCGATGCGGTCGGCGAg                                                                                                                        <  2:334183/148‑1 (MQ=255)                              agCGAACGGCGCCGGCGAAGAAGCGGAACACGTCCGACACGGCGGGCAGTTCATCCTGCAGCACGGCCGTGAACGGCTTGCCGCAGTTGTCCGCCTCCAGGCGGGCGAACACCGCGGCCTGGGCGTCGATGCGGTCGGCGAGGGTGag                                                                                                                  >  3:351368/1‑148 (MQ=255)                                  aaCGGCGCCGGCGAAGAAGCGGAACACGTCCGACACGGCGGGCAGTTCATCCTGCAGCACGGCCGTGAACGGCTTGCCGCAGTTGTCCGCCTCCAGGCGGGCGAACACCGCGGCCTGGGCGTCGATGCGGTCGGCGAGGGTGAgcagca                                                                                                             <  1:78644/149‑1 (MQ=255)                                  aaCGGCGCCGGCGAAGAAGCGGAACACGTCCGACACGGCGGGCAGTTCATCCTGCAGCACGGCCGTGAACGGCTTGCCGCAGTTGTCCGCCTCCAGGCGGGCGAACACCGCGGCCTGGGCGTCGATGCGGTCGGCGAGGGTGAgcagca                                                                                                             <  8:158992/149‑1 (MQ=255)                                                 gaagCGGAACACGTCCGACACGGCGGGCAGTTCATCCTGCAGCACGGCCGTGAACGGCTTGCCGCAGTTGTCCGCCTCCAGGCGGGCGAACACCGCGGCCTGGGCGTCGATGCGGTCGGCGAGGGTGAGCAGCAGTGCCGCACGGTCCt                                                                                              >  8:313031/1‑149 (MQ=255)                                                    gCGGAACACGTCCGACACGGCGGGCAGTTCATCCTGCAGCACGGCCGTGAACGGCTTGCCGCAGTTGTCCGCCTCCAGGCGGGCGAACACCGCGGCCTGGGCGTCGATGCGGTCGGCGAGGGTGAGCAGCAGTGCCGCACGGTCCTTgg                                                                                           <  4:351368/149‑1 (MQ=255)                                                         acacGTCCGACACGGCGGGCAGTTCATCCTGCAGCACGGCCGTGAACGGCTTGCCGCAGTTGTCCGCCTCCAGGCGGGCGAACACCGCGGCCTGGGCGTCGATGCGGTCGGCGAGGGTGAGCAGCAGTGCCGCACGGTCCTTGGGCGgt                                                                                      >  6:163014/1‑149 (MQ=255)                                                          cacGTCCGACACGGCGGGCAGTTCATCCTGCAGCACGGCCGTGAACGGCTTGCCGCAGTTGTCCGCCTCCAGGCGGGCGAACACCGCGGCCTGGGCGTCGATGCGGTCGGCGAGGGTGAGCAGCAGTGCCGCACGGTCCTTGGGCGgtg                                                                                     >  1:228628/1‑149 (MQ=255)                                                             gTCCGACACGGCGGGCAGTTCATCCTGCAGCACGGCCGTGAACGGCTTGCCGCAGTTGTCCGCCTCCAGGCGGGCGAACACCGCGGCCTGGGCGTCGATGCGGTCGGCGAGGGTGAGCAGCAGTGCCGCACGGTCCTTGGGCGGTGttt                                                                                  <  5:163014/149‑1 (MQ=255)                                                                  acacGGCGGGCAGTTCATCCTGCAGCACGGCCGTGAACGGCTTGCCGCAGTTGTCCGCCTCCAGGCGGGCGAACACCGCGGCCTGGGCGTCGATGCGGTCGGCGAGGGTGAGCAGCAGTGCCGCACGGTCCTTGGGCGGTGTTTGCGAc                                                                             >  2:321679/1‑149 (MQ=255)                                                                   cacGGCGGGCAGTTCATACTCCAGCACGGCCGTGAACGGCTTGCCGCAGTTGTCCGCCTACAGGCTGGCGAACACCGCGGCCTGGGCGTCGATGGGGTCGGAGACGGTGAGCAGCAATGCCGCACGGTCCTTGGTCGGGGTTTGCGAc                                                                             >  8:316475/1‑148 (MQ=255)                                                                          ggCAGTTCATCCTGCAGCACGGCCGTGAACGGCTTGCCGCAGTTGTCCGCCTCCAGGCGGGCGAACACCGCGGCCTGGGCGTCGATGCGGTCGGCGAGGGTGAGCAGCAGTGCCGCACGGTCCTTGGGCGGTGTTTGCGACCAGGCTTc                                                                     <  3:1323/149‑1 (MQ=255)                                                                           gCAGTTCATCCTGCAGCACGGCCGTGAACGGCTTGCCGCAGTTGTCCGCCTCCAGGCGGGCGAACACCGCGGCCTGGGCGTCGATGCGGTCGGCGAGGGTGAGCAGCAGTGCCGCACGGTCCTTGGGCGGTGTTTGCGACCAGGCTTCg                                                                    <  7:313031/149‑1 (MQ=255)                                                                              gTTCATCCTGCAGCACGGCCGTGAACGGCTTGCCGCAGTTGTCCGCATCCAGGCGGGCGAACACCGCGGCCTGGGCGTCGATGCGGTCGGCGAGGGTGAGCAGCAGTGCCGCACGGTCCTTGGGCGGTGTTTGCGACCAGGCTTCGAAg                                                                 >  1:79605/1‑149 (MQ=255)                                                                                tCATCCTGCAGCACGGCCGTGAACGGCTTGCCGCAGTTGTCCGCCTCCAGGCGGGCGAACACCGCGGCCTGGGCGTCGATGCGGTCGGCGAGGGTGAGCAGCAGTGCCGCACGGTCCTTGGGCGGTGTTTGCGACCAGGCTTCGAAGGc                                                               >  8:22914/1‑149 (MQ=255)                                                                                 cATCCTGCAGCACGGCCGTGAACGGCTTGCCGCAGTTGTCCGCCTCCAGGCGGGCGAACACCGCGGCCTGGGCGTCGATGCGGTCGGCGAGGGTGAGCAGCAGTGCCGCACGGTCCTTGGGCGGTGTTTGCGACCAGGCTTCGAAGGc                                                               <  1:321679/148‑1 (MQ=255)                                                                                     cTGCAGCACGGCCGTGAACGGCTTGCCGCAGTTGTCCGCCTCCAGGCGGGCGAACACCGCGGCCTGGGCGTCGATGCGGTCGGCGAGGGTGAGCAGCAGTGCCGCACGGTCCTTGGGCGGTGTTTGCGACCAGGCTTCGAAGGCCGCAt                                                          >  5:174321/1‑149 (MQ=255)                                                                                        cagcaCGGCCGTGAACGGCTTGCCGCAGTTGTCCGCCTCCAGGCGGGCGAACACCGCGACCTGGGCGTCGATGCGGTCGGCGAGGGTGAGCAGCAGTGCCGCACGGTCCTTGGGCGGTGTTTGCGACCAGGCTTCGAAGGCCGCATCcg                                                       >  4:30844/1‑149 (MQ=255)                                                                                                     aaCGGCTTGCCGCAGTTGTCCGCCTCCAGGCGGGCGAACACCGCGGCCTGGGCGTCGATGCGGTCGGCGAGGGTGAGCAGCAGTGCCGCACGGTCCTTGGGCGGTGTTTGCGACCAGGCTTCGAAGGCCGCATCCGCGCTTTGCACTGc                                          <  2:228628/149‑1 (MQ=255)                                                                                                      aCGGCTTGCCGCAGTTGTCCGCCTCCAGGCGGGCGAACACCGCGGCCTGGGCGTCGATGCGGTCGGCGAGGGTGAGCAGCAGTGCCGCACGGTCCTTGGGCGGTGTTTGCGACCAGGCTTCGAAGGCCGCATCCGCGCTTTGCACTGCg                                         >  5:317159/1‑149 (MQ=255)                                                                                                         gCTTGCCGCAGTTGTCCGCCTCCAGGCGGGCGAACACCGCGGCCTGGGCGTCGATGCGGTCGGCGAGGGTGAGCAGCAGTGCCGCACGGTCCTTGGGCGGTGTTTGCGACCAGGCTTCGAAGGCCGCATCCGCGCTTTGCACTGCGGCa                                      <  6:174321/149‑1 (MQ=255)                                                                                                           ttGCCGCAGTTGTCCGCCTCCAGGCGGGCGAACACCGCGGCCTGGGCGTCGATGCGGTCGGCGAGGGTGAGCAGCAGTGCCGCACGGTCCTTGGGCGGTGTTTGCGACCAGGCTTCGAAGGCCGCATCCGCGCTTTGCACTGCGGCATg                                    >  2:96217/1‑149 (MQ=255)                                                                                                           aggccgcggtgttcgcccgcctggaggcggacaactgcggcCTGGGCGTCGATGCGGTCGGCGAGGGTGAGCAGCAGTGCCGCACGGTCCTTGGGCGGTGTTTGCGACCAGGCTTCGAAGGCCGCATCCGCGCTTTGCACTGCGGCATg                                    >  6:69987‑M2/42‑149 (MQ=255)                                                                                                              ccgcggtgttcgcccgcctggaggcggacaactgcggcCTGGGCGTCGATGCGGTCGGCGAGGGTGAGCAGCAGTGCCGCACGGTCCTTGGGCGGTGTTTGCGACCAGGCTTCGAAGGCCGCATCCGCGCTTTGCACTGCGGCATGGa                                  >  5:293913‑M2/39‑148 (MQ=255)                                                                                                              ccgcggtgttcgcccgcctggaggcggacaactgcggcCTGGGCGTCGATGCGGTCGGCGAGGGTGAGCAGCAGTGCCGCACGGTCCTTGGGCGGTGTTTGCGACCAGGCTTCGAAGGCCGCATCCGCGCTTTGCACTGCGGCATGGa                                  <  6:293913‑M2/110‑1 (MQ=255)                                                                                                               CGCAGTTGTCCGCCTGCAGGCGGGCGAACACCGCGGCCTGGGCGTCGATGGGGTCGGCGAGGGTGAGCAGGAGTGCCGCACGGTCCTTGGGCTGTGTTTGCGACCAGGCTTCGAAGGCCGCATCCGCGCTTTGCACTGCTGCATGGAc                                 >  7:316475‑M1/1‑148 (MQ=255)                                                                                                                GCAGTTGTCCGCCTCCAGGCGGGCGAACACCGCGGCCTGGGCGTCGATGCGGTCGGCGAGGGTGGGCAGCAGTGCCGCACGGTCCTTGGGCGGTGTTTGCGACCAGGCTTCGAAGGCCGCATCCGCGCTTTGCACTGCGGCATGGACCt                               <  3:198909‑M1/149‑1 (MQ=255)                                                                                                                GCAGTTGTCCGCCTCCAGGCGGGCGAACACCGCGGCCTGGGCGTCGATGCGGTCGGCGAGGGTGAGCAGCAGTGCCGCACGGTCCTTGGGCGGTGTTTGCGACCAGGCTTCGAAGGCCGCATCCGCGCTTTGCACTGCGGCATGGACCt                               <  1:44489‑M1/149‑1 (MQ=255)                                                                                                                          gCCTCCAGGCGGGCGAACACCGCGGCCTGGGCGTCGATGCGGTCGGCGAGGGTGAGCAGCAGTGCCGCACGGTCCTTGGGCGGTGTTTGCGACCAGGCTTCGAAGGCCGCATCCGCGCTTTGCACTGCGGCATGGACCTGTGCATGACt                     >  8:254955/1‑149 (MQ=255)                                                                                                                           ccTCCAGGCGGGCGAACACCGCGACCTGGGCGTCGATGCGGTCGGCGAGGGTGAGCAGCAGTGCCGCACGGTCCTTGGGCGGTGTTTGCGACCAGGCTTCGAAGGCCGCATCCGCGCTTTGCACTGCGGCATGGACCTGTGCATGACTg                    <  3:30844/149‑1 (MQ=255)                                                                                                                                  gcgggcgAACACCGCGGCCTGGGCGTCGATGCGGTCGGCGAGGGTGAGCAGCAGTGCCGCACGGTCCTTGGGCGGTGTTTGCGACCAGGCTTCGAAGGCCGCATCCGCGCTTTGCACTGCGGCATGGACCTGTGCATGACTGGCCTCGt             <  2:44489/149‑1 (MQ=255)                                                                                                                                  gcgggcgAACACCGCGGCCTGGGCGTCGATGCGGTCGGCGAGGGTGAGCAGCAGTGCCGCACGGTCCTTGGGCGGTGTTTGCGACCAGGCTTCGAAGGCCGCATCCGCGCTTTGCACTGCGGCATGGACCTGTGCATGACTGGCCTCGt             >  4:140854/1‑149 (MQ=255)                                                                                                                                  gcgggcgAACACCGCGGCCTGGGCGTCGATGCGGTCGGCGAGGGTGAGCAGCAGTGCCGCACGGTCCTTGGGCGGTGTTTGCGACCAGGCTTCGAAGGCCGCATCCGCGCTTTGCACTGCGGCATGGACCGGTGCATGACTGGCCTCGt             <  3:76424/149‑1 (MQ=255)                                                                                                                                    gggcgAACACCGCGGCCTGGGCGTCGATGCGGTCGGCGAGGGTGGGCAGCAGTGCCGCACGGTCCTTGGGCGGTGTTTGCGACCAGGCTTCGAAGGCCGCATCCGCGCTTTGCACTGCGGCATGGACCTGTGCATGACTGGCCTCGTgg           <  4:198909/149‑1 (MQ=255)                                                                                                                                    gggcgAACACCGCGGCCTGGGCGTCGATGCGGTCGGCGAGGGTGAGCAGCAGTGCCGCACGGTCCTTGGGCGGTGTTTGCGACCAGGCTTCGAAGGCCGCATCCGCGCTTTGCACTGCg                                         <  1:174988/119‑1 (MQ=255)                                                                                                                                    gggcgAACACCGCGGCCTGGGCGTCGATGCGGTCGGCGAGGGTGAGCAGCAGTGCCGCACGGTCCTTGGGCGGTGTTTGCGACCAGGCTTCGAAGGCCGCATCCGCGCTTTGCACTGCg                                         >  2:174988/1‑119 (MQ=255)                                                                                                                                    gggcgAACACCGCGGCCTGGGCGTCGATGCGGTCGGCGAGGGTGAGCAGCAGTGCCGCACGGTCCTTGGGCGGTGTTTGCGACCAGGCTTCGAAGGCCGCATCCGCGCTTTGCACTGCGGCATGGACCTGTGCATGACTGGCCTCGTgg           >  1:351759/1‑149 (MQ=255)                                                                                                                                    gggcgAACACCGCGGCCTGGGCGTCGATGCGGTCGGCGAGGGTGAGCAGCAGTGCCGCACGGTCCTTGGGCGGTGTTTGCGACCAGGCTTCGAAGGCCGCATCCGCGCTTTGCACTGCGGCATGGACCTGTGCATGACTGGCCTCGTgg           <  7:22914/149‑1 (MQ=255)                                                                                                                                       cgAACACCGCGGCCTGGGCGTCGATGCGGTCGGCGAGGGTGAGCAGCAGTGCCGCACGGTCCTTGGGCGGTGTTTGCGACCAGGCTTCGAAGGCCGCATCCGCGCTTTGCACTGCGGCATGGACCTGTGCATGACTGGCCTCGt             <  4:231017/144‑1 (MQ=255)                                                                                                                                       cgAACACCGCGGCCTGGGCGTCGATGCGGTCGGCGAGGGTGAGCAGCAGTGCCGCACGGTCCTTGGGCGGTGTTTGCGACCAGGCTTCGAAGGCCGCATCCGCGCTTTGCACTGCGGCATGGACCTGTGCATGACTGGCCTCGt             >  3:231017/1‑144 (MQ=255)                                                                                                                                        gAACACCGCGGCCTGGGCGTCGATGCGGTCGGCGAGGGTGAGCAGCAGTGCCGCACGGTCCTTGGGCGGTGTTTGCGACCAGGCTTCGAAGGCCGCATCCGCGCTTTGCACTGCGGCATGGACCTGTGCATGACTGGCCTCGTGGGTCt       <  6:317159/149‑1 (MQ=255)                                                                                                                                            acCGCGGCCTGGGCGTCGATGCGGTCGGCGAGGGTGAGCAGCAGTGCCGCACGGTCCTTGGGCGGTGTTTGCGACCAGGCTTCGAAGGCCGCATCcgcg                                                     >  2:223594/1‑99 (MQ=255)                                                                                                                                            acCGCGGCCTGGGCGTCGATGCGGTCGGCGAGGGTGAGCAGCAGTGCCGCACGGTCCTTGGGCGGTGTTTGCGACCAGGCTTCGAAGGCCGCATCcgcg                                                     <  1:223594/99‑1 (MQ=255)                                                                                                                                             ccGCGGCCTGGGCGTCGATGCGGTCGGCGAGGGTGAGCAGCAGTGCCGCACGGTCCTTGGGCGGTGTTTGCGACCAGGCTTCGAAGGCCGCATCCGCGCTTTGCACTGCGGCATGGACCTGTGCATGACTGGCCTCGTGGGTCTGTAcc  <  3:140854/149‑1 (MQ=255)                                                                                                                                                    |                                                                                                                                              GTACTCGCCCGCCGCCGAGCCTTGCAGGCAGCGAACGGCGCCGGCGAAGAAGCGGAACACGTCCGACACGGCGGGCAGTTCATCCTGCAGCACGGCCGTGAACGGCTTGCCGCAGTTGTCCGCCTCCAGGCGGGCGAACACCGCGGCCTGGGCGTCGATGCGGTCGGCGAGGGTGAGCAGCAGTGCCGCACGGTCCTTGGGCGGTGTTTGCGACCAGGCTTCGAAGGCCGCATCCGCGCTTTGCACTGCGGCATGGACCTGTGCATGACTGGCCTCGTGGGTCTGTACC  >  NC_002947/3193083‑3193371

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.