New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913	4226588 =	52 (0.990)	7 (0.130)	6/280	6.1	11.8%	coding (2761/3684 nt)	metH	homocysteine‑N5‑methyltetrahydrofolate transmethylase, B12‑dependent
	?	NC_000913	4226623 =	53 (1.010)					coding (2796/3684 nt)	metH	homocysteine‑N5‑methyltetrahydrofolate transmethylase, B12‑dependent
Rejected: Coverage evenness skew score above cutoff.
Rejected: Frequency below/above cutoff threshold.

CCGCAAGGAGTACGAAACCGTACGTATTCAGCACGGGCGCAAGAAACCGCGCACACCACCGGTCACGCTGGAAGCGGCGCGCGATAACGATTTCGCTTTTGACTGGCAGGCTTACACGCCGCCGGTGGCGCACCGTCTCGGCGTGCAGGAAGTCGAAGCCAGCATCGAAACGCTGCGTAATTACATCGACTGGACACCGTTCTTTATGACCTGGTCGCTGGCCGGGAAGTATCCGCGCATTCTGGAAGATGAAGTGGTGGGCGTTGAGGCGCAGCGGCTGTTTAAAGACGCCAACG  >  NC_000913/4226476‑4226771                                                                                                                                                    |                                                                                                                                                     ccGCAAGGAGTACGAAACCGTACGTATTCAGCACGGGCGCAAGAAACCGCGCACACCACCGGTCACGCTGGAAGCGGCGCGCGATAACGATTTCGCTTTTGACTGGCAGGCTTACACGCCGCCGGTGGCGCACCGTCTCGGCGTGCAgg                                                                                                                                                     >  1:262350/1‑149 (MQ=255) ccGCAAGGAGTACGAAACCGTACGTATTCAGCACGGGCGCAAGAAACCGCGCACACCACCGGTCACGCTGGAAGCGGCGCGCGATAACGATTTCGCTTTTGACTGGCAGGCTTACACGCCGCCGGTGGAGCACCGTCTCGGAGTGAAgg                                                                                                                                                     >  1:34391/1‑149 (MQ=255)   gCAAGGAGTACGAAACCGTACGTATTCAGCACGGGCGCAAGAAACCGCGCACACCACCGGTCACGCTGGAAGCGGCGCGCGATAACGATTTCGCTTTTGACTGGCAGGCTTACACGCCGCCGGTGGCGCACCGTCTCGGCGTGCAGGaa                                                                                                                                                   >  2:596698/1‑149 (MQ=255)           tACGAAACCGTACGTATTCAGCACGGGCGCAAGAAACCGCGCACACCACCGGTCACGCTGGAAGCGGCGCGCGATAACGATTTCGCTTTTGACTGGCAGGCTTACACGCCGCCGGTGGCGCACCGTCTCGGCGTGCAGGAAGTCGAAGc                                                                                                                                           >  2:627063/1‑149 (MQ=255)             cGAAACTGGACGTAGTCAGCCTGTGCGCAAGGAACCGCGCACACCACCGGTCACGCTGGAAGCGGCGCGCGATAACGACTTCGCTTTTGACGGGCAGGCTTACACGCCGCCGGTGGCGCACCGTCTCGGCGTGCAGGAAGTCGAAGCCa                                                                                                                                         <  2:95050/149‑1 (MQ=255)             cGAAACCGTACGTATTCAGCACGGGCGCAAGAAACCGCGCACACCACCGGTCACGCTGGAAGCGGCGCGCGATAACGATTTCGCTTTTGACTGGCAGGCTTACACGCCGCCGGTGGCGCACCGTCTCGGCGTGCAGGAAGTCGAAGCCa                                                                                                                                         >  1:471533/1‑149 (MQ=255)              gAAACCGTACGTATTCAGCACGGGCGCAAGAAACCGCGCACACCACCGGTCACGCTGGAAGCGGCGCGCGATAACGATTTCGCTTTTGACTGGCAGGCGTACACGCCGCCGGTGGCGCACCGTCTCGGCGTGCAGGAAGTCGAAGCCAg                                                                                                                                        >  1:257056/1‑149 (MQ=255)                aaCCGTACGTATTCAGCACGGGCGCAAGAAACCGCGCACACCACCGGTCACGCTGGAAGCGGCGCGCGATAACGATTTCGCTTTTGACTGGCAGGCTTACACGCCGCCGGTGGCGCACCGTCTCGGCGTGCAGGAAGTCGAAGCCAGCa                                                                                                                                      >  2:468238/1‑149 (MQ=255)                    gtacgtaTTCAGCACGGGCGCAAGAAACCGCGCACACCACCGGTCACGCTGGAAGCGGCGCGCGATAACGATTTCGCTTTTGACTGGCAGGCTTACACGCCGCCGGTGGCGCACCGTCTCGGCGTGCAGGAAGTCGAAGCCAGCATCGa                                                                                                                                  >  1:738916/1‑149 (MQ=255)                          aTTCAGAACGGGCGCAAGAAACCGCGCACACCACCGGTCACGCTGGAAGTGGCGCGCGATAACGATTTCGCTTTTGACTGGCAGGCTTACACGCCGCCGGTGGCGCAACGTCTCGGCGTGCAGGAAGTCGAAGCCAGCATCGAAACGCt                                                                                                                            >  2:183837/1‑149 (MQ=255)                           ttCAGCACGGGCGCAAGAAACCGCGCACACCACCGGTCACGCTGGAAGCGGCGCGCGATAACGATTTCGCTTTTGACTGGCAGGCTTACACGCCGCCGGTGGCGCACCGTCTCGGCGTGCAGGAAGTCGAAGCCAGCATCGAAACGCTg                                                                                                                           >  1:884953/1‑149 (MQ=255)                           ttCAGCACGGGCGCAAGAAACCGCGCACACCACCGGTCACGCTGGAAGCGGCGCGCGATAACGATTTCGCTTTTGACTGGCAGGCTTACACGCCGCCGGTGGCGCACCGTCTCGGCGTGCAGGAAGTCGAAGCCAGCATCGAAACGCTg                                                                                                                           <  1:780736/149‑1 (MQ=255)                               gCACGGGCGCAAGAAACCGCGCACACCACCGGTCACGCTGGAAGCGGCGCGCGATAACGATTTCGCTTTTGACTGGCAGGCTTACACGCCGCCGGTGGCGCACCGTCTCGGCGTGCAGGGAGTCGAAGCCAGCATCGAAACGCTGCGTa                                                                                                                       >  1:294824/1‑149 (MQ=255)                               gCACGGGCGCAAGAAACCGCGCACACCACCGGTCACGCTGGAAGCGGCGCGCGATAACGATTTCGCTTTTGACTGGCAGGCTTACACGCCGCCGGTGGCGCACCGTCTCGGCGTGCAGGAAGTCGAAGCCAGCATCGAAACGCTGCGTa                                                                                                                       >  2:556736/1‑149 (MQ=255)                                cACGGGCGCAAGAAACCGCGCACACCACCGGTCACGCTGGAAGCGGCGCGCGATAACGATTTCGCTTTTGACTGGCAGGCTTACACGCCGCCGGTGGCGCACCGTCTCGGCGTGCAGGAAGTCGAAGCCAGCATCGAAACGCTGCGTaa                                                                                                                      <  1:556736/149‑1 (MQ=255)                                 aCGGGCGCAAGAAACCGCGCACACCACCGGTCACGCTGGAAGCGGAGCGCGATAACGATCTGGCTTGGGACTGGCAGGCTTACACGCCGCCGGGGGCGCTCCGGCTCGGCGTGCAGGAAGTCGCAGCCAGCATCGAGACGCTGCGTAAt                                                                                                                     >  1:28184/1‑149 (MQ=255)                                  cGGGCGCAAGAAACCGCGCACACCACCGGTCACGCTGGAAGCGGCGCGCGATAACGATTTCGCTTTTGACTGGCAGGCTTACACGCCGCCGGTGGCGCACCGTCTCGGCGTGCAGGAAGTCGAAGCCAGCATCGAAACGCTGCGTAAtt                                                                                                                    <  2:539560/149‑1 (MQ=255)                                                  cgcACACCACCGGTCACGCTGGAAGCGGCGCGCGATAACGATTTCGCTTTTGACTGGCAGGCTTACACGCCGCCGGTGGCGCACCGTCTCGGCGTGCAGGAAGTCGAAGCCAGCATCGAAACGCTGCGTaa                                                                                                                      <  2:625946/131‑1 (MQ=255)                                                  cgcACACCACCGGTCACGCTGGAAGCGGCGCGCGATAACGATTTCGCTTTTGACTGGCAGGCTTACACGCCGCCGGTGGCGCACCGTCTCGGCGTGCAGGAAGTCGAAGCCAGCATCGAAACGCTGCGTaa                                                                                                                      >  1:625946/1‑131 (MQ=255)                                                       accaccGGTCACGCTGGAAGCGGCGCGCGATAACGATTTCGCTTTTGACTGGCAGGCTTACACGCCGCCGGTGGCGCACCGTCTCGGCGTGCAGGAAGTCGAAGCCAGCATCGAAACGCTGCGTAATTACATCGACTGGACACCGTTCt                                                                                               >  2:130927/1‑149 (MQ=255)                                                         caccGGTCACGCTGGAAGCGGCGCGCGATAACGATTTCGCTTTTGACTGGCAGGCTTACACGCCGCCGGTGGCGCACCGTCTCGGCGTGCAGGAAGTCGAAGCCAGCATTGAAACGCTGCGTAATTACATCGACTTGACACCGTTCttt                                                                                             >  2:780942/1‑149 (MQ=255)                                                         caccGGTCACGCTGGAAGCGGCGCGCGATAACGATTAGTCGTTTTACTGGCAGGCTTACACGCCGCCGGTGGCGCACCGTCTCGGCGTGCAGGAAGTCGAAGCCAGCATCGAAACGCTGCGTAATTACATCGACTGGACACCGTTCttt                                                                                             >  2:617701/1‑149 (MQ=255)                                                           ccGGTCACGCTGGAAGCGGCGCGCGATAACGATTTCGCTTTTGACTGGCAGGCTTACACGCCGCCGGTGGCGCACCGTCTCGGCGTGCAGGAAGTCGAAGCCAGCATCGAAACGCTGCGTAATTACATCGACTGGACACCGTTCTTTAt                                                                                           <  2:738916/149‑1 (MQ=255)                                                            cGGTCACGCTGGAAGCGGCGCGCGATAACGATTTCGCCTTTGACTGGCAGGCTTACACGCCGCCGGTGGCGCACCGTCTCGGCGTGCAGGAAGTCGAAGCCAGCATCGAAACGCTGCGTAATTACATCGACTGGACACCGTTCTTTATg                                                                                          <  2:262350/149‑1 (MQ=255)                                                                cACGCTGGAAGCGGCGCGCGATAACGATTTCGCTTTTGACTGGCAGGCTTACACGCCGCCGGTGGCGCACCGTCTCGGCGTGCAGGAAGTCGAAGCCAGCATCGAAACGCTGCGTAATTACATCGACTGGACACCGTTCTTTATGACCt                                                                                      <  1:21127/149‑1 (MQ=255)                                                                     tGGAAGCGGCGCGCGATAACGATTTCGCTTTTGACTGGCAGGCTTACACGCCGCCGGTGGCGCACCGTCTCGGCGTGCAGGAAGTCGAAGCCAGCATCGAAACGCTGCGTAATTACATCGACTGGACACCGTTCTATATGACCTGGTCg                                                                                 >  1:5154/1‑149 (MQ=255)                                                                                cgcgATAACGATTTCGCTTTTGACTGGCAGGCTTACACGCCGCCGGTGGCGCACCGTCTCGGCGTGCAGGAAGTCGAAGCCAGCATCGAAACGCTGCGTAATTACATCGACTGGACACCGTTCTTTATGACCTGGTCGCTGGCCGGGaa                                                                      <  2:665081/149‑1 (MQ=255)                                                                                    aTAACGATTTCGCTTTTGACTGGCAGGCTTACACGCCGCCGGTGGCGCACCGTCTCGGCGTGCAGGAAGTCGAAGCCAGCATCGAAACGCTGCGTAATTACATCGACTGGACACCGTTCTTTATGACCTGGTCGCTGGCCGGGCAGTAt                                                                  >  1:749442/1‑149 (MQ=255)                                                                                      aaCGATTTCGCTTTTGACTGGCAGGCTTACACGCCGCCGGTGGCGCACCGTCTCGGCGTGCAGGAAGTCGAAGCCAGCATCGAAACGCTGCGTAATTACATCGACTGGACACCg                                                                                                   >  2:687255/1‑114 (MQ=255)                                                                                      aaCGATTTCGCTTTTGACTGGCAGGCTTACACGCCGCCGGTGGCGCACCGTCTCGGCGTGCAGGAAGTCGAAGCCAGCATCGAAACGCTGCGTAATTACATCGACTGGACACCg                                                                                                   <  1:687255/114‑1 (MQ=255)                                                                                           tttCGCTTTTGACTGGCAGGCTTACACGCCGCCGGTGGCGCACCGTCTCGGCGTGCAGGAAGTCGAAGCCAGCATCGAAACGCTGCGTAATTACATCGACTGGACACCGTTCTTTATGACCTGGTCGCTGGCCGGGAAGTATCCGCGCa                                                           >  1:600130/1‑149 (MQ=255)                                                                                                 ttttGACTGGCAGGCTTACAGGGCGCCGGTGGCGCACCGTCTCGGCGTGCAGGAAGTCGAAGCCAGCATCGAAACGCTGCGTAATTACATCGACTGGACACCGTTCTTTATGACCTGGTCGCTGGCCGGGAAGTATCCGCGCATTCTgg                                                     <  2:522342/149‑1 (MQ=255)                                                                                                 ttttGACTGGCAGGCTTACACGCCGCCGGTGGCGCACCGTCTCGGCGTGCAGGAAGTCGAAGCCAGCGTCGAAACGCTGCGTAATTACATCGACTGGACACCGTTCTTTATGACCTGGTCGCTGGCCGGGAAGTGTCCGGCCATTCTgg                                                     >  2:67912/1‑149 (MQ=255)                                                                                                 ttttGACTGGCAGGCTTACACGCCGCCGGTGGCGCACCGTCTCGGCGTGCAGGAAGTCGAAGCCAGCATCGAAACGCTGCGTAATTACATCGACTGGACACCGTTCTTTATGACCTGGTCGCTGGCCGGGAAGTATCCGCGCATTCTgg                                                     >  1:126461/1‑149 (MQ=255)                                                                                                     gACTGGCAGGCTTACACGCCGCCGGTGGCGCACCGTCTCGGCGTGCAGGAAGTCGAAGCCAGCATCGAAACGCTGCGTAATTACATCGACTGGACACCGTTCTTTATGACCTGGTCGCTGGCCGGGAAGTATCCGCGCATTCTGGAAGa                                                 <  2:126461/149‑1 (MQ=255)                                                                                                     gACTGGCAGGCTTACACGCCGCCGGTGGCGCACCGTCTCGGCGTGCAGGAAGTCGAAGCCAGCATCGAAACGCTGCGTAATTACATCGACTGGACACCGTTCTTTATGACCTGGTCGCTGGCCGGGAAGTATCCGCGCACTCCGGAAGa                                                 >  2:716932/1‑149 (MQ=255)                                                                                                      ggtGGCCGGCTTACACGCCGCCGGTGGGGCACCGTCTCGGCGTGAAGGAAGTCGAAGCCAGCATCGAAACGCTGCGTAATTACATCGACTGGACACCGTTCTTTAGGACCTGGTCGCTGGCCGGGAAGTATCCGCGCATTCTGGAAGAt                                                <  2:460112/147‑1 (MQ=255)                                                                                                         ggCACGCTGACTCGCCGCCGGTGGCGCACCGTCTCGGCGTGCAGGAAGTCGAAGCCAGCATCGAACCGCTGCGTAATTACATCGACTGGACACCGTTCTTTATGACCTGGTCGCTGGCCGGGAAGTATCCGCGCATTCTGGAAGATGaa                                             <  2:637917/149‑1 (MQ=255)                                                                                                              gCTTACACGCCGCCGGTGGCGCACCGTCTCGGCGTGCAGGAAGTCGAAGCCAGCATCGAAACGCTGCGTAATTACATCGACTGGACACCGTTCTTTATGACCTGGTCGCTGGCCGGGAAGTATCCGCGCATTCTGGAAGATGAAgtggt                                        <  1:468238/149‑1 (MQ=255)                                                                                                              gCTTACACGCCGCCGGTGGCGCACCGTCTCGGCGTGCAGGAAGTCGAAGCCAGCATCGAAACGCTGCGTAATTACATCGACTGGACACCGTTCTTTATGACCTGGTCGCTGGCCGGGAAGTATCCGCGCATTCTGGAAGATGAAgtggt                                        <  1:777535/149‑1 (MQ=255)                                                                                                                ttATACGCCGCCGGTGGCGCACCGTCTCGGCGTGCAGGAAGTCGAAGCCAGCATCGAAACGCTGCGTAATTACATCGACTGGACACCGTTCTTTATGACCTGGTCGCTGGCCGGGAAGTATCCGCGCATTCTGGAAGATGAAGtggtgg                                      <  1:68191/149‑1 (MQ=255)                                                                                                                    tctccgccgGTGGCGCACGGGGCCGGCGTGAAGTAAGTCGAAGCCAGCATCGAAACGCGGCGTAATTACATCGACTGGACACCGTTCTTTATGACCTGGTCGCTGGCCGGGAAGTATCCGCGCATTCTGGAAGATGAAGTGGTGGGCGt                                  <  2:227241/146‑1 (MQ=255)                                                                                                                    acGCCGCCGGTGGCGCACCGTCTCGGCGTGCAGGAAGTCGAAGCCAGCATCGAAACGCTGCGTAATTACATCGACTGGACACCGTTCTTTATGACCTGGTCGCTGGCCGGGAAGTATCCGCGCATTCTGGAAGATGAAGTGGTGGGCGt                                  >  1:610763/1‑149 (MQ=255)                                                                                                                    acGCCGCCGGTGGCGCACCGTCTCGGCGTGCAGGAAGTCGAAGCCAGCATCGAAACGCTGCGTAATTACATCGACTGGACACCGTTCTTTATGACCTGGTCGCTGGCCGGGAAGTATCCGCGCATTCTGGAAGATGAAGTGGTGGGCGt                                  >  1:756870/1‑149 (MQ=255)                                                                                                                     cgccgcCGGTGGCGCACCGTCTCGGCGTGCAGTAAGTCGAAGCCAGCATCGAAACGCTGCGTAATTACATCGACTGGACACCGTTCTTTATGACCTGGTCGCTGGCCGGGAAGTATCCGCGCATTCTGGAAGATGAAGTGGTGGGCGtt                                 <  2:257056/149‑1 (MQ=255)                                                                                                                                   cACCGTCTCGGCGTGCAGGAAGTCGCAGCCAGCATCGAAACGCTGCGTAATTACATCGACTGGACACCGTTCTTTATGACCTGGTCGCTGGCCGGGAAGTATTCGCGCATTCTGGAAGATGAAGTGGTGGGCGTTGAGGCGCAGCGGCt                   <  1:827473/149‑1 (MQ=255)                                                                                                                                     ccGTCTCGGCGTGCAGGAAGTCGAAGCCAGCATCGAAACGCTGCGTAATTACATCGACTGGACACCGTTCTTTATGACCTGGTCGCTGGCCGGGAAGTATCCGCGCATTCTGGAAGATGAAGTGGTGGGCGTTGAGGCGCAGCGGCTGt                 <  1:63273/149‑1 (MQ=255)                                                                                                                                       gTCACGGAGTGTAGGAAGTCGAAGCCAGCATCGAAAAGCTGCGTAATTGCAGGGACTGGACACCGTTCGTTATGCCCTGGTCGCTGGCCGGGAAGTTTGCGCGCATTGTGGAAGATGATGTGGTGGTCGTTGAGGCGCAGCGGCTGttt               <  2:633064/149‑1 (MQ=255)                                                                                                                                        tctcGGCGTGCAGGAAGTCGAAGCCAGCCTCGCAACGCCGCGTAATTACCTCGACTGGACACCGTTCTTTTTGACCTGGTCGCTGGCCGGGAAGTATCCGCGCATTCTGGAAGATGAAGTGGTGGGCGTTGAGGCGCAGCGGCTGTTTa              <  1:823364/149‑1 (MQ=255)                                                                                                                                            ggCGTGCAGGAAGTCGAAGCCAGCATCGAAACGCTGCGTAATTACATCGACTGGACACCGTTCTTTATGACCTGGTCGCTGGCCGGGAAGTATCCGCGCATTCTGGaa                                                   >  1:444886/1‑108 (MQ=255)                                                                                                                                            ggCGTGCAGGAAGTCGAAGCCAGCATCGAAACGCTGCGTAATTACATCGACTGGACACCGTTCTTTATGACCTGGTCGCTGGCCGGGAAGTATCCGCGCATTCTGGaa                                                   <  2:444886/108‑1 (MQ=255)                                                                                                                                             gagTGCCGGCAGTAGCAGCCAGCATCGCACCGCTGCGTACTTACATCGGCTGGCCCCCGTCCTTTATGACCTGGTCGCTGGCCGGGAAGTATCCGCGCATTCTGGAAGATGAAGTGGTGGGCGTTGAGGCGCAGCGGCTGTTTAAAGAc         <  1:877638/147‑1 (MQ=255)                                                                                                                                             gCGTGAAGGAAGTCGAAGCCAGCATCGAAACGCTGCGTAATTACATCGACTGGACACCGTTCTTTATGACCTGGTCGCTGGCCGGGAAGTATCCGCGCATTCTGGAAGATGAAGTGGTGGGCGTTGAGGCGCAGCGGCTGTTTAAAGAc         <  2:290055/149‑1 (MQ=255)                                                                                                                                                  cAGGAAGTCGAAGCCAGCATCGAAACGCTGCGTAATTACATCGACTGGACACCGTTCTTTATGACCTGGTCGCTGGCCGGGAAGTATCCGCGCATTCTGGAAGATGAAGTGGTGGGCGTTGAGGCGCAGCGGCTGTTTAAAGACGCCaa    >  2:275938/1‑149 (MQ=255)                                                                                                                                                    ggAAGTCGAAGCCAGCATCGAAACGCTGCGTAATTACATCGACTGGACACCGTTCTTTATGACCTGGTAGCTGGCCGGGAAGTATCCGCGCATTCTGGAAGATGCAGTGGTGGGCGTTGAGGCGCAGCGGCTGTTTAAAGACGCCAACg  >  2:224841/1‑149 (MQ=255)                                                                                                                                                    |                                                                                                                                                     CCGCAAGGAGTACGAAACCGTACGTATTCAGCACGGGCGCAAGAAACCGCGCACACCACCGGTCACGCTGGAAGCGGCGCGCGATAACGATTTCGCTTTTGACTGGCAGGCTTACACGCCGCCGGTGGCGCACCGTCTCGGCGTGCAGGAAGTCGAAGCCAGCATCGAAACGCTGCGTAATTACATCGACTGGACACCGTTCTTTATGACCTGGTCGCTGGCCGGGAAGTATCCGCGCATTCTGGAAGATGAAGTGGTGGGCGTTGAGGCGCAGCGGCTGTTTAAAGACGCCAACG  >  NC_000913/4226476‑4226771

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