Predicted mutation
evidence seq id position mutation annotation gene description
JC NC_000913 2,533,399 Δ10 bp intergenic (+21/‑354) cysK → / → ptsH cysteine synthase A, O‑acetylserine sulfhydrolase A subunit/phosphohistidinoprotein‑hexose phosphotransferase component of PTS system (Hpr)

New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913 = 25333982 (0.030)59 (0.900) 48/280 0.3 96.0% intergenic (+20/‑364) cysK/ptsH cysteine synthase A, O‑acetylserine sulfhydrolase A subunit/phosphohistidinoprotein‑hexose phosphotransferase component of PTS system (Hpr)
?NC_000913 2533409 = 3 (0.050)intergenic (+31/‑353) cysK/ptsH cysteine synthase A, O‑acetylserine sulfhydrolase A subunit/phosphohistidinoprotein‑hexose phosphotransferase component of PTS system (Hpr)

GAAGATGAAAGCTTTACCAACAAGAATATTGTGGTTATTCTACCATCATCGGGTGAGCGTTATTTAAGCACCGCATTGTTTGCCGATCTCTTCACTGAGAAAGAATTGCAACAGTAATGCCAGCTTGTTAAAAATGCGTAAAAAAGCACCTTTTTAGGTGCTTTTTTGTGGCCTGCTTCAAACTTTCGCCCCTCCTGGCATTGATTCAGCCTGTCGGAACTGGTATTTAACCAGACTAATTATTTTGATGCGCGAAATTAATCG  >  NC_000913/2533262‑2533525
                                                                                                                                        |                                                                                                                                
gaagaTGAAAGCTTTACCAACAAGAATATTGTGGTTATTCTACCATCATCGGGTGAGCGTTATTTAAGCACCGCATTGTTTGCCGATCTCTTCACTGAGAAAGAATTGCAACAGTAATGCCAGCTTGTTAAAAATGCaccttt                                                                                                                            >  1:1113864‑M1/1‑137 (MQ=255)
    aTGAAAGCTTTACCAACAAGAATATTTAGGTTATGCTACCACCATCGGGTGAGCGTTATTTAAGCACCGCATTGTTTGCCGATCTCTTCACTGAGAAAGAATTGCAACAGTAATGCCAGCTTGTTAACAATGCacctttttag                                                                                                                        >  2:733088‑M1/1‑133 (MQ=255)
    aTGAAAGCTTTACCAACAAGAATATTGTGGTTATTCTACCATCATCGGGTGAGCGTTATTTAAGCACCGCATTGTTTGCCGATCTCTTCCCTGAGAAAGAATTGCAACAGTACTGCCAGCTTGTTAAAAATGCaccttcttag                                                                                                                        >  1:790666‑M1/1‑133 (MQ=255)
    aTGAAAGCTTTACCAACAAGAATATTGTGGTTATTCTACCATCATCGGGTGAGCGTTATTTAAGCACCGCATTGTTTGCCGATCTCTTCACTGAGAAAGAATTGCAACAGTAATGCCAGCTTGTTAAAAATGCacctttttag                                                                                                                        <  2:1090167‑M1/143‑11 (MQ=255)
    aTGAAAGCTTTACCAACAAGAATATTGTGGTTATTCTACCATCATCGGGTGAGCGTTATTTAAGCACCGCATTGTTTGCCGATCTCTTCACTGAGAAAGAATTGCAACAGTAATGCCAGCTTGTTAAAAATGCacctttttag                                                                                                                        <  1:25523‑M1/143‑11 (MQ=255)
    aTGAAAGCGTTACCAACAAGAATATTGTGGTTATTCTACCATCATCGGGTGAGCGTTATTTAAGCACCGCATTGTTTGCCGATCTCTTCACTGAGAAAGAATTGCAACAGTAATGCCAGCTTGTTAAAAATGCacctttttag                                                                                                                        <  2:111615‑M1/143‑11 (MQ=255)
      gAAAGCTTTACCAACAAGAATATTGTGGTTATTCTACCATCATCGGGTGAGCGTTATTTAAGCACCGCATTGTTTGCCGATCTCTTCACTGAGCAAGAATTGCAACAGTAATGCCAGCTTGTTCAATATGCacccttttacgt                                                                                                                      >  2:95021‑M1/1‑131 (MQ=255)
           cTTTACCAACAAGAATATTGTGGTTATTCTACCATCATCGGGTGAGCGTTATTTAAGCACCGCATTGTTTGCCGATCTCTTCACTGAGAAAGAATTGCAACAGTAATGCCAGCTTGTTAAAACTGCaccttttcaggtgcttt                                                                                                                 >  1:222600‑M1/1‑126 (MQ=255)
              tACCAACAAGAATATTGTGGTTATTCTACCATCATCGGGTGAGCGTTATTTAAGCACCGCATTGTTTGCCGATCTCTTCACTGAGAAAGAATTGCAACAGTAATGCCAGCTTGTTAAAAATGCacctttttaggtgctttttt                                                                                                              <  2:304725‑M1/143‑21 (MQ=255)
              tACCAACAAGAATATTGTGGTTATTCTACCATCATCGGGTGAGCGTTATTTAAGCACCGCATTGTTTGCCGATCTCTTCACTGAGAAAGAATTGCAACAGCAATGCCAGCTTGTTAACAATGCacctttttaggtgctttttt                                                                                                              >  1:60085‑M1/1‑123 (MQ=255)
                ccAACAAGAATATTGTGGTTATTCTACCATCATCGGGTGAGCGTTATTTAAGCACCGCATTGTTTGCCGATCTCTTCACTGAGAAAGAATTGCAACAGTAATGCCAGCTTGTTAAAAATGCacctttttaggtgcttttttgt                                                                                                            <  2:761787‑M1/143‑23 (MQ=255)
                ccAACAAGAATATTGTGGTTATTCTACCATCATCGAGTGAGCGTTATTTAAGCACCGCATTGTTTGCCGCTCTCTTCACTGACAAAGAATTGCAACAGTAATGCCAGCTTCTTCAAAATGAacctttttagttgctcatttat                                                                                                            >  1:1015687‑M1/1‑121 (MQ=255)
                ccAACAAGAATATTGTGGTTATTCTACCATCATCAGGTCAACGTGATTTTTTCTCCGCATTGTTTGCCGATCTCTTCACTGAGAAAGAATTGCAACAGTAATGCCAGCTTGTTAAAAATGCacctttttaggtgcttttttgt                                                                                                            <  1:399093‑M1/143‑23 (MQ=255)
                   acaaGAATCGTGTGGTTCTTCTACCATCATCGGGTGAGCGTTATTTAAGCACCGCATTGTTTGCCGATCTCTGCACTGAGAAAGAATTGCAACAGTAATGCCAGCTTGTTAAAAATGCacctttttaggtgcttttttgtggc                                                                                                         <  1:75837‑M1/143‑26 (MQ=255)
                      ggaATATTGTGGTTATTCTACCATCATCGGGTGAGCGTTATTTAAGCACCGCATTGTTTGCCGATCTCTTCACTGAGAAAGAATTGCAACAGTAATGCCAGCTTGTTAAAAATGCacctttttaggtgcttttttgtggcctg                                                                                                      <  1:1050255‑M1/142‑29 (MQ=255)
                      aGAATAGTGTGGTTATTCTACCATCATCGGGTGAGCGTTATTTAAGCACCGCATTGTTTGCCGATCTCTTCACTGAGAAAGAATTGCAACAGTAATGCCAGCTTGTTAAAAATGCacctttttaggtgcttttttgtggcctg                                                                                                      <  2:489113‑M1/143‑29 (MQ=255)
                       gAATATTGTGGTTATTCTACCATCATCGGGTGAGCGTTATTTAAGCACCGCATTGTTTGCCGATCTCTTCACTGAGAAAGAATTGCAACAGTAATGCCAGCTTGTTAAAAATGCacctttggaggtgcttttttgtggactgc                                                                                                     >  2:1134102‑M1/1‑114 (MQ=255)
                             tgtgGTTATTCTACCATCATCGGGTGAGCGTTATTTAAGCACCGCATTGTTTGCCGATCTCTTCACTGAGAAAGAATTGCAACAGTAATGCAAGCTTGTTAAAAATGCacctttttagttgcttttttgtggcctgcttcaaa                                                                                               >  1:798719‑M1/1‑108 (MQ=255)
                                  ttATTCTATCTTCATCGGGTGAGCGTTATTTAAGCACCGCATTGTTTGCCGATCTCTTCACTGAGAAAGAATTGCAACAGTACTGCCAGCTTGTTAAAAATGCacctttttaggtgcttttttgtggcctgcttcaaactttc                                                                                          <  2:1015941‑M1/143‑41 (MQ=255)
                                              catcGGGTGAGCGTTATTTAAGCACCGCATTGTTTGCCGATCTCTTCACTGAGAAAGAATTGCAACAGTAATGCCAGCTTGTTAAAAATGCacctttttaggtgcttttttgtggcctgcttcaaactttcgcccctcctggc                                                                              >  1:1077245‑M1/1‑91 (MQ=255)
                                               atcGGGTGAGCGTTATTTAAGCACCGCATTGTTTGCCGATCTCTTCACTGAGAAAGAATTGCAACAGTAATGCCAGCTTGTTAAAAATGCacctttttaggtgcttttttgtggcctgcttcaaactttcgcccctcctggca                                                                             >  2:947896‑M1/1‑90 (MQ=255)
                                                  gggTGAGCGTTATTTAAGCACCGCATTGTTTGCCGATCTCTTCACTGAGAAAGAATTGCAACAGTAATGCCAGCTTGTTAAAAATGCacctttttaggtgcttttttgtggcctgcttccaactttcgcccctcctcgcattg                                                                          >  2:997707‑M1/1‑87 (MQ=255)
                                                     tGAGCGTTATTTACGCCCCGCATTGTTTGCCGCTCTCTTCACTGAGAAAGAATTGCAACAGTAATGCCAGCTGGTTAAAAATGCacctttgtaggtgcttttttgtggcctgcttcaaaccttcgtccctcctggcattgatt                                                                       <  1:827871‑M1/143‑60 (MQ=255)
                                                     tGAGCGTTATTTAAGCACCGCATTGTTTGCCGATCTCTTCACTGAGAAAGAATTGCAACAGTAATGCCAGCTTGTTAAAAATGCacctttttaggtgcttttttgtggcctgcttcaaactttcgcccctcctggcattgatt                                                                       <  2:335626‑M1/143‑60 (MQ=255)
                                                      gAGCGTTATTTAAGCACCGCATTGTTTGCCGATCTCTTCACTGAGAAAGAATTGCAACAGTAATGCCAGCTTGTTAAAAATGCacctttttaggtgcttttttgtggcctgcttcaaactttcgcccctcctggcattgattc                                                                      >  2:833847‑M1/1‑83 (MQ=255)
                                                             aTTTAAGCACCGCATTGTTTGCCGATCTCTTCACTGAGAAAGAATTGCAACAGTAATGCCAGCTTGTTAAAAACGCaccttttgaggtgcttttttgtggcctgcttcaaactttcgcccctcctggcattgattcagcctgt                                                               <  2:341877‑M1/143‑68 (MQ=255)
                                                              tttAAGCACCGCATTGTTTGCCGATCTCTTCACTGAGAAAGAAGTGCAACAGTAATGCCAGCTTGTTAAAAATGCacctttttaggtgcttttttgtggcctgcttcaaactttcgcccctcctggcattgattcagcctgtc                                                              <  2:720636‑M1/143‑69 (MQ=255)
                                                                tAAGCACCGCATTGTTTGCCGATCTCTTCACTGAGAAAGAATTGCAACAGTAATGCCAGCTTGTTAAAAATGCacctttttaggtgcttttttgtggcctgcttcaaactttcgcccctcctggcattgattcagcctgtcgg                                                            >  1:369154‑M1/1‑73 (MQ=255)
                                                                  aGCACCGCATTGTTTGCCGATCTCTTCACTGAGAAAGAATTGCAACAGTAATGCCAGCTTGTTAAAAATGCacctttttaggtgcttttttgtggcctgcttcaaactttcgcccctcctggcattgattcagcctgtcggaa                                                          <  2:1113864‑M1/143‑73 (MQ=255)
                                                                      ccGCAGTGTGTACCGATCTCTTCACTGAGTAAGAATTGCAACAGTAATGCCAGCTTGTTAAAAATGCacctttttaggtgcttttttgtggcctgcttcaaactttcgcccctcctggcattgattcagcctgtcggaactgg                                                      <  1:748869‑M1/143‑77 (MQ=255)
                                                                         cATTGTTTGCCGATCTCTTCACTGAGAAAGAATTGCAACAGTAATGCCAGCTTGTTAAAAATGCacctttttaggtgcttttttgtggcctgcttcaaactttcgcccctcctggccttgattcagcctgtcggacctggtat                                                   >  2:810561‑M1/1‑64 (MQ=255)
                                                                              tttGCCGATCTCTTCACTGAGAAAGAATTGCAACAGTAATGCCAGCTTGTTAAAAATGCacctttttaggtgcttttttgtggcctgcttcaaactttcgcccctcctggcattgattcagcctgtcggcactggtgtttaac                                              >  2:73086‑M1/1‑59 (MQ=255)
                                                                              tttGCCGAGCGCTTCACTGAGAAAGAATTGCTACAGTAATGCCAGCTTGTTCAAGATGCG‑‑‑‑‑‑‑‑‑‑CCTTTTTAGGTGCTTTTTTGTGGCCTGCATCAAGTTTTCGCTCCTCATGGCATTGAGTCAGCCTGTCGGAACTGGTATTTAAc                                   <  1:467502/143‑1 (MQ=255)
                                                                                tGCCGATCTCTTCACTGAGAAAGAATTGCAGCAGTAATGCCAGCTTGTTAAAAATGCacctttttaggtgcttttttgtggcctgcttctaactttcgcccctcctggcatttattcagcctgtcggaactggtatttaaaca                                            >  2:377873‑M1/1‑57 (MQ=255)
                                                                                     aTCTCTTCACTGAGAAAGAATTGCAACAGTAATGCCAGCTTGTTAAAAATGCacctttttaggtgcttttttgtggcctgcttcaaactttcgcccctcctggcattgattcagcctgtgggaactggtatttaaccagacta                                       >  2:781391‑M1/1‑52 (MQ=255)
                                                                                     aTCTCTTCACTGAGAAAGAATTGCAACAGTAATGCCAGCTTGTTAAAAATGCacctttttaggtgcttttttgtggcctgcttcaaactttcgcccctcctggcattgattcagcctgtcggaactggtatttaaccagacta                                       >  1:696573‑M1/1‑52 (MQ=255)
                                                                                        tctTCACTGAGAAAGAATTGCAACAGTAATGCCAGCTTGTTAAAAAT‑‑‑‑‑‑‑‑‑‑GCACCTTTTTAGGTGCTTTTTTGTGGCCTGCTTACAACTTTCGCCCCCCCTGCCGACCTACGGTTCTGTCGGAACTGGTATTTAACCAGACTAatt                         >  1:1121588/1‑143 (MQ=255)
                                                                                         ctTCACTGAGAAAGAATTGCAACAGTAATGCCAGCTTGTTAAAAATGCacctttttaggtgcttttttgtggcctgcttcaaactttcgcccctcctggcattgattcagcctgtcggaactggtatttaaccagactaatta                                   <  1:91097‑M1/143‑96 (MQ=255)
                                                                                          ttCACTGAGAAAGAATTGCAACAGTAATGCCAGCTTGTTAAAAATGCacctttttaggtgcttttttgtggcctgcttcaaactttcgcccctcctggcattgattcagcctgtcggaactggtatttaaccagactaattat                                  >  1:157509‑M1/1‑47 (MQ=255)
                                                                                          ttCACTGAGAAAGAATTGCAACAGTAATGCCAGCTTGTTAAAAATGCacctttttaggtgcttttgtgttgcctacttcacactttcgcccctcctgtgattgattcagcctgtcggaactggtattgaaccggtatcattag                                  >  2:108371‑M1/1‑47 (MQ=255)
                                                                                               tGAGAAAGAATTGCAACAGTAATGCCAGCTTGTTAAAAATGCacctttttaggtgcttttttgtggcctgcttcaaactttcgcccctcctggcattgattcagcctgtcggaactggtatttaaccagactaattcttttgc                             >  2:520762‑M1/1‑42 (MQ=255)
                                                                                                   aaagaaTTGCAACAGTAATGCCAGCTTGTTAAAAATGCacctttttaggtgcttttttgtggcctgcttcaaactttcgcccctcctggcattgattcagcctgtcggaactggtatttaaccagactaattattttgatgcg                         >  2:53132‑M1/1‑38 (MQ=255)
                                                                                                     agaaTTGCAACAGTAATGCCAGCTTGTTAAAAATGCacctttttaggtgcttttttgtggcctgcttcaaactttcgcccctcctggcattgattcagcctgtcggaactggtatttaaccagactaattattttgatgcgcg                       >  2:517465‑M1/1‑36 (MQ=255)
                                                                                                        gttgCAACAGTAAGGCCAGGTTGTTAAAAAGGCacctttgtaggtccttatttgtggcctgcttcaaacgttcgcccctcctggcattgattcagcctgtcggaactggtatttaagcagactaattattttgatgcgcgaaa                    <  1:74305‑M1/142‑111 (MQ=255)
                                                                                                        aTTGCAACAGTAATGCCAGCTTGTTAAAAATGCacctttttaggtgcttttttgtggcctgcttcaaactttcgcccctcctggcattgattcagcctgtcggaactggtatttaaccagactaattattttgatgcgcgaaa                    <  2:660013‑M1/143‑111 (MQ=255)
                                                                                                         ttGCAACAGTAATGCCAGCTTGTTAAAAATGCacctttttaggtgcttttttgtggcctgcttcaaactttcgcccctcctggcattgattcagcctgtcggaactggtatttaaccagactaattattttgatgcgcgaaat                   >  1:168030‑M1/1‑32 (MQ=255)
                                                                                                         ttGCAACAGTAATGCCAGCTTGTTAAAAATGCacctttttaggtgcttttttgtggcctgcttcaaactttcgcccctcctggcattgattcagcctgtcggaactggtatttaaccagactaattattttgatgcgcgaaat                   <  2:986041‑M1/143‑112 (MQ=255)
                                                                                                            cAACAGTAATGCCAGCTTGTTAAAAATGCaccttttaaggtgctttttagtggcctgcttcaaactttcgcccctcctggcattgattcagcctgtcggaactggtatttaaccagactaattattttgatgcgcgaaattac                >  1:296426‑M1/1‑29 (MQ=255)
                                                                                                             aaCAGTAATGCCAGCTTGTTAAAAATGCacctttttaggtgcttttttgtggcctgcttcaaactttcgcccctcctggcattgattcagcctgtcggaactggtatttaaccagactaattattttgatgcgcgaaattaat               >  1:316550‑M1/1‑28 (MQ=255)
                                                                                                             aaCAGTAATGCCAGCTTGTTAAAAATGCacctttttaggtgcttttttgtggcctgcttcaaacgttcgcccatcctggcattgattcagcctgtcggaactggtatttaaccagactaatgattttgatgcgggaatttaat               >  2:437757‑M1/1‑28 (MQ=255)
                                                                                                               cAGTAATGCCAGCTTGTTAAAAATGCacctttttaggtgcttttttgtggcctgcttcaaactttcgcccctcctggcattgattcagcctgtcggaactggtatttaaccagactaattattttgatgcgcgaaattaatcg             >  2:952870‑M1/1‑26 (MQ=255)
                                                                                                                aGTAATGCCAGCTTGTTAAAAATGCacctttttaggtgcttttttgtggcctgcttcaaactttcgcccctcctggcattgattcagcctgtcggaactggtatttacccagactaattattttgatgcgcgaaattaatcgt            >  1:718210‑M1/1‑25 (MQ=255)
                                                                                                                  tAATGCCAGCTTGTTAAAAATGCacctttttaggtgcttttttgtggcctgcttcaaactttcgcccctcctggcattgattcagcctgtcggaactggtgtttaaccagcctacttattttgatgcgcggaattattcgata          >  1:1063125‑M1/1‑23 (MQ=255)
                                                                                                                   aaTGCCAGCTTGTTAAAAATGCacctttttaggtgcttttttgtggcctgcttcaaactttcgcccctcctggcattgattcagcctgtcggaactggtatttaaccagactaattattttgatgcgcgaaattaatcgttac         >  2:102624‑M1/1‑22 (MQ=255)
                                                                                                                   aaTGCCAGCTGGTTAAAAATGCacctttttaggtgcttttttgtggcctgcttcaaactttcgcccctcctggcattgattcagcctgtcggaactgggatttaaccagactaattattttgatgcgcgaaattaatcgttac         <  2:596288‑M1/143‑122 (MQ=255)
                                                                                                                     tGCTAGCTTGTTAAAAATGCacctttttaggtgcttttttgtggcctgcttcaaactttcggccctcctggctttgatgcagcctgtcggaactggtatttaaccaggctatttgttttgatgtgcgaacttaagcgctacta       >  1:118052‑M1/1‑20 (MQ=255)
                                                                                                                       ccAGCTTGTTAAAAATGCacctttttaggtgcttttttgtggcctgcttcaaactttcgcccctcctggcattgattcagcctgtcggaactggtatttaaccagactaattattttgatgcgcgaaattaatcgttacagga     <  1:947820‑M1/143‑126 (MQ=255)
                                                                                                                        cAGCTTGTTAAAAATGCacctttttaggtgcttttttgtggcctgcttcaaactttcgcccctcctggcattgattcagcctgtcggaactggtatttaaccagactaattattttgatgcgcgaaattaatcgttacaggaa    >  2:1015786‑M1/1‑17 (MQ=255)
                                                                                                                         aGCTTGTTAAAAATGCacctttttaggtgcttttttgtggcctgcttcaaactttcgcccctcctggcattgattcagcctgtcggaactggtatttaaccagactaattattttgatgcgcgaaattaatcgttacaggaaa   <  2:357652‑M1/143‑128 (MQ=255)
                                                                                                                         aGCTGGTTAAAAATGCacctttttaggatcttttttgtagactgattcaaagtttcgcccctcctagaattgattcaggctgtcggaactggtatttaaccagactaattattttgatgcgcgaaattaatcgttacaggaaa   <  1:127973‑M1/143‑128 (MQ=255)
                                                                                                                          gCTTGTTAAAAATGCaccgtgttaggtgcttttttgtggcctgcttcaaactttcgcccctcctggcattgattcagcctgtcggaactggtatttaaccagactaattattttgatgcgcgaaattaatcgttacaggaaaa  <  2:369154‑M1/143‑129 (MQ=255)
                                                                                                                                        |                                                                                                                                
GAAGATGAAAGCTTTACCAACAAGAATATTGTGGTTATTCTACCATCATCGGGTGAGCGTTATTTAAGCACCGCATTGTTTGCCGATCTCTTCACTGAGAAAGAATTGCAACAGTAATGCCAGCTTGTTAAAAATGCGTAAAAAAGCACCTTTTTAGGTGCTTTTTTGTGGCCTGCTTCAAACTTTCGCCCCTCCTGGCATTGATTCAGCCTGTCGGAACTGGTATTTAACCAGACTAATTATTTTGATGCGCGAAATTAATCG  >  NC_000913/2533262‑2533525

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