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

Missing coverage evidence...
   seq id start end size ←reads reads→ gene description
* * ÷ NC_000913 3909312 3910085 774 44 [0] [0] 44 [pstA]–[pstC] [pstA],[pstC]

New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913 = 39093110 (0.000)40 (0.880) 36/268 0.2 100% coding (129/891 nt) pstA phosphate ABC transporter permease
?NC_000913 3910086 = 0 (0.000)coding (313/960 nt) pstC phosphate ABC transporter permease

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

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.