New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_002947	= 65525	43 (0.830)	3 (0.060)	3/258	NT	5.6%	coding (583/1653 nt)	betA‑I	choline dehydrogenase
	?	NC_002947	= 65595	59 (1.170)					coding (653/1653 nt)	betA‑I	choline dehydrogenase

GGGCGAATGGCGGGTCGAGCGGCAACGCTACAGCTGGCCGATCCTTGATGCTTTCCGAGATGCCGCCGAGCAGAGCG‑GTATTGGCAAGGTCGACGATTTCAACACCGGCGACAACCAGGGTTGTG‑GATATTTTCAGGTCAACCAGCGCAGTGGTGTGCGCTGGAATGCTTCCAAAGCTTTCCTGCGGCCGATCAAGGATCGCGCCAACCTCACTGTGCTGACCGGCGTGCAGGTTGACCAAGTACTGCTCGACAACACCCGCGCAAGGGCGGTGAAGGCGCTCTGGCAAGGGGC  >  NC_002947/65380‑65673                                                                                                                                                    |                                                                                                                                                     gggAGAATGGCGGGTCGAGCGGCAACGCTACAGCTGGCCGATCCTTGATGCTTTCCGAGATGCCGCCGAGCAGAGCG‑GTATTGGCAAGGTCGACGATTTCAACACCGGCGACAACCAGGGTTGTG‑GATATTTTCAGGTCAACCAgcg                                                                                                                                                     >  6:214784/1‑147 (MQ=255)    cGAATGGCGGGTCGAGCGGCAACGCTACAGCTGGCCGATCCTTGATGCTTTCCGAGATGCCGCCGAGCAGAGCG‑GTATTGGCAAGGTCGACGATTTCAACACCGGCGACAACCAGGGTTGTG‑GATATTTTCAGGTCAACCAGCGCAgtg                                                                                                                                                >  2:12234/1‑149 (MQ=255)      aaTGGCGGGTCGAGCGGCAACGCTACAGCTGGCCGATCCTTGATGCTTTCCGAGATGCCGCCGAGCAGAGCG‑GTATTGGCAAGGTCGACGATTTCAACACCGGCGACAACCAGGGTTGTG‑GATATTTTCAGGTCAACCAGCGCAgtggt                                                                                                                                              <  5:214784/149‑1 (MQ=255)                    cGGCAACGCTACAGCTGGCCGATCCTTGATGCTTTCCGAGATGCCGCCGAGCAGAGCG‑GTATTGGCAAGGTCGACGATTTCAACACCGGCGACAACCAGGGTTGTG‑GATATTTTCAGGTCAACCAGCGCAGTGGTGTGCGCTGGAATGc                                                                                                                                >  1:282455/1‑149 (MQ=255)                        aaagCTACAGCTGGCCGATCCTTGATGCTTTCCGAGATGCCGCCGAGCAGAGCG‑GTATTGGCAAGGTCGACGATTTCAACACCGGCGACAACCAGGGTTGTG‑GATATTTTCAGGTCAACCAGCGCAGTGGTGTGCGCTGGAATGCTTcc                                                                                                                            <  3:93860/146‑1 (MQ=255)                        aaCGCTACAGCTGGCCGATACTTGATGCTTTCCGAGATGCCGCAGAGCAGAGCGCGTATTGGAAAGGTCGACGATTTCAACACCGGCGACAACCAGGGTTGTG‑GATATTTTCAGGTCAACCAGCGCAGTGGTGTGAGCTGGAATGCTTc                                                                                                                             >  1:89274/1‑149 (MQ=255)                             tACAGCTGGCCGATCCTTGATGCTTTCCGAGATGCCGCCGAGCAGAGCG‑GTATTGGCAAGGTCGACGATTTCAACACCGGCGACAACCAGGGTTGTG‑GATATTTTCAGGTCAACCAGCGCAGTGGTGTGCGCTGGAATGCTTCCAAAGc                                                                                                                       >  8:3156/1‑149 (MQ=255)                                aGCTGGCCGATCCTTGATGCTTTCCGAGATGCCGCCGAGCAGAGCG‑GTATTGGCAAGGTCGACGATTTCAACACCGGCGACAACCAGGGTTGTG‑GATATTTTCAGGTCAACCCGCGCAGTGGTGTGCGCTGGAATGCTTCCAAAGCttt                                                                                                                    <  6:226423/149‑1 (MQ=255)                                 gCTGGCCGATCCTTGATGCTTTCCGAGATGCCGCCGAGCAGAGCG‑GTATTGGCAAGGTCGACGATTTCAACACCGGCGACAACCAGGGTTGTG‑GATATTTTCAGGTCAACCAGCGCAGTGGTGTGCGCTGGAATGCTTCCAAAGCTTTc                                                                                                                   <  8:115407/149‑1 (MQ=255)                                 gCTGGCCGATCCTTGATGCTTTCCGAGATGCCGCCGAGCAGAGCG‑GTATTGGCAAGGTCGACGATTTCAACACCGGCGACAACCAGGGTTGTG‑GATATATTCAGGTCAACCAGCGCAGTGGTGTGCGCTGGAATGCTTCCAAAGCTTTc                                                                                                                   <  8:282549/149‑1 (MQ=255)                                 gCTGGCCGATCCTTGATGCTTTCCGAGATGCCGCCGAGCAGAGCG‑GTATTGGCAAGGGCGACGATTTCAACACCGGCGACAACCAGGGTTGTG‑GATATTTTCAGGTCAACCAGCGCAGTGGTGTGCGCTGGAATGCTTCCAAAGCTTTc                                                                                                                   <  8:148479/149‑1 (MQ=255)                                        gATCCTTGATGCTTTCCGAGATGCCGCCGAGCAGAGCG‑GTATTGGCAAGGTCGACGATTTCAACACCGGCGACAACCAGGGTTGTG‑GATATTTTCAGGTCAACCAGCGCAGTGGTGTGCGCTGGAATGCTTCCAAAGCTTTCCTGCGGc                                                                                                            >  5:143356/1‑149 (MQ=255)                                        gATCCTTGATGCTTTCCGAGATGCCGCCGAGCAGAGCG‑GTATTGGCAAGGTCGACGATTTCAACACCGGCGACAACCAGGGTTGTG‑GATATTTTCAGGTCAACCAGCGCAGTGGTGTGCGCTGGAATGCTTCCAAAGCTTTCCTGCGGc                                                                                                            <  1:32269/149‑1 (MQ=255)                                             ttGATGCTTTCCGAGATGCCGCCGAGCAGAGCG‑GTATTGGCAAGGTCGACGATTTCAACACCGGCGACAACCAGGGTTGTG‑GATATTTTCAGGTCAACCAGCGCAGTGGTGTGCGCTGGAATGCTTCCAAAGCTTTCCTGCGGCCGAt                                                                                                        >  1:327498/1‑148 (MQ=255)                                             ttGATGCTTTCCGAGATGCCGCCGAGCAGAGAG‑GTATTGGCAAGGTCGACGATTTCAACACCGGCGACAACCAGGGTTGTG‑GATATTTTCAGGTCAACCAGAGCAGTGGTGTGCGCTGGAATGCTTCCAAAGCTTTCCTGCGGCCgct                                                                                                        >  2:181360/1‑146 (MQ=255)                                                       ccGAGATGCCGCCGAGCAGAGCG‑GTATTGGCAAGGTCGACGATTTCAACACCGGCGACAACCAGGGTTGTG‑GATATTTTCAGGTCAACCAGCGCAGTGGTGTGCGCTGGAATGCTTCCAAAGCTTTCCTGCGGCCGATCAAGGATcgc                                                                                              <  3:295782/148‑1 (MQ=255)                                                                       cAGAGCG‑GTATTGGCAAGGTCGACGATTTCAACACCGGCGACAACCAGGGTTGTG‑GATATTTTCAGGTCAACCAGCGCAGTGGTGTGCGCTGGAATGCTTCCAAAGCTTTCCTGCGGCCGATCAAGGATCGCGCCAACCTCACTGTgct                                                                             <  1:249794‑M2/149‑4 (MQ=255)                                                                           gCG‑GTATTGGCAAGGTCGACGATTTCAACACCGGCGACAACCAGGGTTGTG‑GATATTTTCAGGTCAACCAGCGCAGTGGTGTGCGCTGGAATGCTTCCAAAGCTTTCCTGCGGCCGAt                                                                                                        <  1:291088/118‑1 (MQ=255)                                                                           gCG‑GTATTGGCAAGGTCGACGATTTCAACACCGGCGACAACCAGGGTTGTG‑GATATTTTCAGGTCAACCAGCGCAGTGGTGTGCGCTGGAATGCTTCCAAAGCTTTCCTGCGGCCGAt                                                                                                        >  2:291088/1‑118 (MQ=255)                                                                               gTATTGGCAAGGTCGACGATTTCAACACCGGCGACAACCAGGGTTGTG‑GATATTTTCAGGTCAACCAGCGCAGTGGTGTGCGCTGGAATGCTTCCAAAGCTTTCCTGCGGCCGATCAAGGATCGCGCCAACCTCACTGTGCTGACCGGc                                                                      <  2:282455/149‑1 (MQ=255)                                                                                 aTTGGCAAGGTCGACGATTTCAACACCGGCGACAACCAGGGTTGTG‑GATATTTTCAGGTCAACCAGCGCAgtggt                                                                                                                                              >  5:96594/1‑75 (MQ=255)                                                                                 aTTGGCAAGGTCGACGATTTCAACACCGGCGACAACCAGGGTTGTG‑GATATTTTCAGGTCAACCAGCGCAgtggt                                                                                                                                              <  6:96594/75‑1 (MQ=255)                                                                                  ttGGCAAGGTCGACGATTTCAACACCGGCGACAACCAGGGTTGTG‑GATATTTTCAGGTCAACCAGCacagtgaggttggcgcgatccttgatcggccgcaggaaagctttggaagcattccagcgcacaccactgcgctggttgacct                                                                    <  1:217313‑M1/148‑83 (MQ=255)                                                                                    ggCAAGGTCGACGATTTCAACACCGGCGACAACCAGGGTTGTG‑GATATTTTCAGGTCAACCAGCacagtgaggttggcgcgatccttgatcggccgcaggaaagctttggaagcattccagcgcacaccactgcgc                                                                              >  5:302465‑M1/1‑64 (MQ=255)                                                                                    ggCAAGGTCGACGATTTCAACACCGGCGACAACCAGGGTTGTG‑GATATTTTCAGGTCAACCAGCacagtgaggttggcgcgatccttgatcggccgcaggaaagctttggaagcattccagcgcacaccactgcgc                                                                              <  6:302465‑M1/136‑73 (MQ=255)                                                                                     gCAAGGTCGACGATTTCAACACCGGCGACAACCAGGGTTGTG‑GATATTTTCAGGCCAACCAGCGCAGTGGTGTGCGCTGGAATGCTTCCAAAGCTTTCCTGCGGCCGATCAAGGATCGCGCCAACCTCACTGTGCTGACCGGCGTGCAg                                                                >  3:114110/1‑149 (MQ=255)                                                                                         ggTCGACGATTTCAACACCGGCGACAACCAGGGTTGTG‑GATATTTTCAGGTCAACCAGCGCAGTGGTGTGCGCTGGAATGCTTCCAAAGCTTTCCTGCGGCCGATCAAGGATCGCGCCAACCTCACTGTGCTGACCGGCGTGCAGGtt                                                             >  5:246270/1‑148 (MQ=255)                                                                                          gTCGACGATTTCAACACCGGCGACAACCAGGGTTGTG‑GATATTTTCAGGTCAACCAGCGCAGTGGTGTGCGCTGGAATGCTTCCAAAGCTTTCCTGCGGCCGATCAAGGATCGCGCCAACCTCACTGTGCTGACCGGCGTGCAGGTTGa                                                           <  3:4734/149‑1 (MQ=255)                                                                                           tccacgaTTTCAACACCGGCGACAACCAGGGTTGTG‑GATATTTTCAGGTCAACCAGCGCAGTGGTCTGCGCTGGAATGCTTCCAAAGCTTTCCTGCTGCCGATCAAGGATCGCGCCAACCTCACTGTGATGACCGGCGTGCAGGTTGAc                                                          >  3:279212/4‑149 (MQ=255)                                                                                            cgacgaTTTCAACACCGGCGACAACCAGGGTTGTG‑GATATTTTCAGGTCAACCAGCGCAGTGGTGTGCGCTGGAATGCTTCCAAAGCTTTCCTGCGGCCGATCAAGGATCGCGCCAACCTCACTGTGCTGACCGGCGTGCAGGTTGAc                                                          >  7:80859/1‑148 (MQ=255)                                                                                              acgaTTTCAACACCGGCGACAACCAGGGTTGTG‑GATATTTTCAGGTCAACCAGCGCAGTGGTGTGCGCTGGAATGCTTCCAAAGCTTTCCTGCGGCCGATCAAGGATCGCGCCAACCTCACTGTGCTGACCGGCGTGCAGGTTGAcc                                                         >  7:5655/1‑147 (MQ=255)                                                                                                 aTTTCAACACCGGCGACAACCAGGGTTGTG‑GATATTTTCAGGTCAACCAGCGCAGTGGTGTGCGCTGGAATGCTTCCAAAGCTTTCCTGCGGCCGATCAAGGATCGCGCCAACCTCACTGTGCTGACCGGCGTGCAGGTTGACCAAGTa                                                    <  1:12234/149‑1 (MQ=255)                                                                                                 aTTTCAACACCGGCGACAACCAGGGTTGTG‑GATATTTTCAGGTCAACCAGCGCAGTGGTGTGCGCTGGAATGCTTCCAAAGCTTTCCTGCGGCCGATCAAGGATCGCGCCAACCTCACTGTGCTGACCGGCGTGCAGGTTGACCAAGTa                                                    <  2:58571/149‑1 (MQ=255)                                                                                                  tttCAACACCGGCGACAACCAGGGTTGTG‑GATATTTTCAGGTCAACCAGCGCAGTGGTGTGCGCTGGAATGCTTCCAAAGCTTTCCTGCGGCCGATCAAGGATCGCGCCAACCTCACTGTGCTGACCGGCGTGCAGGTTGACCAAGTAc                                                   >  7:50964/1‑149 (MQ=255)                                                                                                   ttCAACACCGGCGACAACCAGGGTTGTG‑GATATTTTCAGGTCAACCAGCGCAGTGGTGTGCGCTGGAATGCTTCCAAAGCTTTCCTGCGGCCGATCAAGGATCGCGCCAACCTCACTGTGCTGACCGGCGTGCAGGTTGACCAAGTACt                                                  >  1:255434/1‑149 (MQ=255)                                                                                                   ttCAACACCGGCGACAACCAGGGTTGTG‑GATATTTTCAGGTCAACCAGCGCAGTGGTGTGCGCTGGAATGCTTCCAAAGCTTTCCTGCGGCCGATCAAGGATCGCGCCAACCTCACTGTGCTGACCGGCGTGCAGGTTGACCAAGTACt                                                  >  7:206891/1‑149 (MQ=255)                                                                                                      aaCACCGGCGACAACCAGGGTTGTG‑GATATTTTCAGGTCAACCAGCGCAGTGGTGTGCGCTGGAATGCTTCCAAAGCTTTCCTGCGGCCGATCAAGGATCGCGCCAACCTCACTGTGCTGACCGGCGTGCAGGTTTACCAAGTACTGCt                                               >  1:152920/1‑149 (MQ=255)                                                                                                      aaCACCGGCGACAACCAGGGTTGTG‑GATATTTTCAGGTCAACCAGCGCAGTGGTGTGCGCTGGAATGCTTCCAAAGCTTTCCTGCGGCCGATCAAGGATCGCGCCAACCTCACTGTGCTGACCGGCGTGCAGGTTGACCAAGTACTGCt                                               >  1:282645/1‑149 (MQ=255)                                                                                                         acCTGCGACAACCAGGG‑TGTGTGCTATTTTCAGGTCAACCAGCGCAGTCTGGTGAGCTGGAATGAAGCCAAATCTTTCCTGCGGCCGATCAAGGATCGCGCCAACCTCTCTGTGCTGACCTGCTTGCAGGTTGACCAAGTACTGCTCGa                                            <  8:228444/149‑1 (MQ=255)                                                                                                         acCGGCGACAACCAGGGTTGTG‑GATATTTTCAGGTCAACCAGCGCAGTGGTGTGCGCTGGAATGCTTCCAAAGCTTTCCTGCGGCCGATCAAGGATCGCGCCAACCTCACTGTGCTGACCGGCGTGCAGGTTGACCAAGTACTGCTCGa                                            <  8:50964/149‑1 (MQ=255)                                                                                                              cGACAACCAGGGTTGTG‑GATATTTTCAGGTCAACCAGCGCAGTGGTGTGCGCTGGAATGCTTCCAAAGCTTTCCTGCGGCCGATCAAGGATCGCGCCAACCTCACTGTGCTGACCGGCGTGCAGGTTGACCAAGTACTGCTCGacaaca                                       >  2:94263/1‑149 (MQ=255)                                                                                                               gACAAGGAGGGTTGTG‑GATATTTTCAGGTCAACCAGCGCAGTGGTGTGCGGTGGATTGCTTCCAAAGCTTTGCTGCGGCCGATCAAGGATCGCGCCTACCTCACTGTGCTGACCGGCGTGCAGGTTGACCAAGTACTGCTCGACAacac                                      <  4:279212/149‑1 (MQ=255)                                                                                                                 cAACCAGGGTTGTG‑GATATTTTCAGGTCAACCAGCGCAGTGGTGTGCGCTGGAATGCTTCCAAAGCTTTCCTGCGGCCGATCAAGGATCGCGCCAACCTCACTGTGCTGACCGGCGTGCAGGTTGACCAAGTACTGCTCGACAACAccc                                    <  2:282645/149‑1 (MQ=255)                                                                                                                                gATATTTTCAGGTCAACCAGCGCAGTGGTGTGCGCTGGAATGCTTCCAAAGCTTTCCTGCGGCCGATCAAGGATCGCGCCAACCTCACTGTGCTGACCGGCGTGCAGGTTGACCAAGTACTGCTCGACAACACCCGCGCAAGGGCGGt                       >  1:117957/1‑148 (MQ=255)                                                                                                                                gATATTTTCAGGTCAACCAGCGCAGTGGTGTGCGCTGGAATGCTTCCAAAGCTTTCCTGCGGCCGATCAAGGATCGCGCCAACCTCACTGTGCTGACCGGCGTGCAGGTTGACCAAGTACTGCTCGACAACACCCGCGCAAGGGCGGt                       <  8:206891/148‑1 (MQ=255)                                                                                                                                  tatTTTCAGGTCAACCAGCGCAGTGGTGTGCGCTGGAATGCTTCCAAAGCTTTCCTGCGGCCGATCAAGGATCGCGCCAACCTCACTGTGCTGACCGGCGTGCAGGTTGACCAAGTACTGCTCGACAACACCCGCGCAAGGGCGGTGaa                    <  6:143356/149‑1 (MQ=255)                                                                                                                                         aGGTCAACCAGCGCAGTGGTGTGCGCTGGAATGCTTCCAAAGCTTTCCTGCGGCCGATCAAGGATCGCGCCAACCTCACTGTgctggttgacctgaaaatatccacaaccctggttgtcgccggtgttgaaatcgtcgaccttgccaa              <  1:217313‑M2/148‑67 (MQ=255)                                                                                                                                           gTCAACCAGCGCAGTGGTGTGCGCTGGAATGCTTCCAAAGCTTTCCTGCGGCCGATCAAGGATCGCGCCAACCTCACTGTGCTGACCGGCGTGCAGGTTGACCAAGTACTGCTCGACAACACCCGCGCAAGGGCGGTGAAGGCGCTCTg           >  4:315943/1‑149 (MQ=255)                                                                                                                                            tCAACCAGCGCAGTGGTGTGCGCTGGAATGCTTCCAAAGCTTTCCTGCGGCCGATCAAGGATCGCGCCAACCTCACTGTGCTGACCGGCGTGCAGGTTGACCAAGTACTGCTCGACAACACCCGCGCAAGGGCGGTGAAGGCGCTCTgg          <  2:327498/149‑1 (MQ=255)                                                                                                                                                ccAGCGCAGTGGTGTGCGCTGGAATGCTTCCAAAGCTTTCCTGCGGCCGATCAAGGATCGCGCCAACCTCACTGTGCTGACCGGCGTGCAGGTTGACCAAGTACTGCTCGACAACACCCGCGCAAGGGCGGTGAAGGCGCTCTGGCAAg      <  7:3156/149‑1 (MQ=255)                                                                                                                                                ccAGCGCAGTGGTGTGCGCTGGAATGCTTCCAAAGCTTTCCTGCGGCCGATCAAGGATCGCGCCAACCTCACTGTGCTGACCGGCGTGCAGGTTGACCAAGTACTGCTCGACAACACCCGCGCAAGGGCGGTGAAGGCGCTCTGGCAAg      >  3:72270/1‑149 (MQ=255)                                                                                                                                                  aGCacagtgaggttggcgcgatccttgatcggccgcaggaaagctttggaagcattccagcgcacaccactgcgctggttgacctgaaaatatccacaaccctggttgtcgccggtgttgaaatcgtcgaccttgccaataccgctctg    <  1:249794‑M1/149‑147 (MQ=255)                                                                                                                                                  aGCGCAGTGGTGTGCGCTGGAATGCTTCCAAAGCTTTCCTGCGGCCGATCAAGGATCGCGCCAACCTCACTGTGCTGACCGGCGTGCAGGTTGACCAAGTACTGCTCGACAACACCCGCGCAAGGGCGGTGAAGGCGCTCTGGCAAggg    <  8:5655/149‑1 (MQ=255)                                                                                                                                                   gcgcAGTGGTGTGCGCTGGAATGCTTCCAAAGCTTTCCTGCGGCCGATCAAGGATCGCGCCAACCTCACTGTgctggttgacctgaaaatatccacaaccctggttgtcgccggtgttgaaatcgtcgaccttgcc                <  6:302465‑M2/136‑65 (MQ=255)                                                                                                                                                   gcgcAGTGGTGTGCGCTGGAATGCTTCCAAAGCTTTCCTGCGGCCGATCAAGGATCGCGCCAACCTCACTGTgctggttgacctgaaaatatccacaaccctggttgtcgccggtgttgaaatcgtcgaccttgcc                >  5:302465‑M2/1‑72 (MQ=255)                                                                                                                                                    cgcAGTGGTGTGCGCTGGAATGCTTCCAAAGCTTTCCTGCGGCCGATCAAGGATCGCGCCAACCTCACTGTGCTGACCGGCGTGCAGGTTGACCAAGTACTGCTCGACAACACCCGCGCAAGGGCGGTGAAGGCGCTCTGGCAAGGGGc  <  8:80859/149‑1 (MQ=255)                                                                                                                                                    |                                                                                                                                                     GGGCGAATGGCGGGTCGAGCGGCAACGCTACAGCTGGCCGATCCTTGATGCTTTCCGAGATGCCGCCGAGCAGAGCG‑GTATTGGCAAGGTCGACGATTTCAACACCGGCGACAACCAGGGTTGTG‑GATATTTTCAGGTCAACCAGCGCAGTGGTGTGCGCTGGAATGCTTCCAAAGCTTTCCTGCGGCCGATCAAGGATCGCGCCAACCTCACTGTGCTGACCGGCGTGCAGGTTGACCAAGTACTGCTCGACAACACCCGCGCAAGGGCGGTGAAGGCGCTCTGGCAAGGGGC  >  NC_002947/65380‑65673

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.