Predicted mutation
evidence seq id position mutation annotation gene description
MC JC NC_000913 3,909,311 Δ774 bp [pstA][pstC] [pstA], [pstC]

Missing coverage evidence...
   seq id start end size ←reads reads→ gene description
* * ÷ NC_000913 3909311 3910084 774 40 [0] [0] 39 [pstA]–[pstC] [pstA], [pstC]

New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913 = 39093100 (0.000)33 (0.700) 30/268 0.6 100% coding (129/891 nt) pstA phosphate ABC transporter permease
?NC_000913 3910085 = 0 (0.000)coding (313/960 nt) pstC phosphate ABC transporter permease

AAATTAACAGCCCGCTACCCGCCAGAGCGTTCGCCAGACCACCACCTTCCGTATTGGGCGGCGGCGTCATTTCAGTGAACAGCGCCAGCGACATACCGTCGATACCGCGAGTGATGGTGGACATTAAAATCCAGATCAGCCAGAACAGGCCGAAGGCCATCGTCGCCATCGAGAGCGTCAGCGCAATACGGTTTTTGAGGCGGCGACGCGCCTGCATTTTGCGGCGAGATTCAGCCAGCGCCGCAGTGGTTTGCATTTCAACCATAGCCATTAGCGTGCCCCCT  >  NC_000913/3909168‑3909451
                                                                                                                                              |                                                                                                                                             
aaaGTAGCAGCCCGCTGCCCGCCAGAGCGGTCGCCAGACCACCACCTTCCGTATTGGGCGGCGGCGTCATTTCAGTGAACAGCGCCAGCGACATACCGTCGATACCGCGAGTGATGGTGGACATTAAAATCCAGATCAGCCAG                                                                                                                                               <  1:286023‑M1/143‑1 (MQ=255)
 aaTTAACAGCCCGCTACCCGCCAGAGCGTTCGCCAGACCACCACCTTCCGTATTGGGCGGCGGCGTCATTTCAGTGAACAGCGCCAGCGACATACCGTCGATACCGCGAGTGATGGTGGACATTAAAATCCAGATCAGCCAGg                                                                                                                                              >  2:426694‑M1/1‑142 (MQ=255)
  aTTAACAGCCCGCTACCCGCCAGAGCGTTCGCCAGACCACCACCTTCCGTATTGGGCGGCGGCGTCATTTCAGTGAACAGCGCCAGCGACATACCGTCGATACCGCGAGTGATGGTGGACATTAAAATCCAGATCAGCCAGgc                                                                                                                                             <  1:371994‑M1/143‑3 (MQ=255)
         gCCCGCTACCCGCCAGTGCTTTCGCCAGTCCCCCACCGTCCGTATTGGGCGGCGGCGTCATTTCAGTGAACAGCGCCAGCGACATACCGTCTATACCGCGAGTGATGGTGTACATTAAAATCCAGATCAGCCAGgcgcaagct                                                                                                                                      <  1:155434‑M1/143‑10 (MQ=255)
           ctgCTACTCTCCAGTTCGTTCGCCGGTCCCCCACTTGCCGTATTGTGCGGGTGCGTTATTTTTGTTAATGGCGCCACCACCATACCGTCGATACCGCGAGTGATGGTGGTTATTAAAATCCAGATCAGCCAggcgcaagctca                                                                                                                                    <  1:829093/141‑12 (MQ=255)
            cGCTACCCGCCAGAGCGTTCGCCAGACCACCACCTTCCGTATTGGGCGGCGGCGTCATTTCAGTGAACAGCGCCAGCGACATACCGTCGATACGGCGAGCGATGGTGGACAGTAAAATCCAGATCAGCCAGgcgcaagcgcag                                                                                                                                   >  1:218705‑M1/1‑131 (MQ=255)
                  ccGCCAGAGCGTTCGCCAGACCACCACCTTCCGTCTTGGGCGGCGGCGCCATTTCAGTGAACAGCGCCAGCGACATACCGTCGATACCGCGCGTGATGGTGGACATTAAAATCCAGATCAGCCAGgcgcacgcgcagtcacga                                                                                                                             >  1:192678‑M1/1‑125 (MQ=255)
                   cGCCAGAGCGTTCGCCAGACCACCACCTTCCGTATTGGGCGGCGGCGTCATTTCAGTGAACAGCGCCAGCGACATACCGTCGATACCGCGAGTGATGGTGGACATTAAAATCAAGATCAGCCAGacgcaagctcagtcaggaa                                                                                                                            >  1:544857‑M1/1‑124 (MQ=255)
                             ttCGCCAGACCACCACCTTCCGTATTGGGCGGCGGCGTCATTTCAGTGAGCAGCGCCAGCGACATACCGTCGATACCGCGAGTGATGGTGGACCTTAAAATCCAGATCAGCCAGgctcaagctcagtcaggaacagggcgata                                                                                                                  >  2:460909‑M1/1‑114 (MQ=255)
                                  cAGACCACCACCTTCCGTATTGGGCGGCGGCGTCATTTCAGTGCACAGCGCCAGCGACATACGGTCGATACCGCGAGTGATGGTGGACATTAAGATCCAGATCAGCCAGgcgagagatcagtccggagcagggacatgccgaa                                                                                                             >  1:246196‑M1/1‑109 (MQ=255)
                                         ccaccTTCCGTATTGGGCGGCGGCGTCATTTCAGTGAACAGCGCCAGCGACATACCGTCGATACCGCGAGTGATGGTGGACATTAAAATCCAGATCAGCCAGgcgcaagctcagtcaggaacagggcgataccgaaactcacc                                                                                                      <  2:351291‑M1/143‑42 (MQ=255)
                                                             cggctTCATTTCACTGAAACGGGCCAGCGACATACCGTCTATACCGCGAGTGATGGTGGACATTAAAAGCCAGATCAGCCAGgcgcaagctcagtcaggaacagggcgataccgaaactcaccgggacggcgatcagcagcgc                                                                                  <  1:59259‑M1/143‑62 (MQ=255)
                                                                 gTCATTTCAGTGAACAGCGCCAGCGACATACCGTCGATACCGCGAGTGATGGTGGACATTAAAATCCAGATCAGCCAGgcgcaagctcagtcaggaacagggcgatcacgaaactcaccgggacggcgatcagcagcgcgata                                                                              >  1:29419‑M1/1‑78 (MQ=255)
                                                                          gTGAACAGCGCCAGCGACATACCGTCGATACCGCGAGTGATGGTGGACATTAAAATCCAGATCAGCCAGgcgcaagctcagtcaggaacagggcgataccgaaactcaccgggacggggaacagcagcgccaccaacgcaggc                                                                     >  1:4211‑M1/1‑69 (MQ=255)
                                                                           tGAACAGCGCCAGCGACATACCGTCGATACCGCGAGTGATGGTGGACATTAAAATCCAGATCAGCCAGgcgcaagctcagtcaggaacagggcgataccgaaactcaccgggacggcgatcagcagcgcgatacacgaaggca                                                                    >  2:485751‑M1/1‑68 (MQ=255)
                                                                             aaCAGCGCCAGCGACATACCGTCGATACCGCGAGTGATGCTGGACATTAAAATCCAGATCAGCCAGgcgcaagctcagtcaggaacagggcgataccgaaactcaccgggacggcgatcagcagcgcgctaaacgaagtcacc                                                                  >  1:703554‑M1/1‑66 (MQ=255)
                                                                                aGCGCCAGCGACATACCGTCGATACCGCGAGTGATGGTGGACATTCAAATCCAGATCAGCCAGgcgcaagctcagtcaggaacagggcgataccgaaactcaccgcgacggcgatcagcagcgcgatacacgaagtcaccaac                                                               >  1:67313‑M1/1‑63 (MQ=255)
                                                                                 gcgcCAGCGACATACCGTCGATACCGCGAGTGATGGTGGACATTAAAATCCAGATCAGCCAGgcgcaagctcagtcaggaacagggcgataccgaaactcaccgggacggcgatcagcagcgcgataaacgaagtcaccaacg                                                              >  1:451694‑M1/1‑62 (MQ=255)
                                                                                         gACACACCGTCGATACCGCGAGTGAGGGTGGACAGTAAAATCCAGATCAGCCAGgcgcaagcgcagtcaggaacagggcgataccgaaactcaccgggacggcgatcagcagcgcgataaacgaagtcaccaacgtaccgtag                                                      <  2:53806‑M1/143‑90 (MQ=255)
                                                                                            aTACCGTCGATCCCGCGAGTCATCGTGCACATCAAAATCCAGTTTTGCCAGgcgcaagctcagtcaggaacagggcgataccgaaactcaccgggacggcgatcagcagcgcgataaacgaagtcaccaacgtaccgtagatc                                                   <  2:30115‑M1/143‑93 (MQ=255)
                                                                                               ccggcggggccgcgggggatGGCGGGCATTAAAAGCCAGATCAGCCCGgcgcaagctcagtcaggaacagggcgataccgaaactcaccgggacggcgatcagcagcgcgataaacgaagtcaccaacgtaccgtagatcggc                                                <  2:50579‑M1/134‑96 (MQ=255)
                                                                                                 gTCGATACCGCGAGTGATGGTGGACATTAAAATCCAGATCAGCCAGgcgcaagctcagtcaggaacagggcgataccgaaactcaccgggacggcgatcagcagcgcgataaacgaagtcaccaacgtaccgtagatcggcac                                              >  2:2916‑M1/1‑46 (MQ=255)
                                                                                                     aTACCGCGAGTGATGGTGGACATTAAAATCCAGATCAGCCAGgcgcaagctcagtcaggaacagggcgataccgaaactcaccgggacggcgaccagcagcgcgataaacgcagaccccaacgtaccgcacaacggcaccagc                                          >  1:246177‑M1/1‑42 (MQ=255)
                                                                                                         cgcgAGTGATGGTGGACATTAAAATCCAGATCAGCCAGgcgcaagctcagtcaggaacagggcgataccgaaactcaccgggacggcgatcagcagcgcgataaacgaagtcaccgacgtaccgtagatcggcaccagcgccc                                      >  1:221364‑M1/1‑38 (MQ=255)
                                                                                                           cgAGTGATGGTGGACATTAAAATCCAGATCAGCCAGgcgcaagctcagtcaggaacagggcgataccgaaactcaccgggacggcgatcagcagcgcgataaacgaagtcaccaacgtaccgtagatcggcaccagcggcccg                                    >  2:264819‑M1/1‑36 (MQ=255)
                                                                                                            cagcgatggTGCACATCAAAACCCGGACCAGCCCGgcgcaagctcagccaggaacagggcgataccgatcctcaccgggacggcgatcagcagcgcgataaacgaagtcaccaacgtaccgtagatcggcaccagcgccccgt                                   <  1:193072‑M1/139‑109 (MQ=255)
                                                                                                               tgataatagacattaaaATCAAAAAAAGACAGgctcaagatcagtcaggaacaggcagataccgaaactcaccgggacggcgatcatcagcgcgataaacaaagtcaccaacgtaccgtagatcggcaccagcgccccgtaga                                <  1:254899‑M1/135‑112 (MQ=255)
                                                                                                               tGATGGTGGACGTGAAAATCCAGAGCGGCCAGgcgcaagctcagtcaggaacagggcgataccgaaactcaccgggacggcgagcagcagcgcgataaacgaagtcaccaacgtaccgtagatcggcaccagcgccccgtaga                                <  1:798515‑M1/143‑112 (MQ=255)
                                                                                                                    gtggACATTAAAATCCAGATCAGCCAGgcgcaagctcagtcaggaacagggcgataccgaaactcaccgggacggcgatcagcagcgcgataaacgaagtcaccaacgtaccgtagatcggcaccagcgccccgtagatatcg                           <  2:719922‑M1/143‑117 (MQ=255)
                                                                                                                      ggACATTAAAATCCAGATCAGCCAGgcgcaagctcagtcaggaacagtgcgataccgcaactcaccaagacagcaaacagcagcgcgataaacgaagtcaccaacgtaccgtagatcggcacaagcgacccgtagatatcgtt                         >  2:53882‑M1/1‑25 (MQ=255)
                                                                                                                        aCATTAAAATCCAGATCAGCCAGgcgcaagctcagtcaggaacagggcgataccgaaactcaccgggacggcgatcagcagcgcgataaacgaagtcaccaacgtaccgtagatcggcaccagcgccccgtagatatcgttcg                       >  2:261887‑M1/1‑23 (MQ=255)
                                                                                                                          aTTAAAATCCAGATCAGCCAGgcgcaagctcagtcaggaacagggcgataccgaaactcaccgggacggcgatcagcagcgcgataaacgaagtcaccaacgtaccgtagatcggcaccagcgccccgtagatatcgttcggt                     >  2:635918‑M1/1‑21 (MQ=255)
                                                                                                                                  ccAGATCAGCCAGgcgcaagctcagtcaggaacagggcgataccgaaactcaccgggacggcgatcagcagcgcgataaacgaagtcaccaccgtacaatagaacggcaccagagccccgtagataacgcccgatgcagccca             >  1:374556‑M1/1‑13 (MQ=255)
                                                                                                                                    aGATCAGCCGGgcgcaagctcagtcaggaacagggcgataccgaaactcaccgggacggcgatcagcagcgcgataaacgaagtcaccaacgtaccgtagatcggcaccagcgccccgtagatatcgttcggtgcatcccact           <  1:347780‑M1/143‑133 (MQ=255)
                                                                                                                                    aGATCAGCCAGgcgcaagctcagtcaggaacagggcgataccgaaactcaccgggacggcgatcagcagcgcgataaacgaagtcaccaacgtaccgtagatcggcaccagcgccccgtagatatcgttcggtgcatcccact           <  2:680798‑M1/143‑133 (MQ=255)
                                                                                                                                         gggcaggcgcaagctcagtcaggaacagggcgataccgaaacgcaccgggacggcgatcagcagcgcgataaacgaagtcaccaacgtaccgtagatcggcaccagcgccccgtagatatcgttcggtgcatcccactctttg      <  2:726569‑M1/140‑138 (MQ=255)
                                                                                                                                         aGCCAGgcgcggcctcagccgcgagcagggcgataccgaaactccccgggacggcgatcagcagcgcgataaactaagtcaccaacgcaccgtagatcggcaccagcgccccgtagatatcgttcggtgcatcccactctttg      <  1:489884‑M1/143‑138 (MQ=255)
                                                                                                                                           ccAGgcgcaagctcagtcaggaacagggcgataccgaaactcaccgggacggcgatcagtagcgcgataaacgaagtcaccaacgtaccgtagatcggcaccagcgccccggagatatcgttcggtgcatcccacgctttggt    >  2:722244‑M1/1‑4 (MQ=255)
                                                                                                                                            cAGgcgcaagctcagtcaggaacagggcgataccgaaactcaccgggacggcgatcagcagcgcgataaacgaagtcaccaacgtaccgtagatcggcaccagcgccccgtagatatcgttcggtgcatcccacgctttggtc   >  2:510154‑M1/1‑3 (MQ=255)
                                                                                                                                             aGgcgcaagctcagtcaggaacagggcgataccgaaactcaccgggacggcgatcagcagcgcgataaacgaagtcaccaacgtaccgtagatcggcaccagcgccccgtagatatcgttcggtgcatcccactctttggtcc  >  2:359375‑M1/1‑2 (MQ=255)
                                                                                                                                              |                                                                                                                                             
AAATTAACAGCCCGCTACCCGCCAGAGCGTTCGCCAGACCACCACCTTCCGTATTGGGCGGCGGCGTCATTTCAGTGAACAGCGCCAGCGACATACCGTCGATACCGCGAGTGATGGTGGACATTAAAATCCAGATCAGCCAGAACAGGCCGAAGGCCATCGTCGCCATCGAGAGCGTCAGCGCAATACGGTTTTTGAGGCGGCGACGCGCCTGCATTTTGCGGCGAGATTCAGCCAGCGCCGCAGTGGTTTGCATTTCAACCATAGCCATTAGCGTGCCCCCT  >  NC_000913/3909168‑3909451

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 8 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 41 ≤ 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.