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

GACATGATATTGCCGACCGGCTCCTGAAAGTAAACGGCGAACAGCGGCGCAAAGATAAACAGGCCCCACATGCCGTAAACGATACTTGGAATGGCTGCCAGCAGCTCAATGGCGATACCCAGCGGGCGTTTCAGCCAGCCAGGCGCAAGCTCAGTCAGGAACAGGGCGATACCGAAACTCACCGGGACGGCGATCAGCAGCGCGATAAACGAAGTCACCAACGTACCGTAGATCGGCACCAGCGCCCCGTAGATATCGTTCGGTGCATCCCACTCTTTGGTCC  >  NC_000913/3909943‑3910225
                                                                                                                                              |                                                                                                                                            
aattaacagcccgctacccgccagagcgttcgccagaccaccaccttccgtattgggcggcggcgtcatttcagtgaacagcgccagcgacataccgtcgataccgcgagtgatggtggacattaaaatccagatcagccagg                                                                                                                                              >  2:426694‑M2/143‑143 (MQ=255)
 attaacagcccgctacccgccagagcgttcgccagaccaccaccttccgtattgggcggcggcgtcatttcagtgaacagcgccagcgacataccgtcgataccgcgagtgatggtggacattaaaatccagatcagccaggc                                                                                                                                             <  1:371994‑M2/2‑1 (MQ=255)
        gcccgctacccgccagtgctttcgccagtcccccaccgtccgtattgggcggcggcgtcatttcagtgaacagcgccagcgacataccgtctataccgcgagtgatggtgtacattaaaatccagatcagccagGCGCAAGCt                                                                                                                                      <  1:155434‑M2/9‑1 (MQ=255)
           cgctacccgccagagcgttcgccagaccaccaccttccgtattgggcggcggcgtcatttcagtgaacagcgccagcgacataccgtcgatacggcgagcgatggtggacagtaaaatccagatcagccagGCGCAAGCgcag                                                                                                                                   >  1:218705‑M2/132‑143 (MQ=255)
                 ccgccagagcgttcgccagaccaccaccttccgtcttgggcggcggcgccatttcagtgaacagcgccagcgacataccgtcgataccgcgcgtgatggtggacattaaaatccagatcagccagGCGCACGCGCAGTCACGa                                                                                                                             >  1:192678‑M2/126‑143 (MQ=255)
                  cgccagagcgttcgccagaccaccaccttccgtattgggcggcggcgtcatttcagtgaacagcgccagcgacataccgtcgataccgcgagtgatggtggacattaaaatcaagatcagccagACGCAAGCTCAGTCAGGaa                                                                                                                            >  1:544857‑M2/125‑143 (MQ=255)
                            ttcgccagaccaccaccttccgtattgggcggcggcgtcatttcagtgagcagcgccagcgacataccgtcgataccgcgagtgatggtggaccttaaaatccagatcagccagGCTCAAGCTCAGTCAGGAACAGGGCGATa                                                                                                                  >  2:460909‑M2/115‑143 (MQ=255)
                                 cagaccaccaccttccgtattgggcggcggcgtcatttcagtgcacagcgccagcgacatacggtcgataccgcgagtgatggtggacattaagatccagatcagccagGCGAGAGATCAGTCCGGAGCAGGGACATGCCGaa                                                                                                             >  1:246196‑M2/110‑143 (MQ=255)
                                        ccaccttccgtattgggcggcggcgtcatttcagtgaacagcgccagcgacataccgtcgataccgcgagtgatggtggacattaaaatccagatcagccagGCGCAAGCTCAGTCAGGAACAGGGCGATACCGAAACTCAcc                                                                                                      <  2:351291‑M2/41‑1 (MQ=255)
                                                            cggcttcatttcactgaaacgggccagcgacataccgtctataccgcgagtgatggtggacattaaaagccagatcagccagGCGCAAGCTCAGTCAGGAACAGGGCGATACCGAAACTCACCGGGACGGCGATCAGCAgcgc                                                                                  <  1:59259‑M2/61‑1 (MQ=255)
                                                                gtcatttcagtgaacagcgccagcgacataccgtcgataccgcgagtgatggtggacattaaaatccagatcagccagGCGCAAGCTCAGTCAGGAACAGGGCGATCACGAAACTCACCGGGACGGCGATCAGCAGCGCGATa                                                                              >  1:29419‑M2/79‑143 (MQ=255)
                                                                         gtgaacagcgccagcgacataccgtcgataccgcgagtgatggtggacattaaaatccagatcagccagGCGCAAGCTCAGTCAGGAACAGGGCGATACCGAAACTCACCGGGacggggaacagcagcgccaccaacgcaggc                                                                     >  1:4211‑M2/70‑130 (MQ=255)
                                                                          tgaacagcgccagcgacataccgtcgataccgcgagtgatggtggacattaaaatccagatcagccagGCGCAAGCTCAGTCAGGAACAGGGCGATACCGAAACTCACCGGGACGGCGATCAGCAGCGCGATACACGAAGGCa                                                                    >  2:485751‑M2/69‑143 (MQ=255)
                                                                            aacagcgccagcgacataccgtcgataccgcgagtgatgctggacattaaaatccagatcagccagGCGCAAGCTCAGTCAGGAACAGGGCGATACCGAAACTCACCGGGACGGCGATCAGCAGCGCGCTAAACGAAGTCAcc                                                                  >  1:703554‑M2/67‑143 (MQ=255)
                                                                               agcgccagcgacataccgtcgataccgcgagtgatggtggacattcaaatccagatcagccagGCGCAAGCTCAGTCAGGAACAGGGCGATACCGAAACTCACCGCGACGGCGATCAGCAGCGCGATACACGAAGTCACCAAc                                                               >  1:67313‑M2/64‑143 (MQ=255)
                                                                                gcgccagcgacataccgtcgataccgcgagtgatggtggacattaaaatccagatcagccagGCGCAAGCTCAGTCAGGAACAGGGCGATACCGAAACTCACCGGGACGGCGATCAGCAGCGCGATAAACGAAGTCACCAACg                                                              >  1:451694‑M2/63‑143 (MQ=255)
                                                                                        gacacaccgtcgataccgcgagtgagggtggacagtaaaatccagatcagccagGCGCAAGCGCAGTCAGGAACAGGGCGATACCGAAACTCACCGGGACGGCGATCAGCAGCGCGATAAACGAAGTCACCAACGTACCGTAg                                                      <  2:53806‑M2/89‑1 (MQ=255)
                                                                                           ataccgtcgatcccgcgagtcatcgtgcacatcaaaatccagttttgccagGCGCAAGCTCAGTCAGGAACAGGGCGATACCGAAACTCACCGGGACGGCGATCAGCAGCGCGATAAACGAAGTCACCAACGTACCGTAGATc                                                   <  2:30115‑M2/92‑1 (MQ=255)
                                                                                              ccggcggggccgcgggggatggcgggcattaaaagccagatcagcccgGCGCAAGCTCAGTCAGGAACAGGGCGATACCGAAACTCACCGGGACGGCGATCAGCAGCGCGATAAACGAAGTCACCAACGTACCGTAGATCGGc                                                <  2:50579‑M2/95‑1 (MQ=255)
                                                                                                gtcgataccgcgagtgatggtggacattaaaatccagatcagccagGCGCAAGCTCAGTCAGGAACAGGGCGATACCGAAACTCACCGGGACGGCGATCAGCAGCGCGATAAACGAAGTCACCAACGTACCGTAGATCGGCAc                                              >  2:2916‑M2/47‑143 (MQ=255)
                                                                                                    ataccgcgagtgatggtggacattaaaatccagatcagccagGCGCAAGCTCAGTCAGGAACAGGGCGATACCGAAACTCACCGGGACGGCGACCAGCAGCGCGATAAACGCAGACCCCAACGTACCGCACAACGGCACCAgc                                          >  1:246177‑M2/43‑143 (MQ=255)
                                                                                                        cgcgagtgatggtggacattaaaatccagatcagccagGCGCAAGCTCAGTCAGGAACAGGGCGATACCGAAACTCACCGGGACGGCGATCAGCAGCGCGATAAACGAAGTCACCGACGTACCGTAGATCGGCACCAGCGccc                                      >  1:221364‑M2/39‑143 (MQ=255)
                                                                                                          cgagtgatggtggacattaaaatccagatcagccagGCGCAAGCTCAGTCAGGAACAGGGCGATACCGAAACTCACCGGGACGGCGATCAGCAGCGCGATAAACGAAGTCACCAACGTACCGTAGATCGGCACCAGCGGCCCg                                    >  2:264819‑M2/37‑143 (MQ=255)
                                                                                                           cagcgatggtgcacatcaaaacccggaccagcccgGCGCAAGCTCAGCCAGGAACAGGGCGATACCGATCCTCACCGGGACGGCGATCAGCAGCGCGATAAACGAAGTCACCAACGTACCGTAGATCGGCACCAGCGCCCCGt                                   <  1:193072‑M2/108‑1 (MQ=255)
                                                                                                              tgatggtggacgtgaaaatccagagcggccagGCGCAAGCTCAGTCAGGAACAGGGCGATACCGAAACTCACCGGGACGGCGAGCAGCAGCGCGATAAACGAAGTCACCAACGTACCGTAGATCGGCACCAGCGCCCCGTAGa                                <  1:798515‑M2/111‑1 (MQ=255)
                                                                                                              tgataatagacattaaaatcaaaaaaagacagGCTCAAGATCAGTCAGGAACAGGCAGATACCGAAACTCACCGGGACGGCGATCATCAGCGCGATAAACAAAGTCACCAACGTACCGTAGATCGGCACCAGCGCCCCGTAGa                                <  1:254899‑M2/111‑1 (MQ=255)
                                                                                                                   gtggacattaaaatccagatcagccagGCGCAAGCTCAGTCAGGAACAGGGCGATACCGAAACTCACCGGGACGGCGATCAGCAGCGCGATAAACGAAGTCACCAACGTACCGTAGATCGGCACCAGCGCCCCGTAGATATCg                           <  2:719922‑M2/116‑1 (MQ=255)
                                                                                                                     ggacattaaaatccagatcagccagGCGCAAGCTCAGTCAGGAACAGTGCGATACCGCAACTCACCAAGACAGCAAACAGCAGCGCGATAAACGAAGTCACCAACGTACCGTAGATCGGCACAAGCGACCCGTAGATATCGtt                         >  2:53882‑M2/26‑143 (MQ=255)
                                                                                                                       acattaaaatccagatcagccagGCGCAAGCTCAGTCAGGAACAGGGCGATACCGAAACTCACCGGGACGGCGATCAGCAGCGCGATAAACGAAGTCACCAACGTACCGTAGATCGGCACCAGCGCCCCGTAGATATCGTTCg                       >  2:261887‑M2/24‑143 (MQ=255)
                                                                                                                         attaaaatccagatcagccagGCGCAAGCTCAGTCAGGAACAGGGCGATACCGAAACTCACCGGGACGGCGATCAGCAGCGCGATAAACGAAGTCACCAACGTACCGTAGATCGGCACCAGCGCCCCGTAGATATCGTTCGGt                     >  2:635918‑M2/22‑143 (MQ=255)
                                                                                                                                ccccgaccacctacgcgcAACCGCGGTCAGGAACAGGGCGATACCGAAACTCCCCGGGACGGCGATCACCAGCCCGATAAACGTAGTCACCAACGTACCCTAGATCGGCACCAGCGCTCCGTAGATATCGTTCGGTGCATccc              <  1:674103/129‑1 (MQ=255)
                                                                                                                                 ccagatcagccagGCGCAAGCTCAGTCAGGAACAGGGCGATACCGAAACTCACCGGGACGGCGATCAGCAGCGCGATAAACGAAGTCACCACCGTACAATAGAACGGCACCAGAGCCCCGTAGATAACGCCCGATGCAGCCCa             >  1:374556‑M2/14‑143 (MQ=255)
                                                                                                                                   agatcagccggGCGCAAGCTCAGTCAGGAACAGGGCGATACCGAAACTCACCGGGACGGCGATCAGCAGCGCGATAAACGAAGTCACCAACGTACCGTAGATCGGCACCAGCGCCCCGTAGATATCGTTCGGTGCATCCCAct           <  1:347780‑M2/132‑1 (MQ=255)
                                                                                                                                   agatcagccagGCGCAAGCTCAGTCAGGAACAGGGCGATACCGAAACTCACCGGGACGGCGATCAGCAGCGCGATAAACGAAGTCACCAACGTACCGTAGATCGGCACCAGCGCCCCGTAGATATCGTTCGGTGCATCCCAct           <  2:680798‑M2/132‑1 (MQ=255)
                                                                                                                                        gggcagGCGCAAGCTCAGTCAGGAACAGGGCGATACCGAAACGCACCGGGACGGCGATCAGCAGCGCGATAAACGAAGTCACCAACGTACCGTAGATCGGCACCAGCGCCCCGTAGATATCGTTCGGTGCATCCCACTCTTTg      <  2:726569‑M2/137‑1 (MQ=255)
                                                                                                                                        agccagGCGCGGCCTCAGCCGCGAGCAGGGCGATACCGAAACTCCCCGGGACGGCGATCAGCAGCGCGATAAACTAAGTCACCAACGCACCGTAGATCGGCACCAGCGCCCCGTAGATATCGTTCGGTGCATCCCACTCTTTg      <  1:489884‑M2/137‑1 (MQ=255)
                                                                                                                                          ccagGCGCAAGCTCAGTCAGGAACAGGGCGATACCGAAACTCACCGGGACGGCGATCAGTAGCGCGATAAACGAAGTCACCAACGTACCGTAGATCGGCACCAGCGCCCCGGAGATATCGTTCGGTGCATCCCACGCTTTGGt    >  2:722244‑M2/5‑143 (MQ=255)
                                                                                                                                           cagGCGCAAGCTCAGTCAGGAACAGGGCGATACCGAAACTCACCGGGACGGCGATCAGCAGCGCGATAAACGAAGTCACCAACGTACCGTAGATCGGCACCAGCGCCCCGTAGATATCGTTCGGTGCATCCCACGCTTTGGTc   >  2:510154‑M2/4‑143 (MQ=255)
                                                                                                                                            agGCGCAAGCTCAGTCAGGAACAGGGCGATACCGAAACTCACCGGGACGGCGATCAGCAGCGCGATAAACGAAGTCACCAACGTACCGTAGATCGGCACCAGCGCCCCGTAGATATCGTTCGGTGCATCCCACTCTTTGGTcc  >  2:359375‑M2/3‑143 (MQ=255)
                                                                                                                                              |                                                                                                                                            
GACATGATATTGCCGACCGGCTCCTGAAAGTAAACGGCGAACAGCGGCGCAAAGATAAACAGGCCCCACATGCCGTAAACGATACTTGGAATGGCTGCCAGCAGCTCAATGGCGATACCCAGCGGGCGTTTCAGCCAGCCAGGCGCAAGCTCAGTCAGGAACAGGGCGATACCGAAACTCACCGGGACGGCGATCAGCAGCGCGATAAACGAAGTCACCAACGTACCGTAGATCGGCACCAGCGCCCCGTAGATATCGTTCGGTGCATCCCACTCTTTGGTCC  >  NC_000913/3909943‑3910225

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.