New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_002947	1490726 =	48 (0.920)	3 (0.060)	3/254	NT	5.8%	coding (389/1902 nt)	cysNC	sulfate adenylyltransferase subunit 1/adenylyl‑sulfate kinase
	?	NC_002947	1490781 =	51 (1.030)					coding (444/1902 nt)	cysNC	sulfate adenylyltransferase subunit 1/adenylyl‑sulfate kinase

GGCCTGCAGGCCGAGCGCGAGCAGGGCATCACCATCGATGTCGCCTACCGCTACTTCTCCACCGCCAAGCGCAAGTTCATCATCGCCGACACCCCGGGCCATGAGCAGTACACCCGCAACATGGCCACCGGTGCGTCCACCTGCGACCTGGCGATCATCCTGGTCGATGCCCGCTACGGCGTGCAGACCCAGACCCGTCGCCACAGCTACATCGCGTCGTTGCTGGGTATCAAGCACATCGTGGTCGCGGTCAACAAGATGGACCTGAAGGGCTTCGATGAAGGCGTCTTCGAGTC  >  NC_002947/1490578‑1490873                                                                                                                                                     |                                                                                                                                                    gctgctcggggtcaggttgatggcttcggcgaacttcaggtagtcggccttgatcgactcgaagacgccttcatcgaagcccttcaggtccatcttgttgaccgcgaccacgatgtgcttgatacccagcaacgacgcgatgtagctgt                                                                                                                                                     <  1:207129‑M1/1‑1 (MQ=255) gctgctcggggtcaggttgatggcttcggcgaacttcaggtagtcggccttgatcgactcgaagacgccttcatcgaagcccttcaggtccatcttgttgaccgcgaccacgatgtgcttgatacccagcaacgacgcgatgtagctgt                                                                                                                                                     <  4:302318‑M1/1‑1 (MQ=255)   ccTGCAGGCCGAGCGCGAGCAGGGCATCACCATCGATGTCGCCTACCGCTACTTCTCCACCGCCAAGCGCAAGTTCATCATCGCCGACACCCCGGGCCATGAGCAGTACACCCGCAACATGGCCACCGGTGCGTCCACCTGCGACCTgg                                                                                                                                                   <  8:250523/149‑1 (MQ=255)       cAGGCCGAGCGCGAGCAGGGCATCACCATCGATGTCGCCTACCGCTACTTCTCCACCGCCAAGCGCAAGTTCATCATCGCCGACACCCCGGGCCATGAGTAGTACACCCGCAACATGGCCACCGGTGCGTCCACCTGCGACCTGGCg                                                                                                                                                 <  7:194531/147‑1 (MQ=255)        aGGCCGAGCGCGAGCAGGGCATCACCATCGATGTCGCCTACCGCTACTTCTCCACCGCCAAGCGCAAGTTCATCATCGCCGACACCCCGGGCCATGAGCAGTACACCCGCAACATGGCCACCGGTGCGTCCACCTGCGACCTGGCGatc                                                                                                                                              <  2:118367/149‑1 (MQ=255)        aGGCCGAGCGCGAGCAGGGCATCACCATCGATGTCGCCTACCGCTACTTCTCCACCGCCAAGCGCAAGTTCATCATCGCCGACACCCCGGGCCATGAGCAGTACACCCGCAACATGGCCACCGGTGCGTCCACCTGCGACCTGGCGat                                                                                                                                               >  7:189484/1‑148 (MQ=255)          gCCGAGCGCGAGCAGGGCATCACCATCGATGTCGCCTACCGCTACTTCTCCACCGCCAAGCGCAAGTTCATCATCGCCGACACCCCGGGCCATGAGCAGTACACCCGCAACATGGCCACCGGTGCGTCCACCTGCGACCTGGCGatcat                                                                                                                                            >  4:135420/1‑149 (MQ=255)            cGAGCGCGAGCAGGGCATCACCATCGATGTCGCCTACCGCTACTTCTCCACCGCCAAGCGCAAGTTCATCATCGCCGACACCCCGGGCCATGAGCAGTACACCCGCAACATGGCCACCGGTGCGTCCACCTGCGACCTGGCGATCATcc                                                                                                                                          <  8:189484/149‑1 (MQ=255)               gcgcGAGCAGGGCATCACCATCGATGTCGCCTACCGCTACTTCTCCACCGCCAAGCGCAAGTTCATCATCGCCGACACCCCGGGCCATGAGCAGTACACCCGCAACATGGCCACCGGTGCGTCCACCTGCGACCTGGCGATCATCCTgg                                                                                                                                       <  8:186984/149‑1 (MQ=255)                 gcgAGCAGGGCATCACCATCGATGTCGCCTACCGCTACTTCTCCACCGCCAAGCGCAAGTTCATCATCGCCGACACCCCGGGCCATGAGCAGTACACCCGCAACATGGCCACCGGTGCGTCCACCTGCGACCTGGCGATCATCCTGGTc                                                                                                                                     >  5:159987/1‑149 (MQ=255)                        gggCATCACCATCGATGTCGCCTACCGCTACTTCTCCACCGCCAAGCGCAAGTTCATCATCGCCGACACCCCGGGCCATGAGCAGTACACCCGCAACATGGCCACCGGTGCGTCCACCTGCGACCTGGCGATCATCCTGGTCGATGccc                                                                                                                              >  6:301820/1‑149 (MQ=255)                           cATCACCATCGATGTCGCCTACCGCTACTTCTCCACCGCCAAGCGCAAGTTCATCATCGCCGACACCCCGGGCCATGAGCAGTACACCCGCAACATGGCCACCGGTGCGTCCACCTGCGACCTGGCGATCATCCTGGTCGATGCCCGc                                                                                                                            >  3:60024/1‑148 (MQ=255)                                  aTCGATGTCGCCTACCGCTACTTCTCCACCGCCAAGCGCAAGTTCATCATCGCCGACACCCCGGGCCATGAGCAGTACACCCGCAACATGGCCACCGGTGCGTCCACCTGCGACCTGGCGATCATCCTGGTCGATGCCCGCTACGGc                                                                                                                      <  6:28612/147‑1 (MQ=255)                                  aTCGATGTCGCCTACCGCTACTTCTCCACCGCCAAGCGCAAGTTCATCATCGCCGACACCCCGGGCCATGAGCAGTACACCCGCAACATGGCCACCGGTGCGTCCACCTGCGACCTGGCGATCATCCTGGTCGATGCCCGCTACGGc                                                                                                                      >  5:28612/1‑147 (MQ=255)                                   tCGATGTCGCCTACCGCTACTTCTCCACCGCCAAGCGCAAGTTCATCATCGCCGACACCCCGGGCCATGAGCAGTACACCCGCAACATGGCCACCGGTGCGTCCACCTGCGACCTGGCGATCATCCTGGTCGATGCCCGCTACGGCGTg                                                                                                                   <  4:108502/149‑1 (MQ=255)                                     gATGTCGCCTACCGCTACTTCTCCACCGCCAAGCGCAAGTTCATCATCGCCGACACCCCGGGCCATGAGCAGTACACCCGCAACATGGCCACCGGTGCGTCCACCTGCGACCTGGCGATCATCCTGGTCGATGCCCGCTACGGCGTGca                                                                                                                 <  1:82809/149‑1 (MQ=255)                                             cTACCGCTACTTCTCCACCGCCAAGCGCAAGTTCATCATCGCCGACACCCCGGGCCATGAGCAGTACACCCGCAACATGGCCACCGGTGCGTCCACCTGCGACCTgg                                                                                                                                                   <  2:224885/107‑1 (MQ=255)                                             cTACCGCTACTTCTCCACCGCCAAGCGCAAGTTCATCATCGCCGACACCCCGGGCCATGAGCAGTACACCCGCAACATGGCCACCGGTGCGTCCACCTGCGACCTgg                                                                                                                                                   >  1:224885/1‑107 (MQ=255)                                                 cGCTACTTCTCCACCGCCAAGCGCAAGTTCATCATCGCCGACACCCCGGGCCATGAGTAGTACACCCGCAACATGGCCACCGGTGCGTCCACCTGCGACCTGGCGATCATCCTGGTCGATGCCCGCTACGGCGTGCAGACCCAGACCCg                                                                                                     >  3:302318/1‑149 (MQ=255)                                                   cTACTTCTACACCGCCAAGCGCAAGTTCATCATCGCCGACACCCCGGGCAATGAGCAGTACACCCGCAACATGGCCACCGGTGCGTCCACCTGCGACCTGGCGATCATCCTGGACGATGCCCGCTACGGCGTGCAGACCCAGACCCGTc                                                                                                   >  8:45275/1‑149 (MQ=255)                                                            cgaagacgccttcatcgaagcccttcaggtccatcttgttgaccgcgaccacgatgtgcttgatacccagcaacgacgcgatgtagctgTGGCGATCATCCTGGTCGATGCCCGCTACGGCGTGCAGACACAGACCCGTCGCCACAGCt                                                                                          >  3:28819‑M1/90‑149 (MQ=255)                                                            cACCGCCAAGCGCAAGTTCATCATCGCCGACACCCCGGGCCATGAGCAGTACACCCGCAACATGGCCACCGGTGCGTCCACCTGCGACCTGGCGATCATCCTGGTCGATGCCCGCTACGGCGTGCAGACCCAGACCCGTCGCCACAGCt                                                                                          <  6:159987/149‑1 (MQ=255)                                                               cGCCAAGCGCAAGTTCATCATCGCCGACACCCCGGGCCATGAGCAGTACACCCGCAACATGGCCACCGGTGCGTCCACCTGCGACCTGGCGATCATCCTGGTCGATGCCCGCTACGGCGTGCAGACCCAGACCCGTCGCCACAGCTAc                                                                                        >  4:58159/1‑148 (MQ=255)                                                                      cgcAAGTTCATCATCGCCGACACCCCGGGCCATGAGCAGTACACCCGCAACATGGCCACCGGTGCGTCCACCTGCGACCTGGCGATCATCCTGGTCGATGCCCGCTACGGCGTGCAGACCCAGACCCGTCGCCACAGCTACATCGcgtc                                                                                <  3:58159/149‑1 (MQ=255)                                                                         atcgaagcccttcaggtccatcttgttgaccgcgaccacgatgtgcttgatacccagcaacgacgcgatgtagctgTGGCGATCATCCTGGTCGATGCCCGCTACGGCGTGCAGACCCAGACCCGTCGCCACAGCTACATc                                                                                     >  3:313634‑M1/77‑141 (MQ=255)                                                                         atcgaagcccttcaggtccatcttgttgaccgcgaccacgatgtgcttgatacccagcaacgacgcgatgtagctgTGGCGATCATCCTGGTCGATGCCCGCTACGGCGTGCAGACCCAGACCCGTCGCCACAGCTACATc                                                                                     <  4:313634‑M1/65‑1 (MQ=255)                                                                             tcatcaTCGCCGACACCCCGGGCCATGAGCAGTACACCCGCAACATGGCCACCGGTGCGTCCACCTGCGACCTGGCGATCATCCTGGTCGATGCCCGCTACGGCGTGCAGACCCAGACCCGTCGCCACAGCTACATCGCGTCGTTGCTg                                                                         >  2:291255/1‑149 (MQ=255)                                                                               atcatcGCCGACACCCCGGGCCATGAGCAGTACACCCGCAACATGGCCACCGGTGCGTCCACCTGCGACCTGGCGATCATCCTGGTCGATGCCCGCTACGGCGTGCAGACCCAGACCCGTCGCCACAGCTACATCGCGTCGTTGCTggg                                                                       >  6:251677/1‑149 (MQ=255)                                                                               atcatcGCCGACACCCCGGGCCATGAGCAGTACACCCGCAACATGGCCACCGGTGCGTCCACCTGCGACCTGGCGATCATCCTGGTCGATGCCCGCTACGGCGTGCAGACCCAGACCCGTCGCCACAGCTACATCGCGTCGTTGCTgg                                                                        <  5:301820/148‑1 (MQ=255)                                                                                   tcGCCGACACCCCGGGCCATGAGCAGTACACCCGCAACATGGCCACCGGTGCGTCCACCTGCGACCTGGCGATCATCCTGGTCGATGCCCGCTACGGCGTGCAGACCCAGACCCGTCGCCACAGCTACATCGCGTCGTTGCTGGGTAt                                                                    >  1:70071/1‑148 (MQ=255)                                                                                     gCCGACACCCCGGGCCATGAGCAGTACACCCGCAACATGGCCACCGGTGCGGCCACCTGCGACCTGGCGATCAGCCTGGTCGAGGCCCGCTACGGCGTGCAGACCCAGACCCGTCGCCACAGCTCCATCGCGTCGTTGCTGGGTATCaa                                                                 <  8:212421/149‑1 (MQ=255)                                                                                                   ccATGAGCAGTACACCCGCAACATGGCCACCGGTGCGTCCACCTGCGACCTGGCGATCATCCTGGTCGATGCCCGCTACGGCGTGCAGACCCAGACCCGTCGCCACAGCTACATCGCGTCGTTGCTGGGTATCAAGCACATCGTGGTcg                                                   >  3:13089/1‑149 (MQ=255)                                                                                                     aTGAGCAGTACACCCGCAACATGGCCACCGGTGCGTCCACCTGCGACCTGGCGATCATCCTGGTCGATGCCCGCTACGGCGTGCAGACCCAGACCCGTCGCCACAGCTACATCGCGTCGTTGCTGGGTATCAAGCACATCGTGGTcgcg                                                 <  5:209907/149‑1 (MQ=255)                                                                                                         gCAGTACACCCGCAACATGGCCACCGGTGCGTCCACCTGCGACCTGGCGATCATCCTGGTCGATGCCCGCTACGGCGTGCAGACCCAGACCCGTCGCCACAGCCACATCGCGTCGTTGCTGGGTATCAAGCACATCGTGGTCGCGGTc                                              >  3:51736/1‑148 (MQ=255)                                                                                                            gTACACCCGCAACATGGCCACCGGTGCGTCCACCTGCGACCTGGCGATCATCCTGGTCGATGCCCGCTACGGCGTGCAGACCCAGACCCGTCGCCACAGCTACATCGCGTCGTTGCTGGGTATCAAGCACATCGTGGTCGCGGTcaaca                                          <  3:135420/149‑1 (MQ=255)                                                                                                             tACACCCGCAACATGGCCACCGGTGCGTCCACCTGCGACCTGGCGATCATCCTGGTCGATGCCCGCTACGGCGTGCAGACCCAGACCCGTCGCCACAGCTACATCGCGTCGTTGCTGGGTATCAAGCACATCGTGGTCGCGGTcaac                                           >  2:58548/1‑147 (MQ=255)                                                                                                                acCCGCAACATGGCCACCGGTCCGTCCACCTGCGACCTGGCGATCATCCTGGTCGATGCCCGCTACGGCGTGCAGACCCAGACCCGTCGCCACAGCTACATCGCGTCGTTGCTGGGTATCAAGCACATCGTGGTCGCGGTCAACAAGAt                                      >  2:92446/1‑149 (MQ=255)                                                                                                                 cccGCAACATGGCCACCGGTGCGTCCACCTGCGACCTGGCGATCATCCTGGTCGATGCCCGCTACGGCGTGCAGACCCAGACCCGTCGCCACAGCTACATCGCGTCGTTGCTGGGTATCAAGCACATCGTGGTCGCGGTCAACAAGATg                                     <  4:13089/149‑1 (MQ=255)                                                                                                                  ccGCAACATGGCCACCGGTGCGTCCACCTGCGACCTGGCGATCATCCTGGTCGATGCCCGCTACGGCGTGCAGACCCAGACCCGTCGCCAGAGCTACATCGCGTCGTTGCTGGGTATCAAGCACATCGTGGTCGCGGTCAACAAGATgg                                    <  4:201953/149‑1 (MQ=255)                                                                                                                  ccGCAACATGGCCACCGGTGCGTCCACCTGCGACCTGGCGATCATCCTGGTCGATGCCCGCTACGGCGTGCAGACCCAGACCCGTCGCCACAGCTACATCGCGTCGtt                                                                             >  1:149804/1‑108 (MQ=255)                                                                                                                  ccGCAACATGGCCACCGGTGCGTCCACCTGCGACCTGGCGATCATCCTGGTCGATGCCCGCTACGGCGTGCAGACCCAGACCCGTCGCCACAGCTACATCGCGTCGtt                                                                             <  2:149804/108‑1 (MQ=255)                                                                                                                  ccGCAACATGGCCACCGGTGCGTCCACCTGCGACCTGGCGATCATCCTGGTCGATGCCCGCTACGGCGTGCAGACCCAGACCCGTCGCCACAGCTACATCGCGTCGTTGCTGGGTATCAAGCACATCGTGGTCGCGGTCAACAAGATg                                     >  5:172170/1‑148 (MQ=255)                                                                                                                  ccGCAACATGGCCACCGGTGCGTCCACCTGCGACCTGGCGATCATCCTGGTCGATGCCCGCTACGGCGTGCAGACCCAGACCCGTCGCCACAGCTACATCGCGTCGTTGCTGGGTATCAAGCACATCGTGGTCGCGGTCAACAAGATg                                     >  5:34191/1‑148 (MQ=255)                                                                                                                      aaCATGGCCACCGGTGCGTCCACCTGCGACCTGGCGATCATCCTGGTCGATGCCCGCTACGGCGTGCAGACCCAGACCCGTCGCCACAGCTACATCGCGTCGTTGCTGGGTATCAAGCACATCGTGGTCGCGGTCAACAAGATGGACCt                                <  1:58548/149‑1 (MQ=255)                                                                                                                         aTGGCCACCGGTGCGTCCACCTGCGACCTGGCGATCATCCTGGTCGATGCCCGCTACGGCGTGCAGACCCAGACCCGTCGCCACAGCTACATCGCGTCGTTGCTGGGTATCAAGCACATCGTGGTCGCGGTCAACAAGATGGACCTGaa                             <  6:172170/149‑1 (MQ=255)                                                                                                                                  ggTGCGTCCACCTGCGACCTGGCGATCATCCTGGTCGATGCCCGCTACGGCGTGCAGACCCAGACCCGTCGCCACAGCTACATCGCGTCGTTGCTGGGTATCAAGCACATCGTGGTCGCGGTCAACAAGATGGACCTGAAGGGCTTCGa                    <  4:60024/149‑1 (MQ=255)                                                                                                                                   gTGCGTCCACCTGCGACCTGGCGATCATCCTGGTCGATGCCCGCTACGGCGTGCAGACCCAGACCCGTCGCCACAGCTACATCGCGTCGTTGCTGGGTATCAAGCACATCGTGGTCGCGGTCAACAAGATGGACCTGAAGGGCTTCGAt                   <  6:34191/149‑1 (MQ=255)                                                                                                                                      cGTCCACCTGCGACCTGGCGATCATCCTGGTCGATGCCCGCTACGGCGTGCAGACCCAGACCCGTCGCCACAGCTACATCGCGTCGTTGCTGGGTATCAAGCACATCGTGGTCGCGGTCAACAAGATGGACCTGAAGGGCTTCGATGa                 <  2:70071/148‑1 (MQ=255)                                                                                                                                       gTCCACCTGCGACCTGGCGATCATCCTGGTCGATGCCCGCTACGGCGTGCAGACCCAGACCCGTCGCCACAGCTACATCGCGTCGTTGCTGGGTATCAAGCACATCGTGGTCGCGGTCAACAAGATGGACCTGAAGGGCTTCGATGAAg               <  1:291255/149‑1 (MQ=255)                                                                                                                                         ccACCTGCGACCTGGCGATCATCCTGGTCGATGCCCGCTACGGCGTGCAGACCCAGACCCGTCGCCACAGCTACATCGCGTCGTTGCTGGGTATCAAGCACATCGTGGTCGCGGTCAACAAGATGGACCTGAAGGGCTTCGAt                   >  5:218375/1‑143 (MQ=255)                                                                                                                                         ccACCTGCGACCTGGCGATCATCCTGGTCGATGCCCGCTACGGCGTGCAGACCCAGACCCGTCGCCACAGCTACATCGCGTCGTTGCTGGGTATCAAGCACATCGTGGTCGCGGTCAACAAGATGGACCTGAAGGGCTTCGAt                   <  6:218375/143‑1 (MQ=255)                                                                                                                                         ccACCTGCGACCTGGCGATCATCCTGGTCGATGCCCGCTACGGCGTGCAGACCCAGACCCGTCGCCACAGCTACATCGCGTCGTTGCTGGGTATCAAGCACATCGTGGTCGCGGTCAACAAGATGGACCTGAAGGGCTTCGATGAAGGc             >  7:96415/1‑149 (MQ=255)                                                                                                                                             cTGCGACCTGGCGATCATCCTGGTCGATGCCCGCTACGGCGTGCAGACCCAGACCCGTCGCCACAGCTACATCGCGTCGTTGCTGGGTATCAAGCACATCGTGGTCGCGGTCAACAAGATGGAc                                  >  2:165328/1‑124 (MQ=255)                                                                                                                                             cTGCGACCTGGCGATCATCCTGGTCGATGCCCGCTACGGCGTGCAGACCCAGACCCGTCGCCACAGCTACATCGCGTCGTTGCTGGGTATCAAGCACATCGTGGTCGCGGTCAACAAGATGGAc                                  <  1:165328/124‑1 (MQ=255)                                                                                                                                              tGCGACCTGGCGATCATCCTGGTCGATGCCCGCTACGGCGTGCAGACCCAGACCCGTCGCCACAGCTACATCGCGTCGTTGCTGGGTATCAAGCACATCGTGGTCGCGGTCAACAAGATGGACCTGAAg                            <  4:129039/129‑1 (MQ=255)                                                                                                                                              tGCGACCTGGCGATCATCCTGGTCGATGCCCGCTACGGCGTGCAGACCCAGACCCGTCGCCACAGCTACATCGCGTCGTTGCTGGGTATCAAGCACATCGTGGTCGCGGTCAACAAGATGGACCTGAAg                            >  3:129039/1‑129 (MQ=255)                                                                                                                                               gCGACCTGGCGATCATCCTGGTCGATGCCCGCTACGGCGTGCAGACCCAGACCCGTCGCCACAGCTACATCGCGTCGTTGCTGGGTATCAAGCACATCGTGGTCGCGGTCAACAAGATGGACCTGAAGGGCTTCGATGAAGGCGTCTTc       >  5:62664/1‑149 (MQ=255)                                                                                                                                                cGACCTGGCGATCATCCTGGTCGATGCCCGCTACGGCGTGCAGACCCAGACCCGTCGCCACAGCTACATCGCGTCGTTGCTGGGTATCAAGCACATCGTGGTCGCGGTCAACAAGATGGACCTGAAGGGCTTCGATGAAGGCGTCTTc       <  1:33061/148‑1 (MQ=255)                                                                                                                                                    cTGGCGATCATCCTGGTCGATGCCCGCTACGGCGTGCAGACCCAGACCCGTCGCCACAGCTACATCGCGTCGTTGCTGGGTATCAAGCACATCGTGGTCGCGGTCAACAAGATGGACCTGAAGGGCTTCGATGAAGGCGTCTTCGAGTc  >  6:12531/1‑149 (MQ=255)                                                                                                                                                     |                                                                                                                                                    GGCCTGCAGGCCGAGCGCGAGCAGGGCATCACCATCGATGTCGCCTACCGCTACTTCTCCACCGCCAAGCGCAAGTTCATCATCGCCGACACCCCGGGCCATGAGCAGTACACCCGCAACATGGCCACCGGTGCGTCCACCTGCGACCTGGCGATCATCCTGGTCGATGCCCGCTACGGCGTGCAGACCCAGACCCGTCGCCACAGCTACATCGCGTCGTTGCTGGGTATCAAGCACATCGTGGTCGCGGTCAACAAGATGGACCTGAAGGGCTTCGATGAAGGCGTCTTCGAGTC  >  NC_002947/1490578‑1490873

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.