Predicted mutation
evidence seq id position mutation annotation gene description
MC JC NC_000913 4,139,771 10483 bp→TG [ptsA][eptC] [ptsA], frwC, frwB, pflD, pflC, frwD, yijO, [eptC]

Missing coverage evidence...
   seq id start end size ←reads reads→ gene description
* * ÷ NC_000913 4139771 4150253 10483 36 [1] [0] 34 [ptsA]–[eptC] [ptsA],frwC,frwB,pflD,pflC,frwD,yijO,[eptC]

New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913 = 41397701 (0.020)34 (0.650)
+TG
32/482 0.5 98.6% coding (2452/2502 nt) ptsA putative PTS multiphosphoryl transfer protein PtsA
?NC_000913 4150254 = 0 (0.000)coding (12/1734 nt) eptC phosphoethanolamine transferase EptC

CGGCTCATGGTCTGAGCAAACAGAATGCCCTCATCACCAATCGCTTTTTGCAGACGCGGCAGCACTTCCCATATGGACTCCTTGCTGGCAGCGATAATGCTCGGGTTAGTTGTCACCCCGGCAATGGGGAATATGCGTGCCAGACGTTCGACTTCTGCGACGTTAGCGGTGTCCAGATACAGTTCCATGATGTTTCCTTTTACAGTTCCAGTTCATGTTGCAGCAGGCTGGCGATAGCGTCGGCAGAGGCGGCGTTAACCAGCGCGTTACGGAATTCTTCGTGCATGATGCGGCGAGCGAGGCGCGAGAAAATGCGCATATGCTGATCGCCCGCAGCGTGTTTGTTCAGGGTTAACATAATGATGAATTGCGCTTCATCATCGCCCCAGCGCACCGGCGCTTGCAGACGCGCCACGCTGATGGTGGATTGCTCAATGTGTTCTGATTTGCTGTGTGGAATGGCAAAACTAAAGCCCAGACCGGTAGAGAAAACGGCCTCG  >  NC_000913/4139521‑4140020
                                                                                                                                                                                                                                                         |                                                                                                                                                                                                                                                          
cGGCTCATGGTCTGAGCAAACAGAATGCCCTCATCACCAATCGCTTTTTGCAGACGCGGCAGCACTTCCCATATGGACTCCTTGCTGGCAGCGATAATGCTCGGGTTAGTTGTCACCCCGGCAATGGGGAATATGCGTGCCAGACGTTCGACTTCTGCGACGTTAGCGGTGTCCAGATACAGTTCCATGATGTTTCCTTTTACAGTTCCAGTTCATGTTGCAGCAGGCTTGCGATAGCGTCGGCAGAGGCt                                                                                                                                                                                                                                                           <  1:424954‑M1/251‑2 (MQ=255)
  gCTCATGGTCTGAGCAAACAGAATGCCCTCATCACCAATCGCTTTTTGCAGACGCGGCAGCACTTCCCATATGGACTCCTTGCTGGCAGCGATAATGCTCGGGTTAGTTGTCACCCCGGCAATGGGGAATATGCGTGCCAGACGTTCGACTTCTGCGACGTTAGCGGTGTCCAGATACAGTTCCATGATGTTTCCTTTTACAGTTCCAGTTCATGTTGCAGCAGGCTGGCGATAGCGTCGGCAGAGGCtgt                                                                                                                                                                                                                                                         <  2:314585‑M1/251‑4 (MQ=255)
  gCTCATGGTCTGAGCAAACAGAATGCCCTCATCACCAATCGCTTTTTGCAGACGCGGCAGCACTTCCCATATGGACTCCTTGCTGGCAGCGATAATGCGCGGGTTAGTTGTCACCCCGGCAATGGGGAATATGCGTGCCAGACGTTCGACTTCTGCGACGTTAGCGGTGTCCAGATACAGTTCCATGATGTTTCCTTTTACAGTTCCAGTGCATGTTGCAGCAGGCTGGCGATAGCGTCGGCAGAGGCtgt                                                                                                                                                                                                                                                         <  1:279812‑M1/251‑4 (MQ=255)
     cATGGTCTGAGCAAACAGAATGCCCTCATCACCAATCGCTTTTTGCAGACGCGGCAGCACTTCCCATATGGACTCCTTGCTGGCAGCGATAATGCTCGGGTTAGTTGTCACCCCGGCAATGGGGAATATGCGTGCCAGACGTTCGACTTCTGCGACGTTAGCGGTGTCCAGATACAGTTCCATGATGTTTCCTTTTACAGTTCCAGTTCATGTTGCAGCAGGCTGGCGATAGCGGCGGCAGAGGCtgtgtg                                                                                                                                                                                                                                                      <  1:182328‑M1/251‑7 (MQ=255)
             gAGCAAACAGAATGCCCTCATCACCAATCGCTTTTTGCAGACGCGGCAGCACTTCCCATATGGACTCCTTGCTGGCAGCGATAATGCTCGGGTTAGTTGTCACCCCAGCAATGGGGAATATGCGTGCCAGACGTTCGACTTCTGCGACGTTATCGTTGTCCATATACAGTTCCATGATGTTTCCTTGTACAGTTCCAGTTCATGTTGCAGCAGGCTGGCGATAGCGTCGGCAGAGGCtgtgtggaatg                                                                                                                                                                                                                                                 <  1:297602‑M1/248‑12 (MQ=255)
                cAAACAGAATGCCCTCATCACCAATCGCTTTTTGCAGACGGGGCAGCCCTTCCCATATGGACTCCTTGCTGGCAGCGATAATGCGCGGGTTAGTTGTCACCCCGGCAATGGGGAATATGCGTGCCAGACGTTCGACTTCTGCGACGTTAGCGGTGTCCAGATACAGTTCCATGATGTTTCCTTTTACAGTGCCAGTTCATGTTGCAGCAGGCTGGCGATAGCGTCGGCAGAGGCtgggtggaatgcatagt                                                                                                                                                                                                                                           <  2:212770‑M1/251‑18 (MQ=255)
                      gaatggactcctatccactcgctttttgcagacacggTAGCACTTACCCTATGGACTCCTTGCTGCCAGCGATAATGCTCGGGTTAGTTGTCACCCCGGCAATGGGGAATATGCGTGCCAGCCGTTCGACTTCTGCGACGTTAGCGGTGTCCAGATACAGCTCCATGATGTTTCCTTTTACAGTTCCAGTTCATGTTGCAGCAGGCTGGCGATCGCGTCGGCAGAGGCtgtgtggtatgcatagtaagaac                                                                                                                                                                                                                                     >  2:463470‑M1/19‑228 (MQ=255)
                           cccTCATCACCAATCGCTGTTTGCAGACGCGGCAGCACTTCCCATATGGACTCCTTGCTGGTAGCGATAATGCTCGGGTTAGTTGTCACCCCGGCAATGGGGAAGATGCGTGCCAGACGTTCGACTTCTGCGACGTTAGCGGTGTCCAGATACAGTTCCATGATGTTTCCTTTGGCAGTTCCTGTTCATGTTGCAGCAGGCTGGCGATAGCGGCGGCAGAGGCGGCGTTAACCAGCGCGTGGCGGAAttct                                                                                                                                                                                                                                <  1:396260/251‑1 (MQ=255)
                               catcaCCAATCGCTTTTTGCAGACGCGGCAGCACTTCCCATATGGACTCCTTGCTGGCAGCGATAATGCTCGGGTTAGTTGTCACCCCGGCAATGGGGAATATGCGTGCCAGACGTTCGACTTCTGCGACGTTAGCGGTGTCCAGATACAGTTCCATGATGTTTCCTTTTACAGTTCCAGTTCATGTTGCAGCAGGCTGGCGATAGCGTCGGCAGAGGCtgtgtggaatgcatagtaagaacctgtcttg                                                                                                                                                                                                                             >  1:70314‑M1/1‑219 (MQ=255)
                                               ttGCAGACGCGGCAGCACTTCCCATATGGACTCCTTGCTGGCAGCGATAATGCTCGGGTTAGTTGTCACCCCGGCAATGGGGAATATGCGTGCCAGACGTTCGACTTCTGCGACGTTAGCGGTGTCCAGATACAGTTCCATGATGTTTCCTTTTACAGTTCCAGTTCATGTTGCAGCAGGCTGGCGATAGCGTCGGCAGAGGCtgtgtggaatgcatagtaagaacctgtcttgaaaaaatatcg                                                                                                                                                                                                                  <  2:312919‑M1/245‑43 (MQ=255)
                                               ttGCAGACGCGGCAGCACTTCCCATATGGACTCCTTGCTGGCAGCGATAATGCTCGGGTTAGTTGTCACCCCGGCAATGGGGAATATGCGTGCCAGACGTTCGACTTCTGCGACGTTAGCGGTGTCCAGATACAGTTCCATGATGTTTCCTTTTACAGTTCCAGTTCATGTTGCAGCAGGCTGGCGATAGCGTCGGCAGAGGCtgtgtggaatgcatagtaagaacctgtcttgaaaaaatatcg                                                                                                                                                                                                                  >  1:312919‑M1/1‑203 (MQ=255)
                                                     aCGCGGCAGCACTTCCCATATGGACTCCTTGCTGGCAGCGATAATGCTCGGGTTAGTTGTCACCCCGGCACTGGGGAATATGCGTGCCAGACGTTCGACTTCTGCGACGTTAGCGGTGTCCAGATACAGTTCCATGATGTTTCCTTTTACAGTTCCAGTTCATGTTGCAGCAGGCTGGCGATAGCGTCGGCAGAGGCtgtgtggaatgcatagtaagaccctgtcttgaaaaaatatcgccgaatgtaacg                                                                                                                                                                                                      >  2:430582‑M1/1‑197 (MQ=255)
                                                         ggCAGCACTTCCCATATGGACTCCTTGCTGGCAGCGATAATGCTCGGGTTAGTTGTCACCCCGGCAATGGGGAATATGCGTGCCAGACGTTCGACTTCTGCGACGTTAGCGGTGTCCAGATACAGTTCCATGATGTTTCCTTTTACAGTTCCAGTTCATGTTGCAGCAGGCTGGCGATAGCGTCGGCAGAGGCtgtgtggaatgcatagtaagaacctgtctttaaaaaatatcgccgaatgtaacgacaa                                                                                                                                                                                                  <  2:207264‑M1/251‑59 (MQ=255)
                                                                                      ggCAGCGATAATGCTCGGGTTAGTTGTCACCCCGGCAATGGGGAATATGCGTGCCAGACGTTCGACTTCTGCGACGTTAGCGGTGTCCAGATACAGTTCCATGATGTTTCCTTTTACAGTTCCAGTTCATGTTGCAGCAGGCTGGCGATAGCGTCGGCAGAGGCtgtgtggaatgcatagtaagaacctgtcttgaaaaaatatcgccgaatgtaccgacaattccttaaggatatctgaaggtatattca                                                                                                                                                                     <  1:445935‑M1/251‑88 (MQ=255)
                                                                                      ggCAGCGATAATGCTCGGGTTAGTTGTCACCCCGGCAATGGGGAATATGCGTGCCAGACGTTCGACTTCTGCGACGTTAGCGGTGTCCAGATACAGTTCCATGATGTTTCCTTTTACAGTTCCAGTTCATGTTGCAGCAGGCTGGCGATAGCGTCGGCAGAGGCtgtgtggaatgcatagtaagaacctgtcttgaaaaaatatcgccgaatgtaacgacaattccttaaggatatctgaaggtatattc                                                                                                                                                                      <  1:5963‑M1/250‑87 (MQ=255)
                                                                                           cGATAATGCTCGGGTTAGTTGTCACCCCGGCAATGGGGAATATGCGTGCCAGACGTTCGACTTCTGCGACGTTAGCGGTGGCCAGATACAGTTCCATGATGTTTCCTTTTACAGTTCCAGTTCATGTTGCAGCAGGCTGGGGATAGCGTCGGCAGAGGCtgtgtggaaggcagagtaagaacctgtcttgaaaaaatagcgccgaatgcaacgagaattccttaaggatatctgaaggtatattcagaatt                                                                                                                                                                <  1:463470‑M1/251‑93 (MQ=255)
                                                                                               aaTGCTCGGGTTAGTTGTCACCCCGGCAATGGGGAATATGCGTGCCAGACGTTCGACTTCTGCGACGTTAGCGGTGTCCAGATACAGTTCCATGATGTTTCCTTTTACAGTTCCAGTTCATGTTGCAGCAGGCTGGCGATAGCGTCGGCAGAGGCtgtgtggaatgcatagtaagaacctgtcttgaaaaaatatcgccgaatgtaacgacaattccttaaggatatctgaaggtatattcagaatttga                                                                                                                                                             >  2:307291‑M1/1‑155 (MQ=255)
                                                                                                  gCTCGGGTTAGTTGTCACCCCGGCAATGGGGAATATGCGTGCCAGACGTTCGACTTCTGCGACGTTAGCGGTGTCCAGATACAGTTCCATGATGTTTCCTTTTACAGTGCCAGTTCATGTTGCAGCATTCTGGCGATAGCGTCGGCAGAGGCtgtggggaatgcatagtaagaacctgtcttgaaaaaatatcgccgaatgtaactacaattccttaaggatatctgaaggtatattcagaatttgaata                                                                                                                                                          <  1:444166‑M1/250‑99 (MQ=255)
                                                                                                      gggTTAGTTGTCACCCCGGCAATGGGGAATATGCGTGCCAGACGTTCGACTTCTGCGACGTTAGCGGTGTCCAGATACAGTTCCATGATGTTTCCTTTTACAGTTCCAGTTCATGTTGCAGCAGGCTGGCGATAGCGTCGGCAGAGGCtgtgtggaatgcatagtaagaacctgtcttgaaaaaatatcgccgaatgtaacgacaattccttaaggatatctgaaggtatattcagaatttgaataaaatg                                                                                                                                                     <  2:70314‑M1/251‑104 (MQ=255)
                                                                                                          tAGTTGTCACCCCGGCAATGGGGAATATGCGTGCCAGACGTTCGACTTCTGCGACGTTAGCGGTGTCCAGATACAGTTCCATGATGTTTCCTTTTACAGTTCCAGTTCATGTTGCAGCAGGCTGGCGATAGCGTCGGCAGAGGCtgtgtggaatgcatagtaagaacctgtcttgaaaaaatatcgccgaatgtaacgacaattccttaaggatatctgaaggtatattcagaatttgaataaaatgcag                                                                                                                                                  <  1:375843‑M1/250‑107 (MQ=255)
                                                                                                                   ccccGGCAATGGGGAATATGCGTGCCAGACGTTCGACTTCTGCGACGTTAGCGGTGTCCAGATACAGTTCCATGATGTTTCCTTTTACAGTTCCAGTTCATGTTGCAGCAGGCTGGCGATAGCGTCGGCAGAGGCtgtgtggaatgcatagtaagaacctgtcttgaaaaaatatcgccgaatgtaacgacaattccttaaggatatctgaaggtatattcagaatttgaataaaatgcagacagaaatat                                                                                                                                        >  1:138231‑M1/1‑135 (MQ=255)
                                                                                                                        gCAATGGGGAATATGCGTGCCAGACGTTCGACTTCTGCGACGTTAGCGGTGTCCAGATACAGTTCCATGATGTTTCCTTTTACAGTTCCAGTTCATGTTGCAGCAGGCTGGCGATAGCGTCGGCAGAGGCtgtgtggaatgcatagtaagaac                                                                                                                                                                                                                                     <  1:67647‑M1/153‑24 (MQ=255)
                                                                                                                        gCAATGGGGAATATGCGTGCCAGACGTTCGACTTCTGCGACGTTAGCGGTGTCCAGATACAGTTCCATGATGTTTCCTTTTACAGTTCCAGTTCATGTTGCAGCAGGCTGGCGATAGCGTCGGCAGAGGCtgtgtggaatgcatagtaagaac                                                                                                                                                                                                                                     >  2:67647‑M1/1‑130 (MQ=255)
                                                                                                                         cAATGGGGAATATGCGTGCCAGACGTTCGACTTCTGCGACGTTAGCGGTGTCCAGATACAGTTCCATGATGTTTCCTTTTACAGTTCCAGTTCATGTTGCAGCAGGCTGGCGATAGCGTCGGCAGAGGCtgtgtggaatgcatagtaagaacctgtcttgaaaaaatatcgccgaatgtaacgacaattccttaaggatatctgaaggtatattcagaatttgaataaaatgcagacagaaatatattgaa                                                                                                                                  <  1:307291‑M1/251‑123 (MQ=255)
                                                                                                                                         tGCCAGACGTTCGACTTCTGCGACGTTAGCGGTGTCCATATACAGTTCCATGATGTTTCCTTTTATAGTTCCAGTTCATGTTGCAGCAGGCTGGCGATAGCGGCTGCAGAGGCtgtgtggaatgcatagtaagaacctgtcgtgaaaaaatatcgccgaatgtaacgacaattccttaaggatatctgaaggtatattcagaatttgaataaaatgcagacagaaatatattgaaaacgagggtgttaga                                                                                                                   <  2:258862‑M1/250‑138 (MQ=255)
                                                                                                                                         tGCCAGACGTTCGACTTCTGCGACGTTAGCGGTGTCCAGATACAGTTCCATGATGTTTCCTTTTACAGTTCCAGTTCATGTTGCAGCAGGCTGGCGATAGCGTCGGCAGAGGCtgtgtggaatgcatagtaagaacctgtcttgaaaaaatatcgccgaatgtaacgacaattccttaaggatatctgaaggtatattcagaatttgaataaaatgcagacagaaatatattgaaaacgagggtgttaga                                                                                                                   >  2:112437‑M1/1‑113 (MQ=255)
                                                                                                                                                     gACTTCTGCGACGTTAGCGGTGTCCAGATACAGTTCCATGATGTTTCCTTTTACAGTTCCAGTTCATGTTGCAGCAGGCTGGCGATAGCGTCGGCAGAGGCtgtgtggaatgcatagtaagaacctgtcttgaaaaaatatcgccgaatgtaacgacaattccttaaggatatctgaaggtatattcagaatttgaataaaatgcagacagaaatatattgaaaacgagggtgttagaacagaagtatt                                                                                                        >  1:4512‑M1/1‑101 (MQ=255)
                                                                                                                                                     gACTTCTGCGACGTTAGCGGTGTCCAGATACAGTTCCATGATGTTTCCTTTTACAGTTCCAGTTCATGTTGCAGCAGGCTGGCGATAGCGTCGGCAGAGGCtgtgtggaatgcatagtaagaacctgtcttgaaaaaatatcgccgaatgtaacgacaattccttaaggatatctgaaggtatattcagaatttgaataaaatgcagacagaaatatattgaaaacgagggtgttagaacagaagtatt                                                                                                        <  2:4512‑M1/249‑149 (MQ=255)
                                                                                                                                                          cTGCGACGCTAGCGGAGTCAAGATACAGTTCCATGATGTTTCCTTTTACAGTTCCAGTTCATGTTGCAGCAGGCTGGCGATAGCGTCGGCAAAGGCtgtgtggaatgaaaagtaagaacctgtcttgaaaaaatatcgccgaatgtaacgacaattccttaaggatatctgaaggtatattcagaatttgattaaaatgctgacagaaatatattgaaaacgagggtgttagaacagaagtatttcagaa                                                                                                  >  2:41933‑M1/1‑96 (MQ=255)
                                                                                                                                                                       ggTGTCCAGATACAGTTCCATGATGTTTCCTTTTACAGTTCCAGTTCATGTTGCAGCAGGCTGGCGATAGCGTCGGCAGAGGCtgtgtggaatgcatagtaagaacctgtcttgaaaaaatatcgccgaatgtaacgacaattccttaaggatatctgaaggtatattcagaatttgaataaaatgcagacagaaatatattgaaaacgagggtgttagaacagaagtatttcagaaaaccctcgcgcaa                                                                                     >  2:376125‑M1/1‑83 (MQ=255)
                                                                                                                                                                         tgtCCAGATACAGTTCCCTGATGTTTCCTTTTACAGTTCCAGTTCATGTTGCAGCAGGCTGGCGATAGCGTCGGCAGAGGCtgtgtggaatgcatagtaagaacctgtcttgaaaaaatatcgccgaatgtaacgacaattccttaaggatatctgaaggtatattcagaatttgaataaaatgcagacagacatatattgaaaacgagggtgttagaacagaagtatttcagaaaaccctcgcgcaaa                                                                                    >  2:280751‑M1/1‑81 (MQ=255)
                                                                                                                                                                                      ttCCATGATGTTTCCTTTTACAGTTCCAGTTCATGTTGCAGCAGGCTGGCGATAGCGTCGGCAGAGGCtgtgtggaatgcatagtaagaacctgtcttgaaaaaatatcgccgaatgtaacgacaattc                                                                                                                                                                                               >  1:120598‑M1/1‑68 (MQ=255)
                                                                                                                                                                                      ttCCATGATGTTTCCTTTTACAGTTCCAGTTCATGTTGCAGCAGGCTGGCGATAGCGTCGGCAGAGGCtgtgtggaatgcatagtaagaacctgtcttgaaaaaatatcgccgaatgtaacgacaattc                                                                                                                                                                                               <  2:120598‑M1/129‑62 (MQ=255)
                                                                                                                                                                                                                                 ggCTGGCGATAGCGTCGGAAGAGGCtgtgtggaatgcatagtaagaacctgtcttgaaaaaatatcgccgaatgtaacgacaattccttaaggatatctgaaggtatattcagaatttgaataaaatgcagacagaaatatattgaaaacgagggtgttagaacagaagtatttcagaaaaccctcgcgcaaaagcacgagggtttgcagaagaggaagattagccggtattacgcatacctgccgcaatc                          >  2:322864‑M1/1‑25 (MQ=255)
                                                                                                                                                                                                                                        gATAGCGTCGGCAGAGGCtgtgtggaatgcatagtaagaacctgacttgaaaaaatatcgccgaatggaacaacaattccttaaggatatcggagggtatattcagaatttgaaaaaaatgcagacagaaatatattgaaaacgaggttgttagaacagaagtatttcagaaaaccctcgcgcaaaagcacgagggttttcagaagaggaaaattagccggtattacgcatacctgccgcaatcccggcaa                   <  1:211161‑M1/251‑234 (MQ=255)
                                                                                                                                                                                                                                                         ctgtgtggaatgcatagtaagaacctgtcttgaaaaaatatcgccgaatgtaacgacaattccttaaggatatctgaaggtatattcagaatttgaataaaatgcagacagaaatatattgaaaacgagggtgttagaacagaagtatttcagaaaaccctcgcgcaaaagcacgagggtttgcagaagaggaagattagccggtattacgcatacctgccgcaatcccggcaatagtgaccattaacgct  >  1:399108‑M1/1‑1 (MQ=255)
                                                                                                                                                                                                                                                         |                                                                                                                                                                                                                                                          
CGGCTCATGGTCTGAGCAAACAGAATGCCCTCATCACCAATCGCTTTTTGCAGACGCGGCAGCACTTCCCATATGGACTCCTTGCTGGCAGCGATAATGCTCGGGTTAGTTGTCACCCCGGCAATGGGGAATATGCGTGCCAGACGTTCGACTTCTGCGACGTTAGCGGTGTCCAGATACAGTTCCATGATGTTTCCTTTTACAGTTCCAGTTCATGTTGCAGCAGGCTGGCGATAGCGTCGGCAGAGGCGGCGTTAACCAGCGCGTTACGGAATTCTTCGTGCATGATGCGGCGAGCGAGGCGCGAGAAAATGCGCATATGCTGATCGCCCGCAGCGTGTTTGTTCAGGGTTAACATAATGATGAATTGCGCTTCATCATCGCCCCAGCGCACCGGCGCTTGCAGACGCGCCACGCTGATGGTGGATTGCTCAATGTGTTCTGATTTGCTGTGTGGAATGGCAAAACTAAAGCCCAGACCGGTAGAGAAAACGGCCTCG  >  NC_000913/4139521‑4140020

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 17 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: