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

AGCACAGCGCCGCCAGAGAGGCCGCCCATAGCACCACGCCGATTACTGCGGCAATAATTTTAATCCGCTTCGGAAAGAGGAATACCGGGATCAACCACAGCGAACTGAATAACAGCGAGTCGCGAATGCCGTTAGTGCCACTATAACCACTGATGTAAATAATGGCCTGTAGCAGAGTAGAGAAAAACCAAAAGTAGAGCAGTGCCCAACCCAGGGCTTTCCAGCTAAAAAGAGGTTTAGCCTGGACTTCTGTGGAATGCATAGTAAGAACCTGTCTTGAAAAAATATCGCCGAATGTAACGACAATTCCTTAAGGATATCTGAAGGTATATTCAGAATTTGAATAAAATGCAGACAGAAATATATTGAAAACGAGGGTGTTAGAACAGAAGTATTTCAGAAAACCCTCGCGCAAAAGCACGAGGGTTTGCAGAAGAGGAAGATTAGCCGGTATTACGCATACCTGCCGCAATCCCGGCAATAGTGACCATTAACGCT  >  NC_000913/4150004‑4150501
                                                                                                                                                                                                                                                          |                                                                                                                                                                                                                                                       
gctcatggtctgagcaaacagaatgccctcatcaccaatcgctttttgcagacgcggcagcacttcccatatggactccttgctggcagcgataatgctcgggttagttgtcaccccggcaatggggaatatgcgtgccagacgttcgacttctgcgacgttagcggtgtccagatacagttccatgatgtttccttttacagttccagttcatgttgcagcaggctggcgatagcgtcggcagaggctgt                                                                                                                                                                                                                                                         <  2:314585‑M2/1‑1 (MQ=255)
gctcatggtctgagcaaacagaatgccctcatcaccaatcgctttttgcagacgcggcagcacttcccatatggactccttgctggcagcgataatgcgcgggttagttgtcaccccggcaatggggaatatgcgtgccagacgttcgacttctgcgacgttagcggtgtccagatacagttccatgatgtttccttttacagttccagtgcatgttgcagcaggctggcgatagcgtcggcagaggctgt                                                                                                                                                                                                                                                         <  1:279812‑M2/1‑1 (MQ=255)
   catggtctgagcaaacagaatgccctcatcaccaatcgctttttgcagacgcggcagcacttcccatatggactccttgctggcagcgataatgctcgggttagttgtcaccccggcaatggggaatatgcgtgccagacgttcgacttctgcgacgttagcggtgtccagatacagttccatgatgtttccttttacagttccagttcatgttgcagcaggctggcgatagcggcggcagaggctgtgtg                                                                                                                                                                                                                                                      <  1:182328‑M2/4‑1 (MQ=255)
           gagcaaacagaatgccctcatcaccaatcgctttttgcagacgcggcagcacttcccatatggactccttgctggcagcgataatgctcgggttagttgtcaccccagcaatggggaatatgcgtgccagacgttcgacttctgcgacgttatcgttgtccatatacagttccatgatgtttccttgtacagttccagttcatgttgcagcaggctggcgatagcgtcggcagaggctgTGTGGAATg                                                                                                                                                                                                                                                 <  1:297602‑M2/9‑1 (MQ=255)
              caaacagaatgccctcatcaccaatcgctttttgcagacggggcagcccttcccatatggactccttgctggcagcgataatgcgcgggttagttgtcaccccggcaatggggaatatgcgtgccagacgttcgacttctgcgacgttagcggtgtccagatacagttccatgatgtttccttttacagtgccagttcatgttgcagcaggctggcgatagcgtcggcagaggctgGGTGGAATGCATAGt                                                                                                                                                                                                                                           <  2:212770‑M2/15‑1 (MQ=255)
                    gaatggactcctatccactcgctttttgcagacacggtagcacttaccctatggactccttgctgccagcgataatgctcgggttagttgtcaccccggcaatggggaatatgcgtgccagccgttcgacttctgcgacgttagcggtgtccagatacagctccatgatgtttccttttacagttccagttcatgttgcagcaggctggcgatcgcgtcggcagaggctgTGTGGTATGCATAGTAAGAAc                                                                                                                                                                                                                                     >  2:463470‑M2/231‑251 (MQ=255)
                             catcaccaatcgctttttgcagacgcggcagcacttcccatatggactccttgctggcagcgataatgctcgggttagttgtcaccccggcaatggggaatatgcgtgccagacgttcgacttctgcgacgttagcggtgtccagatacagttccatgatgtttccttttacagttccagttcatgttgcagcaggctggcgatagcgtcggcagaggctgTGTGGAATGCATAGTAAGAACCTGTCTTg                                                                                                                                                                                                                             >  1:70314‑M2/222‑250 (MQ=255)
                                             ttgcagacgcggcagcacttcccatatggactccttgctggcagcgataatgctcgggttagttgtcaccccggcaatggggaatatgcgtgccagacgttcgacttctgcgacgttagcggtgtccagatacagttccatgatgtttccttttacagttccagttcatgttgcagcaggctggcgatagcgtcggcagaggctgTGTGGAATGCATAGTAAGAACCTGTCTTGAAAAAATATCg                                                                                                                                                                                                                  >  1:312919‑M2/206‑245 (MQ=255)
                                             ttgcagacgcggcagcacttcccatatggactccttgctggcagcgataatgctcgggttagttgtcaccccggcaatggggaatatgcgtgccagacgttcgacttctgcgacgttagcggtgtccagatacagttccatgatgtttccttttacagttccagttcatgttgcagcaggctggcgatagcgtcggcagaggctgTGTGGAATGCATAGTAAGAACCTGTCTTGAAAAAATATCg                                                                                                                                                                                                                  <  2:312919‑M2/40‑1 (MQ=255)
                                                   acgcggcagcacttcccatatggactccttgctggcagcgataatgctcgggttagttgtcaccccggcactggggaatatgcgtgccagacgttcgacttctgcgacgttagcggtgtccagatacagttccatgatgtttccttttacagttccagttcatgttgcagcaggctggcgatagcgtcggcagaggctgTGTGGAATGCATAGTAAGACCCTGTCTTGAAAAAATATCGCCGAATGTAACg                                                                                                                                                                                                      >  2:430582‑M2/200‑251 (MQ=255)
                                                       ggcagcacttcccatatggactccttgctggcagcgataatgctcgggttagttgtcaccccggcaatggggaatatgcgtgccagacgttcgacttctgcgacgttagcggtgtccagatacagttccatgatgtttccttttacagttccagttcatgttgcagcaggctggcgatagcgtcggcagaggctgTGTGGAATGCATAGTAAGAACCTGTCTTTAAAAAATATCGCCGAATGTAACGACaa                                                                                                                                                                                                  <  2:207264‑M2/56‑1 (MQ=255)
                                                                                    ggcagcgataatgctcgggttagttgtcaccccggcaatggggaatatgcgtgccagacgttcgacttctgcgacgttagcggtgtccagatacagttccatgatgtttccttttacagttccagttcatgttgcagcaggctggcgatagcgtcggcagaggctgTGTGGAATGCATAGTAAGAACCTGTCTTGAAAAAATATCGCCGAATGTACCGACAATTCCTTAAGGATATCTGAAGGTATATTCa                                                                                                                                                                     <  1:445935‑M2/85‑1 (MQ=255)
                                                                                    ggcagcgataatgctcgggttagttgtcaccccggcaatggggaatatgcgtgccagacgttcgacttctgcgacgttagcggtgtccagatacagttccatgatgtttccttttacagttccagttcatgttgcagcaggctggcgatagcgtcggcagaggctgTGTGGAATGCATAGTAAGAACCTGTCTTGAAAAAATATCGCCGAATGTAACGACAATTCCTTAAGGATATCTGAAGGTATATTc                                                                                                                                                                      <  1:5963‑M2/84‑1 (MQ=255)
                                                                                         cgataatgctcgggttagttgtcaccccggcaatggggaatatgcgtgccagacgttcgacttctgcgacgttagcggtggccagatacagttccatgatgtttccttttacagttccagttcatgttgcagcaggctggggatagcgtcggcagaggctgTGTGGAAGGCAGAGTAAGAACCTGTCTTGAAAAAATAGCGCCGAATGCAACGAGAATTCCTTAAGGATATCTGAAGGTATATTCAGAAtt                                                                                                                                                                <  1:463470‑M2/90‑1 (MQ=255)
                                                                                             aatgctcgggttagttgtcaccccggcaatggggaatatgcgtgccagacgttcgacttctgcgacgttagcggtgtccagatacagttccatgatgtttccttttacagttccagttcatgttgcagcaggctggcgatagcgtcggcagaggctgTGTGGAATGCATAGTAAGAACCTGTCTTGAAAAAATATCGCCGAATGTAACGACAATTCCTTAAGGATATCTGAAGGTATATTCAGAATTTGa                                                                                                                                                             >  2:307291‑M2/158‑250 (MQ=255)
                                                                                                gctcgggttagttgtcaccccggcaatggggaatatgcgtgccagacgttcgacttctgcgacgttagcggtgtccagatacagttccatgatgtttccttttacagtgccagttcatgttgcagcattctggcgatagcgtcggcagaggctgTGGGGAATGCATAGTAAGAACCTGTCTTGAAAAAATATCGCCGAATGTAACTACAATTCCTTAAGGATATCTGAAGGTATATTCAGAATTTGAATa                                                                                                                                                          <  1:444166‑M2/96‑1 (MQ=255)
                                                                                                    gggttagttgtcaccccggcaatggggaatatgcgtgccagacgttcgacttctgcgacgttagcggtgtccagatacagttccatgatgtttccttttacagttccagttcatgttgcagcaggctggcgatagcgtcggcagaggctgTGTGGAATGCATAGTAAGAACCTGTCTTGAAAAAATATCGCCGAATGTAACGACAATTCCTTAAGGATATCTGAAGGTATATTCAGAATTTGAATAAAATg                                                                                                                                                     <  2:70314‑M2/101‑1 (MQ=255)
                                                                                                        tagttgtcaccccggcaatggggaatatgcgtgccagacgttcgacttctgcgacgttagcggtgtccagatacagttccatgatgtttccttttacagttccagttcatgttgcagcaggctggcgatagcgtcggcagaggctgTGTGGAATGCATAGTAAGAACCTGTCTTGAAAAAATATCGCCGAATGTAACGACAATTCCTTAAGGATATCTGAAGGTATATTCAGAATTTGAATAAAATGcag                                                                                                                                                  <  1:375843‑M2/104‑1 (MQ=255)
                                                                                                                 ccccggcaatggggaatatgcgtgccagacgttcgacttctgcgacgttagcggtgtccagatacagttccatgatgtttccttttacagttccagttcatgttgcagcaggctggcgatagcgtcggcagaggctgTGTGGAATGCATAGTAAGAACCTGTCTTGAAAAAATATCGCCGAATGTAACGACAATTCCTTAAGGATATCTGAAGGTATATTCAGAATTTGAATAAAATGCAGACAGAAatat                                                                                                                                        >  1:138231‑M2/138‑251 (MQ=255)
                                                                                                                      gcaatggggaatatgcgtgccagacgttcgacttctgcgacgttagcggtgtccagatacagttccatgatgtttccttttacagttccagttcatgttgcagcaggctggcgatagcgtcggcagaggctgTGTGGAATGCATAGTAAGAAc                                                                                                                                                                                                                                     <  1:67647‑M2/21‑1 (MQ=255)
                                                                                                                      gcaatggggaatatgcgtgccagacgttcgacttctgcgacgttagcggtgtccagatacagttccatgatgtttccttttacagttccagttcatgttgcagcaggctggcgatagcgtcggcagaggctgTGTGGAATGCATAGTAAGAAc                                                                                                                                                                                                                                     >  2:67647‑M2/133‑153 (MQ=255)
                                                                                                                       caatggggaatatgcgtgccagacgttcgacttctgcgacgttagcggtgtccagatacagttccatgatgtttccttttacagttccagttcatgttgcagcaggctggcgatagcgtcggcagaggctgTGTGGAATGCATAGTAAGAACCTGTCTTGAAAAAATATCGCCGAATGTAACGACAATTCCTTAAGGATATCTGAAGGTATATTCAGAATTTGAATAAAATGCAGACAGAAATATATTGaa                                                                                                                                  <  1:307291‑M2/120‑1 (MQ=255)
                                                                                                                                       tgccagacgttcgacttctgcgacgttagcggtgtccatatacagttccatgatgtttccttttatagttccagttcatgttgcagcaggctggcgatagcggctgcagaggctgTGTGGAATGCATAGTAAGAACCTGTCGTGAAAAAATATCGCCGAATGTAACGACAATTCCTTAAGGATATCTGAAGGTATATTCAGAATTTGAATAAAATGCAGACAGAAATATATTGAAAACGAGGGTGTTAGa                                                                                                                   <  2:258862‑M2/135‑1 (MQ=255)
                                                                                                                                       tgccagacgttcgacttctgcgacgttagcggtgtccagatacagttccatgatgtttccttttacagttccagttcatgttgcagcaggctggcgatagcgtcggcagaggctgTGTGGAATGCATAGTAAGAACCTGTCTTGAAAAAATATCGCCGAATGTAACGACAATTCCTTAAGGATATCTGAAGGTATATTCAGAATTTGAATAAAATGCAGACAGAAATATATTGAAAACGAGGGTGTTAGa                                                                                                                   >  2:112437‑M2/116‑250 (MQ=255)
                                                                                                                                                   gacttctgcgacgttagcggtgtccagatacagttccatgatgtttccttttacagttccagttcatgttgcagcaggctggcgatagcgtcggcagaggctgTGTGGAATGCATAGTAAGAACCTGTCTTGAAAAAATATCGCCGAATGTAACGACAATTCCTTAAGGATATCTGAAGGTATATTCAGAATTTGAATAAAATGCAGACAGAAATATATTGAAAACGAGGGTGTTAGAACAGAAGTAtt                                                                                                        >  1:4512‑M2/104‑249 (MQ=255)
                                                                                                                                                   gacttctgcgacgttagcggtgtccagatacagttccatgatgtttccttttacagttccagttcatgttgcagcaggctggcgatagcgtcggcagaggctgTGTGGAATGCATAGTAAGAACCTGTCTTGAAAAAATATCGCCGAATGTAACGACAATTCCTTAAGGATATCTGAAGGTATATTCAGAATTTGAATAAAATGCAGACAGAAATATATTGAAAACGAGGGTGTTAGAACAGAAGTAtt                                                                                                        <  2:4512‑M2/146‑1 (MQ=255)
                                                                                                                                                        ctgcgacgctagcggagtcaagatacagttccatgatgtttccttttacagttccagttcatgttgcagcaggctggcgatagcgtcggcaaaggctgTGTGGAATGAAAAGTAAGAACCTGTCTTGAAAAAATATCGCCGAATGTAACGACAATTCCTTAAGGATATCTGAAGGTATATTCAGAATTTGATTAAAATGCTGACAGAAATATATTGAAAACGAGGGTGTTAGAACAGAAGTATTTCAGaa                                                                                                  >  2:41933‑M2/99‑250 (MQ=255)
                                                                                                                                                                     ggtgtccagatacagttccatgatgtttccttttacagttccagttcatgttgcagcaggctggcgatagcgtcggcagaggctgTGTGGAATGCATAGTAAGAACCTGTCTTGAAAAAATATCGCCGAATGTAACGACAATTCCTTAAGGATATCTGAAGGTATATTCAGAATTTGAATAAAATGCAGACAGAAATATATTGAAAACGAGGGTGTTAGAACAGAAGTATTTCAGAAAACCCTCGCGCaa                                                                                     >  2:376125‑M2/86‑250 (MQ=255)
                                                                                                                                                                       tgtccagatacagttccctgatgtttccttttacagttccagttcatgttgcagcaggctggcgatagcgtcggcagaggctgTGTGGAATGCATAGTAAGAACCTGTCTTGAAAAAATATCGCCGAATGTAACGACAATTCCTTAAGGATATCTGAAGGTATATTCAGAATTTGAATAAAATGCAGACAGACATATATTGAAAACGAGGGTGTTAGAACAGAAGTATTTCAGAAAACCCTCGCGCaaa                                                                                    >  2:280751‑M2/84‑249 (MQ=255)
                                                                                                                                                                                    ttccatgatgtttccttttacagttccagttcatgttgcagcaggctggcgatagcgtcggcagaggctgTGTGGAATGCATAGTAAGAACCTGTCTTGAAAAAATATCGCCGAATGTAACGACAATTc                                                                                                                                                                                               <  2:120598‑M2/59‑1 (MQ=255)
                                                                                                                                                                                    ttccatgatgtttccttttacagttccagttcatgttgcagcaggctggcgatagcgtcggcagaggctgTGTGGAATGCATAGTAAGAACCTGTCTTGAAAAAATATCGCCGAATGTAACGACAATTc                                                                                                                                                                                               >  1:120598‑M2/71‑129 (MQ=255)
                                                                                                                                                                                                                               ggctggcgatagcgtcggaagaggctgTGTGGAATGCATAGTAAGAACCTGTCTTGAAAAAATATCGCCGAATGTAACGACAATTCCTTAAGGATATCTGAAGGTATATTCAGAATTTGAATAAAATGCAGACAGAAATATATTGAAAACGAGGGTGTTAGAACAGAAGTATTTCAGAAAACCCTCGCGCAAAAGCACGAGGGTTTGCAGAAGAGGAAGATTAGCCGGTATTACGCATACCTGCCGCAATc                          >  2:322864‑M2/28‑251 (MQ=255)
                                                                                                                                                                                                                                      gatagcgtcggcagaggctgTGTGGAATGCATAGTAAGAACCTGACTTGAAAAAATATCGCCGAATGGAACAACAATTCCTTAAGGATATCGGAGGGTATATTCAGAATTTGAAAAAAATGCAGACAGAAATATATTGAAAACGAGGTTGTTAGAACAGAAGTATTTCAGAAAACCCTCGCGCAAAAGCACGAGGGTTTTCAGAAGAGGAAAATTAGCCGGTATTACGCATACCTGCCGCAATCCCGGCaa                   <  1:211161‑M2/231‑1 (MQ=255)
                                                                                                                                                                                                                                                       ctgTGTGGAATGCATAGTAAGAACCTGTCTTGAAAAAATATCGCCGAATGTAACGACAATTCCTTAAGGATATCTGAAGGTATATTCAGAATTTGAATAAAATGCAGACAGAAATATATTGAAAACGAGGGTGTTAGAACAGAAGTATTTCAGAAAACCCTCGCGCAAAAGCACGAGGGTTTGCAGAAGAGGAAGATTAGCCGGTATTACGCATACCTGCCGCAATCCCGGCAATAGTGACCATTAACGCt  >  1:399108‑M2/4‑251 (MQ=255)
                                                                                                                                                                                                                                                          |                                                                                                                                                                                                                                                       
AGCACAGCGCCGCCAGAGAGGCCGCCCATAGCACCACGCCGATTACTGCGGCAATAATTTTAATCCGCTTCGGAAAGAGGAATACCGGGATCAACCACAGCGAACTGAATAACAGCGAGTCGCGAATGCCGTTAGTGCCACTATAACCACTGATGTAAATAATGGCCTGTAGCAGAGTAGAGAAAAACCAAAAGTAGAGCAGTGCCCAACCCAGGGCTTTCCAGCTAAAAAGAGGTTTAGCCTGGACTTCTGTGGAATGCATAGTAAGAACCTGTCTTGAAAAAATATCGCCGAATGTAACGACAATTCCTTAAGGATATCTGAAGGTATATTCAGAATTTGAATAAAATGCAGACAGAAATATATTGAAAACGAGGGTGTTAGAACAGAAGTATTTCAGAAAACCCTCGCGCAAAAGCACGAGGGTTTGCAGAAGAGGAAGATTAGCCGGTATTACGCATACCTGCCGCAATCCCGGCAATAGTGACCATTAACGCT  >  NC_000913/4150004‑4150501

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: 
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.