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 40 [0] [0] 40 [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 = 41397700 (0.000)40 (0.660)
+TG
38/490 0.3 100% 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

ACAGCGCCGCCAGAGAGGCCGCCCATAGCACCACGCCGATTACTGCGGCAATAATTTTAATCCGCTTCGGAAAGAGGAATACCGGGATCAACCACAGCGAACTGAATAACAGCGAGTCGCGAATGCCGTTAGTGCCACTATAACCACTGATGTAAATAATGGCCTGTAGCAGAGTAGAGAAAAACCAAAAGTAGAGCAGTGCCCAACCCAGGGCTTTCCAGCTAAAAAGAGGTTTAGCCTGGACTTCTGTGGAATGCATAGTAAGAACCTGTCTTGAAAAAATATCGCCGAATGTAACGACAATTCCTTAAGGATATCTGAAGGTATATTCAGAATTTGAATAAAATGCAGACAGAAATATATTGAAAACGAGGGTGTTAGAACAGAAGTATTTCAGAAAACCCTCGCGCAAAAGCACGAGGGTTTGCAGAAGAGGAAGATTAGCCGGTATTACGCATACCTGCCGCAATCCCGGCAATAGTGACC  >  NC_000913/4150007‑4150492
                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                              
catggtctgagcaaacagaatgccctcatcaccaatcgctttttgcagacgcggcagcacttcccatatggactccttgctggcagcgataatgctcgggttagttgtcaccccggcaatggggaatatgcgtgccagacgttcgacttctgcgacgttagcggtgtccagatacagttccatgatgtttccttttacagttccagttcatgttgcagcaggctggcgatagcgtcggcagaggctgtgt                                                                                                                                                                                                                                              <  1:376929‑M2/3‑1 (MQ=255)
             aacagaatgccctcatcaccaatcgctttttgcagacgcggcagcacttcccatatggactccttgctggcagcgataatgctcgggttagttgtcaccccggcaatggggaatatgcgtgccagacgttcgacttctgcgacgttagcggtgtccagatacagttccatgatgtttccttttacagttccagttcatgttgcagcaggctggcgatagcgtcggcagaggctgTGTGGAATGCATAGTaa                                                                                                                                                                                                                                <  1:29704‑M2/17‑1 (MQ=255)
                           atcaacaatcgcttttagcagacgaggcagcacttcccatatggactcattgttggcagcgataatgctcgggttagttgtcaccccggcaatggggaatatgcgttccagacgttcgacttctgagacgttagcggtgtccagataaagttccatgatgtttccttttacaggtaccgttcatgttgcagcaggagggcgatagcgtaggcagaggctgTGTGGAATGCAAAGGATGAACCTGTCTTgac                                                                                                                                                                                                                  >  2:272448‑M2/221‑250 (MQ=255)
                                           tgcagacgcggcagcacttcccatatggactccttgctggcagcgataatgctcgggttagttgtcaccccggcaatggggaatatgcgtgccagacgttcgacttctgcgacgttagcggtgtccagatacagttccatgatgtttccttttacagttccagttcatgttgcagcaggctggcgatagcgtcggcagaggctgTGTGGAATGCATAGTAAGAACCTGTCTTGAAAAAATATCGCCGAAt                                                                                                                                                                                                   <  1:253525‑M2/46‑1 (MQ=255)
                                                             tcccatatggactccttgctggcagcgataatgctcgggttagttgtcaccccggcaatggggaatatgcgtgccagacgttcgacttctgcgacgttagcggtgtccagatacagttccatgatgtttccttttacagttccagttcatgttgcagcaggctggcgatagcgtcggcagaggctgTGTGGAATGCATAGTAAGAACCTGTCTTGAAAAAATATCGCCGAATGTAACGACAATTCCTTa                                                                                                                                                                                  >  2:68092‑M2/187‑249 (MQ=255)
                                                                 atatggactccttgctggcagcgataatgctcgggttagttgtcaccccggcaatggggaatatgcgtgccagacgttcgacttctgcgacgttagcggtgtccagatacagttccatgatgtttccttttacagttccagttcatgttgcagcaggctggcgatagcgtcggcagaggctgTGTGGAATGCATAGTAAGAACCTGTCTTGAAAAAATATCGCCGAATGTAGCGACAATTCCTTAAGGat                                                                                                                                                                             >  2:331909‑M2/183‑250 (MQ=255)
                                                                     ggactctttgttggctgcgataatgctcgggttagttgtcgccccggcaatggggaatatgcgtgccatacgttcgacttctgcgtcgttagcggtgtccagatacagttccatgatgtttccttttacatttccagttcatgttgcagcaggctggcgatagcgtcggcagaggctgTGTGGAATGCATAGTAATAACCTGTCTTGAAAAAATATCGCCGAATGTAACGACAATTCCTTAAGGATATCt                                                                                                                                                                         <  2:392886‑M2/72‑1 (MQ=255)
                                                                             tgctggcagcgataatgctcgggttagttgtcaccccggcaatggggaatatgcgtgccagacgttcgacttctgcgacgttagcggtgtccagatacagttccatgatgtttccttttacagttccagttcatgttgcagcaggctggcgatagcgtcggcagaggctgTGTGGAATGCATAGTAAGAACCTGTCTTGAAAAAATATCGCCGAATGTAACGACAATTCCTTAAGGATATCTGAAGGtat                                                                                                                                                                 <  1:258930‑M2/80‑1 (MQ=255)
                                                                             ggctggaggccataaggctcgggttaggtgtcaccccggcaatggggaataggcgtgccagacgttcgacttctgcgacgttagcggtgtccagatacagttccatgatgtttccttttacagttccagttcattgtgcagcaggctggcgatagcgtctgcagaggctgTGTGGAATGCATAGTAAGAACCTGTCTTGAAAAAATATCGCCGAATGTAACGGTAATTCCTTAAGTATATCTGAAGGtat                                                                                                                                                                 <  2:128439‑M2/80‑1 (MQ=255)
                                                                                   cagcgataatgctcgggttagttgtcaccccggcaatggggaatatgcgtgccagacgttcgacttctgcgacgttagcggtgtccagatacagttccatgatgtttccttttacagttccagttcatgttgcagcaggctggcgatagcgtcggcagaggctgTGTGGAATGCATAGTAAGAACCTGTCTTGAAAAAATATCGCCGAATGTAACGACAATTCCTTAAGGATATCTGAAGGTATATTc                                                                                                                                                             >  1:318598‑M2/165‑248 (MQ=255)
                                                                                       gataatgctcgggttagttgtcaccccggcaatggggaatatgcgtgccagacgttcgacttctgcgacgttagcggtgtccagatacagttccatgatgtttccttttacagttccagttcatgttgcagcaggctggcgatagcgtcggcagaggctgTGTGGAATGCATAGTAATAACCTGTCTTGAAAAAATATCGCCGAATGTAACGACAATTCCTTAAGGATATCTGAAGGTATATTCAGAAttt                                                                                                                                                      <  1:354710‑M2/91‑1 (MQ=255)
                                                                                           atgctcgggttagttgtcaccccggcaatggggaatatgcgtgccagacgttcgacttctgcgacgttagcggtgtccagatacagttccatgatgtttccttttacagttccagttcatgttgcagcaggctggcgatagcgtcggcagaggctgTGTGGAATGCATAGTAAGAACCTGTCTTGAAAAAATATCGCCGAATGTAACGACAATTCCTTAAGGATATCTGAAGGTATATTCAGAATTTGAAt                                                                                                                                                  <  1:129514‑M2/95‑1 (MQ=255)
                                                                                                 gggttagttgtcaccccggcaatggggaatatgcgtgccagacgttcgacttctgcgacgttagcggtgtccagatacagttccatgatgtttccttttacagttccagttcatgttgcagcaggctggcgatagcgtcggcagaggctgTGTGGAATGCATAGTAAGAACCTGTCTTGAAAAAATATCGCCGAATGTAACGACAATTCCTTAAGGATATCTGAAGGTATATTCAGAATTTGAATAAAATg                                                                                                                                            >  1:77714‑M2/151‑251 (MQ=255)
                                                                                                       gttgtcaccccggcaatggggaatatgcgtgccagacgttcggcttctgcgacgttagcggtgtccagatacagttccatgatgtttccttttacagttccagttcatgttgcagcaggctggcgatagcgtcggcagaggctgTGTGGAATGCATAGTAAGATCCTGTCTTGAAAAAATATCGCCGAATGTAACGACAAATCCGTAAGGATATCTGAAGGTATATTCAGAATTTGAAGATAATGcagaca                                                                                                                                      <  2:144082‑M2/107‑1 (MQ=255)
                                                                                                          gtctccccggcaatggggaatatgcgtgccagacgtgcgacttctgcgacgttagcggtgtccagatacagttccaggatgtttccttttacagttccagttcatgttgcagcaggctggcgatagcgtcggcagaggctgTGTGGAATGCATAGTAAGAACCTGTCGTGAAAAAATATCGCCGAATGTAACGACAATTCCTTAAGGATATCTGAAGGTATATTCAGAATTTGAATAAAATGCAGCCAGaa                                                                                                                                   <  2:489106‑M2/110‑1 (MQ=255)
                                                                                                          gtcaccccggcaatggggaatatgcgtgccagacgttcgacttctgcgacgttagcggtgtccagatacagttccatgatgtttccttttacagttccagttcatgttgcagcaggctggcgatagcgtcggcagaggctgTGTGGAATGCATAGTAAGAACCTGTCTTGAAAAAATATCGCCGAATGTAACGACAATTCCTTAAGGATATCTGAAGGTATATTCAGAATTTGAATAAAATGCAGACAGaa                                                                                                                                   <  1:447809‑M2/110‑1 (MQ=255)
                                                                                                              ccccggcaatggggaatatgcgtgccagacgttcgacttctgcgacgttagcggtgtccagatacagttccatgatgtttccttttacagttccagttcatgttgcagcaggctggcgatagcgtcggcagaggctgTGTGGAATGCATAGTAAGAACCTGTCTTGAAAAAATATCGCCGAATGTAACGACAATTCCTTAAGGATATCTGAAGGTATATTCAGAATTTGAATAAAATGCAGACAGAAatat                                                                                                                               <  1:198925‑M2/114‑1 (MQ=255)
                                                                                                                           gaatatgcgtgacagacgttcgacttctgcgaagttagcggtgtccagatacagttccatgatgtttccttttacagttccagttcatgttgcagcaggctggcgatagcgtcggcagaggctgTGTGGAATGCATAGTAAGAACCTGTCTTGAAAAAATATCGCCGAATGTAACGACAATTACTTAAGGATATCTGAAGGTATATTCAGAATTTGAATAAAATACAGACAGAAATATATTGAAAAAGAgg                                                                                                                  >  1:406269‑M2/125‑251 (MQ=255)
                                                                                                                                  cgtgccagacgttcgacttctgcgacgttagcggtgtccagatacagttccatgatgtttccttttacagttccagttcatgttgcagcaggctggcgatagcgtcggcagaggctgTGTGGAATGCATAGTAAGAACCTGTCTTGAAAAAATATCGCAGAATGTAACGACAATTCCTTAAGGATATCTGAAGGTATATTCAGAATTTGAATAAAATGCAGACAGAAATATATTGAAAACGAGGGTGTTAg                                                                                                           <  2:427626‑M2/134‑1 (MQ=255)
                                                                                                                                           cgttcgacttctgcgacgttagcggtgtccagatacagttccatgatgtttccttttacagttccagttcatgttgcagcaggctggcgatagcgtcggcagaggctgTGTGGAATGCATAGTAAGAACCTGTCTTGAAAAAATATCGCCGAATGTAACGACAATTCATTAAGGATATCTGAAGGTATATTTAGAATTTGAATAAAATGCAGACAGAAATATATTGAAAACGAGGGTGTTAGAACAGAAGt                                                                                                  >  1:298078‑M2/109‑251 (MQ=255)
                                                                                                                                                gacttctgcgacgttagcggtgtccagatacagttccatgatgtttccttttacagttccagttcatgttgcagcaggctggcgatagcgtcggcagaggctgTGTGGAATGCATAGTAAGAACCTGTCTTGAAAAAATATCGCCGAATGTAACGACAATTCCTTAAGGATATCTGAAGGTATATTCAGAATTTGAATAAAATGCAGACAGAAATATATTGAAAACGAGGGTGTTAGAACAGAAGTATTTc                                                                                             <  1:474606‑M2/148‑1 (MQ=255)
                                                                                                                                                       gcgacgttagcggtgtccagatacagttccatgatgtttccttttacagttccagttcatgttgcagcaggctggcgatagcgtcggcagaggctgTGTGGAATGCATAGTAAGAACCTGTCTTGAAAAAATATCGCCGAATGTAACGACAATTCCTTAAGGATATCTGAAGGTATATTCAGAATTTGAATAAAATGCAGACAGAAATATATTGAAAACGAGGGTGTTAGAACAGAAGTATTTCAGaaaa                                                                                       <  1:327887‑M2/154‑1 (MQ=255)
                                                                                                                                                            gttagcggtgtccagatacagttccatgatgtttccttttacagttccagttcatgttgcagcaggctggcgatagcgtcggcagaggctgTGTGGAATGCATa                                                                                                                                                                                                                                    <  1:134912‑M2/13‑1 (MQ=255)
                                                                                                                                                            gttagcggtgtccagatacagttccatgatgtttccttttacagttccagttcatgttgcagcaggctggcgatagcgtcggcagaggctgTGTGGAATGCATa                                                                                                                                                                                                                                    >  2:134912‑M2/92‑104 (MQ=255)
                                                                                                                                                              tagcggtgtccagataaagttccatgatgattccttttaccgttccagttcatgttgcagcaggctggcgatagcgtctgcagaggctgTGTGGAATGCATAGTAAGAAGCTGTCTTGAAAAAATATCGCCGAATGTAACGACAATTCCTTAAGGATATCTGAAGGTATCTTCAGAATTTGAATAAAATGCAGACAGAAATATATTGAAAACGCGGGTGTTAGAACAGAAGTATTTCAGAAAACCCtcgcc                                                                               >  2:121340‑M2/90‑250 (MQ=255)
                                                                                                                                                                 cggtgtccagatacagttccatgatgtttccttttacagttccagttcatgttgcagcaggctggcgatagcgtcggcagaggctgTGTGGAATGCATAGTAAGAACCTGTCTTGAAAAAATATCGCCGAATGTAACGACAATTCCTTAAGGATATCTGACGGTATATTCAGAATTTGAATAAAATGCCGACAGAAATATATTGAAAACGAGGGGGTTAGAACAGAAGTATTTCAGAAAACCATCGCGCa                                                                             >  2:237091‑M2/87‑250 (MQ=255)
                                                                                                                                                                   gtgtccagatacagttccatgatgtttccttttacagttccagttcatgttgcagcaggctggcgatagcgtcggcagaggctgTGTGGAATGCATAGTAAGAACCTGTCTTGAAAAAATATCGCCGAATGTAACGACAATTCCTTAAGGATATCTGAAGGTATATTCAGAATTTGAATAAAATGCAGACAGAAATATATTGAAAACGAGGGTGTTAGAACAGAAGTATTTCAGAAAACCCTCGCGCaaa                                                                           >  1:254049‑M2/85‑250 (MQ=255)
                                                                                                                                                                     gtccagatgcagttccatgatgtttccttttacagttccatttcatgttgcagcaggctggcgatagcgtcggcagaggctgTGTGGAATGCATAGTAAGAACCTGTCTTGAAAAAATATCGCCGAATGTAACGACAATTCCTTAAGTATATCTGAAGGTATATTCAGAATTTGAATAAAATGCAGACAGAAATATATTGAAAACGAGGGTGTTAGAACAGAAGTATTTCAGAAAACCCTCGCGCAAAAGc                                                                        <  2:318598‑M2/169‑1 (MQ=255)
                                                                                                                                                                               agttccatgatgtttccttttacagttccagttcatgttgcagcaggctggcgatagcgtcggaagaggctgTGTGGAATGCATAGTAAGAACCTGTCTTGAAAAAATATCGCCGAATGTAACGACAATTCCTTAAGGATATCTGAAGGTATATTCAGAATTTGAATAAAATGCAGACAGAAATATATTGAAAACGCGGGTGTTAGAACAGAAGTATTTCAGAAAACCCTCGCGAAAAAGCACGAGGGttt                                                              >  2:228201‑M2/73‑251 (MQ=255)
                                                                                                                                                                                gttccatgatgtttccttttacagttccagttcatgttgcagcaggctggcgatagcgtcggcagaggctgTGTGGAATGCATAGTAAGAACCTGTCTTGAAAAAATATCGCCGAATGTAACGACAATTCCTTAAGGATATCTGAAGGTATATTCAGAATTTGAATAAAATGCAGACAGAAATATATTGAAAACGAGGGTGTTAGAACAGAAGTATTTCAGAAAACCCTCGCGCAAAAGCACGAGGGTTTg                                                             >  1:500545‑M2/72‑251 (MQ=255)
                                                                                                                                                                                              ccttttacagttccagttcatgttgcagcaggctggcgatagcgtcggcagaggctgTGTGGAATGCATAGTAAGAACCTGTCTTGAAAAAATATCGCCGAATGTAACGACAATTCCTTAAGGATATCTGAAGGTATATTCAGAATTTGAATAAAATGCAGACAGAAATATATTGAAAACGAGGGTGTTAGAACAGAAGTATTTCAGAAAACCCTCGCGCAAAAGCACGAGTTTTTGCAGAAGAGGAAGAt                                               <  2:438720‑M2/194‑1 (MQ=255)
                                                                                                                                                                                                        ttccagttcatgttgcagcaggctggcgatagcgtcggcagaggctgTGTGGAATGCATAGTAAGAACCTGTCTTGAAAAAATATCGCCGAATGTAACGACa                                                                                                                                                                                          >  1:67613‑M2/48‑102 (MQ=255)
                                                                                                                                                                                                        ttccagttcatgttgcagcaggctggcgatagcgtcggcagaggctgTGTGGAATGCATAGTAAGAACCTGTCTTGAAAAAATATCGCCGAATGTAACGACa                                                                                                                                                                                          <  2:67613‑M2/55‑1 (MQ=255)
                                                                                                                                                                                                           cagttcatgttgcagcaggctggcgatagcgtcggcagaggctgTGTGGAATGCATAGTAAGAACCTGTCTTGAAAAAATATCGCCGAATGTAACGACAATTCCTTAAGGATATCTGAAGGTATATTCAGAATTTGAATAAAATGCAGACAGAAATATATTGAAAACGAGGGTGTTAGAACAGAAGTATTTCAGAAAACCCTCGCGCCAAAGCACGAGGGTTTGCAGAAGAGGAAGATTAGCCGGTATTAc                                  >  1:159058‑M2/45‑251 (MQ=255)
                                                                                                                                                                                                                 atgttgcagcaggctggcgatagcgtcggcagaggctgTGTGGAATGCATAGTAAGAACCTGTCTTGAAAAAATATCGCCGAATGTAACGACAATTCCTTAAGGATATCTGAAGGTATATTCAGAATTTGAATAAAATGCAGACAGAAATATATTGAAAACGAGGGTGTTAGAACAGAAGTATTTCAGAAAACCCTCGCGCAAAAGCACGAGGGTTTGCAGAAGAGGAAGATTAGCCGGTATTACGcatcc                            >  2:400299‑M2/39‑249 (MQ=255)
                                                                                                                                                                                                                            ggctggcgatagggtcggcagaggctgTGTGGAATGCATAGTAAGAACCTGTCTTGAAAAAATATCGCCGAATGTAACGACAATTCCTTAAGGATATCTGAAGGTATATTCAGAATTTGAATAAAATGCAGACAGAAATATATTGAAAACGAGGGTGTTAGAACAGAAGTATTTCAAAAACCCCTCGCGCAAAAGCACGAGGGTTTGCAGAAGAGGAAGATTAGCCGTTATTACGCATACCTGCCGCAATc                 <  2:77714‑M2/224‑1 (MQ=255)
                                                                                                                                                                                                                              ctggcgatagcgtcggcagaggctgTGTGGAATGCATAGTAAGAACCTGTCTTGAAAAAATATCGCCGAATGTAACGACAATTCCTTAAGGATATCTGAAGGTATATTCAGAATTTGAATAAAATGCAGACAGAAATATATTGAAAACGAGGGTGTTAGAACAGAAGTATTTCAGAAAACCCTCGCGCAAAAGCACGAGGGTTTGCAGAAGAGGAAGATTAGCCGGTATTACGCATACCTGCCGCAATccc               >  1:205695‑M2/26‑251 (MQ=255)
                                                                                                                                                                                                                                      agcgtcggcagaggctgTGTGGAATGCATAGTAAGAACCTGTCTTGAAAAAATATCGCCGAATGTAACGACAATTCCTTAAGGATATCTGAAGGTATATTCAGAATTTGAATAAAATGCAGACAGAAATATATTGAAAACGAGGGTGTTAGAACAGAAGTATTTCAGAAAACCCTCGCGCAAAAGCACGAGGGTTTGCAGAAGAGGAAGATTAGCCGGTATTACGCATACCTGCCGCAATCCCGGCAAt         >  1:11870‑M2/18‑249 (MQ=255)
                                                                                                                                                                                                                                         gtcggcagaggctgTGTGGAATGCATAGTAAGAACCTGTCTTGAAAAAATATCGCCGAATGTAACGACAATTCCTTAAGGATATCTGAAGGTATATTCAGAATTTGAATAAAATGCAGACAGAAATATATTGAAAACGAGGGTGTTAGAACAGAAGTATTTCAGAAAACCCTCGCGCAAAAGCACGAGGGTTTGCAGAAGAGGAAGATTAGCCGGTATTACGCATACCTGCCGCAATCCCGGCAATAGTg     >  1:256080‑M2/15‑250 (MQ=255)
                                                                                                                                                                                                                                           cggcagaggctgTGTGGAATGCATAGTAAGAACCTGTCTTGAAAAAATATCGCCGAATGTAACGACAATTCCTTAAGGATATCTGAAGGTATATTCAGAATTTGAATAAAATGCAGACAGAAATATATTGAAAACGAGGGTGTTAGAACAGAAGTATTTCAGAAAACCCTCGCGCAAAAGCACGAGGGTTTGCAGAAGAGGAAGATTAGCCGGTATTACGCATACCTGCCGCAATCCCGGCAATAGTGAcc  <  1:253841‑M2/239‑1 (MQ=255)
                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                              
ACAGCGCCGCCAGAGAGGCCGCCCATAGCACCACGCCGATTACTGCGGCAATAATTTTAATCCGCTTCGGAAAGAGGAATACCGGGATCAACCACAGCGAACTGAATAACAGCGAGTCGCGAATGCCGTTAGTGCCACTATAACCACTGATGTAAATAATGGCCTGTAGCAGAGTAGAGAAAAACCAAAAGTAGAGCAGTGCCCAACCCAGGGCTTTCCAGCTAAAAAGAGGTTTAGCCTGGACTTCTGTGGAATGCATAGTAAGAACCTGTCTTGAAAAAATATCGCCGAATGTAACGACAATTCCTTAAGGATATCTGAAGGTATATTCAGAATTTGAATAAAATGCAGACAGAAATATATTGAAAACGAGGGTGTTAGAACAGAAGTATTTCAGAAAACCCTCGCGCAAAAGCACGAGGGTTTGCAGAAGAGGAAGATTAGCCGGTATTACGCATACCTGCCGCAATCCCGGCAATAGTGACC  >  NC_000913/4150007‑4150492

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