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 57 [0] [0] 57 [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)54 (0.860) 48/268 0.3 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

GATATTGCCGACCGGCTCCTGAAAGTAAACGGCGAACAGCGGCGCAAAGATAAACAGGCCCCACATGCCGTAAACGATACTTGGAATGGCTGCCAGCAGCTCAATGGCGATACCCAGCGGGCGTTTCAGCCAGCCAGGCGCAAGCTCAGTCAGGAACAGGGCGATACCGAAACTCACCGGGACGGCGATCAGCAGCGCGATAAACGAAGTCACCAACGTACCGTAGATCGGCACCAGCGCCCCGTAGATATCGTTCGGTGCATCCCACTCTTTGGTCC  >  NC_000913/3909949‑3910226
                                                                                                                                         |                                                                                                                                            
ccaggccgcgacccgccggcgcgttcgccagaccgccgccgtccgtattgggcggcggcgtcatttcagtgaacagcgccagcgacataccggcgataccgcgagtgatggtggacattaaaatccagatcagccagGCGCaa                                                                                                                                         <  2:535385‑M2/6‑1 (MQ=255)
acagcccgctacccgccagagccttcgccagaccaccaccttccgtattgggcggcggcgtcatttcagtgaacagcgccagcgacattccgtcgataccgcgagtgatggtggacattaaaatccagatcagccagGCGCaa                                                                                                                                         <  1:619210‑M2/6‑1 (MQ=255)
  agcccgctacccgccagagcgttcgccagaccaccaccttccgtattgggcggcggcgtcatttcagtgaacagcgccagcgacctaccgtcgataccgcgagtgatggtggacattacaatccagatcagccagGCGCAAGc                                                                                                                                       >  1:72506‑M2/136‑143 (MQ=255)
    cccgctacccgccagagcgttcgccagaccacccccttccgtattgggcggcggcgtcatttcagtgaacagcgccagcgacataccgtcgataccgcgagtgatggtggacattaaaatccagatcagccagGCGCAAGCtc                                                                                                                                     <  2:965955‑M2/10‑1 (MQ=255)
    cccgctacccgccagagcgttcgccagaccaccaccttccgtactgggcggcggcgtcatttcagtgaacagcgccagcgacataccgtcgataccgcgagtgatggtggacattaaaatccagatcagccagGCGCAAGCtc                                                                                                                                     >  2:1042648‑M2/134‑143 (MQ=255)
         tacccgccagagcgttcgccataccaccacctgccgtattgggcggctgcgtcatttcagtgaacagcgccagcgacataccgtcgataccgcgagtgatggtggacattaaaatccagatcagccagGCGCAAGCtcagtca                                                                                                                                <  2:131317‑M2/15‑1 (MQ=255)
             cgccagagcgttcgccagaccaccaccttccgtattgggcggcggcgtcatttcagtgaacagcgccagcgacataccgtcgataccgcgagtgatggtggacattaaaatccagatcagccagGCGCAAGCTCAGTCAGGaa                                                                                                                            <  1:1045702‑M2/19‑1 (MQ=255)
             cgccagagcgttcgccagaccaccaccttccgtattgggcggcggcgtcatttcagtgaacagcgccagcgacataccgtcgataccgcgagtgatggtggacattaaaatccagatcagccagGCGCAAGCTCAGTCAGGaa                                                                                                                            >  2:913357‑M2/125‑143 (MQ=255)
              gcaagggcgttcgccagaccaccaccttccggattgggcggcggcgtcatttcagtgaacagcgccagcgagatgccgtcgataccgcgagtgatggtggacattaaaatccagatcagccagGCGCAAGCTCAGTCAGGAAc                                                                                                                           <  2:196400‑M2/20‑1 (MQ=255)
                     cgttcgccagaccaccaccttccgtattgggcggcggcgtcatttcagtgaacagcgccagcgacatacggtcgataccgcgagtgatggtggacattaagatccagatcagccagGCGCGAGCTCAGTCAGGAGCAGGGCGa                                                                                                                    >  1:474154‑M2/117‑143 (MQ=255)
                          gccagaccaccaccttccgtattgggcggcggcgtcatttcagtgaacagcgccagcgacataccgtcgataccgcgagtgatggtggacattaaaatccagatcagccagGCGCAAGCTCAGTCAGGAACAGGGCGATACCg                                                                                                               <  2:898663‑M2/32‑1 (MQ=255)
                           ccagaccaccaccttccgtattgggcggcggcgtcatttcagtgaacagcgccagctacatcccgtcgataccgcgagtgatggtggacattaaaatccagatcagccagGCGCAAGCTCAGTCAGGAACAGGGCGATACCGa                                                                                                              <  2:1027373‑M2/33‑1 (MQ=255)
                             agaccgccgccgcccggagtgggcggcgccgtcatttcagtgaacagcgccagcgacataccgtcgataccgcgagtgatggtggacattaaaatccagatcagccagGCGCAAGCTCAGTCAGGAACAGGGCGATACCGaaa                                                                                                            <  2:401286‑M2/35‑1 (MQ=255)
                              gaccaccaccttccgtattgggcggcggcgtcatttcagtgaacagcgccagcgacataccgtcgataccgcgagtgatggtggacattaaaacccagatcagccagGCGCAAGCTCAGTCAGGAACAGGGCGATACCGAAAc                                                                                                           >  2:412771‑M2/108‑143 (MQ=255)
                               accaccaccttccgtattgggcggcggcgtcatttcagtgaacagcgccagcgacataccgtcgataccgcgagtgatggtggacattaaaatccagatcagccagGCGCAAGCTCAGTCAGGAACAGGGCGCTACcgaaaat                                                                                                          >  1:59425‑M2/107‑141 (MQ=255)
                                 caccaccttccgtattgggcggcggcgtcatttcagtgaacagcgccagcgacataccgtcgataccgcgagtgatggtggacattaaaatccatatcagccagGAGCAAGCTCAGTCAGGAACAGGGCGATACCGAAACTCa                                                                                                        >  1:1052711‑M2/105‑143 (MQ=255)
                                  accaccttccgtattgggcggcggcgtcatttcagtgaacagcgccagcgacatacagtcgattccgcgagtgagggtggacattaaaatccagatcagccagACGCAAGCTCAGTCAGGAACAGGGCGATACCGAAACTCAc                                                                                                       >  1:1063936‑M2/104‑143 (MQ=255)
                                   ccaccttccgtattgggcggcggcgtcatttcagtgaacagcgccagcgacataccgtcgataccgcgagtgatggtggacattaaaatccagatcagccagGCGCAAGCTCAGTCAGGAACAGGGCGATACCGAAACTCAcc                                                                                                      >  1:953137‑M2/103‑143 (MQ=255)
                                    caccttccgtattgggcggcggcgtcatttcagtgaacagcgccagcgacataccgtcgataccgcgagtgatggtggacattaaaatccagatcagccagGCGCAAGCTCAGTCAGGAACAGGGCGATACCGAAACTCACCg                                                                                                     >  1:365760‑M2/102‑143 (MQ=255)
                                         tccgtattgggcggcggcgtcatttcagtgaacagcgacagcgacataccgtcgataccgcgagtgatggtggacattaaaatccagatcagccagGCGCAAGCTCAGTCATGAACCGCGCGATACCGACACTCACCGGGGCg                                                                                                >  2:254328‑M2/97‑143 (MQ=255)
                                             tattgggcggcggcgtcatttcagtgaacagcgccagcgacataccgtcgataccgcgagtgatggtggacattaaaatccagatcagccagGCGCAAGCTCAGTCAGGAACAGGGCGATACCGAAACTCACCGGGACGGCGa                                                                                            <  1:788295‑M2/51‑1 (MQ=255)
                                                  ggcggcggcgtcatttcagtgaacagcgccagcgacataccgtcgataccgcgcgtgatggtggacattaaaatccagatcagccagGCCCAAGCCCAGGCAGGAACAGGGCGATACCGAACCTCACCCGGACGGCGagcacg                                                                                       >  2:458317‑M2/88‑141 (MQ=255)
                                                       cggcgtcatttcagtgaacagcgccagcgacataccgtcgataccgcgagtgatggtggacattaaaatccagatcagccagGCGCAAGCTCAGTCAGGAACAGGGCGATACCGAAACTCACCGGGACAGCGGCCAGCAgcgc                                                                                  >  1:1113453‑M2/83‑143 (MQ=255)
                                                        ggcgtcatttcagtgatcagcgccagcgacataccgtcgataccgcgagtgatggtggacattaaaatccagatcagccagGCGCAAGCTCAGTCAGGAACAGGGCGATACCGAAACTCACCGGGACGGCGATCAGCAGcgcg                                                                                 <  1:1148792‑M2/62‑1 (MQ=255)
                                                        ggcgtcatttcagtgaacagcgccagcgacataccgtcgataccgcgagtgatggtggacattaaaatccagatcagccagGCGCAAGCTCAGTCAGGAACAGGGCGATACCGAAACTCACCGGGACGGCGATCAGCAGcgcg                                                                                 >  1:1014228‑M2/82‑143 (MQ=255)
                                                              accacagtgaacagcgccagcgacataccggcgataccgcgagtgatggtggacattaaaatccagatcagccagGCGCAAGCTCAGTCAGGAACAGGGCGATACCGAAACTCACCGGGACGGCGATCAGCAGCGCGATAAAc                                                                           <  1:728025‑M2/68‑1 (MQ=255)
                                                               tttcagtgaacagcgccagcgacataccgtcgataccgcgagtgatggtggacattaaaatccagatcagccagGCGCAAGCTACGTCACGATCAGGGCGATACCGAAACCCACCGGGCCGGCGATCAGCAGCGAGAGCAACg                                                                          >  2:894440‑M2/75‑143 (MQ=255)
                                                                 tgagtgaccagcgccagctgcataccgtcgataccgcgagtgatggtggagattaaaatccagatcagccagGCGCAAGGTCAGTCAGGAACAGGGCGATACCGAAGCTCACCGGGACGGCGATCAGCAGCGCGATAAACGaa                                                                        <  2:1096868‑M2/71‑1 (MQ=255)
                                                                     cgaacagcgccagcgacataccgtcgataccgcgagtgatggtggacattaaaatccagatcagccagGCGCAAGCTCAGTCAGGACCAGGGCCATCCCGGATCTCACCGGGACGGCGATCAGCAGCGCGATAAACGAAGTCa                                                                    <  1:612623‑M2/75‑1 (MQ=255)
                                                                          agcgccagcgacgtaccgtcgataccgcgagtgatggtggacattaaaagccagatcagccagGCGCAAGCTCAGTCAGGAACAGGGCGATACCGAAACTCACCGGGACGGCGATCAGCAGCGCGATAAACGAAGTCACCAAc                                                               <  1:679095‑M2/80‑1 (MQ=255)
                                                                                        accgtcgataccgcgagtgatggtggacattaaaatccagatcagccagGCGCAAGCTCAGTCAGGAACAGGGCGATACCGAAACTCACCGGGACGGCGATCAGCAGCGCGATAAACGCAGTCACGAACGTAccgtagcgcgg                                                 >  2:84436‑M2/50‑138 (MQ=255)
                                                                                        accgccgataccgcgagtgatggtggacattaaaatccagctcagccagGCGCAAGCTCAGTCAGGAACAGGGCGATACCGAAACTCACCGGGACGGCGATCAGCAGCGCGATAAACGAAGTCACCAACGTACCGTAGATCgg                                                 <  1:789890‑M2/94‑1 (MQ=255)
                                                                                        accgccgataccgcgagggatgggggacattaaaatccagagcagccagGCGCAAGCTCAGTCAGGAACAGGGCGATACCGAAACTCACCGGGACGGCGATCAGCAGCGCGATAAACGAAGTCACCAACGTACCGTAGATCgg                                                 <  1:992089‑M2/94‑1 (MQ=255)
                                                                                         ccgtcgataccgcgagtgatggtggacattaaaatccagatcagccagGCGCAAGCTCAGTCAGGAACAGGGCGATACCGAAACTCACCGGGACGGAGAACAGCAGCGCGACAAACGAAGTCACCAACGTACCGTAGACCGGc                                                >  2:468181‑M2/49‑143 (MQ=255)
                                                                                           gtcgataccgcgagtgatggtggacattaaaatccagatcagccagGCGCAAGCTCAGTCAGGAACAGGGCGATACCGAAACTCACCGGGACGGCGATCAGCAGCGCGATAAACGAAGTCACCAACGTACCGTAGATCGGCAc                                              <  1:1042648‑M2/97‑1 (MQ=255)
                                                                                             cgataccgcgtgttatgggggacattaaaatccagatcagccagGCGCAAGCTCAGTCAGGAACAGGGCGATACCGAAACTCACCGGGACGGCGATCAGCAGCGCGATAAACGAAGTCACCAACGTACCGTAGATCGGCACCa                                            <  2:1052711‑M2/99‑1 (MQ=255)
                                                                                                  gcgcgagtgacggtggcccttaaacaccagatcagccagGCGCAAGCTCAGTCAGGAACAGGGCGATACCGAAACTCACCGGGACGGCGATCAGCAGCGCGATAAACGAAGTCACCAACGTACCGTAGATCGGCACCAGCGcc                                       <  1:488571‑M2/104‑1 (MQ=255)
                                                                                                  ccgcgagtgatggtggacattaaaatccaaatcagccagGCGCAAGCTCAGTCAGGAACAGGGCGATACCGAAACTAACCGGGACGACGATCAGCAGCGCGAAAAACGCAGTCACCAACGTACGGTAGACCGGCACCAACGcc                                       >  2:212873‑M2/40‑143 (MQ=255)
                                                                                                    gcgagtgatggtggacattaaaatccagatcagccagGCGCAAGCTCAGTCAGGAACAGGGCGATACCGAAACTCACCGGGACGGCGATCAGCAGCGCGATAAACGAAGTCACCAACGTACCGTAGATCGGCACCAGCgcccg                                     >  1:158711‑M2/38‑142 (MQ=255)
                                                                                                       agtgatggtggacattaaaatccagatcagccagGCGCAAGCTCAGTCAGGAACAGGGCGATACCGAAACTCACCGGGACGGCGATTAGCAGCGCGATAAACGAAGTCACCAACGTACCGTAGATCGGCACCAGCGCCCCGTa                                  <  1:510507‑M2/109‑1 (MQ=255)
                                                                                                       agtgatggtggacattaaaacccagatcagccagGCGCAAGCTCAGTCAGGAACAGGGCGATACCGAAACTCACCGGGACGGCGATCAGCAGCGCGATAAACGAAGTCACCAACGTACCGTAGATCGGCACCAGCGCCCCGTa                                  <  1:412771‑M2/109‑1 (MQ=255)
                                                                                                         tgatggtggacattaaaatccagatcagccagGCGCAAGCTCAGTCAGGAACAGGGCGATACCGAAACTCACCGGGACGGCGATCAGCAGCGCGATAAACGAAGTCACCAACGTACCGTAGATCGGCACCAGCGCCCCGTAGa                                >  1:956421‑M2/33‑143 (MQ=255)
                                                                                                              gtggacattaaaatccagatcagccagGCGCAAGCTCAGTCAGGAACAGGGCGATACCGAAACTCACCGGGACGGCGATCAGCAGCGCGATAAACGAAGTCACCAACGTGCCGTAGATCGGCACCAGCGCCCCGTAGATATCg                           >  1:148184‑M2/28‑143 (MQ=255)
                                                                                                              gtggacattaaaatccagatcagccagGCGCAAGCTCAGTCAGGAACAGGGCGATACCGAAACTCACCGGGACGGCGATCAGCAGCGCGATAAACGAAGTCACCAACGTACCGTAGATCGGCACCAGCGCCCCGTAGATATCg                           >  2:939356‑M2/28‑143 (MQ=255)
                                                                                                              gtggacattaaaatccagatcagccagGCGCAAGCTCAGTCAGGAACAGGGCGATACCGAAACTCACCGGGACGGCGATCAGCAGCGCGATAAACGAAGGCACCAACGGACCGTAGATCGGCACCAGCGCCCCGGAGATGTCg                           >  2:784617‑M2/28‑143 (MQ=255)
                                                                                                                ggacattaaaatccagatcagccagGCGCAAGCTCAGTCAGGAACAGGGCGATACCGAAACTCACCGGGACGGCGATCGGCAGCGCGATAAACGAAGTCACCAACGTACCGCAGATAGGCACCAGAGCCCCGTAGACAGCgta                         >  2:92620‑M2/26‑142 (MQ=255)
                                                                                                                 gacattaaaatccagatcaaccagGCGCAAGCTCAGTCAGGAACAGGGCGATACCGAAACTCACCGGGACGGCGATCAGCAGCGCGATAAACGAAGTCACCAACGTACCGTAGATCGGCACCAGCGCCCCGTAGATATCGTTc                        >  1:250987‑M2/25‑143 (MQ=255)
                                                                                                                    attaaaatccagatcagccagGCGCAAGCTCAGTCAGGAACAGGGCGATACCGAAACTCACCGGGACGGCGCTCAGCAGCCCGATAAACGAAGTCACCAACGTACCGTAGCGAGGCACCCGCGCCCCGTAGATATCGTTCGGt                     >  1:599835‑M2/22‑143 (MQ=255)
                                                                                                                    attaaaatccagatcagccagGCGCAAGCTCAGTCAGGAACAGGGCGATACCGAAACTCACCGGGACGGCGATCAGCAGCGCGATAAACGAAGTCACCAACGTACCGTAGATCGGCACCAGCGCCCCGTAGATATCGTTCGGt                     >  1:398369‑M2/22‑143 (MQ=255)
                                                                                                                          gtccagatcagcaagGCGCAAGCTCAGGCAGGAACAGGGCGATACCGAAACGCACCGGGACGGCGATCAGCAGCGCGATAAACGAAGTCACCAACGTACCGTAGATCGGCACCAGCGCCCCGTAGATATCGTTCGGTGCATcc               <  1:982752‑M2/128‑1 (MQ=255)
                                                                                                                            ccagatcagccagGCGCAAGCTCAGTCAGGAACAGGGCGATACCGAAACTCCCCGGGCCAGCAAATAGCAGCGCGATAACCGAAGTCACCAACGTACCGTAGATCGGCGCCAGCGCCCCGTAAATAACGCtcggagccagccc             >  2:780566‑M2/14‑137 (MQ=255)
                                                                                                                             cagatcagccagGCGCAAGCTCAGTCAGGAACAGGGCGATACCGAAACTCACCGGGACGGCGATCAGCAGCGCGATAAACGAAGTCACCAACGTACCGTAGATCGGCACCAGCGCCCCGTAGATATCGTTCGGTGCATCCCAc            <  1:772486‑M2/131‑1 (MQ=255)
                                                                                                                               gatcagccagGCGCAAGCTCAGTCAGGAACAGGGCGATACCGAAACTCACCGGGACGGCGATCAGCAGCGCGATAAACGAAGTCACCAACGTACCGTAGATCGGCACCAGCGCCCCGTAGATATCGTTCGGTGCATCCCActc          >  1:953673‑M2/11‑143 (MQ=255)
                                                                                                                                 tcagccagGCGCAAGCTCAGTCAGGAACAGGGCGATACCGAAACTCACCGGGACGGCGATCAGCAGCGCGATAAACGAAGTCACCAACGTACCGTAGATCGGCACCAGCGCCCCGTAGATATCGTTCGGTGCATCCCACTCtt        <  1:458454‑M2/135‑1 (MQ=255)
                                                                                                                                   agccagGCGCAAGCTCAGTCCGGGACAGGGCGATACCGAAGCTCACCGGGACGGCGATCAGCAGCGCGATAAACGAAGTCACCAACGTACCGTAGATCGGCACCAGCGCCCCGTAGATATCGTTCGGTGCATCCCACTCTTTg      <  2:1183175‑M2/137‑1 (MQ=255)
                                                                                                                                      cagGCGCAAGCTCAGTCAGGAACAGGGCGATACCGAAACTCACCGGGACGGCCATCAGCAGCGCGATAAACGAAGTCACCAACGTACCGTAGATCGGCACCAGCGCCCCGTAGAGATCGTTCGGTGCATCGCACTCTTTGGTc   >  1:251356‑M2/4‑143 (MQ=255)
                                                                                                                                       agGAGCAAGCTCAGTCAGGAACAGGGCGATACCGAAACTCACCGGGACGGCGATTAGCAGCGCGATAAACGAAGTCACCAACGTACCGTAGATCGGCACCAGAGCTCCGTAGATATCGTTCGGTGCATCCGACTCTTTggtac  >  1:619033‑M2/3‑141 (MQ=255)
                                                                                                                                         |                                                                                                                                            
GATATTGCCGACCGGCTCCTGAAAGTAAACGGCGAACAGCGGCGCAAAGATAAACAGGCCCCACATGCCGTAAACGATACTTGGAATGGCTGCCAGCAGCTCAATGGCGATACCCAGCGGGCGTTTCAGCCAGCCAGGCGCAAGCTCAGTCAGGAACAGGGCGATACCGAAACTCACCGGGACGGCGATCAGCAGCGCGATAAACGAAGTCACCAACGTACCGTAGATCGGCACCAGCGCCCCGTAGATATCGTTCGGTGCATCCCACTCTTTGGTCC  >  NC_000913/3909949‑3910226

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: