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.

CCGATACCCAGCGTGATGATTTTGTCGCTCGTACCCGCAAGGAGTACGAAACCGTACGTATTCAGCACGGGCGCAAGAAACCGCGCACACCACCGGTCACGCTGGAAGCGGCGCGCGATAACGATTTCGCTTTTGACTGGCAGGCTTACACGCCGCCGGTGGCGCACCGTCTCGGCGTGCAGGAAGTCGAAGCCAGCATCGAAACGCTGCGTAATTACATCGACTGGACACCGTTCTTTATGACCTGGTCGCTGGCCGGGAAGTATCCGCGCATTCTGGAAGATGAAGTGGTGG  >  NC_000913/4226442‑4226735                                                                                                                                                   |                                                                                                                                                    ccGATACCCAGCGTGATGATTTTGTCGCTCGTACCCGCAAGGAGTACGAAACCGTACGTATTCAGCACGGGCGCAAGAAACCGCGCACACCACCGGTCACGCTGGAAGCGGCGCGCGATAACGATTTCGCTTTTGACTGGCAGGCTTac                                                                                                                                                   <  1:757037/149‑1 (MQ=255)  cGATACCCAGCGTGATGATTTTGTCGCTCGTACCCGCAAGGAGTACGAAACCGTACGTATTCAGCACGGGCGCAAGAAACCGCGCACACCACCGGTCACGCTGGAAGCGGCGCGCGATAACGATTTCGCTTTTGACTGGCAGGCTTaca                                                                                                                                                  <  2:795058/149‑1 (MQ=255)   gATACCCAGCGTGATGATTTTGTCGCTCGTACCCGCAAGGAGTACGAAACCGTACGTATTCAGCACGGGCGCAAGAAACCGCGCACACCACCGGTCACGCTGGAAGCGGCGCGCGATAACGATTTCGCTTTTGACTGGCAGGCTTacac                                                                                                                                                 >  2:797161/1‑149 (MQ=255)       cccAGCGTGATGATTTTGTCGCTCGTACCCGCAAGGAGTACGAAACCGTACGTATTCAGCACGGGCGCAAGAAACCGCGCACACCACCGGTCACGCTGGAAGCGGCGCGCGATAACGATTTCGCTTTTGACTGCCAGGCTTACAcgccg                                                                                                                                             >  2:201955/1‑149 (MQ=255)          aGCGTGATGATTTTGTCGCTGGTACCCGCAAGGAGTATGGAACGCTGCGGATGCAGCACGGGCGCAAGAAACCGCGCACACCACCGGTCACGCTGGAAGCGGCGCGCGATAACGATTTCGCTTTTGACTGGCAGGCTTACACGccgccg                                                                                                                                          <  2:329015/149‑1 (MQ=255)          aGCGTGATGATTTTGTCGCTCGTACCCGCAAGGAGTACGAAACCGTACGTATTCAGCACGGGCGCAAGAAACCGCGCACACCACCGGTCACGCTGGAAGCGGCGCGCGATAACGATTTCGCTTTTGACTGGCAGGCTTACACGccgccg                                                                                                                                          >  1:539560/1‑149 (MQ=255)           gCGTGATGATTTTGTCGCTCGTACCCGCAAGGAGTACGAAACCGTACGTATTCAGCAAGGGCGCAAGAAACCGCGCACACCACCGGTCACGCTGGAAGCGGCGCGCGATAACGATTTCGCTTTTGACTGGCCGGCCTCCACgccccccg                                                                                                                                         >  2:282108/1‑147 (MQ=255)             gTGATGATTTTGTCGCTCGTACCCGCAAGGAGTACGAAACCGTACGTATTCAGCACGGGCGCAAGAAACCGCGCACACCACCTGTCACGCTGGAAGCGGCGCGCGATAACGATTTCGCTTTTGACTGGCAGGCTTACACGCCGCCGGTg                                                                                                                                       <  2:49578/149‑1 (MQ=255)                  gatTTTGTCGCTCGTACCCGCAAGGAGTACGAAACCGTACGTATTCAGCACGGGCGCAAGAAACCGCGCACACCACCGGTCACGCTGGAAGCGGCGCGCGATAACGATTTCGCTTTTCACTGCCAGGCTTACACGCCGCCGGTGGCGCa                                                                                                                                  >  2:850317/1‑149 (MQ=255)                    ttttGTCGCGCGTGCCCGCAAGGAGTACGAAACCGTACGTATTCAGCACGGGCGCAAGAAACCGCGCACACCACCGGTCACGCTGGAAGCGGCGCGCGATAACGATTTCGCTTTTGACTGGCAGGCTTACACGCCGCCGGTGGCGCAcc                                                                                                                                <  1:76567/149‑1 (MQ=255)                                tACCCGCAAGGAGTACGAAACCGTACGTATTCAGCACGGGCGCAAGAAACCGCGCACACCACCGGTCACGCTGGAAGCGGCGCGCGATAACGATTTCGCTTTTGACTGGCAGGCTTACACGCCGCCGGTGGCGCACCGTCTCGGCGTGc                                                                                                                    <  1:797161/149‑1 (MQ=255)                                   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:780736/149‑1 (MQ=255)                                                             ttCAGCACGGGCGCAAGAAACCGCGCACACCACCGGTCACGCTGGAAGCGGCGCGCGATAACGATTTCGCTTTTGACTGGCAGGCTTACACGCCGCCGGTGGCGCACCGTCTCGGCGTGCAGGAAGTCGAAGCCAGCATCGAAACGCTg                                                                                       >  1:884953/1‑149 (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)                                                                                                                                                   |                                                                                                                                                    CCGATACCCAGCGTGATGATTTTGTCGCTCGTACCCGCAAGGAGTACGAAACCGTACGTATTCAGCACGGGCGCAAGAAACCGCGCACACCACCGGTCACGCTGGAAGCGGCGCGCGATAACGATTTCGCTTTTGACTGGCAGGCTTACACGCCGCCGGTGGCGCACCGTCTCGGCGTGCAGGAAGTCGAAGCCAGCATCGAAACGCTGCGTAATTACATCGACTGGACACCGTTCTTTATGACCTGGTCGCTGGCCGGGAAGTATCCGCGCATTCTGGAAGATGAAGTGGTGG  >  NC_000913/4226442‑4226735

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