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

ACAGCCCGCTACCCGCCAGAGCGTTCGCCAGACCACCACCTTCCGTATTGGGCGGCGGCGTCATTTCAGTGAACAGCGCCAGCGACATACCGTCGATACCGCGAGTGATGGTGGACATTAAAATCCAGATCAGCCAGAACAGGCCGAAGGCCATCGTCGCCATCGAGAGCGTCAGCGCAATACGGTTTTTGAGGCGGCGACGCGCCTGCATTTTGCGGCGAGATTCAGCCAGCGCCGCAGTGGTTTGCATTTCAACCATAGCCATTAGCGTGCCCCCT  >  NC_000913/3909175‑3909452
                                                                                                                                        |                                                                                                                                             
ccaGGCCGCGACCCGCCGGCGCGTTCGCCAGACCGCCGCCGTCCGTATTGGGCGGCGGCGTCATTTCAGTGAACAGCGCCAGCGACATACCGGCGATACCGCGAGTGATGGTGGACATTAAAATCCAGATCAGCCAGgcgcaa                                                                                                                                         <  2:535385‑M1/142‑7 (MQ=255)
aCAGCCCGCTACCCGCCAGAGCCTTCGCCAGACCACCACCTTCCGTATTGGGCGGCGGCGTCATTTCAGTGAACAGCGCCAGCGACATTCCGTCGATACCGCGAGTGATGGTGGACATTAAAATCCAGATCAGCCAGgcgcaa                                                                                                                                         <  1:619210‑M1/143‑7 (MQ=255)
  aGCCCGCTACCCGCCAGAGCGTTCGCCAGACCACCACCTTCCGTATTGGGCGGCGGCGTCATTTCAGTGAACAGCGCCAGCGACCTACCGTCGATACCGCGAGTGATGGTGGACATTACAATCCAGATCAGCCAGgcgcaagc                                                                                                                                       >  1:72506‑M1/1‑135 (MQ=255)
    cccGCTACCCGCCAGAGCGTTCGCCAGACCACCCCCTTCCGTATTGGGCGGCGGCGTCATTTCAGTGAACAGCGCCAGCGACATACCGTCGATACCGCGAGTGATGGTGGACATTAAAATCCAGATCAGCCAGgcgcaagctc                                                                                                                                     <  2:965955‑M1/143‑11 (MQ=255)
    cccGCTACCCGCCAGAGCGTTCGCCAGACCACCACCTTCCGTACTGGGCGGCGGCGTCATTTCAGTGAACAGCGCCAGCGACATACCGTCGATACCGCGAGTGATGGTGGACATTAAAATCCAGATCAGCCAGgcgcaagctc                                                                                                                                     >  2:1042648‑M1/1‑133 (MQ=255)
         tACCCGCCAGAGCGTTCGCCATACCACCACCTGCCGTATTGGGCGGCTGCGTCATTTCAGTGAACAGCGCCAGCGACATACCGTCGATACCGCGAGTGATGGTGGACATTAAAATCCAGATCAGCCAGgcgcaagctcagtca                                                                                                                                <  2:131317‑M1/143‑16 (MQ=255)
             cGCCAGAGCGTTCGCCAGACCACCACCTTCCGTATTGGGCGGCGGCGTCATTTCAGTGAACAGCGCCAGCGACATACCGTCGATACCGCGAGTGATGGTGGACATTAAAATCCAGATCAGCCAGgcgcaagctcagtcaggaa                                                                                                                            <  1:1045702‑M1/143‑20 (MQ=255)
             cGCCAGAGCGTTCGCCAGACCACCACCTTCCGTATTGGGCGGCGGCGTCATTTCAGTGAACAGCGCCAGCGACATACCGTCGATACCGCGAGTGATGGTGGACATTAAAATCCAGATCAGCCAGgcgcaagctcagtcaggaa                                                                                                                            >  2:913357‑M1/1‑124 (MQ=255)
              gcaagggcgttcgCCAGACCACCACCTTCCGGATTGGGCGGCGGCGTCATTTCAGTGAACAGCGCCAGCGAGATGCCGTCGATACCGCGAGTGATGGTGGACATTAAAATCCAGATCAGCCAGgcgcaagctcagtcaggaac                                                                                                                           <  2:196400‑M1/137‑21 (MQ=255)
                     cGTTCGCCAGACCACCACCTTCCGTATTGGGCGGCGGCGTCATTTCAGTGAACAGCGCCAGCGACATACGGTCGATACCGCGAGTGATGGTGGACATTAAGATCCAGATCAGCCAGgcgcgagctcagtcaggagcagggcga                                                                                                                    >  1:474154‑M1/1‑116 (MQ=255)
                          gCCAGACCACCACCTTCCGTATTGGGCGGCGGCGTCATTTCAGTGAACAGCGCCAGCGACATACCGTCGATACCGCGAGTGATGGTGGACATTAAAATCCAGATCAGCCAGgcgcaagctcagtcaggaacagggcgataccg                                                                                                               <  2:898663‑M1/143‑33 (MQ=255)
                           ccAGACCACCACCTTCCGTATTGGGCGGCGGCGTCATTTCAGTGAACAGCGCCAGCTACATCCCGTCGATACCGCGAGTGATGGTGGACATTAAAATCCAGATCAGCCAGgcgcaagctcagtcaggaacagggcgataccga                                                                                                              <  2:1027373‑M1/143‑34 (MQ=255)
                             aGACCGCCGCCGCCCGGAGTGGGCGGCGCCGTCATTTCAGTGAACAGCGCCAGCGACATACCGTCGATACCGCGAGTGATGGTGGACATTAAAATCCAGATCAGCCAGgcgcaagctcagtcaggaacagggcgataccgaaa                                                                                                            <  2:401286‑M1/143‑36 (MQ=255)
                              gACCACCACCTTCCGTATTGGGCGGCGGCGTCATTTCAGTGAACAGCGCCAGCGACATACCGTCGATACCGCGAGTGATGGTGGACATTAAAACCCAGATCAGCCAGgcgcaagctcagtcaggaacagggcgataccgaaac                                                                                                           >  2:412771‑M1/1‑107 (MQ=255)
                               accaccaccTTCCGTATTGGGCGGCGGCGTCATTTCAGTGAACAGCGCCAGCGACATACCGTCGATACCGCGAGTGATGGTGGACATTAAAATCCAGATCAGCCAGgcgcaagctcagtcaggaacagggcgctaccgaaaat                                                                                                          >  1:59425‑M1/1‑106 (MQ=255)
                                 caccaccTTCCGTATTGGGCGGCGGCGTCATTTCAGTGAACAGCGCCAGCGACATACCGTCGATACCGCGAGTGATGGTGGACATTAAAATCCATATCAGCCAGgagcaagctcagtcaggaacagggcgataccgaaactca                                                                                                        >  1:1052711‑M1/1‑104 (MQ=255)
                                  accaccTTCCGTATTGGGCGGCGGCGTCATTTCAGTGAACAGCGCCAGCGACATACAGTCGATTCCGCGAGTGAGGGTGGACATTAAAATCCAGATCAGCCAGacgcaagctcagtcaggaacagggcgataccgaaactcac                                                                                                       >  1:1063936‑M1/1‑103 (MQ=255)
                                   ccaccTTCCGTATTGGGCGGCGGCGTCATTTCAGTGAACAGCGCCAGCGACATACCGTCGATACCGCGAGTGATGGTGGACATTAAAATCCAGATCAGCCAGgcgcaagctcagtcaggaacagggcgataccgaaactcacc                                                                                                      >  1:953137‑M1/1‑102 (MQ=255)
                                    caccTTCCGTATTGGGCGGCGGCGTCATTTCAGTGAACAGCGCCAGCGACATACCGTCGATACCGCGAGTGATGGTGGACATTAAAATCCAGATCAGCCAGgcgcaagctcagtcaggaacagggcgataccgaaactcaccg                                                                                                     >  1:365760‑M1/1‑101 (MQ=255)
                                         tCCGTATTGGGCGGCGGCGTCATTTCAGTGAACAGCGACAGCGACATACCGTCGATACCGCGAGTGATGGTGGACATTAAAATCCAGATCAGCCAGgcgcaagctcagtcatgaaccgcgcgataccgacactcaccggggcg                                                                                                >  2:254328‑M1/1‑96 (MQ=255)
                                             tATTGGGCGGCGGCGTCATTTCAGTGAACAGCGCCAGCGACATACCGTCGATACCGCGAGTGATGGTGGACATTAAAATCCAGATCAGCCAGgcgcaagctcagtcaggaacagggcgataccgaaactcaccgggacggcga                                                                                            <  1:788295‑M1/143‑52 (MQ=255)
                                                  ggcggcggcGTCATTTCAGTGAACAGCGCCAGCGACATACCGTCGATACCGCGCGTGATGGTGGACATTAAAATCCAGATCAGCCAGgcccaagcccaggcaggaacagggcgataccgaacctcacccggacggcgagcacg                                                                                       >  2:458317‑M1/1‑87 (MQ=255)
                                                       cggcgTCATTTCAGTGAACAGCGCCAGCGACATACCGTCGATACCGCGAGTGATGGTGGACATTAAAATCCAGATCAGCCAGgcgcaagctcagtcaggaacagggcgataccgaaactcaccgggacagcggccagcagcgc                                                                                  >  1:1113453‑M1/1‑82 (MQ=255)
                                                        ggcgTCATTTCAGTGATCAGCGCCAGCGACATACCGTCGATACCGCGAGTGATGGTGGACATTAAAATCCAGATCAGCCAGgcgcaagctcagtcaggaacagggcgataccgaaactcaccgggacggcgatcagcagcgcg                                                                                 <  1:1148792‑M1/143‑63 (MQ=255)
                                                        ggcgTCATTTCAGTGAACAGCGCCAGCGACATACCGTCGATACCGCGAGTGATGGTGGACATTAAAATCCAGATCAGCCAGgcgcaagctcagtcaggaacagggcgataccgaaactcaccgggacggcgatcagcagcgcg                                                                                 >  1:1014228‑M1/1‑81 (MQ=255)
                                                              accacagtgAACAGCGCCAGCGACATACCGGCGATACCGCGAGTGATGGTGGACATTAAAATCCAGATCAGCCAGgcgcaagctcagtcaggaacagggcgataccgaaactcaccgggacggcgatcagcagcgcgataaac                                                                           <  1:728025‑M1/139‑69 (MQ=255)
                                                               tttCAGTGAACAGCGCCAGCGACATACCGTCGATACCGCGAGTGATGGTGGACATTAAAATCCAGATCAGCCAGgcgcaagctacgtcacgatcagggcgataccgaaacccaccgggccggcgatcagcagcgagagcaacg                                                                          >  2:894440‑M1/1‑74 (MQ=255)
                                                                 tgagtGACCAGCGCCAGCTGCATACCGTCGATACCGCGAGTGATGGTGGAGATTAAAATCCAGATCAGCCAGgcgcaaggtcagtcaggaacagggcgataccgaagctcaccgggacggcgatcagcagcgcgataaacgaa                                                                        <  2:1096868‑M1/141‑72 (MQ=255)
                                                                     cgaACAGCGCCAGCGACATACCGTCGATACCGCGAGTGATGGTGGACATTAAAATCCAGATCAGCCAGgcgcaagctcagtcaggaccagggccatcccggatctcaccgggacggcgatcagcagcgcgataaacgaagtca                                                                    <  1:612623‑M1/142‑76 (MQ=255)
                                                                          aGCGCCAGCGACGTACCGTCGATACCGCGAGTGATGGTGGACATTAAAAGCCAGATCAGCCAGgcgcaagctcagtcaggaacagggcgataccgaaactcaccgggacggcgatcagcagcgcgataaacgaagtcaccaac                                                               <  1:679095‑M1/143‑81 (MQ=255)
                                                                                        aCCGTCGATACCGCGAGTGATGGTGGACATTAAAATCCAGATCAGCCAGgcgcaagctcagtcaggaacagggcgataccgaaactcaccgggacggcgatcagcagcgcgataaacgcagtcacgaacgtaccgtagcgcgg                                                 >  2:84436‑M1/1‑49 (MQ=255)
                                                                                        aCCGCCGATACCGCGAGTGATGGTGGACATTAAAATCCAGCTCAGCCAGgcgcaagctcagtcaggaacagggcgataccgaaactcaccgggacggcgatcagcagcgcgataaacgaagtcaccaacgtaccgtagatcgg                                                 <  1:789890‑M1/143‑95 (MQ=255)
                                                                                        aCCGCCGATACCGCGAGGGATGGGGGACATTAAAATCCAGAGCAGCCAGgcgcaagctcagtcaggaacagggcgataccgaaactcaccgggacggcgatcagcagcgcgataaacgaagtcaccaacgtaccgtagatcgg                                                 <  1:992089‑M1/143‑95 (MQ=255)
                                                                                         ccGTCGATACCGCGAGTGATGGTGGACATTAAAATCCAGATCAGCCAGgcgcaagctcagtcaggaacagggcgataccgaaactcaccgggacggagaacagcagcgcgacaaacgaagtcaccaacgtaccgtagaccggc                                                >  2:468181‑M1/1‑48 (MQ=255)
                                                                                           gTCGATACCGCGAGTGATGGTGGACATTAAAATCCAGATCAGCCAGgcgcaagctcagtcaggaacagggcgataccgaaactcaccgggacggcgatcagcagcgcgataaacgaagtcaccaacgtaccgtagatcggcac                                              <  1:1042648‑M1/143‑98 (MQ=255)
                                                                                             cGATACCGCGTGTTATGGGGGACATTAAAATCCAGATCAGCCAGgcgcaagctcagtcaggaacagggcgataccgaaactcaccgggacggcgatcagcagcgcgataaacgaagtcaccaacgtaccgtagatcggcacca                                            <  2:1052711‑M1/143‑100 (MQ=255)
                                                                                                  gcgcgaGTGACGGTGGCCCTTAAACACCAGATCAGCCAGgcgcaagctcagtcaggaacagggcgataccgaaactcaccgggacggcgatcagcagcgcgataaacgaagtcaccaacgtaccgtagatcggcaccagcgcc                                       <  1:488571‑M1/142‑105 (MQ=255)
                                                                                                  ccGCGAGTGATGGTGGACATTAAAATCCAAATCAGCCAGgcgcaagctcagtcaggaacagggcgataccgaaactaaccgggacgacgatcagcagcgcgaaaaacgcagtcaccaacgtacggtagaccggcaccaacgcc                                       >  2:212873‑M1/1‑39 (MQ=255)
                                                                                                    gcgAGTGATGGTGGACATTAAAATCCAGATCAGCCAGgcgcaagctcagtcaggaacagggcgataccgaaactcaccgggacggcgatcagcagcgcgataaacgaagtcaccaacgtaccgtagatcggcaccagcgcccg                                     >  1:158711‑M1/1‑37 (MQ=255)
                                                                                                       aGTGATGGTGGACATTAAAATCCAGATCAGCCAGgcgcaagctcagtcaggaacagggcgataccgaaactcaccgggacggcgattagcagcgcgataaacgaagtcaccaacgtaccgtagatcggcaccagcgccccgta                                  <  1:510507‑M1/143‑110 (MQ=255)
                                                                                                       aGTGATGGTGGACATTAAAACCCAGATCAGCCAGgcgcaagctcagtcaggaacagggcgataccgaaactcaccgggacggcgatcagcagcgcgataaacgaagtcaccaacgtaccgtagatcggcaccagcgccccgta                                  <  1:412771‑M1/143‑110 (MQ=255)
                                                                                                         tGATGGTGGACATTAAAATCCAGATCAGCCAGgcgcaagctcagtcaggaacagggcgataccgaaactcaccgggacggcgatcagcagcgcgataaacgaagtcaccaacgtaccgtagatcggcaccagcgccccgtaga                                >  1:956421‑M1/1‑32 (MQ=255)
                                                                                                              gtggACATTAAAATCCAGATCAGCCAGgcgcaagctcagtcaggaacagggcgataccgaaactcaccgggacggcgatcagcagcgcgataaacgaagtcaccaacgtgccgtagatcggcaccagcgccccgtagatatcg                           >  1:148184‑M1/1‑27 (MQ=255)
                                                                                                              gtggACATTAAAATCCAGATCAGCCAGgcgcaagctcagtcaggaacagggcgataccgaaactcaccgggacggcgatcagcagcgcgataaacgaagtcaccaacgtaccgtagatcggcaccagcgccccgtagatatcg                           >  2:939356‑M1/1‑27 (MQ=255)
                                                                                                              gtggACATTAAAATCCAGATCAGCCAGgcgcaagctcagtcaggaacagggcgataccgaaactcaccgggacggcgatcagcagcgcgataaacgaaggcaccaacggaccgtagatcggcaccagcgccccggagatgtcg                           >  2:784617‑M1/1‑27 (MQ=255)
                                                                                                                ggACATTAAAATCCAGATCAGCCAGgcgcaagctcagtcaggaacagggcgataccgaaactcaccgggacggcgatcggcagcgcgataaacgaagtcaccaacgtaccgcagataggcaccagagccccgtagacagcgta                         >  2:92620‑M1/1‑25 (MQ=255)
                                                                                                                 gACATTAAAATCCAGATCAACCAGgcgcaagctcagtcaggaacagggcgataccgaaactcaccgggacggcgatcagcagcgcgataaacgaagtcaccaacgtaccgtagatcggcaccagcgccccgtagatatcgttc                        >  1:250987‑M1/1‑24 (MQ=255)
                                                                                                                    aTTAAAATCCAGATCAGCCAGgcgcaagctcagtcaggaacagggcgataccgaaactcaccgggacggcgctcagcagcccgataaacgaagtcaccaacgtaccgtagcgaggcacccgcgccccgtagatatcgttcggt                     >  1:599835‑M1/1‑21 (MQ=255)
                                                                                                                    aTTAAAATCCAGATCAGCCAGgcgcaagctcagtcaggaacagggcgataccgaaactcaccgggacggcgatcagcagcgcgataaacgaagtcaccaacgtaccgtagatcggcaccagcgccccgtagatatcgttcggt                     >  1:398369‑M1/1‑21 (MQ=255)
                                                                                                                          gtcCAGATCAGCAAGgcgcaagctcaggcaggaacagggcgataccgaaacgcaccgggacggcgatcagcagcgcgataaacgaagtcaccaacgtaccgtagatcggcaccagcgccccgtagatatcgttcggtgcatcc               <  1:982752‑M1/142‑129 (MQ=255)
                                                                                                                            ccAGATCAGCCAGgcgcaagctcagtcaggaacagggcgataccgaaactccccgggccagcaaatagcagcgcgataaccgaagtcaccaacgtaccgtagatcggcgccagcgccccgtaaataacgctcggagccagccc             >  2:780566‑M1/1‑13 (MQ=255)
                                                                                                                             cAGATCAGCCAGgcgcaagctcagtcaggaacagggcgataccgaaactcaccgggacggcgatcagcagcgcgataaacgaagtcaccaacgtaccgtagatcggcaccagcgccccgtagatatcgttcggtgcatcccac            <  1:772486‑M1/143‑132 (MQ=255)
                                                                                                                               gATCAGCCAGgcgcaagctcagtcaggaacagggcgataccgaaactcaccgggacggcgatcagcagcgcgataaacgaagtcaccaacgtaccgtagatcggcaccagcgccccgtagatatcgttcggtgcatcccactc          >  1:953673‑M1/1‑10 (MQ=255)
                                                                                                                                 tCAGCCAGgcgcaagctcagtcaggaacagggcgataccgaaactcaccgggacggcgatcagcagcgcgataaacgaagtcaccaacgtaccgtagatcggcaccagcgccccgtagatatcgttcggtgcatcccactctt        <  1:458454‑M1/143‑136 (MQ=255)
                                                                                                                                   aGCCAGgcgcaagctcagtccgggacagggcgataccgaagctcaccgggacggcgatcagcagcgcgataaacgaagtcaccaacgtaccgtagatcggcaccagcgccccgtagatatcgttcggtgcatcccactctttg      <  2:1183175‑M1/143‑138 (MQ=255)
                                                                                                                                      cAGgcgcaagctcagtcaggaacagggcgataccgaaactcaccgggacggccatcagcagcgcgataaacgaagtcaccaacgtaccgtagatcggcaccagcgccccgtagagatcgttcggtgcatcgcactctttggtc   >  1:251356‑M1/1‑3 (MQ=255)
                                                                                                                                       aGgagcaagctcagtcaggaacagggcgataccgaaactcaccgggacggcgattagcagcgcgataaacgaagtcaccaacgtaccgtagatcggcaccagagctccgtagatatcgttcggtgcatccgactctttggtac  >  1:619033‑M1/1‑2 (MQ=255)
                                                                                                                                        |                                                                                                                                             
ACAGCCCGCTACCCGCCAGAGCGTTCGCCAGACCACCACCTTCCGTATTGGGCGGCGGCGTCATTTCAGTGAACAGCGCCAGCGACATACCGTCGATACCGCGAGTGATGGTGGACATTAAAATCCAGATCAGCCAGAACAGGCCGAAGGCCATCGTCGCCATCGAGAGCGTCAGCGCAATACGGTTTTTGAGGCGGCGACGCGCCTGCATTTTGCGGCGAGATTCAGCCAGCGCCGCAGTGGTTTGCATTTCAACCATAGCCATTAGCGTGCCCCCT  >  NC_000913/3909175‑3909452

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: