New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913 4343130 =37 (0.930)3 (0.080) 3/252 NT 7.5% coding (1220/1356 nt) melA alpha‑galactosidase
?NC_000913 4343158 = 39 (1.020)coding (1248/1356 nt) melA alpha‑galactosidase

TGGCATTCAGCCGACCAAAGTCGGTACGCTACCTTCGCATCTGGCCGCCCTGATGCAAACCAACATCAACGTACAGACGCTGCTGACCGAAGCTATTCTTACGGAAAATCGCGACCGTGTTTACCACGCCGCGATGATGGACCCGCATACTGCCGCCGTGCTGGGCATTGACGAAATATATGCTCTTGTTGACGACCTGATTGCCGCCCACGGCGACTGGCTGCCAGGCTGGTTGCACCGTTAAAACGCGACTAAACGCTACTGCGCCGGGGGATTT  >  NC_000913/4343023‑4343299
                                                                                                                                       |                                                                                                                                             
tGGCATTCAGCCGACAAAAGTCGGTACGCTACCTTCGCATCTGGCCGCCCTGATGCAAACCAACATCAACGTACAGACGCTGCTGACCGAAGCTATTCTTACGGAAAATCGCGACCGTGTTTACCACGCCGCGATGATGGAcc                                                                                                                                        >  1:785730/1‑143 (MQ=255)
          ccGACCAAAGTCGGTACGCTACCTTCGCATCTGGCCGCCCTGATGCAAACCAACATCAACGTACAGACGCTGCTGACCGAAGCTATTCTTACGGCAAATCGCGACCGTGTTTACCACGCCGCGCTGATGGACCCGCATACTgc                                                                                                                              >  2:12673/1‑143 (MQ=255)
             gcctggcagccagtcgccgtgggcggcaatcaggtcgtcaacaagagcatatatttcgtcaatgcccagcacggcggcagtatgcgggtccatcATCGCGACCGTGTTTACCACGCCGCgatgat                                                                                                                                             <  2:664998‑M1/31‑1 (MQ=255)
             gcctggcagccagtcgccgtgggcggcaatcaggtcgtcaacaagagcatatatttcgtcaatgcccagcacggcggcagtatgcgggtccatcATCGCGACCGTGTTTACCACGCCGCgatgat                                                                                                                                             >  1:664998‑M1/95‑125 (MQ=255)
               cAAAGTCGGTACGCTACCTTCGCATCTGGCCGCCCTGATGCAAACCAACATCAACGTACAGACGCTGCTGACCGAAGCTATTCTTACGGAAAATCGCGACCGTGTTTACCACGCCGCGATGATGGACCCGCATACTGccgccg                                                                                                                         >  1:625477/1‑143 (MQ=255)
               cAAAGTCGGTACGCTACCTTCGCATCTGGCCGCCCTGATGCAAACCAACATCAACGTACAGACGCTGCTGACCGAAGCTATTCTTACGGAAAATCGCGACCATGTTTAGCACGCCGCGATGATGGACCCGCATACTGccggcg                                                                                                                         >  1:48732/1‑143 (MQ=255)
                     cGGTACGCTACCTTCGCATCTGGCCGCCCTGATGCAAACCAACATCAACGTACAGACGCTGCTGACCGAAGCTATTCTTACGGAAAATCGCGACCGTGTTTACCACGCCGCGATGATGGACCCGCATACTGCCGCCGTGCTgg                                                                                                                   <  2:329494/143‑1 (MQ=255)
                       gTACGCTACCTTCGCATCTGGCCGCCCTGATGCAAACCAACATCAACGTCAAGACGCTGCTGACCGAAGATACTCATACGGAAAATCGCGACCGGGTTTACCACGCCGCCGTGATGGACCCGCATACTGCCGGCGTGCTggac                                                                                                                 >  2:278595/1‑141 (MQ=255)
                            cTACCTTCGCATCTGGCCTCCCTGATGCAAAGCAACATCAACGTACAGACGCTGCTGACCGAGGCTATTCTTACGGAAAATCGCGACCGTGTTTACCACGCCGCGATGATGGACCCGCATACTGCCGCCGTGCTGGGCATTGa                                                                                                            <  2:633606/143‑1 (MQ=255)
                               ccTTCGCATCTGGCCGCCCTGATGCAAACCAACATCAACGTACAGACGCTGCTGACCGAAGCTATTCTTACGGACAATCGCGACCGTGTTTACCACGCCGCGATGATGGACCCGCCTACTGCCGCCGTGCCGGGCATTCACgc                                                                                                         >  2:537665/1‑142 (MQ=255)
                                  tCGCATCTGGCCGCCCTGATGCAAACCAACATCAACGTACAGACGCTGCTGACCGAAGCTATTCTTACGGAAAATCGCGACCGTGTTTACCACGCCGCGATGATGGACCCGCATACTGCCGCCGTGCTGGGCATTGACGAAat                                                                                                      <  1:189491/143‑1 (MQ=255)
                                  tCGCATCTGGCCGCCCTGATGCAAACCAACATCAACGTACAGACGCTGCTGACCGAAGCTATTCTTACGGAAAATCGCGACCGTGTTTACCACGCCGCGATGATGGACCCGCATACTGCCGCCGTGCTGGGCATTGACGAAat                                                                                                      >  1:464031/1‑143 (MQ=255)
                                         tGGCCGCCCTGATGCAAACCAACATCAACGTACAGACGCTGCTGACCGAAGCTATTCTTACGGAAAATCGCGACCGTGTTTACCACGCCGCGATGATGGACCCGCATACTGCCGCCGTGCTGGGCATTGACGAAATATATGct                                                                                               >  1:699898/1‑143 (MQ=255)
                                               cccTGATGCAAACCAACATCAACGTACAGACGCTGCTGACCGAAGCTATTCTTACGGAAAATCGCGACCGTGTTTACCACGCCGCGATGATGGACCCGCATACTGCCGCCGTGCTGGGCATTGACGAAATATATGCTCttgtt                                                                                         <  1:468877/143‑1 (MQ=255)
                                                    aTGCAAACCAACATCAACGTACAGACGCTGCTGACCGAAGCTATTCTTACGGAAAATCGCGACCGTGTTTACCACGCCGCGATGATGGACCCGCATACTGCCGCCGTGCTGGGCATTGACGAAATATATGCTCTTGTTgacga                                                                                    <  2:35643/143‑1 (MQ=255)
                                                     tGCAAACCAACATCAACGTACAGACGCTGCTGACCGAAGCTATTCTTACGGAAAATCGCGACCGTGTTTACCACGCCGCGATGATGGACCCGCATACTGCCGCCGTGCTGGGCATTGACGAAATATATGCTCTTGTTgacgac                                                                                   <  1:481235/143‑1 (MQ=255)
                                                         aaCCAACATCAACGTACAGACGCTGCGGACCGAAGCTATTCTTACGGAAAATCGCGACCGTGTTTACCACGCCGCGATGATGGACCCGCATACTGCCGCCGTGCTGGGCATTGACGAAATATATGCTCTTGTTGACGACCTGa                                                                               <  1:676370/143‑1 (MQ=255)
                                                            caccATCTACGTACAGACGCTGCTGACCGCAGCTATTCTTACGGAAAATCGCGACCGTGTTTACCACGCCGCGATGATGGATCCGCATACTGCCGCCGTGCTTGGCATTGACGAAATATATGCTCTTGTTGACGACCTGATTg                                                                            <  2:102484/140‑1 (MQ=255)
                                                            cAACATCAACGTACAGACGCTGCTGACCGAAGCTATTCTTACGGAAAATCGCGACCGTGTTTACCACGCCGCGATGATGGACCCGCATACTGCCGCCGGGCTGGGCATTGACGAAATATATGCTCGTGTTGACGACCTGATTg                                                                            >  1:394733/1‑143 (MQ=255)
                                                               cATCAACGTACAGACGCTGCTGACCGAAGCTATTCTTACGGAAAATCGCGACCGTGTTTACCACGCCGCGATGATGGACCCGCATACTGCCGCCGTGCTGGGCATTGACGAAATATATGCTCTTGTTGACGACCTGATTgccg                                                                         <  1:39608/143‑1 (MQ=255)
                                                                aTCAACGTACAGACGCTGCTGACCGAAGCTATTCTTACGGAAAATCGCGACCGTGTTTACCACGCCGCGATGATGGACCCGCATACTGCCGCCGTGCTGGGCAtt                                                                                                              >  2:440756/1‑105 (MQ=255)
                                                                aTCAACGTACAGACGCTGCTGACCGAAGCTATTCTTACGGAAAATCGCGACCGTGTTTACCACGCCGCGATGATGGACCCGCATACTGCCGCCGTGCTGGGCAtt                                                                                                              <  1:440756/105‑1 (MQ=255)
                                                                  cAACGTACAGACGCTGCTGACCGAAGCTATTCTTACGGAAAATCGCGACCGTGTTTACCACGCCGCGATGATGGACCCGCATACTGCCGCCGTGCTGGGCATTGACGAAATATATGCTCTTGTTGACGACCTGATTGCCGccc                                                                      <  2:430534/143‑1 (MQ=255)
                                                                   aaCGGACAGACGCTGCTGACCGAGGCTGTTCTTACGGAAAAGCGCGGCCGTGTGGACCACGCCGCGATGATGGACCCGCATACTGCCGCCGTGCTGGGCATTGACGAAATATATGCTCTTGTTGACGACCTGATTGCCGCCCa                                                                     <  2:224665/143‑1 (MQ=255)
                                                                       tACAGCCGCTGCTGACCGAAGCTATTCTTACGGAAAATCGCGACCGTGTTTACCACGCCGCGATGATGGACCCGCATACTGCCGCCGTGCTGGGCATTGACGAAATATATGCTCTTGTTGACGACCTGATTGCCGCCCACGGc                                                                 <  1:294032/143‑1 (MQ=255)
                                                                            aCGCTGCTGACCGAAGCTATTCTTACGGAAAATCGCGACCGTGTTTACCACGCCGCGATGATGGACCCGCATACTGCCGCCGTGCTGGGCATTGACGAAATATATGCTCTTGTTGACGACCTGATTGCCGCCCACGGCGAcgg                                                            >  2:152551/1‑141 (MQ=255)
                                                                              gctgctGACCGAAGCTATTCTTACGGAAAATCGCGACCGTGTTTACCACGCCGCGATGATGGACCCCCATACTGCCGCAGCGCTGGGCATTGACGAAATATATGCTCTTGTTGACGACCTGATTGCCCCCCACCGCGCCTGGc                                                          >  2:472419/1‑143 (MQ=255)
                                                                               ctgctgCCCGAAGCTAGTCTTACGGAAAATCGCGACCGTGTTTACCACGCCGCGATGATGGACCCGCATACTGCCGCCGTGCTGGGCATTGACGAAATATATGCTCTTGTTGACGACCTGATTGCCGCCCACGGCGACTGGCt                                                         <  2:674592/143‑1 (MQ=255)
                                                                                tgctgACCGAAGCTATTCTTACGGAAAATCGCGACCGTGTTTACCACGCCGCGATGATGGACCCGCATACTTCTGCCGTGCTGGGCATTGACGAAATATATGCTCTTGTTGACGACCTGATTGCCGCCCACGGCGACTGGCTg                                                        <  1:537665/143‑1 (MQ=255)
                                                                                     aCCGAAGCTATTCTTACGGAAAATCGCGACCGTGTTTACCACGCCGCGATGATGGACCCGCATACTGCCGCCGTGCTGGGCATTGACGAAATATATGCTCTTGTTGACGACCTGATTGCCGCCCACGGCGACTGGCTGCCAgg                                                   >  1:24649/1‑143 (MQ=255)
                                                                                        gAAGCTATTCTTACGGAAAATCGCGACCGTGTTTACCACGCCGCGATGATGGACCCGCATACTGCCGCCGTGCTGGGCATTGACGAAATATATGCTCTTGTTGACGACCTGATTTCCGCCCACGGCGACTGGCTGCCGGGCTg                                                >  2:507922/1‑143 (MQ=255)
                                                                                        gAAGCTATTCTTACGGAAAATCGCGACCGTGTTTACCACGCCGCGATGATGGACCCGCATACTGCCGCCGTGCTGGGCATTGACGAAATATATGCTCTTGTTGACGACCTGATTGCCACCCAGCGCGACTGGCTGCGAGGCTg                                                >  2:214556/1‑143 (MQ=255)
                                                                                        gAAGCTAGTCTTACGGAAAATCGCGACCGTGTTTACCACGCCGCGATGATGGACCCGCATACTGCCGCCGTGCTCCGCATTGACGAAATATATGCTCTTGTTGACGACCTGATTGCCGCCCACGGCGACTGGCTGCCAGGCTg                                                <  2:516917/143‑1 (MQ=255)
                                                                                               ttCTTACGGAAACTCGCGTCCGTCTTTACCACGCCGCGATGATGGACCCGCATACTGCCGCCGTGCTGGGCATTGACTAAATATATGCTCTTGTTGGCGACCTGATTGCCGCCCACGGCGACTGGCTGCCAGGCTGGTTGCAc                                         <  2:488823/143‑1 (MQ=255)
                                                                                                      ggAAAATCGCGACCGTGTTTACCACGCCGCGATGATGGACCCGCATACTGCCGCCGTGCTGGGCATTGACGAAATATATGCTCTTGTTGACGACCTGATTGCCGCCCACGGCGGCTGGCAGCCGGGCTGTGTGCACCGTaaac                                  >  2:261470/1‑142 (MQ=255)
                                                                                                        atcatcgcggcgtggtaaacacggtcgcgatGATGGACCCGCATACTGCCGCCGTGCTGGGCATTGACGAAATATATGCTCTTGTTGACGACCTGATTGCCGCCCACGGCGACTGGCTGCCAGGc                                                  >  1:664998‑M2/32‑125 (MQ=255)
                                                                                                        atcatcgcggcgtggtaaacacggtcgcgatGATGGACCCGCATACTGCCGCCGTGCTGGGCATTGACGAAATATATGCTCTTGTTGACGACCTGATTGCCGCCCACGGCGACTGGCTGCCAGGc                                                  <  2:664998‑M2/94‑1 (MQ=255)
                                                                                                        aaaaTCGCGACCGTGTTTACCACGCCGCGATGATGGACCCGCATACTGCGGCCGTGCTGGGCATTGACGAAGTAGATGCTCTTGTTGACGACCTGATTGCCGCCCACGGCCATGGGCTGCCAGGCTGGTTGCACCggggaaac                                >  2:442344/1‑136 (MQ=255)
                                                                                                          aaTCGCGACCGTGTTTACCACGCCGCGATGATTGACCCGCATACTGCCGCCGTGCTGGGCATTGACGAAATATATGCTCTTGTTGACGACCTGATTGCCGCCCACGGCGACTGGCTGCCAGGCTGGTTGCACCGTTAAAAcgc                              <  1:672488/143‑1 (MQ=255)
                                                                                                          aaTCGCGACCGTGTTTACCACGCCGCGATGATGGACCCGCATACTGCCGCCGTGCTGGTCATTGACGAAGTATACGCTCTTGTTGACGACCTGATTGCCGCCCACGGCGACTGGCTGCCAGGCTGGTTGCACCGTTAAAAcgc                              <  2:707303/143‑1 (MQ=255)
                                                                                                           atcgcggcgtggtaaacacggtcgcgatGATGGACCCGCATACTGCCGCCGTGCTGGGCCTGGACGAAATAGATGCTCTTGTTGACGACCGGATTGCCGCCCACGGCGAGTGGATGCCCGGCTGGTTGCACCGTTAGACcgcg                             >  2:197744‑M2/29‑143 (MQ=255)
                                                                                                              GCGACCGTGTTTACCACGCCGCGATGATGGACCCGCATACTGCCGCCGTGCTGGGCATTGACGAAATATATGCTCTTGTTGACGACCTGATTGTCGCCCACGGCGAGTGGCTGCCAGGCTGGTTGCACCGTTAAAACGCGACt                          <  1:761264‑M1/143‑1 (MQ=255)
                                                                                                              GCGACCGTGTTTACCACGCCGCGATGATGGACCCGCATACTGCCGCCGTGCTGGGCATTGACGAAATATATGCTCTTGTTGACGACCTGATTGCCGCCCACGGCGACTGGCTGCCAGGCTGGTTGCACCGTTAAAACGCGACt                          <  2:737917‑M1/143‑1 (MQ=255)
                                                                                                              GCGACCGTGTTTACCACGCCGCGATGATGGACCCGCATACTGCCGCCGTGCTGGGCATTGACGAAATATATGCTCTTGTTGACGACCTGATTGCCGCCCACGGCGACTGGCTGCCACGCTGGTTGCACCGTTCAAACGCGACt                          <  1:497415‑M1/143‑1 (MQ=255)
                                                                                                               CGACCGTGTTTACCACGCCGCGATGATGGACCCGCATACTGCCGCCGTGCTGGGCATTGACGAAATATATGCTCTTGTTGACGACCTGATTGGCGCCCACGGCGACTGGCTGCCCGGCTGGTTGCACCGgtaaaacgatacta                         <  2:255516‑M1/143‑7 (MQ=255)
                                                                                                                 aCCGTGTTTACCACGCCGCGATGATGGACCCGCATACTGCCGCCGTGCTGGGCATTGACGAAATATATGCTCTTGTTGACGACCTGATTGCCGCCCACGGCGACTGGCTGCCAGGCTGGTTGCACCGTTAAAACGCGACTaaa                       >  1:328452/1‑143 (MQ=255)
                                                                                                                       tttACCACGCCGCGATGATGGACCCGCATACTGCCGCCGTGCTGGGCATTGACGAAATATATGCTCTTGTTGACGACCTGATTGCCGCCCACGGCGACTGGCTGCCAGGCTGGTTGCACCGTTaaaa                                 <  1:103086/127‑1 (MQ=255)
                                                                                                                       tttACCACGCCGCGATGATGGACCCGCATACTGCCGCCGTGCTGGGCATTGACGAAATATATGCTCTTGTTGACGACCTGATTGCCGCCCACGGCGACTGGCTGCCAGGCTGGTTGCACCGTTaaaa                                 >  2:103086/1‑127 (MQ=255)
                                                                                                                                 cgcgatGATGGACCCGCATACTGCCGCCGTGCTGGGCATTGACGAAATATATGCTCTTGTTGACGACCTGATTGCCGCCCACGGCGACTGGCTGCCAGGCTGGTTGCACCGTTAAAACGCGACTAAACGCTACTGCGCCgggg       <  1:307766‑M2/137‑1 (MQ=255)
                                                                                                                                 cgcgatGATGGACCCGCATACTGCCGCCGTGCTGGGCATTGACGAAATATATGCTCTTGTTGACGACCTGATTGCCGCCCACGGCGACTGGCTGCCAGGCTGGTTGCACCGTTAAAACGCGACTAAACGCTACTGCGCCgggg       <  2:699898‑M2/137‑1 (MQ=255)
                                                                                                                                    gatGATGGACCCGCATACTGCCGCCGTGCTGGGCATTGACGAAATATATGCTCTTGTTGACGACATGATTGCCGCCCACGGCGACTGGCTGCCAGGCTGGTTGCACCGTTAAAACGCGACTAAACGCTACTGCGCCGGGGGAt    >  2:93812‑M2/4‑143 (MQ=255)
                                                                                                                                      tGATGGACCCGCATACTGCCGCCGTGCTGGGCCTCCCCGAACTATTTGCTCTTGTTGACGACCTGATTGCCGCCCACGGCGACTGGCTGCCAGGCTGGTTGCACCGTTAAAACGCGACTAAACGCTACTGCGCCGGGGGAttt  <  1:27425‑M2/142‑1 (MQ=255)
                                                                                                                                       |                                                                                                                                             
TGGCATTCAGCCGACCAAAGTCGGTACGCTACCTTCGCATCTGGCCGCCCTGATGCAAACCAACATCAACGTACAGACGCTGCTGACCGAAGCTATTCTTACGGAAAATCGCGACCGTGTTTACCACGCCGCGATGATGGACCCGCATACTGCCGCCGTGCTGGGCATTGACGAAATATATGCTCTTGTTGACGACCTGATTGCCGCCCACGGCGACTGGCTGCCAGGCTGGTTGCACCGTTAAAACGCGACTAAACGCTACTGCGCCGGGGGATTT  >  NC_000913/4343023‑4343299

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.