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

AAATTAACAGCCCGCTACCCGCCAGAGCGTTCGCCAGACCACCACCTTCCGTATTGGGCGGCGGCGTCATTTCAGTGAACAGCGCCAGCGACATACCGTCGATACCGCGAGTGATGGTGGACATTAAAATCCAGATCAGCCAGAACAGGCCGAAGGCCATCGTCGCCATCGAGAGCGTCAGCGCAATACGGTTTTTGAGGCGGCGACGCGCCTGCATTTTGCGGCGAGATTCAGCCAGCGCCGCAGTGGTTTGCATTTCAACCATAGCCATTAGCGTGCCCCC  >  NC_000913/3909168‑3909450
                                                                                                                                              |                                                                                                                                            
acattAACACCCCGCTACCCGCCAGAGCGTTCGCCAGACCACCACCTTCCGTATTGGGCGGCGGCGTCATTTCAGTGAACAGCGCCAGCGACATACCGTCGATACCGCGAGTGATGGTGGACATTAAAATCCAGATCAGCCAG                                                                                                                                              <  1:248133‑M1/141‑1 (MQ=255)
      ccctcccGCTTCCTTCCAGAGCGTCCGCCCGACCACCTCCTCCCGTAATAGGCGGCTGCGTCATTTCACTGCACAGTGCCAGCGACATACCGTCGATACCCCGAGTGATGGTGCACATTAAAATCCAGTTCAGCCAggcgcaa                                                                                                                                        <  1:75107/139‑7 (MQ=255)
       cAGCCCGCTACCCGCCAGAGCGTTCGCCAGACCACCACCGTCCGTATTGGGCGGCGGCGTCATTTCAGTGAACAGCGCCAGCGACATACCGTCGATACCGCGAGTGATGGTGGACATTAAAAGCGAGCTCAGCCAGgcgcaag                                                                                                                                       >  1:872884‑M1/1‑136 (MQ=255)
         gCCCTCTGCGCGCCAGGGCGTTCGCCAGACCACCACCTTCCGTATTGAGCGGCGGCTTCATTTCAGTGAGGAGCGCCAGCGACGTGTCGTCGACACCGCGAGTGATGGTGGACATTAAAATCCAGATCAGCCAGgcgcaagct                                                                                                                                     <  1:810586‑M1/143‑10 (MQ=255)
          cccGCTACCCGCCAGAGCGTTCGCCAGACCACCACCTTCCGTATTGGGCGGCGGCGTCATTTCAGTGAACAGCGCCAGCGACATACCGTCGATACCGCGAGTGATGGTGGACATTAAAATCCAGATCAGCCACgcgcaagctc                                                                                                                                    >  2:693273‑M1/1‑133 (MQ=255)
           ccGCTACCCGCCAGAGCGTTCGCCAGACCACCACCTTCCGTCTTGGGCGGCGGCGTCACTTCAGTGAACAGCGCCAGCGACATACCGTCGATACCGCGAGTGAAGGTGGACATTAAAATCCAGATCAGCCAGgcgcaagctca                                                                                                                                   >  1:195743‑M1/1‑132 (MQ=255)
           ccGCTACCCGCCAGAGCGTTCGCCAGACCACCACCTTCCGTATTGGGCGGCGGCGTCATTTCAGTGAACAGCGCCAGCGACATACCGTCGATACCGCGAGTGATGGTGGACATTAAAATCCAGATCAGCCCGgcgcaagctcg                                                                                                                                   >  1:194787‑M1/1‑132 (MQ=255)
                  ccGCCATAGCGTTCGCCGGACCACCCCCTTCCGTATTGGGCGGCGGCGTCAGTTCAGTGAACAGCGCCAGCGACATACCGTCGATACCGCGAGTGATGGTGGACATTAAAATCCAGATCAGCCAGgcgcaagctcagtcagga                                                                                                                            <  2:254219‑M1/143‑19 (MQ=255)
                    gCCAGAGCGTTCGCCAGACCACCACCTTCCGTATTGGGCGGCGGCGTCATTTCAGTGAACAGCGCCAGCGACATACCGTCGATACCGCGAGTGATGGTGGACATTCACATCCAGATCAGCCGGgcgcaggcccggtcagtaac                                                                                                                          >  2:780978‑M1/1‑123 (MQ=255)
                       agagCGTTCGCCAGACCACCACCTTCCGTATTGGGCGGCGGCGTCATTTCAGTGAACAGCGCCAGCGACATACCGTCGATACCGCGAGTGATGGTGGACATTAAAATCCAGATCAGCCAGgcgcaagctccgtcaggagcagg                                                                                                                       >  2:577802‑M1/1‑120 (MQ=255)
                       agagCGTTCGCCAGACCACCACCTTCCGTATTGGGCGGCGGCGTCATTTCAGTGAACAGCGCCAGCGACATACCGTCGATACCGCGAGTGATGGTGGACATTAAAATACAGATCAGCCAGgcgcaagctcaggcaggaacaag                                                                                                                       >  1:337308‑M1/1‑120 (MQ=255)
                                gCCAGACCACCACCTGCGGTAGTGGGCGGCGGCGTCATTGCAGTGAACAGCGCCAGCGACATACCGTCGATACCGCGAGTGATGGTGGACATTAAAATCCAGATCAGCCAGgcgcaagctcagtcaggaacagggcgataccg                                                                                                              <  2:312515‑M1/143‑33 (MQ=255)
                                     accaccaccTTCCGTATTGGGCGGCGGCGTCCTTACAGTGAACAGCGCCCGCGACATACCGTCGATACCGCGAGTGAGGGTGGACAGTAACATCCAGATCAGCCCGccgcaagctcagtcaggcaaacggcaacaccgaacga                                                                                                         >  1:42827‑M1/1‑106 (MQ=255)
                                     accaccaccTTCCGTATTGGGCGGCGGCGTCATTTCAGTGAACAGCGCCAGCGACATACCGTCGATACCGCGAGTGATGGTGGACATTAAAATCCAGATCAGCCAGgcgcaagctcagtcaggaacagggcgataccgaaact                                                                                                         <  1:337156‑M1/143‑38 (MQ=255)
                                       caccaccTTCCGTATTGGGCGGCGGCGTCATTTCAGTGAACAGCGCCAGCGACATACCGTCGATACCGCGAGTGATGGTGGACATTAAAATCCAGATCAGCCAGgcgcaagctcagtcaggaacagggcgataccgaaactca                                                                                                       <  2:577333‑M1/143‑40 (MQ=255)
                                       caccaccTTCCGTATTGGGCGGCGGCGTCATTTCAGTGAACAGCGCCAGCGACATACCGTCGATACCGCGAGTGATGGTGGACATTAAAATCCAGATCAGCCAGgcgcaagctcagtcaggaacagggcgataccgaaactca                                                                                                       >  2:32088‑M1/1‑104 (MQ=255)
                                        accaccTTCCGTATCGGGCGGCGGCGTCATTTCAGTGAACAGCGCCAGCGACATACCGTCGATACCGCGAGTGATGGTGGACATTAAAATCCAGATCAGCCAGgcgcaagctcagtcaggaacagggcgataccgaaactcac                                                                                                      <  1:359371‑M1/143‑41 (MQ=255)
                                        accaccGTCCGTATTGGGCGGCGGCGTCATTTCAGTGAACAGCGCCAGCGACATACCGTCGATACCGCGAGTGATGGTGGACATTAAAATCCAGATCAGCCAGgcgcaagctcagtcaggaacagggcgataccgaaactcac                                                                                                      <  1:656171‑M1/143‑41 (MQ=255)
                                                    cttgGGCGGCGGCGTCATTTGAGTGAACAGCGCCAGCGACATACCGTCGATACCGCGAGTGATGGTGGACATTAAAATCCAGATCAGCCAGgcgcaagctcagtcaggaacagggcgataccgaaactcaccgggacggcgat                                                                                          <  1:150179‑M1/142‑53 (MQ=255)
                                                           ggcggcgTCATTTCAGTGAACAGCGCCAGCGACATACCGTCGATACCGCGAGTGATGGTGGACATTACAATCCAGATCACCCAGgcgcaagctcagccagcaccagcgccataccgaacctccccgggcccgcgatcacccgc                                                                                   >  1:204804‑M1/1‑84 (MQ=255)
                                                             cggcgTCATTTCAGTGAACAGCGCCAGCGACATACCGTCGATACCGCGAGTGATGGTGGACATTAAAATCCAGATCAGCCAGgcgcaagctcagtcaggaacagggcgataccgaaacacaccgggacggcgatcagcagcgc                                                                                 >  1:514018‑M1/1‑82 (MQ=255)
                                                             cggcgTCATTTCAGTGAACAGCGCCAGCGACATACCGTCGATACCGCGAGTGATGGTGGACATTAAAATCCAGAGCAGCCAGgcgcaagctcagtcaggaacagggcgataccgaaactcaccgggacggcgatcagcagcgc                                                                                 <  2:416319‑M1/143‑62 (MQ=255)
                                                              ggcgTCATTTCAGTGAACAGCGCCAGCGACATACCGTCGATACCGCGAGTGATGGGGGACATTAAAATCCAGATCAGCCAGgcgcaagctcagtcaggaacagggcgataccgaaactcaccgggacggcgatcagcagcgcg                                                                                <  2:860405‑M1/143‑63 (MQ=255)
                                                               gcgTCACTTCAGTGAACAGCGCCAGCGACATACCGTCGATACCGCGAGTGATGGTGGACATTAAAATCCAGATCAGCCAGgcgcaagctcagtcaggaacagggcgataccgaaactcaccgggacggcgatcagcagcgcga                                                                               <  1:218328‑M1/143‑64 (MQ=255)
                                                                  tCATTTCAGTGAACAGCGCCAGCCACACACCTTTGCTACCGCGAGTGATGGTGGACATTAAAGTCCAGATCAGCCAGgcgcaagctcagtcaggaacagggcgataccgaaactcaccgggacggcgatcagcagcgcgataa                                                                            <  2:594210‑M1/143‑67 (MQ=255)
                                                                   gacttcagTGCACAGCGCCAGCGACATACCGTCGATACCGCGAGTCATGGTGGACATGAAAATCCAGATCAGCCAGgcgcaagctcagtcaggaacagggcgataccgaaactcaccgggacggcgatcagcagcgcgataaa                                                                           <  1:688166‑M1/140‑68 (MQ=255)
                                                                       tCAGTGAACAGCGCCAGCGACATACCGTCGATACCGCGAGTGATGGTGGACATTAAAATCCAGATCAGCCAGgcgcaagctcagtcaggaacagggcgataccgaaactcaccgggacggcgatcagcagcgcgataaacgaa                                                                       <  2:2266‑M1/143‑72 (MQ=255)
                                                                           tGAACAGCGCCAGCGACATACCGTCGATACCGCGAGTGATGGTGGACATTAAAATCCAGATCAGCCAGgcgcaagctcagtcaggaacagggcgataccgaaactcaccgggacggcgatcagcagcgcgataaacgaagtca                                                                   >  2:286708‑M1/1‑68 (MQ=255)
                                                                           tGAACAGCGCCAGCGACAGCCCGTCGATACCGCTAGTGATGGTGTACATTAAAATCCAGATCAGCCAGgcgcaagctcagtcaggaacagggcgataccgaaactcaccgggacggcgatcagcagcgcgataaacgaagtca                                                                   <  1:47542‑M1/143‑76 (MQ=255)
                                                                                aGCGCCAGCGACATCCCGGCGATACCGCGAGTGATGGTGGACATTAAAATCCAGATCAGCCAGgcgcaagctcagtcaggaacagggcgataccgaaactcaccgggacggcgatcagcagcgcgataaacgaagtcaccaac                                                              <  2:186780‑M1/143‑81 (MQ=255)
                                                                                 gcgcCAGCGACATCCCGGCGATACCGGGAGTGATGGTGGACATTAAAATCCAGATCAGCCGGgcgcgaggtcagtcaggaacagggcgataccgaaactcaccgggacggcgatcagcagcgcgataaacgaagtcaccaacg                                                             <  2:405481‑M1/143‑82 (MQ=255)
                                                                                    ccatccacttaccttcgatACCGCGAGTGATGGTGGACATTAAAATCCAGATCAGCCAGgcgcaagctcagtcaggaacagggcgataccgaaactcaccgggacggcgatcagcagcgcgataaacgaagtcaccaacgtac                                                          <  1:566677‑M1/134‑85 (MQ=255)
                                                                                       gCGACATACCGTCGATACCGCGAGTGATGGTGGACATTAAAATCCAGATCAGCCAGgcgcaagctcagtcaggaacagggcgataccgaaactcaccgggacggcgatcagcagcgcgataaacgaagtcaccaacgtaccgt                                                       >  2:472198‑M1/1‑56 (MQ=255)
                                                                                         gACATACCGTCGATACCGCGAGTGATGGTGGACATTAAAATCCAGATCAGCCAGgcgcaagctcagtcaggaacagggcgataccgaaactcaccgggacggcgatcagcagcgcgataaacgaagtcaccaacgtaccgtag                                                     <  2:414313‑M1/143‑90 (MQ=255)
                                                                                          aCATACCGTCGATACCGCGAGTGATGGTGGACATTAAAATCCAGATCAGCCAGgcgcaagctcagtcaggaacagggcgataccgaaactcaccgggacggcgatcagcagcgcgataaacgaagtcaccaacgtaccgtaga                                                    >  2:394493‑M1/1‑53 (MQ=255)
                                                                                             tACCGTCGATACCGCGAGTGATGGTGGACATTAAAATCCAGATCAGCCAGgcgcaagctaagtcaggaacagggcgataccgacactcaccgggacggcgatcagcagcgcgataaacgaagtccccaacgtaccgcaaatcc                                                 >  2:395304‑M1/1‑50 (MQ=255)
                                                                                                 gTCGATACCGCGAGTGATGGTGGACATTAAAATCCAGATCAGCCAGgcgcaagctcagtcaggaacagggcgataccgaaactcaccgggacggcgatcagcagcgcgataaacgaagtcaccaacgtaccgtagataggcac                                             >  1:727138‑M1/1‑46 (MQ=255)
                                                                                                          gcgAGTGATGGTGGACATTAAAATCCAGATCAGCCAGgcgcaagctcagtcaggaacagggcgataccgaaactcaccgggacggcgatcagcagcgcgatacacgaagtcacccacgtaccgtagttcgtcacccgcgcccc                                    >  2:344007‑M1/1‑37 (MQ=255)
                                                                                                              gTGATGGTGGACATTAAAATCCAGATCAGCCAGgcgcaagctcagtcaggaacagggcgataccgaaactcaccgggacggcgatcagcagcgcgataaacgaagtcaccaacgtaccgtagatcggcaccagatcggaagag                                <  1:833301‑M1/143‑111 (MQ=255)
                                                                                                                       gACATTAAAATCCAGATCAGCCAGgcgcaagctcagtcaggaacagggcgataccgaaactccccgggacggcgatcagcagcgcgataaacgaagtcaccaacgtaccgtggatcggcacgagcgccccgtaggtgtcggtc                       >  2:178679‑M1/1‑24 (MQ=255)
                                                                                                                       gACATTAAAATCCAGATCAGCCAGgcgcaagctcagtcaggaacagggcgataccgaaactcaccgggacggcgatcagcagcgcgataaacgaaggcaccaacgtaccgtagatcggcaccagcgccccggagatatcgttc                       >  2:178678‑M1/1‑24 (MQ=255)
                                                                                                                          aTTAAAGTCCAGATCAGCCAGgcgcaagctcagtcaggaacagggcgataccgaaactcaccgggacggcgatccgcaacgcgataaacgaaataacacacgtaccgcagatcagcaccagcgcccagtaggactcgtagcga                    >  1:292383‑M1/1‑21 (MQ=255)
                                                                                                                          aTTAAAATCCAGATCCGCCAGgcgcaagctcagtcaggaacagggcgataccgaaactcaccgggacggcgatcagcagcgcgataaacgaagtcaccaacgtaccgtagatcggcaccagcgccccgtagatatcgttcggt                    <  1:245648‑M1/143‑123 (MQ=255)
                                                                                                                          aTTAAAATCCAGATCAGCCAGgcgcaagctcagtcaggaacagggcgataccgaaactcaccgggacggcgatcagcagcgcgataaacgaagtcaccaacgtaccgtagatcggcaccagcgccccgtagatatcgttcggt                    >  1:450416‑M1/1‑21 (MQ=255)
                                                                                                                           tcaaaatcCAGATCACCCAGgcgcaagctcagtcaggaacagggcgataccgacactccccgtgacggcgatcagcagcgcgatcaacgtagtcaccaacgtaccgtagatcggcaccagcgccccgtagatatcgttcggtg                   <  2:109‑M1/141‑124 (MQ=255)
                                                                                                                                           ccAGgcgcaagctcagtcaggaacagggcgataccgaaactcaccgggacggcgatcagcagcgcgataaacgaagtcaccaacgtaccgtagatcggcaccagcgccccgtagatatcgttcggtgcatcccactctttggt   <  2:727138‑M1/143‑140 (MQ=255)
                                                                                                                                           ccAGgcgcaagctcagtcaggaacagggcgataccgaaactcaccgggacggcgatcagcagcgcgataaacgaagtcaccaacgtaccgtagatcggcaacagcgccccgtagatatcgttcggtgcatcccactctttagc   >  1:475177‑M1/1‑4 (MQ=255)
                                                                                                                                            cAGgcacaagctcagtcaggaacagggcgataccgaaactcaccgggacggcgatcagcagcgcgataaacgaagtcaccaacgtaccgtagatcggcaccagcgccccgtagatatcgttcggtgcatcccactctttggtc  <  2:425111‑M1/143‑141 (MQ=255)
                                                                                                                                              |                                                                                                                                            
AAATTAACAGCCCGCTACCCGCCAGAGCGTTCGCCAGACCACCACCTTCCGTATTGGGCGGCGGCGTCATTTCAGTGAACAGCGCCAGCGACATACCGTCGATACCGCGAGTGATGGTGGACATTAAAATCCAGATCAGCCAGAACAGGCCGAAGGCCATCGTCGCCATCGAGAGCGTCAGCGCAATACGGTTTTTGAGGCGGCGACGCGCCTGCATTTTGCGGCGAGATTCAGCCAGCGCCGCAGTGGTTTGCATTTCAACCATAGCCATTAGCGTGCCCCC  >  NC_000913/3909168‑3909450

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.