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 31 [1] [0] 29 [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)29 (0.640)
+TG
26/480 0.6 98.3% 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

GCCAGAGAGGCCGCCCATAGCACCACGCCGATTACTGCGGCAATAATTTTAATCCGCTTCGGAAAGAGGAATACCGGGATCAACCACAGCGAACTGAATAACAGCGAGTCGCGAATGCCGTTAGTGCCACTATAACCACTGATGTAAATAATGGCCTGTAGCAGAGTAGAGAAAAACCAAAAGTAGAGCAGTGCCCAACCCAGGGCTTTCCAGCTAAAAAGAGGTTTAGCCTGGACTTCTGTGGAATGCATAGTAAGAACCTGTCTTGAAAAAATATCGCCGAATGTAACGACAATTCCTTAAGGATATCTGAAGGTATATTCAGAATTTGAATAAAATGCAGACAGAAATATATTGAAAACGAGGGTGTTAGAACAGAAGTATTTCAGAAAACCCTCGCGCAAAAGCACGAGGGTTTGCAGAAGAGGAAGATTAGCCGGTATTACGCATACCTGCCGCAATCCCGGCAATAGTGACCATTAA  >  NC_000913/4150015‑4150497
                                                                                                                                                                                                                                               |                                                                                                                                                                                                                                                   
gagcaaacagaatgccctcatcaccaatcgctttttgcagacgcggcagcacttcccatatggactccttgctggcagcgataatgctcgggttagttgtcaccccggcaatggggaatatgcgtgccagacgttcgacttctgcgacgttagcggtttccagatacagttccatgatgtgtccttttacaggtccagttcatgttgcagcaggctggcgatagcggcggcagaggctgTGTGGAATGCAt                                                                                                                                                                                                                                          <  2:443398‑M2/12‑1 (MQ=255)
gagcaaacagaatgccctcatcaccaatcgctttttgcagacgcggcagcacttcccatatggactccttgctggcagcgataatgctcgggttagttgtcaccccggcaatggggaatatgcgtgccagacgttcgacttctgcgacgttagcggtgtccagatacagttccatgatgtttccttttacagttccagttcatgttgcagcaggctggcgatagcgtcggcagatgctgTTTTGAATGCAt                                                                                                                                                                                                                                          <  2:61488‑M2/12‑1 (MQ=255)
         gaatgccctcatcgccaatcgctttttgcagacgctgcagcacttcccatatggactccttgctggcagcgataatgctcgggttagttgtcaccccggcaatggggaatatgcgtgccagacgttcgacttctgcgacgttagcggtgtccagatacagttccatgatgtttccttttacagtgccagttcatgttgcagcaggctggcgatagcgtcggcagaggctgTGTGGAATGCATAGTAAGAAc                                                                                                                                                                                                                                 <  2:256370‑M2/21‑1 (MQ=255)
         gaatgccctcatcaccaatcgctttttgcagacgcggcagcacttcccatatggactccttgctggcagcgataatgctcgggttagttgtcaccccggcaatggggaatatgcgtgccagacgttcgacttctgcgacgttagcggtgtccagatacagttccatgatgtttccttttacagttccagttcatgttgcagcatgctggcgatagcgtcggcagaggctgTGTGGAATGCATAGTAAGAAc                                                                                                                                                                                                                                 <  1:297189‑M2/21‑1 (MQ=255)
                  catcacccatagctttttgcagacgcggcagcacttcccatatggactccttgctggcagcgataatgctcgggttagttgtcaccccggcaatggggaatatgcgtgccagacgttcgacttctgcgacgttagcggtgtccagatacagttccatgatgtttccttttacagttccagttcatgttgcagcaggctggcgatagcgtcggcagaggctgTGTGGAATGCCTAGTAAGAACCTGTCTTg                                                                                                                                                                                                                         >  2:298036‑M2/222‑250 (MQ=255)
                                    gcagacgcggcagcacttcccatatggactccttgctggcagcgataatgctcgggttagttgtcaccccggcaatggggaatatgcgtgccagacgttcgacttctgcgacgttagcggtgtccagatacagttccatgatgtttccgtttacagttccagttcatgttgcagcaggctggcgatagcgtcggcagaggctgTGTGGAATGCATAGTAAGAACCTGTCTTGAAAAAATATCGCCGAATg                                                                                                                                                                                                       <  1:187604‑M2/47‑1 (MQ=255)
                                        actcggcagcacttcccatatggactccttgctggcagcgataatgctcgggttagttgtcaccccggcaatggggaatatgcgtgccagacgtgcgacttctgcgacgttagcggtgtccagatacagttccatgatgtttccttttacagttccagttcatgttgcatcaggctggcgatagcgtcggcagaggcggTGTGGAATGCATAGTAAGAACCTTTCGTTAAAAAAGATCGCCGAATGTAACg                                                                                                                                                                                                  <  2:271542‑M2/52‑1 (MQ=255)
                                               agcacttcccatatggactccttgctggcagcgataatgctcgggttagttgtcaccccggcaatggggaatatgcgtgccagacgttcgacttctgcgacgttagcggtgtccagatacagttccatgatgtttccttttacagttccagttcatgttgcagcaggctggcgatagcgtcggcagaggctgTGTGGAATGCATAg                                                                                                                                                                                                                                        <  1:184518‑M2/14‑1 (MQ=255)
                                               agcacttcccatatggactccttgctggcagcgataatgctcgggttagttgtcaccccggcaatggggaatatgcgtgccagacgttcgacttctgcgacgttagcggtgtccagatacagttccatgatgtttccttttacagttccagttcatgttgcagcaggctggcgatagcgtcggcagaggctgTGTGGAATGCATAg                                                                                                                                                                                                                                        >  2:184518‑M2/193‑206 (MQ=255)
                                                     tcccatatggactccttgctggcagcgataatgctcgggttagttgtcaccccggcaatggggaatatgcgtgccagacgttcgacttctgcgacgttagcggtgtccagatacagttccatgatgtttccttttacagttccagttcatgttgcagcaggctggcgatagcgtcggcagaggctgTGTGGAATGCATAGTAAGAACCTGTCTTGAAAAAATATCGCCGAATGTAACGACAATTCCTTaa                                                                                                                                                                                      <  1:10211‑M2/64‑1 (MQ=255)
                                                                  ccttgctggcagcgataatgctcgggttagttgtcaccccggcaatggggaatatgcgtgccagccgttcgacttctgcgacgttagcggtgtccagatacagttccatgatgtttccttttacagttccagttcatgttgcagcaggctggcgatagcgtcggcagaggctgTGTGGAATGCCTAGTAAGAACCTGTCTTGAAAAAATATCGCCGAATGTAACGACAATTCCTTAAGGATATCTGAAGGt                                                                                                                                                                        >  1:469716‑M2/174‑251 (MQ=255)
                                                                         ggaagcgatactgctcgggttagttggcaccccgggaagggggaatatgcgtgccagacgttcgacttctgcgacgttagcggtgtccagatacaggtccatgaggtttccttttacagttccagttcatgttgcagcaggctggcgatagcgtcggcagaggctgGGTGGAATGCATAGTAAGAACCTGTCTTGAAAAAATATCGCCGAATGTAACGACAATTCCTTAAGGAGATCTGAAGgtatag                                                                                                                                                                    <  2:77701‑M2/82‑2 (MQ=255)
                                                                         ggaagcgataatgctcgggttagttgtcaccccggcaatggggaatatgcgtgccagacggtcgacgtcgtcgacgtgagcggtgtccagatacagttccatgatgtttccttttacagtgccagttcatgttgcagcaggctgtcgatagcgtcggcagaggctgTGGGGAATGCATAGTAAGAACCGGTCATGAAAAAATATCGCCGAATGTAACGACAATTCCTGAAGGATATCTGAAGGTATATTc                                                                                                                                                                  <  2:64630‑M2/84‑1 (MQ=255)
                                                                              cgataatgctcgggttagttgtcaccccggcaatggggaatatgcgtgccagccgttcgacttctgcgacgttagcggtgtccagatacagttccatgatgtttccttttacagttccagttcatgttgcagcaggctggcgatagcgtcggcagaggctgTGTGGAATGCATAGTAAGAACCTGTCTTGAAAAAATATCGCCGAATGTAACGACAATTCCTTTAGGATATCTGAAGGTACATTCAGAAtt                                                                                                                                                            >  1:407413‑M2/162‑251 (MQ=255)
                                                                               gataatgctcgggttagttgtcaccccggcaatggggaatatgcgtgccagacgttcgacttctgcgacgttagcggtgtccagatacagttccatgatgtttccttttacagttccagttcatgttgcagcaggctggcgatagcgtcggcagaggctgTGTGGAATGCATAGTAAGAACCTGTCTTGAAAAAATATCGCCGAATGTAACGACAATTCCTTAAGGATATCTGAAGGTATATTCAGAAttt                                                                                                                                                           <  1:442973‑M2/91‑1 (MQ=255)
                                                                                         gggttagttgtcacctcggcaatggggaatatgcgtgccagacgttcgacttctgcgacgttagcggtgtccagatacagttccatgatgtttccttttacagttccagttcatgttgcagcaggctggcgatagcgtcggcagaggctgTGTGGAATGCATAGTAAGAACCTGTCTTGAAAAAATATCGCCGAATGTAACGACAATTCCTTAAGGATATCTGAAGGTATATTCAGAATTTGAATAAAATg                                                                                                                                                 <  1:429589‑M2/101‑1 (MQ=255)
                                                                                                   tcaccccggcaatggggaatatgcgtgccagacgttcgacttctgcgacgttagcggtgtccagatacagttccatgatgtttccttttacagttccagttcatgttgcagcaggctggcgatagcgtcggcagaggctgTGTGGAATGCATAGTAAGAACCTGTCTTGAAAAAATATCGCCGAATGTAACGACa                                                                                                                                                                                               >  1:24651‑M2/141‑195 (MQ=255)
                                                                                                   tcaccccggcaatggggaatatgcgtgccagacgttcgacttctgcgacgttagcggtgtccagatacagttccatgatgtttccttttacagttccagttcatgttgcagcaggctggcgatagcggcggcagaggctgTGTGGAATGCATAGTAAGAACCTGTCTTGAAAAAATATCGCCGAATGTAACGACa                                                                                                                                                                                               <  2:24651‑M2/55‑1 (MQ=255)
                                                                                                    caccccggcaatggggaatatgcgtgccagacgttcgacttctgcgacgttagcggtgtccagatacagttccatgatgtttccttttacagttccagttcatgttgcagcaggctggcgatagcgtcggcagaggctgTGTGGAATGCATAGTAAGAACCTGTCTTGAAAAAATATCGCCGAATGGAACGACAATTCCTTAAGGATATCTGAAGGTATATTCAGAATTTGAATAAAATGCAGACAGaaa                                                                                                                                       <  1:298036‑M2/111‑1 (MQ=255)
                                                                                                                   gaatatgcgttccagacgttcgacttctgcgacgttagcggtgtccagatacagttccatgatgtttccttttacagttccagttcatgttgcagcaggctggcgatagcgtcggcagaggctgTGTGGAATGCATAGTAAGAACCTGTCTTGAAAAAATATCGCCGAATGTAACGACAATTCCTTAAGGATATCTGAAGGTATATTCAGAATTTGAATAAAATGCAGACAGAAATATATTGAAAACGAgg                                                                                                                       <  1:214602‑M2/127‑1 (MQ=255)
                                                                                                                                          cttctgcgacgttagcggtgtccagatacagttccatgatgtttccttttacagttccagttcatgttgcagcaggctggcgatagcgtcggcagaggctgTGTGGAATGCATAGTAAGAACCTGTCTTGAAAAAATATCGCCGAATGTAACGACAATTCCTTAAGGATATCTGAAGGTATATTCAGAATTTGAATAAAATGCAGACAGAAATATATTGAAAACGAGGGTGTTAGAACAGAAGTATTTc                                                                                                  >  1:49943‑M2/102‑249 (MQ=255)
                                                                                                                                                   cgttagcggtgtccagatacagttccatgatgtttccttttacagttccagttcatgttgcagcaggctggcgatagcgtcggcagaggctgTGTGGAATGCATAGTAAGAACCTGTCTTGAAAAAATATCGCCGAATGTAACGACAATTCCTTAAGGATATCTGAAGGTATATTCAGAATTTGAATAAAATGCAGACAGAAATATATTGAAAACGAGGGTGTTAGAACAGAAGTATTTCAGAAAACCCTc                                                                                       <  2:468109‑M2/159‑1 (MQ=255)
                                                                                                                                                       agcggtgtccagatacagttccatgatgtttccttttacagttccagttcatgttgcagcaggctggagatagcgtcggcagaggctgTGTGGAATGACTAGTAAGAACCTGTCTTGACAAAATATCGCCGAATGTAACGACAATTCCTTAAGGATATCTGAAGGTATATTCAGAATTTGAATAAAATGCAGACAGAAATATATTGAAAACGAGGGTGTTAGAACAGAAGTATTTCAGAAAACCCTAgcgc                                                                                   >  1:311668‑M2/89‑251 (MQ=255)
                                                                                                                                                               ccagatacaggtccatgatgttgccttttacagttccagttcatgttgcagcaagctggcgatagcgtcggcagaggctgTGTGGAATGCAGAGGAAGAACCTGTCTTGAAAAAATAGCGCCGAATGTAACGACAATTTCTGAAGGATATCTGAAGGTATATTCAGAATTTGAATAAAATGGAGACAGAAAGATATTTAAAACGAGGGTGTTAGAACAGAAGTATTTCAGAAAACCCTCTCGCAAAAGCAc                                                                           <  2:469716‑M2/171‑1 (MQ=255)
                                                                                                                                                                                      ccttttacagttccagttcatgttgcagcaggctggcgatagcgtcggcagaggctgTGTGGAATGCATAGTAAGAACCTGTCTTGAAAAAATATCGCCGAATGTAACGACAATTCCTTAAGGATATCTGAAGGTATATTCAGAATTTGAATAAAATGCAGACAGAAATATATTGAAAACGAGGGTGTTAGAACAGAAGTATTTCAGAAAACCCTCGCGCAAAAGCACGAGGGTTTGCAGAAGAGGAAGa                                                     <  1:295972‑M2/193‑1 (MQ=255)
                                                                                                                                                                                                                   aggctggcgatagcgtcggcagatgctgTGTGGGCTGCATAGTAAGAACCTGTCTTGACAAAATATCGCCGAATGTAACGACAATTCCTTAAGGATATCTGAAGGTATATTCAGAATTTGAATAAAATGCAGACAGAACTATATTGAAAACGAGGGTGTGAGAACAGAAGTATTTCAGAAAACCCTCGCGCAAAAGCACGAGGTTTTGCAGAAGAGGAAGATAAGCCGGTATTACGCATACCTGCCGCAAt                       >  2:408846‑M2/29‑251 (MQ=255)
                                                                                                                                                                                                                        ggcgatggcgtcggcggaggctgTGTGGAATGCATAGAAAGAACCGGTCTTGAAAAAATTTCGCCGAATGTAACGACAATTTCTTAAGGATATCTGAAGGTATATTCAGAATTTGAATAAAATGCAGACAGAAATATATTGAAAACGAGGGTGTTAGAACAGAAGTATTTCATAAAACCCTCGCGCAAAAGCACGCGGGTTTGCAGAAGAGGAAGATTAGCCGGTATTACGCATACCTGCCGCAATccc                    <  2:452830‑M2/226‑1 (MQ=255)
                                                                                                                                                                                                                        ggcgatcgcgtcggcagaggatgTGTGGAATGCATAGTACGAACCTGTCTTGAAACAATATCGCCGAATGTAACGACAATTCCTTAAGGATCTCTGAAGGTATATTCAGAATTTGAATAAAATGCAGACAGAAATATATTGAAAACGAGGGTGTTAGAACAGAAGTATTTCAGAAAACCCTCGCGCAAAAGCACGAGGGTTTGCAGAAGAGGAAGATTAGCCGGTATTACGCATACCTGCAGCAATCCCgg                  >  2:301565‑M2/24‑251 (MQ=255)
                                                                                                                                                                                                                                        gaggctgTGTGGAATGCATAGTAAGAACCTGTCTTGAAAAAAGATCGCCGAATGTAACGACAATTCCTTAAGGATATCTGAAGGTATATTCAGAATTTGAATAAAATGCAGACAGAAATATATTGAAAACGAGGGTGTTAGAACAGAAGTTTTTCAGAAAACCCTCGCGCAAAAGCACGAGGGTTTGCAGAAGAGGAAGATTATCCGGTGTTACGCATACCTGCCGCAATCCCGGCAATAGTGACCATTaa  <  1:279334‑M2/244‑1 (MQ=255)
                                                                                                                                                                                                                                               |                                                                                                                                                                                                                                                   
GCCAGAGAGGCCGCCCATAGCACCACGCCGATTACTGCGGCAATAATTTTAATCCGCTTCGGAAAGAGGAATACCGGGATCAACCACAGCGAACTGAATAACAGCGAGTCGCGAATGCCGTTAGTGCCACTATAACCACTGATGTAAATAATGGCCTGTAGCAGAGTAGAGAAAAACCAAAAGTAGAGCAGTGCCCAACCCAGGGCTTTCCAGCTAAAAAGAGGTTTAGCCTGGACTTCTGTGGAATGCATAGTAAGAACCTGTCTTGAAAAAATATCGCCGAATGTAACGACAATTCCTTAAGGATATCTGAAGGTATATTCAGAATTTGAATAAAATGCAGACAGAAATATATTGAAAACGAGGGTGTTAGAACAGAAGTATTTCAGAAAACCCTCGCGCAAAAGCACGAGGGTTTGCAGAAGAGGAAGATTAGCCGGTATTACGCATACCTGCCGCAATCCCGGCAATAGTGACCATTAA  >  NC_000913/4150015‑4150497

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