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 26 [0] [0] 25 [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)25 (0.820)
+TG
25/486 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

CAGAGAGGCCGCCCATAGCACCACGCCGATTACTGCGGCAATAATTTTAATCCGCTTCGGAAAGAGGAATACCGGGATCAACCACAGCGAACTGAATAACAGCGAGTCGCGAATGCCGTTAGTGCCACTATAACCACTGATGTAAATAATGGCCTGTAGCAGAGTAGAGAAAAACCAAAAGTAGAGCAGTGCCCAACCCAGGGCTTTCCAGCTAAAAAGAGGTTTAGCCTGGACTTCTGTGGAATGCATAGTAAGAACCTGTCTTGAAAAAATATCGCCGAATGTAACGACAATTCCTTAAGGATATCTGAAGGTATATTCAGAATTTGAATAAAATGCAGACAGAAATATATTGAAAACGAGGGTGTTAGAACAGAAGTATTTCAGAAAACCCTCGCGCAAAAGCACGAGGGTTTGCAGAAGAGGAAGATTAGCCGGTATTACGCATACCTGCCGCAATCCCGGCAATAGTGACCATT  >  NC_000913/4150017‑4150495
                                                                                                                                                                                                                                             |                                                                                                                                                                                                                                                 
gcaaacagaatgccctcatcaccaatcgctttttgcagacgcggcagcacttcccatatggactccttgctggcagcgataatgctcgggttagttgtcaccccggcaatggggaatatgcgtgccagacgttcgacttctgcgacgttagcggtgtccagatacagttccatgatgtttccttttacagttccagttcatgttgcagcaggctggcgatagcgtcggcagaggctgTGTGGAATGCAt                                                                                                                                                                                                                                        >  1:170438‑M2/238‑249 (MQ=255)
           gccctcatcaccaatcgctttttgcagacgcggcagcacttcccatatggactccttgctggcagcgataatgctcgggttagttgtccccccggcaatggggaatatgcgtgccagacgttcgacttctgcgacgttagcggtgtccagatacagttccatgatgtttccttttacagttccagttcatgttgcagcaggctggcgatagcgtcggcagaggctgTGTGGAATGCATAGTAAGAACCTGt                                                                                                                                                                                                                           <  1:174302‑M2/25‑1 (MQ=255)
                                   cagacgcggcagcacttcccatatggactccttgctggcagcgataatgctcgggttagttggcaccccggcaatggggaatatgcgtgccagacgttcgacttctgcgacgttagcggtgtccagatacagttccatgatgtttccttttacagttccagttcatgttgcagcaggctggcgatagcgtcggcagaggctgTGTGGAATGCATAGTAAGAACCTGTCTTGAAAAAATATCGCCGAATGTa                                                                                                                                                                                                   <  1:170456‑M2/49‑1 (MQ=255)
                                         cggcagcacttcccatatggactcattgctggcagcgataatgctcgggttcgttgtcaccccggcaatggggaatatgcgtgccagccgttcgacttctgcgacgttagcggtgtccagatacagttccatgatgtttccttttacagttcccgttcatgttgcagcaggctggcgctagcgtcggcagaggctgTGTGGAATGCATAGTAAGAACCTGTCTTGCAACAATATCGCAGAATGTAACGAc                                                                                                                                                                                              >  2:301096‑M2/197‑250 (MQ=255)
                                            cagcacttcccatatggactccttgctggcagcgataatgctcgggttagttgtcaccccggcaatggggaatatgcgtgccagacgttcgacttctgcgacgttagcggtgtccagatacagttccatgatgtttccttttacagttccagttcatgttgcagcaggctggcgatagcgtcggcagaggctgTGTGGAATGCATAGTAAGAACCTGTCTTGAAAAAATATCGCCGAATGTAACGACAAtt                                                                                                                                                                                          <  1:278578‑M2/58‑1 (MQ=255)
                                                      catatggactccttgctggcagcgataatgctcgggttagttgtcaccccggcaatggggaatatgcgtgccagacgttcgacttctgcgacgttagcggtgtccagatacagttccatgatgtgtccttttagagtgccatttcatgtttcagcaggctggcgatagcgtcagcagagggtgTGTGGGATGCAGAAGAAGAACCTGTCTGGAAGAAATATCTCCTGATGTAACGAGAATTGCTTAAGGat                                                                                                                                                                                <  2:285213‑M2/68‑1 (MQ=255)
                                                                 cttgctggcagcgataatgctcgggttagttgtcaccccggcaatggggaatatgcgtgccagacgttcgacttctgcgacgttagcggtgtccagatacagttccatgatgtttccttttacagttccagttcatgttgcagcaggctggcgatagcgtcggcagaggctgTGTGGAATGCATAGTAAGAACCTGTCTTGAAAAAATATCGCCGAATGTAACGACAATTCCTTAAGGATATCTGAagga                                                                                                                                                                      >  1:23482‑M2/173‑249 (MQ=255)
                                                                        gcagcgataatgctcgggttagttgtcaccccggcaatggggaatatgcgtgccaggcgttcgacttctgcgacgttagcggtgtccagatacagttccatgatgtttccttttacagttccagtgcatgttgcagcaggctggcgatagcgtcggcagaggcggTGTGTAATGCATAGTAAGAACCTGTCTTGAAAAAATATCGCCGAATGTATCGACAATTCCTTAATGGTATCTGAAGGTATATTCAg                                                                                                                                                              <  2:300991‑M2/86‑1 (MQ=255)
                                                                              ataatgctcgggttagttgtcaccccggcaatggggaatatgcgtgccagacgttcgacttctgcgacgttagcggtgtccagatacagttccatgatgtttccttttacagttccagttcatgttgcagcaggctggcgatagcgtcggcagaggctgTGTGGAATGCATAGTAAGAACCTGTCTTGAAAAAATATCGCCGAATGTAACGACAATTCCTTAAGGATATCTGAAGGTATATTCAGAATTTg                                                                                                                                                        >  1:182266‑M2/160‑251 (MQ=255)
                                                                                   gctcgggttagttgtcaccccggcaatggggaatatgcgtgccagacgttcgacttctgcgacgttagcggtgtccagatacagttccatgatgtttccttttacagttccagttcatgttgcagcaggctggcgatagcgtcggcagaggctgTGTGGAATGCATAGTAAGAACCTGTCTTGAAAAAATATCGCCGAATGTAACGACAATTCCTTAAGGATATCTGAAGGTATATTCAGAATTTGAATaa                                                                                                                                                   <  1:295489‑M2/97‑1 (MQ=255)
                                                                                            agttgtcaccccggcaatggggaatatgcgtgccagacgttcgacttctgcgacgttagcggtgtccagatacagttccatgatgtttccttttacagttccagttcatgttgcagcaggctggcgatagcgtcggcagaggctgTGTGGAATGCATAGTAAGAACCTGTCTTGAAAAAATATCGCCGAATGTAACGACa                                                                                                                                                                                             >  1:228526‑M2/146‑200 (MQ=255)
                                                                                            agttgtcaccccggcaatggggaatatgcgtgccagacgttcgacttctgcgacgttagcggtgtccagatacagttccatgatgtttccttttacagttccagttcatgttgcagcaggctggcgatagcgtcggcagaggctgTGTGGAATGCATAGTAAGAACCTGTCTTGAAAAAATATCGCCGAATGTAACGACa                                                                                                                                                                                             <  2:228526‑M2/55‑1 (MQ=255)
                                                                                                                       gcgtgccagacgttcgacttctgcgacgttagcggtgtccagatacagttccatgatgtttccttttacagttccagttcatgttgcagcaggctggcgatagcgtcggcagaggctgTGTGGAATGCATAGTAAGAACCTGTCTTGAAAAAATATCGCCGAATGTAACGACAATTCCTTAAGGATATCTGAAGGTATATTCAGAATTTGAATAAAATGCAGACAGAAATATATTGAAAACGAGGGTGTTa                                                                                                               <  2:199921‑M2/133‑1 (MQ=255)
                                                                                                                           gccagacgttcgacttctgcgacgttagcggtgtccagatacagttccatgatgtttcctttgacagttccagtgcatgatgcagcaggctggcgatagcgtcggcagaggctgTGTGGAATGCATAGTAAGAACCTGTCTTGAAAAAATATCGCCGAATGGAACGGCAATTCCTTAAGGATATCTGAAGGTATATTCAGAATGTGAATGAAATGCAGACAGAAATATATTGAAAACGAGGGTGTTAGAAc                                                                                                           <  1:301096‑M2/137‑1 (MQ=255)
                                                                                                                              agacgttcgacttctgcgacgttagcggtgtccagctacagttccatgatgtttccttttacagttccagttcatgttgcagcaggctggcgatagcgtcggcagaggctgTGTGGAATGCATAGTAAGAACCTGTCTTGAAAAAATATCGCCGAATGTAACGACAATTCCTTAAGGATATCTGAAGGTATATTCAGAATTTGAATAAAATGCAGACAGAAATATATTGAAAACGAGGGTGTTAGAAc                                                                                                           >  2:121990‑M2/112‑248 (MQ=255)
                                                                                                                                         ttctgcgacgttagcggtgtccagatacagttccatgatgtttccttttacagttccagttcatgttgcagcaggctggcgatagcgtcggcagaggctgTGTGGAATGCATAGTAAGAACCTGTCTTGAAAAAATATCGCCGAATGTAACGACAATTCCTTAAGGATATCTGAAGGTATATTCAGAATTTGAATAAAATGCAGACAGAAATATATTGAAAACGAGGGTGTTAGAACAGAAGTATTTc                                                                                                >  1:117254‑M2/101‑248 (MQ=255)
                                                                                                                                              cgacgttagcggtgtccagatacagttccatgatgtttccttttacagttccagttcatgttgcagcaggctggcgatagcgtcggcagaggctgTGTGGAATGCATAGTAAGAACCTGTCTTGAAAAAATATCGCCGAATGTAACGACAATTCCTTAAGGATATCTGAAGGTATATTCAGAATTTGAATAAAATGCAGACAGAAATATATTGAAAACGAGGGTGTTAGAACAGAAGTATTTCAGAAAAcc                                                                                        >  1:272784‑M2/96‑251 (MQ=255)
                                                                                                                                                             ccagatacagttccatgatgtttccttttacagttccagttcatgttgcagcaggctggcgatagcgtcggcagaggctgTGTGGAATGCATAGTAAGAACCTGTCTTGAAAAAATATCGCCGAATGTAACGACAATTCCTTAAGGATATCTGAAGGTATATTCAGAATTTGAATAAAATGCAGACAGAAATATATTGAAAACGAGGGTGTTAGAACAGAAGTATTTCAGAAAACCCTCGCGCAAAAGCAc                                                                         <  1:218488‑M2/171‑1 (MQ=255)
                                                                                                                                                                      gttccatgatgtttccttttacagttccaggtcatgttgcagcaggctggcgatagcgtcggcagaggctgTGTGGAATGCATAGTAAGAACCTGTCTTGAAAAAATATCGCCGAATGTAACGACAATTCCTTAAGGATATCTGAAGGTATATTCAGAATTTGAATAAAATGCAGACAGAAATATATTGAAAACGAGGGTGTTAGAACAGAAGTATTTCAGAAAACCCTCGCTCAAAAGCACGAGGGTTTg                                                                >  2:116529‑M2/72‑251 (MQ=255)
                                                                                                                                                                                 tttccttttacagttccagttcatgttgcagcaggctggcgatagcgtcggcagaggctgTGTGGAATGCATAGTAAGAACCTGTCTTGAAAAAATATCGCCGAATGTAACGACAATTCCTTAAGGATATCTGAAGGTATATTCAGAATTTGAATAAAATGCAGACAGAAATATATTGAAAACGAGGGTGTTAGAACAGAAGTATTTCAGAAAACCCTCGCGCAAAAGCACGAGGGTTTGCAGAAGAGGa                                                      >  2:29502‑M2/61‑250 (MQ=255)
                                                                                                                                                                                          acagttccagttcatgttgcagcaggctggcgatagcgtcggcagaggctgTGTGGAATGCATAGTAAGAACCTGTCTTGAAAAAATATCGCCGAATGTAACGACAATTCCTTAAGGATATCTGAAGGTATATTCAGAATTTGAATAAAATGCAGACAGAAATATATTGAAAACGAGGGTGTTAGAACAGAAGTATTTCAGAAAACCCTCGCGCAAAAGCACGAGGGTTTGCAGAAGAGGAAGATTAGcc                                             <  2:54276‑M2/199‑1 (MQ=255)
                                                                                                                                                                                           cagttccagttcatgttgcagcaggctggcgatagcgtcggcagaggctgTGTGGAATGCATAGTAAGAACCTGTCTTGAAAAAATATCGCCGAATGTAACGACAATTCCTTCAGGATATCTGAAGGTATATTCAGAATTTGAATAAAATGCAGACAGAAATATATTGAAAACGAGGGTGTTAGAACAGAAGTATTTCAGAAAACCCTCGCGCAAAAGCACGAGGGTTTGCAGAAGAGGAAGATTAGCCgg                                           >  2:142945‑M2/51‑251 (MQ=255)
                                                                                                                                                                                                                      ggcgatagcgtcggcagaggctgTGTGGAATGCATAGTAAGAGCCTGTCTTGAAAAAATATCGCCGAATGTAACGACAATTCCTTAAGGATATCTGAAGGTATATTCAGAATTTGAATAAAATGCAGACAGAAATATATTGAAACCGAGGGTGTTAGAACAGAAGTATTTCAGAAACCCCTCGCCCAAAAGCACGAGGGTTTGCAGAAGAGGAAGATTAGCCGGTATTACGCATACCTGCCGCAATCCCgg                <  2:75370‑M2/228‑1 (MQ=255)
                                                                                                                                                                                                                      ggcgatagcgtcggcagaggctgTGTGGAATGCATAGAAAGAACCTGTCTTGAAAAAATATCGCCGAATGTAACGACAATCCCTTAAGGATATCGGAAGGTATATTCAGAATTTGAATAAAATGCAGACAGAAATATATTGAAAACGAGGGTGTTAGAACAGAAGTATTTCAGAAACCCCTCCCGCAAAAGCACGAGGGTTTGCAGAAGAGGAAGATTAGCCGGTATTACCCATACCTGCCGCAATccc                  <  2:263622‑M2/226‑1 (MQ=255)
                                                                                                                                                                                                                                    cagaggctgTGTGGAATTCATAGTAAGAACCTGTCTTGAAAAAATTTCGCCGAATGTAACGACAATTCTTTAAGGATATCTGAAGGTATATTCAGAATTTGAATAAAATGCAGACAGAAATATATTGAAAACGAGGGTGTTAGAACTGAAGTTTTTCAGAAAACCGTCGCGCTAAAGCACGGGGGGTTGCAGAAGAGGAAGATTAGCAGTGAGTGCGCATACCTGCCGCAATCCCGGCAATAGTGACCAtt  <  1:292356‑M2/242‑1 (MQ=255)
                                                                                                                                                                                                                                             |                                                                                                                                                                                                                                                 
CAGAGAGGCCGCCCATAGCACCACGCCGATTACTGCGGCAATAATTTTAATCCGCTTCGGAAAGAGGAATACCGGGATCAACCACAGCGAACTGAATAACAGCGAGTCGCGAATGCCGTTAGTGCCACTATAACCACTGATGTAAATAATGGCCTGTAGCAGAGTAGAGAAAAACCAAAAGTAGAGCAGTGCCCAACCCAGGGCTTTCCAGCTAAAAAGAGGTTTAGCCTGGACTTCTGTGGAATGCATAGTAAGAACCTGTCTTGAAAAAATATCGCCGAATGTAACGACAATTCCTTAAGGATATCTGAAGGTATATTCAGAATTTGAATAAAATGCAGACAGAAATATATTGAAAACGAGGGTGTTAGAACAGAAGTATTTCAGAAAACCCTCGCGCAAAAGCACGAGGGTTTGCAGAAGAGGAAGATTAGCCGGTATTACGCATACCTGCCGCAATCCCGGCAATAGTGACCATT  >  NC_000913/4150017‑4150495

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