Predicted mutation
evidence seq id position mutation freq annotation gene description
JC NC_000913 3,548,550 24 bp→TTTTTTTTTAATGATACGGCGACC 6.5% coding (1498‑1521/2085 nt) malQ ← 4‑alpha‑glucanotransferase

New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913 = 354854958 (0.820)3 (0.060)
+TTTTTTTTTAATGATACGGCGACC
3/128 NT 6.5% coding (1522/2085 nt) malQ 4‑alpha‑glucanotransferase
?NC_000913 3548574 = 62 (0.870)coding (1497/2085 nt) malQ 4‑alpha‑glucanotransferase

ACGGAACGTCTTCTCGTGGTCGTTTTCGAAATAGAGCACTTTGTAAGAGTACACACCGCTGCTGCGCAGCTTACCGACAATCTCTACCGGTACGGTACCGAGATCTTCACCAATCACCATACAGCGATGACGTTTACTTTCGAGTGCCAGAATCGAGAGCAGATCATCCACCGGATAGTGAACATACGC  >  NC_000913/3548478‑3548666
                                                                                                |                                                                                            
aCGGAACGTCTTCTCGTGGTCGTTTTCGAAATAGAGCACTTTGTAAGAGTACACACCGCTGCTGCGCAGCTTACCGACAATCTCTACCGGTACGGTAcc                                                                                            <  1:435801/99‑1 (MQ=255)
   gAACGTCTTCTCGTGGTCGTTTTCGAAATAGAGCACTTTGTAAGAGTACACACCGCTGCTGCGCAGCTTACCGACAATCTCTACCGGTACGGTACCgag                                                                                         >  1:1223795/1‑99 (MQ=255)
   gAACGTCTTCTCGTGGTCGTTTTCGAAATAGAGCACTTTGTAAGAGTACACACCGCTGCTGCGCAGCTTACCGACAATCTCTACCGGTACGGTACCgag                                                                                         <  2:124448/99‑1 (MQ=255)
    aaCGTCTTCTCGTGGTCGTTTTCGAAATAGAGCACTTTGTAAGAGTACACACCGCTGCTGCGCAGCTTACCGACAATCTCTACCGGTACGGTACCgaga                                                                                        <  2:863082/99‑1 (MQ=255)
    aaCGTCTTCTCGTGGTCGTTTTCGAAATAGAGCACTTTGTAAGAGTACACACCGCTGCTGCGCAGCTTACCGACAATCTCTACCGGTACGGTACCgaga                                                                                        >  2:674502/1‑99 (MQ=255)
    aaCGTCTTCTCGTGGTCGTTTTCGAAATAGAGCACTTTGTAAGAGTACACACCGCTGCTGCGCAGCTTACCGACAATCTCTACCGGTACGGTACCgaga                                                                                        <  1:1551842/99‑1 (MQ=255)
       gTCTTCTCGTGGTCGTTTTCGAAATAGAGCACTTTGTAAGAGTACACACCGCTGCTGCGCAGCTTACCGACAATCTCTACCGGTACGGTACCGAGATCt                                                                                     <  2:1198915/99‑1 (MQ=255)
        tcttctCGTGGTCGTTTTCGAAATAGAGCACTTTGTAAGAGTACACACCGCTGCTGCGCAGCTTACCGACAATCTCTACCGGTACGGTACCGAGATCtt                                                                                    <  1:333365/99‑1 (MQ=255)
        tcttctCGTGGTCGTTTTCGAAATAGAGCACTTTGTAAGAGTACACACCGCTGCTGCGCAGCTTACCGACAATCTCTACCGGTACGGTACCGAGATCtt                                                                                    >  1:1276682/1‑99 (MQ=255)
        tcttctCGTGGTCGTTTTCGAAATAGAGCACTTTGTAAGAGTACACACCGCTGCTGCGCAGCTTACCGACAATCTCTACCGGTACGGTACCGAGATCtt                                                                                    >  1:787646/1‑99 (MQ=255)
         cttctCGTGGTCGTTTTCGAAATAGAGCACTTTGTAAGAGTACACACCGCTGCTGCGCAGCTTACCGACAATCTCTACCGGTACGGTACCGAGATCTTc                                                                                   <  2:575614/99‑1 (MQ=255)
          ttctCGTGGTCGTTTTCGAAATAGAGCACTTTGTAAGAGTACACACCGCTGCTGCGCAGCTTACCGACAATCTCTACCGGTACGGTACCGAGATCTTCa                                                                                  >  1:1926042/1‑99 (MQ=255)
                   tCGTTTTCGAAATAGAGCACTTTGTAAGAGTACACACCGCTGCTGCGCAGCTTACCGACAATCTCTACCGGTACGGTACCGAGATCTTCACCAATCAcc                                                                         >  1:1675720/1‑99 (MQ=255)
                    cGTTTTCGAAATAGAGCACTTTGTAAGAGTACACACCGCTGCTGCGCAGCTTACCGACAATCTCTACCGGTACGGTACCGAGATCTTCACCAATCACCa                                                                        <  1:357685/99‑1 (MQ=255)
                    cGTTTTCGAAATAGAGCACTTTGTAAGAGTACACACCGCTGCTGCGCAGCTTACCGACAATCTCTACCGGTACGGTACCGAGATCTTCACCAATCACCa                                                                        <  1:22431/99‑1 (MQ=255)
                     gTTTTCGAAATAGAGCACTTTGTAAGAGTACACACCGCTGCTGCGCAGCTTACCGACAATCTCTACCGGTACGGTACCGAGATCTTCACCAATCACCAt                                                                       >  2:1924533/1‑99 (MQ=255)
                      ttttCGAAATAGAGCACTTTGTAAGAGTACACACCGCTGCTGCGCAGCTTACCGACAATCTCTACCGGTACGGTACCGAGATCTTCACCAATCACCATa                                                                      >  1:574737/1‑99 (MQ=255)
                      ttttCGAAATAGAGCACTTTGTAAGAGTACACACCGCTGCTGCGCAGCTTACCGACAATCTCTACCGGTACGGTACCGAGATCTTCACCAATCACCATa                                                                      >  2:1772206/1‑99 (MQ=255)
                       tttCGAAATAGAGCACTTTGTAAGAGTACACACCGCTGCTGCGCAGCTTACCGACAATCTCTACCGGTACGGTACCGAGATCTTCACCAATCACCATAc                                                                     <  1:1375984/99‑1 (MQ=255)
                            aaaTAGAGCACTTTGTAAGAGTACACACCGCTGCTGCGCAGCTTACCGACAATCTCTACCGGTACGGTACCGAGATCTTCACCAATCACCATACAGCGa                                                                >  1:1115778/1‑99 (MQ=255)
                             aaTAGAGCACTTTGTAAGAGTACACACCGCTGCTGCGCAGCTTACCGACAATCTCTACCGGTACGGTACCGAGATCTTCACCAATCACCATACAGCGAt                                                               >  2:210668/1‑99 (MQ=255)
                              aTAGAGCACTTTGTAAGAGTACACACCGCTGCTGCGCAGCTTACCGACAATCTCTACCGGTACGGTACCGAGATCTTCACCAATCACCATACAGCGATg                                                              >  1:611899/1‑99 (MQ=255)
                               tAGAGCACTTTGTAAGAGTACACACCGCTGCTGCGCAGCTTACCGACAATCTCTACCGGTACGGTACCGAGATCTTCACCAATCACCATACAGCGATGa                                                             <  1:1678009/99‑1 (MQ=255)
                                   gcactttgtaagagtacacaccgctgctgcgcagctttttttttttaatgatacggcgaccACCGAGATCTACACCAATCACCATACAGCGATGACGtt                                                         <  2:230520‑M2/38‑1 (MQ=255)
                                   gCACTTTGTAAGAGTACACACCGCTGCTGCGCAGCTTACCGACAATCTCTACCGGTACGGTACCGAGATCTTCACCAATCACCATACAGCGATGACGtt                                                         <  1:1330714/99‑1 (MQ=255)
                                   gCACTTTGTAAGAGTACACACCGCTGCTGCGCAGCTTACCGACAATCTCTACCGGTACGGTACCGAGATCTTCACCAATCACCATACAGCGATGACGtt                                                         >  2:826600/1‑99 (MQ=255)
                                   gCACTTTGTAAGAGTACACACCGCTGCTGCGCAGCTTACCGACAATCTCTACCGGTACGGTACCGAGATCTTCACCAATCACCATACAGCGATGACGtt                                                         <  1:485840/99‑1 (MQ=255)
                                     aCTTTGTAAGAGTACACACCGCTGCTGCGCAGCTTACCGACAATCTCTACCGGTACGGTACCGAGATCTTCACCAATCACCATACAGCGATGACGTTTa                                                       <  1:256609/99‑1 (MQ=255)
                                      ctttgtaagagtacacaccgctgctgcgcagctttttttttttaatgatacggcgaccACCGAGATCTACACCAATCACCATACAGCGATGACGTTTAc                                                      <  2:1208720‑M2/41‑1 (MQ=255)
                                        ttGTAAGAGTACACACCGCTGCTGCGCAGCTTACCGATATTATCTACCGGTACGCTACCGAGATCTACACCAATCACCATACAGCGATGACGTTTACtt                                                    <  2:233947/99‑1 (MQ=255)
                                        ttGTAAGAGTACACACCGCTGCTGCGCAGCTTACCGACAATCTCTACCGGTACGGTACCGAGATCTTCACCAATCACCATACAGCGATGACGTTTACtt                                                    <  1:1773363/99‑1 (MQ=255)
                                          gTAAGAGTACACACCGCTGCTGCGCAGCTTACCGACAATCTCTACCGGTACGGTACCGAGATCTTCACCAATCACCATACAGCGATGACGTTTACTTTc                                                  <  2:317101/99‑1 (MQ=255)
                                             agagTACACACCGCTGCTGCGCAGCTTACCGATTATCTCTACCCGTACACTACCGAGATCTACACCAATCACCATACAGCGATGACGTTTACTTTCGAg                                               <  2:971504/99‑1 (MQ=255)
                                             agagTACACACCGCTGCTGCGCAGCTTACCGACAATCTCTACCGGTACGGTACCGAGATCTTCACCAATCACCATACAGCGATGACGTTTACTTTCGAg                                               >  2:451886/1‑99 (MQ=255)
                                                gTACACACCGCTGCTGCGCAGCTTACCGACAATCTCTACCGGTACGGTACCGAGATCTTCACCAATCACCATACAGCGATGACGTTTACTTTCGAGTGc                                            >  2:197074/1‑99 (MQ=255)
                                                gTACACACCGCTGCTGCGCAGCTTACCGACAATCTCTACCGGTACGGTACCGAGATCTTCACCAATCACCATACAGCGATGACGTTTACTTTCGAGTGc                                            >  1:1427986/1‑99 (MQ=255)
                                                     cacCGCTGCTGCGCAGCTTACCTACTTTCTCGACCCGTCCGCCACCGAGATCTACACCAATCACCATACAGCGATGACGTTTACTTTCGAGTGCCAGaa                                       <  2:1803610/99‑1 (MQ=255)
                                                      acCGCTGCTGCGCAGCTTACCGACATTCTCTACCCGTACGCTACCGAGATCTTCACCAATCACCATACAGCGATGACGTTTACTTTCGAGTGCCAGAAt                                      <  2:271919/99‑1 (MQ=255)
                                                      acCGCTGCTGCGCAGCTTACCGACAATCTCTACCGGTACGGTACCGAGATCTTCACCAATCACCATACAGCGATGACGTTTACTTTCGAGTGCCAGAAt                                      >  1:1701405/1‑99 (MQ=255)
                                                       ccGCTGCTGCGCAGCTTACCGACAATCTCTACCGGTACGGTACCGAGATCTTCACCAATCACCATACAGCGATGACGTTTTCTTTCGAGTGCCAGAATc                                     <  1:1809358/99‑1 (MQ=255)
                                                       ccGCTGCTGCGCAGCTTACCGACAATCTCTACCGGTACGGTACCGAGATCTTCACCAATCACCATACAGCGATGACGTTTACTTTCGAGTGCCAGAATc                                     >  1:1876923/1‑99 (MQ=255)
                                                         gctgctGCGCAGCTTACCGACAATCTCTACCGGTACGGTACCGAGATCTTCACCAATCACCATACAGCGATGACGTTTACTTTCGAGTGCCAGAATCga                                   <  2:1263170/99‑1 (MQ=255)
                                                          ctgctgCGCAGCTTACCGACAATCTCTACCGGTACGGTACCGAGATCTTCACCAATCACCATACAGCGATGACGTTTACTTTCGAGTGCCAGAATCgag                                  >  1:1274950/1‑99 (MQ=255)
                                                          ctgctgCGCAGCTTACCGACAATCTCTACCGGTACGGTACCGAGATCTTCACCAATCACCATACAGCGATGACGTTTACTTTCGAGTGCCAGAATCgag                                  >  1:1019094/1‑99 (MQ=255)
                                                          ctgctgCGCAGCTTACCGACAATCTCTACCGGTACGGTACCGAGATCTTCACCAATCACCATACAGCGATGACGTTTACTTTCGAGTGCCAGAATCgag                                  >  1:806474/1‑99 (MQ=255)
                                                              tgcGCAGCTTACCGACAATCTCTACCGGTACGGTACCGAGATCTTCACCAATCACCATACAGCGATGACGTTTACTTTCGAGTGCCAGAATCGAGAGCa                              <  2:902503/99‑1 (MQ=255)
                                                                 gcagcttatcttctttcatgatccgtccaccACCGAGATCTACACCAATCACCATACAGCGATGACGTTTACTTTCGAGTGCCAGAATCGAGAGCAGat                           <  2:1485308‑M2/68‑1 (MQ=255)
                                                                 gcAGCTTACCGACAATCTCTACCGGTACGGTACCGAGATCTTCACCAATCACCATACAGCGATGACGTTTACTTTCGAGTGCCAGAATCGAGAGCAGat                           <  1:848695/99‑1 (MQ=255)
                                                                  cAGCTTACCGACAATCTCTACCGGTACGGTACCGAGATCTTCACCAATCACCATACAGCGATGACGTTTACTTTCGAGTGCCAGAATCGAGAGCAGatc                          >  2:1528360/1‑99 (MQ=255)
                                                                  cAGCTTACCGACAATCTCTACCGGTACGGTACCGAGATCTTCACCAATCACCATACAGCGATGACGTTTACTTTCGAGTGCCAGAATCGAGAGCAGatc                          <  2:1514195/99‑1 (MQ=255)
                                                                    gCTTACCGACAATCTCTACCGGTACGGTACCGAGATCTTCACCAATCACCATACAGCGATGACGTTTACTTTCGAGTGCCAGAATCGAGAGCAGatcat                        >  1:632378/1‑99 (MQ=255)
                                                                    gCTTACCGACAATCTCTACCGGTACGGTACCGAGATCTTCACCAATCACCATACAGCGATGACGTTTACTTTCGAGTGCCAGAATCGAGAGCAGatcat                        >  1:1276612/1‑99 (MQ=255)
                                                                     cTTACCGACAATCTCTACCGGTACGGTACCGAGATCTTCACCAATCACCATACAGCGATGACGTTTACTTTCGAGTGCCAGAATCGAGAGCAGatcatc                       >  1:1384667/1‑99 (MQ=255)
                                                                     cTTACCGACAATCTCTACCGGTACGGTACCGAGATCTTCACCAATCACCATACAGCGATGACGTTTACTTTCGAGTGCCAGAATCGAGAGCAGatcatc                       <  1:241791/99‑1 (MQ=255)
                                                                        aCCGACAATCTCTACCGGTACGGTACCGAGATCTTCACCAATCACCATACAGCGATGACGTTTACTTTCGAGTGCCAGAATCGAGAGCAGATCATCCAc                    <  1:721586/99‑1 (MQ=255)
                                                                             cAATCTCTACCGGTACGGTACCGAGATCTTCACCAATCACCATACAGCGATGACGTTTACTTTCGAGTGCCAGAATCGAAAGCAGATCATCCACCGGAt               >  2:1338127/1‑99 (MQ=255)
                                                                              aaTCTCTACCGGTACGGTACCGAGATCTTCACCAATCACCATACAGCGATGACGTTTACTTTCGAGTGCCAGAATCGAGAGCAGATCATCCACCGGATa              <  2:1621461/99‑1 (MQ=255)
                                                                               aTCTCTACCGGTACGGTACCGAGATCTTCACCAATCACCATACAGCGATGACGTTTACTTTCGAGTGCCAGAATCGAGAGCAGATCATCCACCGGATAg             >  2:997440/1‑99 (MQ=255)
                                                                                    tACCGGTACGGTACCGAGATCTTCACCAATCACCATACAGCGATGACGTTTACTTTCGAGTGCCAGAATCGAGAGCAGATCATCCACCGGATAGTGAAc        <  2:8820/99‑1 (MQ=255)
                                                                                     aCCGGTACGGTACCGAGATCTTCACCAATCACCATACAGCGATGACGTTTACTTTCGAGTGCCAGAATCGAGAGCAGATCATCCACCGGATAGTGAACa       <  2:1933590/99‑1 (MQ=255)
                                                                                     aCCGGTACGGTACCGAGATCTTCACCAATCACCATACAGCGATGACGTTTACTTTCGAGTGCCAGAATCGAGAGCAGATCATCCACCGGATAGTGAACa       <  1:603014/99‑1 (MQ=255)
                                                                                     aCCGGTACGGTACCGAGATCTTCACCAATCACCATACAGCGATGACGTTTACTTTCGAGTGCCAGAATCGAGAGCAGATCATCCACCGGATAGTGAACa       <  2:696498/99‑1 (MQ=255)
                                                                                     aCCGGTACGGTACCGAGATCTTCACCAATCACCATACAGCGATGACGTTTACTTTCGAGTGCCAGAATCGAGAGCAGATCATCCACCGGATAGTGAACa       <  1:270378/99‑1 (MQ=255)
                                                                                         gtacggtacCGAGATCTTCACCAATCACCATACAGCGATGACGTTTACTTTCGAGTGCCAGAATCGAGAGCAGATCATCCACCGGATAGTGAACATAcg   >  1:301823/1‑99 (MQ=255)
                                                                                         gtacggtacCGAGATCTTCACCAATCACCATACAGCGATGACGTTTACTTTCGAGTGCCAGAATCGAGAGCAGATCATCCACCGGATAGTGAACATAcg   <  1:132392/99‑1 (MQ=255)
                                                                                          tacggtacCGAGATCTTCACCAATCACCATACAGCGATGACGTTTACTTTCGAGTGCCAGAATCGAGAGCAGATCATCCACCGGATAGTGAACATAcgc  >  1:1402727/1‑99 (MQ=255)
                                                                                                |                                                                                            
ACGGAACGTCTTCTCGTGGTCGTTTTCGAAATAGAGCACTTTGTAAGAGTACACACCGCTGCTGCGCAGCTTACCGACAATCTCTACCGGTACGGTACCGAGATCTTCACCAATCACCATACAGCGATGACGTTTACTTTCGAGTGCCAGAATCGAGAGCAGATCATCCACCGGATAGTGAACATACGC  >  NC_000913/3548478‑3548666

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