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

Missing coverage evidence...
   seq id start end size ←reads reads→ gene description
* * ÷ NC_000913 3909311 3910084 774 34 [0] [0] 34 [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)30 (0.910) 28/484 NT 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

CATCATCACCGGGGTTTGTTCGAACACATCACGCATTACCGCCGCAATGTACGGAATAATCATGATGGCGAGGATCACGCCTGCCGCGAGGATACCGATACCAAATGCGGGGCCAGAGAACAGCGCGCCAACAATCGGGATATTCGACATGATATTGCCGACCGGCTCCTGAAAGTAAACGGCGAACAGCGGCGCAAAGATAAACAGGCCCCACATGCCGTAAACGATACTTGGAATGGCTGCCAGCAGCTCAATGGCGATACCCAGCGGGCGTTTCAGCCAGCCAGGCGCAAGCTCAGTCAGGAACAGGGCGATACCGAAACTCACCGGGACGGCGATCAGCAGCGCGATAAACGAAGTCACCAACGTACCGTAGATCGGCACCAGCGCCCCGTAGATATCGTTCGGTGCATCCCACTCTTTGGTCCATAGGAAAGCCAGACCAAATTTCTGAATGCTCGGCCAGGAGGAGATGATCAGAGAGACAATAATGCCACCCAACATCAATAGCACAATCAGCGCCGCCAGTTTTACCAGCACGCTGAAAATTATGTCGCCCTTTTTAC  >  NC_000913/3909798‑3910363
                                                                                                                                                                                                                                                                                               |                                                                                                                                                                                                                                                                                      
acagaccaaccacaatcgacggcgcagagagcagaatgtcgttaatgaagcgaatcacttctgccagccaggatttacgaccatattccgccagataaatccccgccataatgcccagcggcgtaccgaataccgtggcccacaaaattaacagcccgctacccgccagagcgttcgccagaccaccaccttccgtattgggcggcggcgtcatttcagtgaacagcgccagcgacataccgtcgataccgcgagtgatggtggacattaaaatccagatcagccaggc                                                                                                                                                                                                                                                                                       >  1:27100‑M2/288‑289 (MQ=255)
   gaccaaccacaatcgacggcgcagagagcagaatgtcgttaatgaagcgaatcacttctgccagccaggatttacgaccatattccgccagataaatccccgccataatgcccagcggcgtaccgaataccgtggcccacaaaattaacagcccgctacccgccagagcgttcgccagaccaccaccttccgtattgggcggcggcgtcatttcagtgaacagcgccagcgacataccgtcgataccgcgagtgatggtggacattaaaatccagatcagccagGCGCa                                                                                                                                                                                                                                                                                    >  2:266206‑M2/285‑289 (MQ=255)
     caaacaccaatggccggcgcagagagcggaaggtcgttaatgaagcgaatcactcctgccagccaggatttccgaccatattccgccagataaacccccgccataatgcccagcggcgtaccgaataccgtggcccacaaaattaacagcccgctacccgccagagcgttcgccagaccaccaccttccgtattgggcggcggcgtcatttcagtgaacagcgccagcgacataccgtcgataccgcgagtgatggtggacattaaaatccagatcagccagGCGCAAg                                                                                                                                                                                                                                                                                  <  2:61463‑M2/7‑1 (MQ=255)
              atcgacggcgcagagagcagaatgtcgttaatgaagcgaatcacttctgccagccaggatttacgaccatattccgccagataaatccccgccataatgcccagcggcgtaccgaataccgtggcccacaaaattaacagcccgctacccgccagagcgttcgccagaccaccaccttccgtattgggcggcggcgtcatttcagtgaacagcgccagcgacataccgtcgataccgcgagtgatggtggacattaaaatccagatcagccagGCGCAAGCtcagtaag                                                                                                                                                                                                                                                                         >  1:317866‑M2/274‑289 (MQ=255)
                       caagaaggcagaaggtcgtaagtgaagcgaatcccttctgccagccaggatttacgaccatattccgccagataaatccccgccataatgcccagcggcgtaccgaataccgtggcccacaaaattaacagcccgctacccgccagagcgttcgccagaccaccaccttccgtattgggcggcggcgtcatttcagtgaacagcgccagcgacataccgtcgataccgcgagtgatggtggacattaaaatccagatcagccagGCGCAAGCTCAGTCAGGAACAggg                                                                                                                                                                                                                                                                 <  2:317866‑M2/24‑1 (MQ=255)
                                     gtcgttaatgaagcgaatcacttctgccagccaggatttacgaccatattccgccagataaatccccgccataatgcccagcggcgtaccgaataccgtggcccacaaaattaacagcccgctacccgccagagcgttcgccagaccaccaccttccgtattgggcggcggcgtcatttcagtgaacagcgccagcgacataccgtcgataccgcgagtgatggtggacattaaaatccagatcagccagGCGCAAGCTCAGTCAGGAACAGGGCGATACCGAAACTc                                                                                                                                                                                                                                                   >  2:304820‑M2/251‑288 (MQ=255)
                                                                                          ccagataaatccccgccataatgcccagcggcgtaccgaataccgtggcccacaaaattaacagcccgctacccgccagagcgttcgccagaccaccaccttccgtattgggcggcggcgtcatttcagtgaacagcgccagcgacataccgtcgataccgcgagtgatggtggacattaaaatccagatcagccagGCGCAAGCTCAGTCAGGAACAGGGCGATACCGAAACTCACCGGGACGGCGATCAGCAGCGCGATAAACGAAGTCACCAACGTACCGTAGATc                                                                                                                                                                                             >  1:165986‑M2/198‑289 (MQ=255)
                                                                                                     cccgccataatgcccagcggcgtaccgaataccgtggcccacaaaattaacagcccgctacccgccagagcgttcgccagaccaccaccttccgtattgggcggcggcgtcatttcagtgaacagcgccagcgacataccgtcgataccgcgagtgatggtggacattaaaatccagatcagccagGCGCAAGCTCAGTCAGGAACAGGGCGATACCGAAACTCACCGGGACGGCGATCAGCAGCGCGATAAACGAAGTCACCAACGTACCGTAGATCGGCACCAgcgc                                                                                                                                                                                  <  1:266206‑M2/103‑1 (MQ=255)
                                                                                                               tgcccagcggcgtaccgaataccgtggcccacaaaattaacagcccgctacccgccagagcgttcgccagaccaccaccttccgtattgggcggcggcgtcatttcagtgaacagcgccagcgacataccgtcgataccgcgagtgatggtggacattaaaatccagatcagccagGCGCAAGCTCAGTCAGGAACAGGGCGATACCGAAACTCACCGGGACGGCGATcagcag                                                                                                                                                                                                                               >  2:93433‑M2/177‑234 (MQ=255)
                                                                                                               tgcccagcggcgtaccgaataccgtggcccacaaaattaacagcccgctacccgccagagcgttcgccagaccaccaccttccgtattgggcggcggcgtcatttcagtgaacagcgccagcgacataccgtcgataccgcgagtgatggtggacattaaaatccagatcagccagGCGCAAGCTCAGTCAGGAACAGGGCGATACCGAAACTCACCGGGACGGCGATcagcag                                                                                                                                                                                                                               <  1:93433‑M2/58‑1 (MQ=255)
                                                                                                                      cggcgtaccgaataccgtggcccacaaaattaacagcccgctacccgccagagcgttcgccagaccaccaccttccgtattgggcggcggcgtcatttcagtgaacagcgccagcgacataccgtcgataccgcgagtgatggtggacattaaaatccagatcagccagGCGCAAGCTCAGTCAGGAACAGGGCGATACCGAAACTCACCGGGACGGCGATCAGCAGCGCGATAAACGAAGTCACCAACGTACCGTAGATCGGCACCAGCGCCCCGTAGATATCGTTCg                                                                                                                                                                 >  2:43365‑M2/170‑289 (MQ=255)
                                                                                                                          gtaccgaataccgtggcccacaaaattaacagcccgctacccgccagagcgttcgccagaccaccaccttccgtattgggcggcggcgtcatttcagtgaacagcgccagcgacataccgtcgataccgcgagtgatggtggacattaaaatccagatcagccagGCGCAAGCTCAGTCAGGAACAGGGCGATACCGAAACTCACCGGGACGGCGATCAGCAGCGCGATAAACGAAGTCACCAACGTACCGTAGATCGGCACCAGCGCCCCGTAGATATCGTTCGGTg                                                                                                                                                              <  1:275053‑M2/123‑1 (MQ=255)
                                                                                                                          gcaccaaaaaccgtggcccacaaaattaacagcccccgccccgccaaagcgttcgccagaccaccaccttccgtattgggcgtcggcgtcatttcagtgaacagcgccagcgacataccgtcgataccgcgagtgatggtggacagtaaaatccagatcagccagGCGCAAGCTCAGTCAGGAACAGGGCGATACCGAAACTCACCGGGACGGCGATCAGCAGCGCGATAAACGAAGTCACCAACGTACCGTAGATCGGCACCAGCGCCCCGTAGATATCGTTCGGTg                                                                                                                                                              <  2:165986‑M2/123‑1 (MQ=255)
                                                                                                                            accgaataccgtcgcccacaaaataaccagccccctacacgccagggcgttcgccagaccaccaccttccgtattgggcggcggcgtcatttcagtgaacagcgccagcgacataccgtcgataccgcgagtgatggtggacattaaaatccagatcagccagGCGCAAGCTCAGTCAGGAACAGGGCGATACCGAAACTCACCGGGACGGCGATCAGCAGCGCGATAAACGAAGTCACCAACGTACCGTAGATCGGCACCAGCGCCCCGTAGATATCGTTCGGTGCa                                                                                                                                                            <  2:27100‑M2/125‑1 (MQ=255)
                                                                                                                                       tggcccacaaaattaacagcccgctacccgccagagcgttcgccagaccaccaccttccgtattgggcggcggcgtcatttcagtgaacagcgccagcgacataccgtcgataccgcgagtgatggtggacattaaaatccagatcagccagGCGCAAGCTCAGTCAGGAACAGGGCGATACCGAAACTCACCGGGACGGCGATCAGCAGCGCGATAAACGAAGTCACCAACGTACCGTAGATCGGCACCAGCGCCCCGTAGATATCGTTCGGt                                                                                                                                                               <  1:21935‑M2/122‑1 (MQ=255)
                                                                                                                                       tggcccacaaaattaacagcccgctacccgccagagcgttcgccagaccaccaccttccgtattgggcggcggcgtcatttcagtgaacagcgccagcgacataccgtcgataccgcgagtgatggtggacattaaaatccagatcagccagGCGCAAGCTCAGTCAGGAACAGGGCGATACCGAAACTCACCGGGACGGCGATCAGCAGCGCGATAAACGAAGTCACCAACGTACCGTAGATCGGCACCAGCGCCCCGTAGATATCGTTCGGt                                                                                                                                                               >  2:21935‑M2/153‑274 (MQ=255)
                                                                                                                                       tggcccacaaaattaacagcccgctacccgccagagcgttcgccagaccaccaccttccgtattgggcggcggcgtcatttcagtgaacagcgccagcgacataccgtcgataccgcgagtgatggtggacattaaaatccagatcagccagGCGCAAGCTCAGTCAGGAACAGGGCGATACCGAAACTCACCGGGACGGCGATCAGCAGCGCGATAAACGAAGTCACCAACGTACCGTAGATCGGCACCAGCGCCCCGTAGACATCGTTCGGTGCATCCCACTCTTTg                                                                                                                                                >  2:318121‑M2/153‑289 (MQ=255)
                                                                                                                                        ggcccacaaaattaacagcccgctacccgccagagcgttcgccagaccaccaccttccgtattgggcggcggcgtcatttcagtgaacagcgccagcgacataccgtcgataccgcgagtgatggtggccattaaaatccaaatcagccagGCGCAAGCTCAGTCAGGAACAGGGCGATACCGAAACTCACCGGGACGGCGATCAGCAGCGCGATAAACGAAGTCCCCAAAGTACCGTAGCTCGGCAACAGCGCCCAGTAGATATCGTTCGGTGCATCCCACTCTTTgg                                                                                                                                               >  2:213695‑M2/152‑289 (MQ=255)
                                                                                                                                          cccacaaaattaacagcccgctacccgccagagcgttcgccagaccaccaccttccgtattgggcggcggcgtcatttcagtgaccagcgccagcgacataccgtcgataccgcgagtgatggtggacattaaaatccagatcagccagGCGCAAGCTCAGTCAGGAACAGGGCGATACCGAAACTCACCGGGACGGCGATGAGCAGCGCGATAAACTAAGTCACCAACGCACCGTAGATCGCCACCAGCGCCACGTAGATATCGATCGGTGCATCCACCTCTTTGGt                                                                                                                                              >  1:298530‑M2/150‑288 (MQ=255)
                                                                                                                                                         gcccgctacccgccagagcgttcgccagaccaccaccttccgtattgggcggcggcgtcatttcagtgaacagcgccagcgacataccgtcgataccgcgagtgatggtggacattaaaatccagatcagccagGCGCAAGCTCAGTCAGGAACAGGGCGATACCGAAACTCACCGGGACGGCGATCAGCAGCGCGATAACCGAAGt                                                                                                                                                                                                                >  1:276534‑M2/135‑207 (MQ=255)
                                                                                                                                                         gcccgctacccgccagagcgttcgccagaccaccaccttccgtattgggcggcggcgtcatttcagtgaacagcgccagcgacataccgtcgataccgcgagtgatggtggacattaaaatccagatcagccagGCGCAAGCTCAGTCAGGAACAGGGCGATACCGAAACTCACCGGGACGGCGATCAGCAGCGCGATAAACGAAGt                                                                                                                                                                                                                <  2:276534‑M2/73‑1 (MQ=255)
                                                                                                                                                            cgctacccgccagagcgttcgccagaccaccaccttccgtattgggcggcggcgtcatttcagtgaacagcgccagcgacataccgtcgataccgcgagtgatggtggacattaaaatccagatcagccagGCGCAAGCTCAGTCAGGAACAGGGCGATACCGAAACTCACCGGGACGGCGATCAGCAGCGCGATAAACGAAGTCACCAACGTACCGTAGATCGGAACCAGCGCCCCGTAGATATCGTTCGGTTCATCCCACTCTTTGGTCCATAGGAAAGCACGACCa                                                                                                                           >  2:89616‑M2/132‑289 (MQ=255)
                                                                                                                                                                         agcgttcgccagaccaccaccttccgtattgggcggcggcgtcatttcagtgaacagcgccagcgacataccgtcgataccgcgagtgatggtggacattaaaatccagatcagccagGCGCAAGCTCAGTCAGGAACAGGGCGATACCGAAACTCACCGGGACGGCGATCAGCAGCGCGATAAACGAAGTCACCAACGTACCGTAGATCGGCACCAGCGCCCCGTAGATATCGTTCGGTGCATCCCACTCTTTGGTCCATAGGAAAGCCAGACCAAATTTCTGAATg                                                                                                               >  2:182459‑M2/119‑288 (MQ=255)
                                                                                                                                                                                    taccaccaccttccgtattgggcggcggcgtcatttcagtgaacagcgccagcgacataccgtcgataccgcgagtgatggtggacattaaaatccagatcagccagGCGCAAGCTCAGTCAGGAACAGGGCGATACCGAAACTCACCGGGACGGCGATCAGCAGCGCGATAAACGAAGTCACCAACGTACCGTAGATCGGCACCAGCGCCCCGTAGATATCGTTCGGTGCATCCCACTCTTTGGTCCATAGGAAAGCCAGACCAAATTTCTGAATGCTCGGCCaggag                                                                                                   <  1:304820‑M2/182‑1 (MQ=255)
                                                                                                                                                                                      ccaccacctgccgtattgggcggcggcgtcatgtcagtgaacagcgccagcgacataccgtcgataccgcgagtgatggtggacattaaaatccagatcagccagGCGCAAGCTCAGTCAGGAACAGGGCGATACCGAAACTCACCGGGACGGCGATCAGCAGCGCGATAAACGAAGTCACCAACGTACCGTAGATCGGCACCAGCGCCCCGTAGATATCGTTCGGTGCATCCCACTCTTTGGTCCATAGGAAAGCCAGACCAAATTTCTGAATGCTCGGCCAggagga                                                                                                 <  1:89616‑M2/184‑1 (MQ=255)
                                                                                                                                                                                                        ggcggcggcgtcatttcagtgaacagcgccagcgacataccgtcgataccgcgagtgatggtggacattaaaatccagatcagccagGCGCAAGCTCAGTCAGGAACAGGGCGATACCGAAACTCACCGGGACGGCGATCAGCAGCGCGATAAACGAAGTCACCAACGTACCGTAGATCGGCACCAGCGCCCCGTAGATATCGTTCGGTGCATCCCACTCTTTGGTCCATAGGAAAGCCAGACCAAATTTCTGAATGCTCGGCCAGGAGGAGATGATCAGAGAGACaat                                                                               >  2:83483‑M2/88‑289 (MQ=255)
                                                                                                                                                                                                             cggcgtcatttcagtgaacagcgccagcgacataccgtcgataccgcgagtgatggtggacattaaaatccagatcagccagGCGCAAGCTCAGTCAGGAACAGGGCGATACCGAAACTCACCGGGACGGCGATCAGCAGCGCGATAAACGAAGTCACCAACGTACCGTAGATCGGCACCAGCGCCCCGTAGATATCGTTCGGTGCATCCCACTCTTTGGTCCATAGGAAAGCCAGACCAAATTTCTGAATGCTCGGCCAGGAGGAGATGATCag                                                                                        >  2:71859‑M2/83‑275 (MQ=255)
                                                                                                                                                                                                             cggcgtcatttcagtgaacagcgccagcgacataccgtcgataccgcgagtgatggtggacattaaaatccagatcagccagGCGCAAGCTCAGTCAGGAACAGGGCGATACCGAAACTCACCGGGACGGCGATCAGCAGCGCGATAAACGAAGTCACCAACGTACCGTAGATCGGCACCAGCGCCCCGTAGATATCGTTCGGTGCATCCCACTCTTTGGTCCATAGGAAAGCCAGACCAAATTTCTGAATGCTCGGCCAGGAGGAGATGATCag                                                                                        <  1:71859‑M2/193‑1 (MQ=255)
                                                                                                                                                                                                                                       gcgacataccgtcgataccgcgagtgatggtggacattaaaatccagatcagccagGCGCAAGCTCAGTCAGGAACAGGGCGATACCGAAACTCACCGGGACGGCGATCAGCAGCGCGATAAACGAAGTCACCAACGTACCGTAGATCGGCACCAGCGCCCCGTAGATATCGTTCGGTGCATCCCACTCTTTGGTCCATAGGAAAGCCAGACCAAATTTCTGAATGCTCGGCCAGGAGGAGATGATCAGAGAGACAATAATGCCACCCAACATCAATAGCACAATCAgc                                                >  1:6826‑M2/57‑289 (MQ=255)
                                                                                                                                                                                                                                              accgtcgataccgcgagtgatggtggacattaaaatccagatcagccagGCGCAAGCTCAGTCAGGAACAGGGCGATACCGAAACTCACCGGGACGGCGATCAGCAGCGCGATAAACGAAGTCACCAACGTACCGTAGATCGGCACCAGCGCCCCGTAGATATCGTTCGGTGCATCCCACTCTTTGGTCCATAGGAAAGCCAGACCAAATTTCTGAATGCTCGGCCAGGAGGAGATGATCAGAGAGACAATAATGCCACCCAACATCAATAGCACAATCAGCGCCGCCa                                         <  1:43365‑M2/240‑1 (MQ=255)
                                                                                                                                                                                                                                                                gatggtggacattaaaatccagatcagccagGCGCAAGCTCAGTCAGGAACAGGGCGATACCGAAACTCACCGGGACGGCGATCAGCAGCGCGATAAACGAAGTCACCAACGTACCGTAGATCGGCACCAGCGCCCCGTAGATATCGTTCGGTGCATCCCACTCTTTGGTCCATAGGAAAGCCAGACCAAATTTCTGAATGCTCGGCCAGGAGGAGATGATCAGAGAGACAATAATGCCACCCAACATCAATAGCACAATCAGCGCCGCCAGTTTTACCAGCACGCt                         <  1:182459‑M2/256‑1 (MQ=255)
                                                                                                                                                                                                                                                                                 tccagatcagccagGCGCAAGCTCAGTCAGGAACAGGGCGATACCGAAACTCACCGGGACGGCGATCAGCAGCGCGATAAACGAAGTCACCAACGTACCGTAGATCGGCACCAGCGCCCCGTAGATATCGTTCGGTGCATCCCACTCTTTGGTCCATAGGAAAGCCAGACCAAATTTCTGAATGCTCGGCCAGGAGGAGATGATCAGAGAGACACTAATGCCACCCAACATCAATAGCACAATCAGCCCCGCCAGTTTTACCAGCCAGCTGAACATTCTGTCGACCttt      >  2:26743‑M2/15‑289 (MQ=255)
                                                                                                                                                                                                                                                                                      atcagccagGCGCAAGCTCAGTCAGGAACAGGGCGATACCGAAACTCACCGGGACGGCGATCAGCAGCGCGATAAACGAAGTCACCAACGTACCGTAGATCGGCACCAGCGCCCAGTAGATATCGTTCGGTGCATCCCACTCTTTGGTCCATAGGAAAGCCAGACCAAATTTCTGAATGCTCGGCCCGGAGGAGATGATCAGAGAGACAATAATGCCACCCAACATCAATAGCACAATCAGCCCCGCCAGTTTTACCAGCACGCTGAAAATTAAGTCCCCATTTTTAc  >  2:109736‑M2/10‑288 (MQ=255)
                                                                                                                                                                                                                                                                                      atcagccagGCGCAAGCTCAGTCAGGAACAGGGCGATACCGAAACTAACCGGGACGGCGATCAGCAGCGCGATAAACGAAGTCACCAACGTACCGTAGATCGGCACCAGCGCCCCGTAGATATCGTTCGGTGCATCCCACTCTTTGGTCCACAGGAAAGCCAGACCAAATTTCTGAATGCTCGGCCAGGAGGAGATGATCAGAGAGACAATAATGCCACACAACATCAATAGCCCAATCAGCGCCGCAAGTTTTACCAGAACGCTGAAAATAATGTAGCCCTTTTTAc  >  1:244222‑M2/10‑288 (MQ=255)
                                                                                                                                                                                                                                                                                               |                                                                                                                                                                                                                                                                                      
CATCATCACCGGGGTTTGTTCGAACACATCACGCATTACCGCCGCAATGTACGGAATAATCATGATGGCGAGGATCACGCCTGCCGCGAGGATACCGATACCAAATGCGGGGCCAGAGAACAGCGCGCCAACAATCGGGATATTCGACATGATATTGCCGACCGGCTCCTGAAAGTAAACGGCGAACAGCGGCGCAAAGATAAACAGGCCCCACATGCCGTAAACGATACTTGGAATGGCTGCCAGCAGCTCAATGGCGATACCCAGCGGGCGTTTCAGCCAGCCAGGCGCAAGCTCAGTCAGGAACAGGGCGATACCGAAACTCACCGGGACGGCGATCAGCAGCGCGATAAACGAAGTCACCAACGTACCGTAGATCGGCACCAGCGCCCCGTAGATATCGTTCGGTGCATCCCACTCTTTGGTCCATAGGAAAGCCAGACCAAATTTCTGAATGCTCGGCCAGGAGGAGATGATCAGAGAGACAATAATGCCACCCAACATCAATAGCACAATCAGCGCCGCCAGTTTTACCAGCACGCTGAAAATTATGTCGCCCTTTTTAC  >  NC_000913/3909798‑3910363

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 9 ≤ ATCG/ATCG < 20 ≤ ATCG/ATCG < 35 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: