New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913 4227233 =19 (0.940)3 (0.160) 3/280 2.1 20.1% coding (3406/3684 nt) metH homocysteine‑N5‑methyltetrahydrofolate transmethylase, B12‑dependent
?NC_000913 4227352 = 6 (0.310)coding (3525/3684 nt) metH homocysteine‑N5‑methyltetrahydrofolate transmethylase, B12‑dependent

ATGGTGAAAGCGCTTGCCGACCGTTTAGCCGAAGCCTTTGCGGAGTATCTCCATGAGCGTGTGCGTAAAGTCTACTGGGGCTATGCGCCGAACGAGAACCTCAGCAACGAAGAGCTGATCCGCGAAAACTACCAGGGCATCCGTCCGGCACCGGGCTATCCGGCCTGCCCGGAACATACGGAAAAAGCCACCATCTGGGAGCTGCTGGAAGTGGAAAAACACACTGGCATGAAACTCACAGAATCTTTCGCCATGTGGCCCGGTGCATCGGTTTCGGGTTGGTA  >  NC_000913/4227086‑4227369
                                                                                                                                                   |                                                                                                                                        
taccccagattcggtgccagccagcgctcaacttcggtaacgctcatacctttacggcgggcataatcttcaacctgatcgcgctgaatttgtgctacagcgtagtacttgctgtccgggtggctgaagtaccaacccgaaaccgatGCAc                                                                                                                                       >  2:77152‑M1/148‑151 (MQ=255)
aTGGTGAAAGCGCTTGCCGACCGTTTAGCCGAAGCCTTTGCGGAGTATCTCCATGAGCGTGTGCGTAAAGTCTACTGGGGCTATGCGCCGAACGAGAACCTCAGCAACGAAGAGCTGATCCGCGAAAACTACCAGGGCATCCGTCCGGCAc                                                                                                                                       <  2:39401/151‑1 (MQ=255)
  ccccagattcggtgtcagccagcgctcaacttcggtaacgctcatacctttacggcgggcctaattttcaacctgatcgcgctgaatttgtgctacagcggagtacttgctgttcgggtgggtgaagtaacaaaccgaaaacgaaGCAcc                                                                                                                                      <  1:133569‑M1/5‑1 (MQ=255)
                        ttAGCCGAAGCCTTTGCGGAGTATCTCCATGAGCGTGTGCGTAAAGTCTACTGGGGCTATGCGCCGAACGAGAACCTCAGCAACGAAGAGCTGATCCGCGAAAACTACCAGGGCATCCGTCCGGCACCGGGCTATCCGGCCTGCCCGGAAc                                                                                                               <  1:70637/151‑1 (MQ=255)
                                   cTTTGCGGAGTATCTCCATGAGCGTGTGCGTAAAGTCTACTGGGGCTATGCGCCGAACGAGAACCTCAGCAACGAAGAGCTGATCCGCGAAAACTACCAGGGCATCCGTCCGGCACCGGGCTATCCGGCCTGCCCGGAACATACGGaaaaa                                                                                                    >  1:276608/1‑151 (MQ=255)
                                            gTATCTCCATGAGCGTGTGCGTAAAGTCTACTGGGGCTATGCGCCGAACGAGAACCTCAGCAACGAAGAGCTGATCCGCGAAAACTACCAGGGCATCCGTCCGGCACCGGGCTATCCGGCCTGCCCGGAACATACGGAAAAAGCCACCAt                                                                                            >  2:304252/1‑150 (MQ=255)
                                             tATCTCCATGAGCGTGTGCGTAAAGTCTACTGGGGCTATGCGCCGAACGAGAACCTCAGCAACGAAGAGCTGATCCGCGAAAACTACCAGGGCATCCGTCCGGCACCGGGCTATCCGGCCTGCCCGGAACATACGGAAAAAGCCACCATc                                                                                           <  2:36002/150‑1 (MQ=255)
                                              aTCTCCATGAGCGTGTGCGTAAAGTCTACTGGGGCTATGCGCCGAACGAGAACCTCAGCAACGAAGAGCTGATCCGCGAAAACTACCAGGGCATCCGTCCGGCACCGGGCTATCCGGCCTGCCCGGAACATACGGAAAAAGCCACCAt                                                                                            >  1:75679/1‑148 (MQ=255)
                                                cctttacggcgggcataatcttcaacctgatcgcgctgaatttgtgctacagcgtagtacttgctgtccgggtggctgaagtaccaacccgaaaccgatGCACCGGGCTATCCGGCCTGCCCGGAACATACGGAAAAAGCCACCATCTggg                                                                                       >  2:38880‑M1/100‑151 (MQ=255)
                                                 tcCATGAGCGTGTGCGTAAAGTCTACTGGGGCTATGCGCCGAACGAGAACCTCAGCAACGAAGAGCTGATCCGCGAAAACTACCAGGGCATCCGTCCGGCACCGGGCTATCCGGCCTGCCCGGAACATACGGAAAAAGCCACCATCTGGGa                                                                                      <  2:101777/151‑1 (MQ=255)
                                                  ccATGAGCGTGTGCGTAAAGTCTACTGGGGCTATGCGCCGAACGAGAACCTCAGCAACGAAGAGCTGATCCGCGAAAACTACCAGGGCATCCGTCCGGCACCGGGCTATCCGGCCTGCCCGGAACATACGGAAAAAGCCACCATCTGGGAg                                                                                     <  1:307874/151‑1 (MQ=255)
                                                       aGCGTGTGCGTAAAGTCTACTGGGGCTATGCGCCGAACGAGAACCTCAGCAACGAAGAGCTGATCCGCGAAAACTACCAGGGCATCCGTCCGGCACCGGGCTATCCGGCCTGCCCGGAACATACGGAAAAAGCCACCATCTGGGAgctg                                                                                  >  2:166998/1‑149 (MQ=255)
                                                             tgCGTAAAGTCTACTGGGGCTATGCGCCGAACGAGAACCTCAGCAACGAAGAGCTGATCCGCGAAAACTACCAGGGCATCCGTCCGGCACCGGGCTATCCGGCCTGCCCGGAACATACGGAAAAAGCCACCATCTGGGAGCTGCTGGAAGt                                                                          <  2:149804/151‑1 (MQ=255)
                                                                     gTCTACTGGGGCTATGCGCCGAACGAGAACCTCAGCAACGAAGAGCTGATCCGCGAAAACTACCAGGGCATCCGTCCGGCACCGGGCTATCCGGCCTGCCCGGAACATACGGAAAAAGCCACCATCTGGGAGCTGCTGGAAGTGGaaaaa                                                                   >  2:318294/1‑150 (MQ=255)
                                                                      tCTACTGGGGCTATGCGCCGAACGAGAACCTCAGCAACGAAGAGCTGATCCGCGAAAACTACCAGGGCATCCGTCCGGCACCGGGCTATCCGGCCTGCCCGGAACATACGGAAAAAGCCACCATCTGGGAGCTGCTGGAAGTGGAAAAac                                                                  <  1:72370/150‑1 (MQ=255)
                                                                        acctgatcgcgctgaatttgtgctacagcgtagtacttgctgtccgggtggctgaagtaccaacccgaaaccgctGCACCGGGCTATCCGGCCTGCCCGGAACATACGGAAAAAGACACCATCTGGGAGCTGCTGGCAGTGGAAAAacaca                                                               <  2:42906‑M1/76‑1 (MQ=255)
                                                                                                        cAACGAAGAGCTGATCCGCGAAAACTACCAGGGCATCCGTCCGGCACCGGGCTATCCGGCCTGCCCGGAACATACGGAAAAAGCCACCATCTGGGAGCTGCTGGAAGTGGAAAAACACACTGGCATGAAACTCACAGAATCTTTCGCCAtg                               <  1:98438/151‑1 (MQ=255)
                                                                                                          aCGAAGAGCTGATCCGCGAAAACTACCAGGGCATCCGTCCGGCACCGGGCTATCCGGCCTGCCCGGAACATACGGAAAAAGCCACCATCTGGGAGCTGCTGGAAGTGGAAAAACACACTGGCATGAAACTCACAGAATCTTTCGCCAtgtg                             <  1:57692/151‑1 (MQ=255)
                                                                                                            gAAGAGCTGATCCGCGAAAACTACCAGGGCATCCGTCCGGCACCGGGCTATCCGGCCTGCCCGGAACATACGGAAAAAGCCACCATCTGGGAGCTGCTGGAAGTGGAAAAACACACTGGCATGAAACTCACAGAATCTTTCGCCATGTgg                            >  2:22167/1‑150 (MQ=255)
                                                                                                                  cTGATCCGCGAAAACTACCAGGGCATCCGTCCGGCACCGGGCTATCCGGCCTGCCCGGAACATACGGAAAAAGCCACCATCTGGGAGCTGCTGGAAGTGGAAAAACACACTGGCATGAAACTCACAGAATCTTTCGCCATGTGGCCCGGTg                     <  2:99273/151‑1 (MQ=255)
                                                                                                                    gATCCGCGAAAACTACCAGGGCATCCGTCCGGCACCGGGCTATCCGGCCTGCCCGGAACATACGGAAAAAGCCACCATCTGGGAGCTGCTGGAAGTGGAAAAACACACTGGCATGAAACTCACAGAATCTTTCGCCATGTGGCCCGGTGCa                   >  2:259511/1‑151 (MQ=255)
                                                                                                                      tCCGCGAAAACTACCAGGGCATCCGTCCGGCACCGGGCTATCCGGCCTGCCCGGAACATACGGAAAAAGCCACCATCTGGGAGCTGCTGGAAGTGGAAAAACACACTGGCATGAAACTCACAGAATCTTTCGCCATGTGGCCCGGTGCATc                 <  2:192143/151‑1 (MQ=255)
                                                                                                                                cTACCAGGGCATCCGTCCGGCACCGGGCTATCCGGCCTGCCCGGAACATACGGAAAAAGCCACCATCTGGGAGCTGCTGGAAGTGGAAAAACACACTGGCATGAAACTCACAGAATCTTTCGCCATGTGGCCCGGTGCATCGGTTTCGGGt       >  1:229076/1‑151 (MQ=255)
                                                                                                                                       cccgaaaccgatGCACCGGGCTATCCGGCCTGCCCGGAACATACGGAAAAAGCCACCATCTGGGAGCTGCTGGAAGTGGAAAAACACACTGGCATGAAACTCACAGAATCTTTCGCCATGTGGCCCGGTGCATCGGTTTCGGGTTGGTa  <  2:264566‑M1/137‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                        
ATGGTGAAAGCGCTTGCCGACCGTTTAGCCGAAGCCTTTGCGGAGTATCTCCATGAGCGTGTGCGTAAAGTCTACTGGGGCTATGCGCCGAACGAGAACCTCAGCAACGAAGAGCTGATCCGCGAAAACTACCAGGGCATCCGTCCGGCACCGGGCTATCCGGCCTGCCCGGAACATACGGAAAAAGCCACCATCTGGGAGCTGCTGGAAGTGGAAAAACACACTGGCATGAAACTCACAGAATCTTTCGCCATGTGGCCCGGTGCATCGGTTTCGGGTTGGTA  >  NC_000913/4227086‑4227369

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