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 48 [0] [0] 47 [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)43 (0.900) 37/268 0.2 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

ATATTGCCGACCGGCTCCTGAAAGTAAACGGCGAACAGCGGCGCAAAGATAAACAGGCCCCACATGCCGTAAACGATACTTGGAATGGCTGCCAGCAGCTCAATGGCGATACCCAGCGGGCGTTTCAGCCAGCCAGGCGCAAGCTCAGTCAGGAACAGGGCGATACCGAAACTCACCGGGACGGCGATCAGCAGCGCGATAAACGAAGTCACCAACGTACCGTAGATCGGCACCAGCGCCCCGTAGATATCGTTCGGTGCATCCCACTCTTTGGTCCAT  >  NC_000913/3909949‑3910227
                                                                                                                                        |                                                                                                                                              
cagcccgctacccgccagagcgttcgccagaccaccaccgtccgtattgggcggcggcgtcatttcagtgaacagcgccagcgacataccgtcgataccgcgagtgatggtggacattaaaagcgagctcagccagGCGCAAg                                                                                                                                          >  1:872884‑M2/137‑143 (MQ=255)
  gccctctgcgcgccagggcgttcgccagaccaccaccttccgtattgagcggcggcttcatttcagtgaggagcgccagcgacgtgtcgtcgacaccgcgagtgatggtggacattaaaatccagatcagccagGCGCAAGCt                                                                                                                                        <  1:810586‑M2/9‑1 (MQ=255)
   cccgctacccgccagagcgttcgccagaccaccaccttccgtattgggcggcggcgtcatttcagtgaacagcgccagcgacataccgtcgataccgcgagtgatggtggacattaaaatccagatcagccacGCGCAAGCtc                                                                                                                                       >  2:693273‑M2/134‑143 (MQ=255)
    ccgctacccgccagagcgttcgccagaccaccaccttccgtcttgggcggcggcgtcacttcagtgaacagcgccagcgacataccgtcgataccgcgagtgaaggtggacattaaaatccagatcagccagGCGCAAGCtca                                                                                                                                      >  1:195743‑M2/133‑143 (MQ=255)
    ccgctacccgccagagcgttcgccagaccaccaccttccgtattgggcggcggcgtcatttcagtgaacagcgccagcgacataccgtcgataccgcgagtgatggtggacattaaaatccagatcagcccgGCGCAAGctcg                                                                                                                                      >  1:194787‑M2/133‑142 (MQ=255)
           ccgccatagcgttcgccggaccacccccttccgtattgggcggcggcgtcagttcagtgaacagcgccagcgacataccgtcgataccgcgagtgatggtggacattaaaatccagatcagccagGCGCAAGCTCAGTCAGGa                                                                                                                               <  2:254219‑M2/18‑1 (MQ=255)
             gccagagcgttcgccagaccaccaccttccgtattgggcggcggcgtcatttcagtgaacagcgccagcgacataccgtcgataccgcgagtgatggtggacattcacatccagatcagccggGCGCAGGCCCGGTCAGTAAc                                                                                                                             >  2:780978‑M2/124‑143 (MQ=255)
                agagcgttcgccagaccaccaccttccgtattgggcggcggcgtcatttcagtgaacagcgccagcgacataccgtcgataccgcgagtgatggtggacattaaaatccagatcagccagGCGCAAGCTCCGTCAGGAGCAgg                                                                                                                          >  2:577802‑M2/121‑143 (MQ=255)
                agagcgttcgccagaccaccaccttccgtattgggcggcggcgtcatttcagtgaacagcgccagcgacataccgtcgataccgcgagtgatggtggacattaaaatacagatcagccagGCGCAAGCTCAGGCAGGAacaag                                                                                                                          >  1:337308‑M2/121‑141 (MQ=255)
                         gccagaccaccacctgcggtagtgggcggcggcgtcattgcagtgaacagcgccagcgacataccgtcgataccgcgagtgatggtggacattaaaatccagatcagccagGCGCAAGCTCAGTCAGGAACAGGGCGATACCg                                                                                                                 <  2:312515‑M2/32‑1 (MQ=255)
                              accaccaccttccgtattgggcggcggcgtccttacagtgaacagcgcccgcgacataccgtcgataccgcgagtgagggtggacagtaacatccagatcagcccgCCGCAAGCTCAGTCAGGCAAACGGCAACAccgaacga                                                                                                            >  1:42827‑M2/107‑140 (MQ=255)
                              accaccaccttccgtattgggcggcggcgtcatttcagtgaacagcgccagcgacataccgtcgataccgcgagtgatggtggacattaaaatccagatcagccagGCGCAAGCTCAGTCAGGAACAGGGCGATACCGAAACt                                                                                                            <  1:337156‑M2/37‑1 (MQ=255)
                                caccaccttccgtattgggcggcggcgtcatttcagtgaacagcgccagcgacataccgtcgataccgcgagtgatggtggacattaaaatccagatcagccagGCGCAAGCTCAGTCAGGAACAGGGCGATACCGAAACTCa                                                                                                          <  2:577333‑M2/39‑1 (MQ=255)
                                caccaccttccgtattgggcggcggcgtcatttcagtgaacagcgccagcgacataccgtcgataccgcgagtgatggtggacattaaaatccagatcagccagGCGCAAGCTCAGTCAGGAACAGGGCGATACCGAAACTCa                                                                                                          >  2:32088‑M2/105‑143 (MQ=255)
                                 accaccttccgtatcgggcggcggcgtcatttcagtgaacagcgccagcgacataccgtcgataccgcgagtgatggtggacattaaaatccagatcagccagGCGCAAGCTCAGTCAGGAACAGGGCGATACCGAAACTCAc                                                                                                         <  1:359371‑M2/40‑1 (MQ=255)
                                 accaccgtccgtattgggcggcggcgtcatttcagtgaacagcgccagcgacataccgtcgataccgcgagtgatggtggacattaaaatccagatcagccagGCGCAAGCTCAGTCAGGAACAGGGCGATACCGAAACTCAc                                                                                                         <  1:656171‑M2/40‑1 (MQ=255)
                                             cttgggcggcggcgtcatttgagtgaacagcgccagcgacataccgtcgataccgcgagtgatggtggacattaaaatccagatcagccagGCGCAAGCTCAGTCAGGAACAGGGCGATACCGAAACTCACCGGGACGGCGAt                                                                                             <  1:150179‑M2/52‑1 (MQ=255)
                                                    ggcggcgtcatttcagtgaacagcgccagcgacataccgtcgataccgcgagtgatggtggacattacaatccagatcacccagGCGCAAGCTCAGCCAGCACCAGCGCCATACCGAACCTCCCCGGGCCCgcgatcacccgc                                                                                      >  1:204804‑M2/85‑138 (MQ=255)
                                                      cggcgtcatttcagtgaacagcgccagcgacataccgtcgataccgcgagtgatggtggacattaaaatccagatcagccagGCGCAAGCTCAGTCAGGAACAGGGCGATACCGAAACACACCGGGACGGCGATCAGCAgcgc                                                                                    >  1:514018‑M2/83‑143 (MQ=255)
                                                      cggcgtcatttcagtgaacagcgccagcgacataccgtcgataccgcgagtgatggtggacattaaaatccagagcagccagGCGCAAGCTCAGTCAGGAACAGGGCGATACCGAAACTCACCGGGACGGCGATCAGCAgcgc                                                                                    <  2:416319‑M2/61‑1 (MQ=255)
                                                       ggcgtcatttcagtgaacagcgccagcgacataccgtcgataccgcgagtgatgggggacattaaaatccagatcagccagGCGCAAGCTCAGTCAGGAACAGGGCGATACCGAAACTCACCGGGACGGCGATCAGCAGcgcg                                                                                   <  2:860405‑M2/62‑1 (MQ=255)
                                                        gcgtcacttcagtgaacagcgccagcgacataccgtcgataccgcgagtgatggtggacattaaaatccagatcagccagGCGCAAGCTCAGTCAGGAACAGGGCGATACCGAAACTCACCGGGACGGCGATCAGCAGCGCGa                                                                                  <  1:218328‑M2/63‑1 (MQ=255)
                                                           tcatttcagtgaacagcgccagccacacacctttgctaccgcgagtgatggtggacattaaagtccagatcagccagGCGCAAGCTCAGTCAGGAACAGGGCGATACCGAAACTCACCGGGACGGCGATCAGCAGCGCGATaa                                                                               <  2:594210‑M2/66‑1 (MQ=255)
                                                            gacttcagtgcacagcgccagcgacataccgtcgataccgcgagtcatggtggacatgaaaatccagatcagccagGCGCAAGCTCAGTCAGGAACAGGGCGATACCGAAACTCACCGGGACGGCGATCAGCAGCGCGATaaa                                                                              <  1:688166‑M2/67‑1 (MQ=255)
                                                                tcagtgaacagcgccagcgacataccgtcgataccgcgagtgatggtggacattaaaatccagatcagccagGCGCAAGCTCAGTCAGGAACAGGGCGATACCGAAACTCACCGGGACGGCGATCAGCAGCGCGATAAACGaa                                                                          <  2:2266‑M2/71‑1 (MQ=255)
                                                                    tgaacagcgccagcgacataccgtcgataccgcgagtgatggtggacattaaaatccagatcagccagGCGCAAGCTCAGTCAGGAACAGGGCGATACCGAAACTCACCGGGACGGCGATCAGCAGCGCGATAAACGAAGTCa                                                                      >  2:286708‑M2/69‑143 (MQ=255)
                                                                    tgaacagcgccagcgacagcccgtcgataccgctagtgatggtgtacattaaaatccagatcagccagGCGCAAGCTCAGTCAGGAACAGGGCGATACCGAAACTCACCGGGACGGCGATCAGCAGCGCGATAAACGAAGTCa                                                                      <  1:47542‑M2/75‑1 (MQ=255)
                                                                         agcgccagcgacatcccggcgataccgcgagtgatggtggacattaaaatccagatcagccagGCGCAAGCTCAGTCAGGAACAGGGCGATACCGAAACTCACCGGGACGGCGATCAGCAGCGCGATAAACGAAGTCACCAAc                                                                 <  2:186780‑M2/80‑1 (MQ=255)
                                                                          gcgccagcgacatcccggcgataccgggagtgatggtggacattaaaatccagatcagccggGCGCGAGGTCAGTCAGGAACAGGGCGATACCGAAACTCACCGGGACGGCGATCAGCAGCGCGATAAACGAAGTCACCAACg                                                                <  2:405481‑M2/81‑1 (MQ=255)
                                                                             ccatccacttaccttcgataccgcgagtgatggtggacattaaaatccagatcagccagGCGCAAGCTCAGTCAGGAACAGGGCGATACCGAAACTCACCGGGACGGCGATCAGCAGCGCGATAAACGAAGTCACCAACGTAc                                                             <  1:566677‑M2/84‑1 (MQ=255)
                                                                                gcgacataccgtcgataccgcgagtgatggtggacattaaaatccagatcagccagGCGCAAGCTCAGTCAGGAACAGGGCGATACCGAAACTCACCGGGACGGCGATCAGCAGCGCGATAAACGAAGTCACCAACGTACCGt                                                          >  2:472198‑M2/57‑143 (MQ=255)
                                                                                  gacataccgtcgataccgcgagtgatggtggacattaaaatccagatcagccagGCGCAAGCTCAGTCAGGAACAGGGCGATACCGAAACTCACCGGGACGGCGATCAGCAGCGCGATAAACGAAGTCACCAACGTACCGTAg                                                        <  2:414313‑M2/89‑1 (MQ=255)
                                                                                   acataccgtcgataccgcgagtgatggtggacattaaaatccagatcagccagGCGCAAGCTCAGTCAGGAACAGGGCGATACCGAAACTCACCGGGACGGCGATCAGCAGCGCGATAAACGAAGTCACCAACGTACCGTAGa                                                       >  2:394493‑M2/54‑143 (MQ=255)
                                                                                      taccgtcgataccgcgagtgatggtggacattaaaatccagatcagccagGCGCAAGCTAAGTCAGGAACAGGGCGATACCGACACTCACCGGGACGGCGATCAGCAGCGCGATAAACGAAGTCCCCAACGTACCGCAAAtcc                                                    >  2:395304‑M2/51‑142 (MQ=255)
                                                                                          gtcgataccgcgagtgatggtggacattaaaatccagatcagccagGCGCAAGCTCAGTCAGGAACAGGGCGATACCGAAACTCACCGGGACGGCGATCAGCAGCGCGATAAACGAAGTCACCAACGTACCGTAGATAGGCAc                                                >  1:727138‑M2/47‑143 (MQ=255)
                                                                                                   gcgagtgatggtggacattaaaatccagatcagccagGCGCAAGCTCAGTCAGGAACAGGGCGATACCGAAACTCACCGGGACGGCGATCAGCAGCGCGATACACGAAGTCACCCACGTACCGTAGTTCGTCACCCGCGcccc                                       >  2:344007‑M2/38‑143 (MQ=255)
                                                                                                       gtgatggtggacattaaaatccagatcagccagGCGCAAGCTCAGTCAGGAACAGGGCGATACCGAAACTCACCGGGACGGCGATCAGCAGCGCGATAAACGAAGTCACCAACGTACCGTAGatcggcaccagatcggaagag                                   <  1:833301‑M2/110‑11 (MQ=255)
                                                                                                                gacattaaaatccagatcagccagGCGCAAGCTCAGTCAGGAACAGGGCGATACCGAAACTCCCCGGGACGGCGATCAGCAGCGCGATAAACGAAGTCACCAACGTACCGTGGATCGGCACGAGCGCCCCGTAGGTGTCGGTc                          >  2:178679‑M2/25‑143 (MQ=255)
                                                                                                                gacattaaaatccagatcagccagGCGCAAGCTCAGTCAGGAACAGGGCGATACCGAAACTCACCGGGACGGCGATCAGCAGCGCGATAAACGAAGGCACCAACGTACCGTAGATCGGCACCAGCGCCCCGGAGATATCGTTc                          >  2:178678‑M2/25‑143 (MQ=255)
                                                                                                                   attaaagtccagatcagccagGCGCAAGCTCAGTCAGGAACAGGGCGATACCGAAACTCACCGGGACGGCGATCCGCAACGCGATAAACGAAATAACACACGTACCGCAGATCAGCACcagcgcccagtaggactcgtagcga                       >  1:292383‑M2/22‑131 (MQ=255)
                                                                                                                   attaaaatccagatccgccagGCGCAAGCTCAGTCAGGAACAGGGCGATACCGAAACTCACCGGGACGGCGATCAGCAGCGCGATAAACGAAGTCACCAACGTACCGTAGATCGGCACCAGCGCCCCGTAGATATCGTTCGGt                       <  1:245648‑M2/122‑1 (MQ=255)
                                                                                                                   attaaaatccagatcagccagGCGCAAGCTCAGTCAGGAACAGGGCGATACCGAAACTCACCGGGACGGCGATCAGCAGCGCGATAAACGAAGTCACCAACGTACCGTAGATCGGCACCAGCGCCCCGTAGATATCGTTCGGt                       >  1:450416‑M2/22‑143 (MQ=255)
                                                                                                                    tcaaaatccagatcacccagGCGCAAGCTCAGTCAGGAACAGGGCGATACCGACACTCCCCGTGACGGCGATCAGCAGCGCGATCAACGTAGTCACCAACGTACCGTAGATCGGCACCAGCGCCCCGTAGATATCGTTCGGTg                      <  2:109‑M2/123‑1 (MQ=255)
                                                                                                                                    ccagGCGCAAGCTCAGTCAGGAACAGGGCGATACCGAAACTCACCGGGACGGCGATCAGCAGCGCGATAAACGAAGTCACCAACGTACCGTAGATCGGCACCAGCGCCCCGTAGATATCGTTCGGTGCATCCCACTCTTTGGt      <  2:727138‑M2/139‑1 (MQ=255)
                                                                                                                                    ccagGCGCAAGCTCAGTCAGGAACAGGGCGATACCGAAACTCACCGGGACGGCGATCAGCAGCGCGATAAACGAAGTCACCAACGTACCGTAGATCGGCAACAGCGCCCCGTAGATATCGTTCGGTGCATCCCActctttagc      >  1:475177‑M2/5‑140 (MQ=255)
                                                                                                                                     cagGCACAAGCTCAGTCAGGAACAGGGCGATACCGAAACTCACCGGGACGGCGATCAGCAGCGCGATAAACGAAGTCACCAACGTACCGTAGATCGGCACCAGCGCCCCGTAGATATCGTTCGGTGCATCCCACTCTTTGGTc     <  2:425111‑M2/140‑1 (MQ=255)
                                                                                                                                        gcgcAAGCTCAGTCAGGACCAGGGCGATACCGAAACTCACCGGGACGGCGATCAGCAGCGCGATAAACGAAGTCACCAACGTACCGTAGATCGGCGCCAGCGCCCCGTAGATATCGTTCGGTGCATCCCACTCTTTGGTCCAt  >  2:32404/1‑143 (MQ=255)
                                                                                                                                        |                                                                                                                                              
ATATTGCCGACCGGCTCCTGAAAGTAAACGGCGAACAGCGGCGCAAAGATAAACAGGCCCCACATGCCGTAAACGATACTTGGAATGGCTGCCAGCAGCTCAATGGCGATACCCAGCGGGCGTTTCAGCCAGCCAGGCGCAAGCTCAGTCAGGAACAGGGCGATACCGAAACTCACCGGGACGGCGATCAGCAGCGCGATAAACGAAGTCACCAACGTACCGTAGATCGGCACCAGCGCCCCGTAGATATCGTTCGGTGCATCCCACTCTTTGGTCCAT  >  NC_000913/3909949‑3910227

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.