New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_002947	3859566 =	44 (0.870)	5 (0.110)	5/242	NT	10.6%	coding (400/753 nt)	cobM	precorrin‑4 C(11)‑methyltransferase
	?	NC_002947	3859608 =	43 (0.910)					coding (442/753 nt)	cobM	precorrin‑4 C(11)‑methyltransferase

TACGGCGCCATCGGTGAGCAGATTCGCCACCTGCAAGCGCTGGGCATCGATTACCAGATCATCCCCGGCGTTACCGCCACGGCAGCGAGCGCCGCCCTGCTCGGCTGCGAGCTGACGTTGCCTCAAGTGGCGCAGACCGTGATCCTCACCCGCTATGGCGACAGCTCGCCCATGCCGCCGGGTGAGCAGTTGGGTGACCTGGCGCGCCATGGCAGCACCCTGGCAATTCATCTGGGGGTCAAGCACCTGTCGCGCATTGTCGCTGAGTTGATGCCGCACTACGGCGCGCACTG  >  NC_002947/3859422‑3859714                                                                                                                                                 |                                                                                                                                                     tACGGCGCCATCGGTGAGCAGATTCGCCACCTGCAAGCGCTGGGCATCGATTACCAGATCATCCCCGGCGTTACCGCCACGGCAGCGAGCGCCGCCCTGCTCGGCTGCGAGCTGACGTTGCCTCAAGTGGCGCAGACCGTGATCCTCAc                                                                                                                                                  >  8:104062/1‑149 (MQ=255)       gcCATCGGTGAGCAGATTCGCCACCTGCAAGCGCTGGGCATCGATTACCAGATCATCCCCGGCGTTACCGCCACGGCAGCGAGCGCCGCCCTGCTCGGCTGCGAGCTGACGTTGCCTCAAGTGGCGCAGACCGTGATcc                                                                                                                                                      >  1:316913/1‑139 (MQ=255)       gcCATCGGTGAGCAGATTCGCCACCTGCAAGCGCTGGGCATCGATTACCAGATCATCCCCGGCGTTACCGCCACGGCAGCGAGCGCCGCCCTGCTCGGCTGCGAGCTGACGTTGCCTCAAGTGGCGCAGACCGTGATcc                                                                                                                                                      <  2:316913/139‑1 (MQ=255)           tCGGTGAGCAGATTCGCCACCTGCAAGCGCTGGGCATCGATTACCAGATCATCCCCGGCGTTACCGCCACGGCAGCGAGCGCCGCCCTGCTCGGCTGCGAGCTGACGTTGCCTCAAGTGGCGCAGACCGTGATCCTCACCCGCTATGGc                                                                                                                                        <  4:149366/149‑1 (MQ=255)            tggcccgatgcaccaccgcaaccgggcagtgcgcgccgtagtgcggcatcaactcagcgacaatgcgcgacaggtgcttgacccccagatgaattgccagggtgctgccatggcgcgccaggtcacccaactgCTCACCCGCTATGGc                                                                                                                                        <  6:206690‑M1/15‑1 (MQ=255)             ggTGAGCAGATTCGCCACCTGCAAGCGCTGGGCATCGATTACCAGATCATCCCCGGCGTTACCGCCACGGCAGCGAGCGCCGCCCTGCTCGGCTGCGAGCTGACGTTGCCTCAAGTGGCGCAGACCGTGATCCTCACCCGCTATGGCGa                                                                                                                                      <  5:329785/149‑1 (MQ=255)               tGAGCAGATTCGCCACCTGCAAGCGCTGGGCATCGATTACCAGATCATCCCCGGCGTTACCGCCACGGCAGCGAGCGCCGCCCTGCTCGGCTGCGAGCTGACGTTGCCTCAAGTGGCGCAGACCGTGATCCTCACCCGCTATGGCGACa                                                                                                                                    <  2:49003/149‑1 (MQ=255)                      aTTCGCCACCTGCAAGCGCTGGGCATCGATTACCAGATCATCCCCGGCGTTACCGCCACGGCAGCGAGCGCCGCCCTGCTCGGCTGCGAGCTGACGTTGCCTCAAGTGGCGCAGACCGTGATCCTCACCCGCTATGGCGACAGCTCGcc                                                                                                                             <  5:23863/149‑1 (MQ=255)                        tCGCCACCTGCAAGCGCTGGGCATCGATTACCAGATCATCCCCGGCGTTACCGCCACGGCAGCGAGCGCCGCCCTGCTCGGCTGCGAGCTGACGTTGCCTCAAGTGGCGCAGACCGTGATCCTCACCCGCTATGGCGACAGCTCGCCCa                                                                                                                           >  4:241855/1‑149 (MQ=255)                          gCCACCTGCAAGCGCTGGGCATCGATTACCAGATCATCCCCGGCGTTACCGCCACGGCAGCGAGCGCCGCCCTGCTCGGCTGCGAGCTGACGTTGCCTCAAGTGGCGCAGACCGTGATCCTCACCCGCTATGGCGACAGCTCGCCCAt                                                                                                                          <  7:104062/148‑1 (MQ=255)                          gCCACCTGCAAGCGCTGGGCATCGATTACCAGATCATCCCCGGCGTTACCGCCACGGCAGCGAGCGCCGCCCTGCTCGGCTGCGAGCTGACGTTGCCTCAAGTGGCGCAGACCGTGATCCTCACCCGCTATGGCGACAGCTCGCCCATg                                                                                                                         <  6:293843/149‑1 (MQ=255)                            cACCTGCAAGCGCTGGGCATCGATTACCAGATCATCCCCGGCGTTACCGCCACGGCAGCGAGCGCCGCCCTGCTCGGCTGCGAGCTGACGTTGCCTCAAGTGGCGCAGACCGTGATCCTCACCCGCTATGGCGACAGCTCGCCCATgcc                                                                                                                       <  7:132341/149‑1 (MQ=255)                              ccTGCAAGCGCTGGGCATCGATTACCAGATCATCCCCGGCGTTACCGCCACGGCAGCGAGCGCCGCCCTGCTCGGCTGCGAGCTGACGTTGCCTCAAGTGGCGCAGACCGTGATCCTCACCCGCTATGGCGACAGCTCGCCCATgccgc                                                                                                                     >  5:177497/1‑149 (MQ=255)                                      cagtgcgcgccgtagtgcggcatcaactcagcgacaatgcgcgacaggtgcttgacccccagatgaattgccagggtgctgccatggcgcgccaggtcacccaactgCTCACCCGCTATGGCGACAGCTCGCCCATGCCGCCGGGTGAg                                                                                                             >  5:294415‑M1/108‑149 (MQ=255)                                       ccTGGGCATCGATTACCAGATCATCCCCGGCGTTACCGCCACGGCAGCGAGCGCCGCCCTGCTCGGCTGCGAGCTGACGTTGCCTCAAGTGGCGCAGACCGTGATCCTCACCCGCTATGGCGACAGCTCGCCCATGCCGCCGGGTGAGc                                                                                                            <  6:189307/148‑1 (MQ=255)                                        gtgcgcgccgtagtgcggcatcaactcagcgacaatgcgcgacaggtgcttgacccccagatgaattgccagggtgctgccatggcgcgccaggtcacccaactgCTCACCCGCTATGGCGACAGCTCGCCCATGCCGCCGGGTGAGCa                                                                                                           >  5:206690‑M1/106‑149 (MQ=255)                                           ggCATCGATTACCAGATCATCCCCGGCGTTACCGCCACGGCAGCGAGCGCCGCCCTGCTCGGCTGCGAGCTGACGTTGCCTCAAGTGGCGCAGACCGTGATCCTCACCCGCTATGGCGACAGCTCGCCCATGCCGCCGGGTGAGCAGt                                                                                                         <  2:210067/148‑1 (MQ=255)                                            gCATCGATTACCAGATCATCCCCGGCGTTACCGCCACGGCAGCGAGCGCCGCCCTGCTCGGCTGCGAGCTGACGTTGCCTCAAGTGGCGCAGACCGTGATCCTCACCCGCTATGGCGACAGCTCGCCCATGCCGCCGGGTGAGCAGTTg                                                                                                       >  4:45124/1‑149 (MQ=255)                                             cATCGATTACCAGATCATCCCCGGCGTTACCGCCACGGCAGCGAGCGCCGCCCTGCTCGGCTGCGAGCTGACGTTGCCTCAAGTGGCGCAGACCGTGATCCTCACCCGCTATGGCGACAGCTCGccc                                                                                                                            >  3:8815/1‑127 (MQ=255)                                             cATCGATTACCAGATCATCCCCGGCGTTACCGCCACGGCAGCGAGCGCCGCCCTGCTCGGCTGCGAGCTGACGTTGCCTCAAGTGGCGCAGACCGTGATCCTCACCCGCTATGGCGACAGCTCGccc                                                                                                                            <  4:8815/127‑1 (MQ=255)                                             cATCGATTACCAGATCATCCCCGGCGTTACCGCCACGGCAGCGAGCGCCGCCCTGCTCGGCTGCGAGCTGACGTTGCCTCAAGTGGCGCAGACCGTGATCCTCACCCGCTATGGCGACAGCTCGCCCATGCCGCCGGGTGAGCAGTTgg                                                                                                      >  5:199839/1‑149 (MQ=255)                                             cATCGATTACCAGATCATCCCCGGCGTTACCGCCACGGCAGCGAGCGCCGCCCTGCTCGGCTGCGAGCTGACGTTGCCTCAAGTGGCGCAGACCGTGATCCTCACCCGCTATGGCGACAGCTCGCCCATGCCGCCGGGTGAGCAGTTgg                                                                                                      >  8:221467/1‑149 (MQ=255)                                               tCGATTACCAGATCATCCCCGGCGTTACCGCCACGGCAGCGAGCGCCGCCCTGCTCGGCTGCGAGCTGACGTTGCCTCAAGTGGCGCAGACCGTGATCCTCACCCGCTATGGCGACAGCTCGCCCATGCCGCCGGGTGAGCAGTTGGGt                                                                                                    <  3:45124/149‑1 (MQ=255)                                                      ccAGATCATCCCCGGCGTTACCGCCACGGCAGCGAGCGCCGCCCTGCTCGGCTGCGAGCTGACGTTGCCTCAAGTGGCGCAGACCGTGATCCTCACCCGCTATGGCGACAGCTCGCCCATGCCGCCGGGTGAGCAGTTGGGTGACCTgg                                                                                             >  7:49993/1‑149 (MQ=255)                                                             atcCCCGGCGTTACCGCCACGGCAGCGAGCGCCGCCCTGCTCGGCTGCGAGCTGACGTTGCCTCAAGTGGCGCAGACCGTGATCCTCACCCGCTATGGCGACAGCTCGCCCATGCCGCCGGGTGAGCAGTTGGGTGACCTGGCGCGcc                                                                                       <  8:240686/148‑1 (MQ=255)                                                              aactcagcgacaatgcgcgacaggtgcttgacccccagatgaattgccagggtgctgccatggcgcgccaggtcacccaactgCTCACCCGCTATGGCGACAGCTCGCCCATGCCGCCGGGTGAGCAGTTGGGt                                                                                                    <  4:139353‑M1/51‑1 (MQ=255)                                                              aactcagcgacaatgcgcgacaggtgcttgacccccagatgaattgccagggtgctgccatggcgcgccaggtcacccaactgCTCACCCGCTATGGCGACAGCTCGCCCATGCCGCCGGGTGAGCAGTTGGGt                                                                                                    >  3:139353‑M1/84‑134 (MQ=255)                                                                cccGGCGTTACCGCCACGGCAGCGAGCGCCGCCCTGCTCGGCTGCGAGCTGACGTTGCCTCAAGTGGCGCAGACCGTGATCCTCACCCGCTATGGCGACAGCTCGCCCATGCCGCCGGGTGAGCAGTTGGGTGACCTGGCGCGCCATgg                                                                                   <  6:199839/149‑1 (MQ=255)                                                                cccGGCGTTACCGCCACGGCAGCGAGCGCCGCCCTGCTCGGCTGCGAGCTGACGTTGCCTCAAGTGGCGCAGACCGTGATCCTCACCCGCTATGGCGACAGCTCGCCCATGCCGCCGGGTGAGCAGTTGGGTGACCTGGCGCGCCATgg                                                                                   <  3:241855/149‑1 (MQ=255)                                                                  cGGCGTTACCGCCACGGCAGCGAGCGCCGCCCTGCTCGGCTGCGAGCTGACGTTGCCTCAAGTGGCGCAGACCGTGATCCTCACCCGCTATGGCGACAGCTCGCCCATGCCGCCGGGTGAGCAGTTGGGTGACCTGgcgcgc                                                                                        <  6:44985/142‑1 (MQ=255)                                                                  cGGCGTTACCGCCACGGCAGCGAGCGCCGCCCTGCTCGGCTGCGAGCTGACGTTGCCTCAAGTGGCGCAGACCGTGATCCTCACCCGCTATGGCGACAGCTCGCCCATGCCGCCGGGTGAGCAGTTGGGTGACCTGgcgcgc                                                                                        >  5:44985/1‑142 (MQ=255)                                                                  cGGCGTTACCGCCACGGCAGCGAGCGCCGCCCTGCTCGGCTGCGAGCTGACGTTGCCTCAAGTGGCGCAGACCGTGATCCTCACCCGCTATGGCGACAGCTCGCCCATGCCGCCGGGTGAGCAGTTGGGTGACCTGGCGCGCCATgg                                                                                   >  2:17466/1‑147 (MQ=255)                                                                      gTTACCGCCACGGCAGCGAGCGCCGCCCTGCTCGGCTGCGAGCTGACGTTGCCTCAAGTGGCGCAGACCGTGATCCTCACCCGCTATGGCGACAGCTCGCCCATGCCGCCGGGTGAGCAGTTGGGTGACCTGGCGCGCCATGGCAGCAc                                                                             <  8:49993/149‑1 (MQ=255)                                                                      gTTACCGCCACGGCAGCGAGCGCCGCCCTGCTCGGCTGCGAGCTGACGTTGCCTCAAGTGGCGCAGACCGTGATCCTCACCCGCTATGGCGACAGCTCGCCCATGCCGCCGGGTGAGCAGTTGGGTGACCTGGCGCGCCATGGCAGCAc                                                                             >  6:132525/1‑149 (MQ=255)                                                                      gTTACCGCCACGGCAGCGAGCGCCGCCCTGCTCGGCTGCGAGCTGACGTTGCCTCAAGTGGCGCAGACCGTGATCCTCACCCGCTATGGCGACAGCTCGCCCATGCCGCCGGGTGAGCAGTTGGGTGACCTGGCGCGCCATGGCAGCAc                                                                             <  2:50087/149‑1 (MQ=255)                                                                             ccACGGCAGCGAGCGCCGCCCTGCTCGGCTGCGAGCTGACGTTGCCTCAAGTGGCGCAGACCGTGATCCTCACCCGCTATGGCGACAGCTCGCGCATGCCGCCGGGTGAGCAGTTgg                                                                                                      >  5:51729/1‑117 (MQ=255)                                                                             ccACGGCAGCGAGCGCCGCCCTGCTCGGCTGCGAGCTGACGTTGCCTCAAGTGGCGCAGACCGTGATCCTCACCCGCTATGGCGACAGCTCGCCCATGCCGCCGGGTGAGCAGTTgg                                                                                                      <  6:51729/117‑1 (MQ=255)                                                                             ccACGGCAGCGAGCGCCGCCCTGCTCGGCTGCGAGCTGACGTTGCCTCAAGTGGCGCAGACCGTGATCCTCACCCGCTATGGCGACAGCTCGCCCATGCCGCCGGGTGAGCAGTTGGGTGACCTGGCGCGCCATGGCAGCACCCTGGc                                                                       >  2:262928/1‑148 (MQ=255)                                                                              cACGGCAGCGAGCGCCGCCCTGCTCGGCTGCGAGCTGACGTTGCCTCAAGTGGCGCAGACCGTGATCCTCACCCGCTATGGCGACAGCTCGCCCATGCCGCCGGGTGAGCAGTTGGGTGACCTGGCGCGCCATGGCAGCACCCTGGCaa                                                                     <  8:181462/149‑1 (MQ=255)                                                                                 ggCAGCGAGCGCCGCCCTGCTCGGCTGCGAGCTGACGTTGCCTCAAGTGGCGCAGACCGTGATCCTCACCCGCTATGGCGACAGCTCGCCCATGCCGCCGGGTGAGCAGTTGGGTGACCTGGCGCGCCATGGCAGCACCCTGGCAATTc                                                                  <  5:94263/149‑1 (MQ=255)                                                                                     gtgcttgacccccagatgaattgccagggtgctgccatggcgcgccaggtcacccaactgc                                                                                                                                                      <  2:21171‑M1/1‑1 (MQ=255)                                                                                     gtgcttgacccccagatgaattgccagggtgctgccatggcgcgccaggtcacccaactgc                                                                                                                                                      >  1:21171‑M1/61‑61 (MQ=255)                                                                                         gcgcCGCCCTGCTCGGCTGCGAGCTGACGTTGCCTCAAGTGGCGCAGACCGTGATCCTCACCCGCTATGGCGACAGCTCGCCCATGCCGCCTGGTGAGCAGTTGGGTGACCTGGCGCGCCATGGCAGCACCCTGGCAATTCATCTgggg                                                          <  3:303913/149‑1 (MQ=255)                                                                                         gcgcCGCCCTGCTCGGCTGCGAGCTGACGTTGCCTCAAGTGGCGCAGACCGTGATCCTCACCCGCTATGGCGACAGCTCGCCCATGCCGCCGGGTGAGCAGTTGGGTGACCTGGCGCGCCATGGCAGCACCCTGGCAATTCATCTgggg                                                          <  1:126156/149‑1 (MQ=255)                                                                                         gcgcCGCCCTGCTCGGCTGCGAGCTGACGTTGCCTCAAGTGGCGCAGACCGTGATCCTCACCCGCTATGGCGACAGCTCGCCCATGCCGCCGGGTGAGCAGTTGGGTGACCTGGCGCGCCATGGCAGCACCCTGGCAATTCATCTgggg                                                          <  1:262928/149‑1 (MQ=255)                                                                                                ccTGCTCGGCTGCGAGCTGACGTTGCCTCAAGTGGCGCAGACCGTGATCCTCACCCGCTATGGCGACAGCTCGCCCATGCCGCCGGGTGAGCAGTTGGGTGACCTGGCGCGCCATGGCAGCACCCTGGCAAtt                                                                   >  6:300126/1‑133 (MQ=255)                                                                                                ccTGCTCGGCTGCGAGCTGACGTTGCCTCAAGTGGCGCAGACCGTGATCCTCACCCGCTATGGCGACAGCTCGCCCATGCCGCCGGGTGAGCAGTTGGGTGACCTGGCGCGCCATGGCAGCACCCTGGCAAtt                                                                   <  5:300126/133‑1 (MQ=255)                                                                                                 cTGCTCGGCTGCGAGCTGACGTTGCCTCAAGTGGCGCAGACCGTGATCCTCACCCGCTATGGCGACAGCTCGCCCATGCCGCCGGGTGAGCAGTTGGGTGACCTGGCGCGCCATGGCAGCACCCTGGCAATTCATCTGGGGGTCAAGCa                                                  >  8:304959/1‑149 (MQ=255)                                                                                                          tGCGAGCTGACGTTGCCTCAAGTGGCGCAGACCGTGATCCTCACCCGCTATGGCGACAGCTCGCCCATGCCGCCGGGTGAGCAGTTGGGTGACCTGGCGCGCCATGGCAGCACCCTGGCAATTCATCTGGGGGTCAAGCACCTGTcgcg                                         >  3:175753/1‑149 (MQ=255)                                                                                                               gCTGACGTTGCCTCAAGTGGCGCAGACCGTGATCCTCACCCGCTATGGCGACAGCTCGCCCATGCCGCCGGGTGAGCAGTTGGGTGACCTGGCGCGCCATGGCAGCACCCTGGCAATTCATCTGGGGGTCAAGCACCTGTCGCGCATTg                                    >  7:158203/1‑149 (MQ=255)                                                                                                                   aCGTTGCCTCAAGTGGCGCAGACCGTGATCCTCACCCGCTATGGCGACAGCTCGCCCATGCCGCCGGGTGAGCAGTTGGGTGACCTGGCGCGCCATGGCAGCACCCTGGCAATTCATCTGGGGGTCAAGCACCTGTCGCGCATTGTCGc                                >  7:123450/1‑149 (MQ=255)                                                                                                                         ccTCAAGTGGCGCAGACCGTGATCCTCACCCGCTATGGCGACAGCTCGCCCATGCCGCCGGGTGAGCAGTTGGGTGACCTGGCGCGCCATGGCAGCACCCTGGCAATTCATCTGGGGGTCAAGCACCTGTCGCGCATTGTCGCTGAGt                           <  7:304959/148‑1 (MQ=255)                                                                                                                                    gcAGACCGTGATCCTCACCCGCTATGGCGACAGCTCGCCCATGCCGCCGGGTGAGCAGTTGGGTGACCTGGCGCGCCATGGCAGCACCCTGGCAATTCATCTGGGGGTCAAGCACCTGTCGCGCATTGTCGCTGAGTTGATGCCGCAc                >  1:105925/1‑148 (MQ=255)                                                                                                                                        aCCGTGATCCTCACCCGCTATGGCGACAGCTCGCCCATGCCGCCGGGTGAGCAGTTGGGTGACCTGGCGCGCCATGGCAGCACCCTGGCAATTCATCTGGGGGTCAAGCACCTGTCGCGCATTGTCGCTGAGTTGATGCCGCACTACgg           <  4:175753/149‑1 (MQ=255)                                                                                                                                            tGATCCTCACCCGCTATGGCGACAGCTCGCCCATGCCGCCGGGTGAGCAGTTGGGTGACCTGGCGCGCCATGGCAGCACCCTGGCAATTCATCTGGGGGTCAAGCACCTGTCGCGCATTGTCGCTGAGTTGATGCCGCACTACGgcgcg       <  3:282673/149‑1 (MQ=255)                                                                                                                                                 CTCACCCGCTATGGCGACAGCTCGCCCATGCCGCCGGGTGAGCAGTTGGGTGACCTGGCGCGCCATGGCAGCACCCTGGCAATTCATCTGGGGGTCAAGCACCTGTCGCGCATTGTCGCTGAGTTGATGCCGCACTACGGCGCGCACTg  <  7:21234‑M1/149‑1 (MQ=255)                                                                                                                                                 CTCACCCGCTATGGCGACAGCTCGCCCATGCCGCCGGGTGAGCAGTTGGGTGACCTGGCGCGCCATGGCAGCACCCTGGCAATTCATCTGGGGGTCAAGCACCTGTCGCGCATTGTCGCTGAGGTGATGCCGCACTACGGCGCGCACTg  <  5:331555‑M1/149‑1 (MQ=255)                                                                                                                                                 |                                                                                                                                                     TACGGCGCCATCGGTGAGCAGATTCGCCACCTGCAAGCGCTGGGCATCGATTACCAGATCATCCCCGGCGTTACCGCCACGGCAGCGAGCGCCGCCCTGCTCGGCTGCGAGCTGACGTTGCCTCAAGTGGCGCAGACCGTGATCCTCACCCGCTATGGCGACAGCTCGCCCATGCCGCCGGGTGAGCAGTTGGGTGACCTGGCGCGCCATGGCAGCACCCTGGCAATTCATCTGGGGGTCAAGCACCTGTCGCGCATTGTCGCTGAGTTGATGCCGCACTACGGCGCGCACTG  >  NC_002947/3859422‑3859714

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.