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 50 [0] [0] 50 [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)48 (0.950) 42/268 0.1 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

GACATGATATTGCCGACCGGCTCCTGAAAGTAAACGGCGAACAGCGGCGCAAAGATAAACAGGCCCCACATGCCGTAAACGATACTTGGAATGGCTGCCAGCAGCTCAATGGCGATACCCAGCGGGCGTTTCAGCCAGCCAGGCGCAAGCTCAGTCAGGAACAGGGCGATACCGAAACTCACCGGGACGGCGATCAGCAGCGCGATAAACGAAGTCACCAACGTACCGTAGATCGGCACCAGCGCCCCGTAGATATCGTTCGGTGCATCCCACTCTTTGGTC  >  NC_000913/3909943‑3910224
                                                                                                                                              |                                                                                                                                           
aattaacagcccgctacccgccagagcgttcgccagaccaccaccttccgtattgggcggcggcgtcatttcagtgaacagcgccagcgacataccgtcgataccgcgagtgatgacgggcatcaaactcccgatcagccagg                                                                                                                                             >  1:713767‑M2/143‑143 (MQ=255)
   ttacagcccgctacccgccatagcgttcgccagaccagcgccttccggattgggcggcggcgtcatttcagtgaacagcgccagcgacataccgtcgataccgcgagtgatggtggacattaaaatccagatcagccaggcgc                                                                                                                                          <  1:609451‑M2/4‑1 (MQ=255)
          ccgctacccgccagagcgttcgccagaccaccaccttccgtattgggcggcggcgtcatttcagtgaacagcgccagcgatataccgtcgataccgcgagtgatggtggacattaaaatccagatcagccaaGCGAAAGctcg                                                                                                                                   >  1:466183‑M2/133‑142 (MQ=255)
                       gagcgttcgccagaccaccaccttccgtattgggcggcggcgtcatttcagtgaacagcgccagcgacataccgtcgataccgcgagtgatggtggacattaaaatccagatcagccagGCGCAAGCTCAGTCAGGAACAggg                                                                                                                      <  1:753444‑M2/24‑1 (MQ=255)
                               gccagaccaccaccttccgtattgggcggcggcgtcatttcagtgaacagcgccagcgacataccgtcgataccgcgagtgatggtggacattaaaatccagatcagccagGCGCAAGCTCAGTCAGGAACAGGGCGATACCg                                                                                                              >  1:628390‑M2/112‑143 (MQ=255)
                                       accccctcccgtatcgggcggcggcgtcattgcagtgaccagcgccagcgacataccgtcgataccgcgagtgatggtggacattaaaatccagatcagccagGCGCAAGCTCAGTCAGGAACAGGGCGATACCGAAACTCAc                                                                                                      <  2:658207‑M2/40‑1 (MQ=255)
                                          accttccgtattgggcggcggcgtcatttcagtgaacagcgccagcgacataccgtcgataccgcgagtgatggtggacattaaaatccagatcagccagGCGCAAGCTCAGTCAGGAACAGGGCGATACCGAAACTCAc                                                                                                      >  2:456732‑M2/101‑140 (MQ=255)
                                          accttccgtattgggcggcggcgtcatttcagtgaacagcgccagcgacataccgtcgataccgcgagtgatggtggacattaaaatccagatcagccagGCGCAAGCTCAGTCAGGAACAGGGCGATACCGAAACTCAc                                                                                                      <  1:456732‑M2/40‑1 (MQ=255)
                                            cttccgtattgggcggcggcgtcagttcagggagcggctccagcgacataccggcgataccgcgagtgatggtggacattaaaagccggatcagccagGCGCAAGCTCAGTCAGGAACAGGGCGATACCGGAACTCGCCGGGa                                                                                                 <  2:621938‑M2/45‑1 (MQ=255)
                                             ttccgtattgggcggcggcgtcatttcagtgaacagcgccagcgacataccgtcgataccgcgagtgatggtggacattaaaatccagatcagccagGCGCAAGCTCAGTCAGGAACAGGGCGATACCGAAACTCACCGGGAc                                                                                                >  1:93202‑M2/98‑143 (MQ=255)
                                             ttccgtattgggcggcggcgtcatttcagtgaacagcgccagcgacataccgtcgataccgcgagtgatggtggacattaaaatccagatcagccagGCGCAAGCGCAGTCAGGACCAGGACGATACCGCACCTCACCGCGAc                                                                                                >  2:649191‑M2/98‑143 (MQ=255)
                                               tcctattgggcggtggcgtttttttagtgagctgcgccagcgacctaccgtcgataccgcgagtgatgttggacattaaaatccagatcagccagGCGCAAGCTCAGTAAAGAACAGGCCTATACCTAAACTCACCGGGACgg                                                                                              <  2:332218‑M2/48‑1 (MQ=255)
                                                 gtattgggcggcggcgtcatttcagtgaacagcgccagcgacataccgtcgataccgcgagtgatggtggacattaaaatccagatcagccagGCGCAAGCTCAGTCAGGAACAGGGCGATACCGAAACTCACCGGGACGGCg                                                                                            >  1:909641‑M2/94‑143 (MQ=255)
                                                  tattgggcggcggcgtcatttcagtgaacagcgccagcgacataccgtcgataccgcgagtgatggtggacattaaaatccagatcagccagGCGCAAGCTCAGTCAGGAACAGGGCGATACCGAAACTCACCGGGACGGCGa                                                                                           >  1:101646‑M2/93‑143 (MQ=255)
                                                  tagcgggcggcggcgtcatttcagtgaacagcgccagcgacataccggcgataccgcgagtgatggtggacattaaaatccagatcagccagGCGCAAGCTCAGTCAGGAACAGGGCGATACCGAAACTCACCGGGACGGCGa                                                                                           <  1:407049‑M2/51‑1 (MQ=255)
                                                   attgggcggcggcgtcatttcagtgaacagcgccagcgacataccgtcgataccgcgagtgatggtggacattaaaatccagatcagccagGCGCAAGCTCAGTCAGGAACAGGGCGATACCGAAACTCACCGGGACGGCGAt                                                                                          <  1:623220‑M2/52‑1 (MQ=255)
                                                   atggggcggcggcgtcattccagtgaacagcgccagcgacataccgtcgataccgcgagtgatggtggacattaaaatccagatcagctagGCGCAAGCTCAGTCAGGAACAGGGCGATACCGAAGCTCACCGGGACGGCGAt                                                                                          <  1:484574‑M2/52‑1 (MQ=255)
                                                      gggcggcggcgtcatttcagtgaacagcgccagcgacataccgtcgataccgcgagtgatggtggacattaaaatccagatcagccagGCGCAAGCTCAGTCAGGAACAGGGCGATACCGAAACTCACCGGGACGGCGATcag                                                                                       <  1:561916‑M2/55‑1 (MQ=255)
                                                           gcggcgtcatttcagtgaacagcgccagcgacataccgtcgataccgcgagtgatggtggacattaaaatccagatcagccagGCGCAAGCTCAGTCAGGAACAGGGCGATACCGAAACTCACCGGGACGGCGATCAGCAgcg                                                                                  <  1:812032‑M2/60‑1 (MQ=255)
                                                             ggcgtcatttcagtgaacagcgccagcgacataccgtcgataccgcgagtgatggtggacattaaaatccagatcagccagGCGCAAGCTCAGTCAGGAACAGGGCGATACCGAAACTCACCGGGACGGCGATCAGCAGcgcg                                                                                <  1:473383‑M2/62‑1 (MQ=255)
                                                                    tggcagtgagcagcgccagcgacataccggcgatgccgcgagtgatggtggacattaaaatccagatcagcaagGCGCAAGCTCAGTCAGGAACAGGGCGATACCGAAACTAACCGGGACGGCGATCAGCAGCGCGATAAACg                                                                         <  2:253131‑M2/69‑1 (MQ=255)
                                                                      tcagtgaacagcgccagcgacagaccgtcgataccgcgagtgatggtggacattaaaatccagatcagccagGCGCAAGCTCAGTCAGGAACAGGGCGATACCGAGACTCACCGGGACGGCGATCAGCAGCGCGATAAACGaa                                                                       <  1:818682‑M2/71‑1 (MQ=255)
                                                                         gtgaacagcgccagcgacataccgtcgataccgcgagtgatggtggacattaaaatccagatcagccagGCGCAAACTCAGTCAGGAACAGGGCGATACCGAAACTCACCGGGACGGCGATCAGCAGCGCGATCGACGAAGTc                                                                    >  2:756123‑M2/70‑143 (MQ=255)
                                                                          ttacccgcgccagcgacagaccggcgataccgcgatcgagtgtggacattaaaatccagatcagccagGCGCAAGCTCAGTCAGGTACAGGGCGATATCGAAACTCACCGTGACGGCGATCAGCAGCGCGATAAACGAAGTCa                                                                   <  1:103759‑M2/75‑1 (MQ=255)
                                                                          tgaacagcgccagcgacataccgtcgataccgcgagtgatggtggacattaaaatccagatcagccagGCGCAAGCTCAGTCAGGAACAGGGCGATACCGAAACTCACCGGGACGGCGATCAGCAGCGCGATAAACGAAGTCa                                                                   <  1:248148‑M2/75‑1 (MQ=255)
                                                                            aacagcgccagcgacataccgtcgataccgcgagtgatggtggacattaaaatccagatcagccagGCGCAAGCTCAGTCAGGAACAGGGCGATACCGAAACTCACCCGGAAGGCGTTCAGCCTCGCGATAAACGAAGTCAcc                                                                 >  1:722260‑M2/67‑143 (MQ=255)
                                                                              cagcgccagcgacataccgtcgataccgcgagtgatggtggacattaaaatccagatcagccagGCGCAAGCTCAGTCAGGAACAGGGCGATACCGAAACTCACCGGGACGGCGATCAGCAGCGCGATAAACGAAGTCACCaa                                                               >  2:527728‑M2/65‑143 (MQ=255)
                                                                               agcgccagcgacataccgtcgataccgcgagtgatggtggacattaaaatccagatcagccagGCGCAAGCTCAGTCAGGAACAAGGCGATAATTAAACTCACCGGCACGGAGCACAGCAGCGCAAGCAAAGAAGTAACCAAc                                                              >  2:107004‑M2/64‑143 (MQ=255)
                                                                                   ccagcgacataccgtcgagaccgcgagtgagggtggacattaaaagccagatcagccagGCGCAAGCTCAGTCAGGAACAGGGCGATACCGAAACTCACCGGGACGACGCACCGCCGCGCGATAAACGAAGTCACCAACGTAc                                                          <  2:566187‑M2/84‑1 (MQ=255)
                                                                                      gcgacataccgccgataccgcgagtgacggtggacattaaaatccagttcagccagCCGCAAGCTCAGTCAGGAACAGGGCGATACCGACACTCACCGGGACGGCGATCAGCAGCGCGATAAACGAAGTCACCAACGTACCGt                                                       <  1:110940‑M2/87‑1 (MQ=255)
                                                                                            taccgtcgataccgcgagtgatggtggacattaaaatccagatcagccagGCGCAAGCGCAGTCAGGAACAGTGCGATACCGAAACTCACCGGGACGGCGATCAGCATCGCGATAAACGAAGTCACCGACGTACCGTAgatag                                                 >  1:444306‑M2/51‑141 (MQ=255)
                                                                                             accgtcgataccgcgagtgatggtggacattaaaatccagatcagccagGCGCAAGCTCAGTCAGGAACAGGGCGATACCGAAACTCACCGGGACGGCGATCAGCAGCGCGATAAACGAAGTCACCAACGTACCGTAGATCgg                                                <  2:101646‑M2/94‑1 (MQ=255)
                                                                                              ccgtcgataccgcgagtgatggtggacattaaaatccagatcagccagGCGCAAGCTCAGTCAGGAACAGGGCTATACCGAAACTCACCGGGACGGCGATCAGCAGCGCGATAAACGAAGTCACCAACGTACCGTAGATCGGc                                               >  1:567397‑M2/49‑143 (MQ=255)
                                                                                                    ataccgcgagtgatggcggacattaaaatccagaacagccagGCGCAAGCTCAGTCAGGAACAGGGCGATACCGAAACTCACCGGGACGGCGATCAGCAGCGCGATAAACGAAGTCACCAACGTACCGTAGATCGGCACCAgc                                         >  1:708962‑M2/43‑143 (MQ=255)
                                                                                                        cgcgggggatggtggacattaaaatccagatcagccagGCGCAAGCTCAGTCAGGAACAGGGCGATACCGAAACTCACCGGGACGGCGATCAGCAGCGCGATAAACGAAGTCACCAACGTACCGTAGATCGGCACCAGCGccc                                     <  2:470108‑M2/105‑1 (MQ=255)
                                                                                                           gagtcatgagggacagtacacgccagaccagccacGCGCAAGCTCAGTCAGGAACAGCGCGGTACCGAACGTCACCGGGACGGCGATCAGCAGCGCGATAAACGAAGTCACCAACGTACCGTAGATCGGCACCAGCGCCCCGt                                  <  2:505363‑M2/108‑1 (MQ=255)
                                                                                                            agtgatggtggacattaaaatccagatcagccagGCGCAAGCTCAGTCAGGAACAGGGCGATACCGAAACTCACCGGGACGGCGATCAGCAGCGCGATAAACGAAGTCACCAACGTACCGTAGATCGGCACCAGCGCCCCGTa                                 <  1:102536‑M2/109‑1 (MQ=255)
                                                                                                              tgatggtggacattaaaatccagatcagccagGCGCAAGCTCAGTCAGGAACAGGGCGATACCGAAACGCACCGGGACGGCGATCAGCAGCGCGATAAACGAAGTCACCAACGTACGGTACATCGGCACCAGCGCCCGGTGGa                               >  2:675023‑M2/33‑143 (MQ=255)
                                                                                                                  tgtgtcccttacaatccagctccgccacGCGCAAGCTCCGTCAGGAACAGGGCGATACCGAAACTCACCGGGACGGCGATCAGCAGCGCGATAAACGAAGTCACCAACGTACCGTAGATCGGCACCAGCGCCCCGTAGATATc                           <  1:878354‑M2/115‑1 (MQ=255)
                                                                                                                    tggacattaaaatccagatcagccagGCGCAAGCTCAGTCAGGAACAGGGCGATACCGAAACTCACCGGGACGGCGATCAGCAGCGCGATAAACGAAGTCACCAACGTACCGTCGATCGGCACCAGCGCCCCGTAGCTATCGt                         >  1:181914‑M2/27‑143 (MQ=255)
                                                                                                                    tggacattaaaatccagatcagccagGCGCAAGCTCAGTCAGGAACAGGGCAATACCAAAACACACCGGGACGGCGATCAGCAGCGCGATAAACGAAGTCACCAACGTACCGTAGATCGGCACCAGCGCCCCGTAGATATCGt                         >  1:264276‑M2/27‑143 (MQ=255)
                                                                                                                       acattagaatccagatcagccagGCGCAAGCTCAGTCAGGAACAGGGCGATACCGAAACTCACCGGGACGGCGATCAGCAGCGCGATAAACGAAGTCACCAACGTACCGTAGATCGGCACCAGCGCCCCGTAGATATCGTTCg                      <  2:901907‑M2/120‑1 (MQ=255)
                                                                                                                        ccttaaaatccagatcgcccagGCGCAGGCGCAGTCAGGAACAGGGCGATACCGAAACTCACCGGGACGGCGATCAGCAGCGCGATAAACGAAGTCACCAACGTACCGTAGATCGGCACCAGCGCCCCGTAGATATCGTTCgg                     <  2:131351‑M2/121‑1 (MQ=255)
                                                                                                                         attaaaatccagatcagccagGCGCAAGCTCAGTCAGGAACAGGGCGATACCGAAACTCACCGGGACGGCGATCAGCAGCGCGATAAACGaa                                                                       >  2:335040‑M2/22‑92 (MQ=255)
                                                                                                                         attaaaatccagatcagccagGCGCAAGCTCAGTCAGGAACAGGGCGATACCGAAACTCACCGGGACGGCGATCAGCAGCGCGATAAACGaa                                                                       <  1:335040‑M2/71‑1 (MQ=255)
                                                                                                                         attaaaacccagatcagccagGCGCAAGCTCAGTCAGGAACAGGGCGATACCGAAACTCACCGGGACGGCGATCAGCAGCGCGATAAACGAAGTCACCAACGTACCGTAGATCGGCACCAGCGCCCCGTAGATATCGTTCGGt                    <  1:536190‑M2/122‑1 (MQ=255)
                                                                                                                                 ccagatcagccagGCGCAAGCTCAGTCAGGAACAGGGCGATACCGAAACTCACCGGGACGGCGATCAGCAGCGCGATAAACGAAGTCACCAACGTACCGTAGATCGGCACCAGCGCCCCGTAGATATCGTTCGGTGCATCCCa            >  1:365145‑M2/14‑143 (MQ=255)
                                                                                                                                   agatcagccagGCGCAAGCTCAGTCAGGAACAGGGCGATACCGAAACTCACCGGGACGGCGATCAGCAGCGCGATAAACGCAGTCACCAACGTCGCGTAGATGGGCACCAGCGCCCCGTAGAAATTATTTGGtgcatcacgcg          >  1:465988‑M2/12‑138 (MQ=255)
                                                                                                                                          ccagGCGCACGCCCGGCCAGGAACAGAGCGATACCGTAACTCACCGGGACGGCGATCAGCAGCGCGATAAACGACGTCACCAACGTACCGTAGATCGGCACCAGCGCCCCGTAGATATCGTTCGGTGCATCCCACTCTTTGGt   <  2:592993‑M2/139‑1 (MQ=255)
                                                                                                                                           cagGCGCAAGCTCAGTCAGGAACAGGGCGATACCGAACCTCACCCGGGCGGCGATCAGCAGCGCGATAAACGAAGTCACCAACGTACCGTAGATCGGCACCAGCGCCCCGTAGATATCGTTCGGTGCATCCCACTCTTTGGTc  <  1:253601‑M2/140‑1 (MQ=255)
                                                                                                                                              |                                                                                                                                           
GACATGATATTGCCGACCGGCTCCTGAAAGTAAACGGCGAACAGCGGCGCAAAGATAAACAGGCCCCACATGCCGTAAACGATACTTGGAATGGCTGCCAGCAGCTCAATGGCGATACCCAGCGGGCGTTTCAGCCAGCCAGGCGCAAGCTCAGTCAGGAACAGGGCGATACCGAAACTCACCGGGACGGCGATCAGCAGCGCGATAAACGAAGTCACCAACGTACCGTAGATCGGCACCAGCGCCCCGTAGATATCGTTCGGTGCATCCCACTCTTTGGTC  >  NC_000913/3909943‑3910224

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.