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 27 [0] [0] 28 [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)27 (0.670) 24/270 0.7 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

ACATGATATTGCCGACCGGCTCCTGAAAGTAAACGGCGAACAGCGGCGCAAAGATAAACAGGCCCCACATGCCGTAAACGATACTTGGAATGGCTGCCAGCAGCTCAATGGCGATACCCAGCGGGCGTTTCAGCCAGCCAGGCGCAAGCTCAGTCAGGAACAGGGCGATACCGAAACTCACCGGGACGGCGATCAGCAGCGCGATAAACGAAGTCACCAACGTACCGTAGATCGGCACCAGCGCCCCGTAGATATCGTTCGGTGCATCCCACTCTTTGGTCCAT  >  NC_000913/3909944‑3910227
                                                                                                                                             |                                                                                                                                              
cttaaccgcccgcgacccgccagagcgttcgccagaccaccaccttccgtattgggcggcggcgtcatttcagtgaaccgcgccagcgacataccgtcgataccgcgagtgatggtggacattaaaatccagatcagccaggc                                                                                                                                               <  2:349932‑M2/2‑1 (MQ=255)
   accagcccgccacccgccagagcgttcgccagaccaccaccttccgtattgggcggcggcgtcatttcagtgaacagcgccagcgacataccgtcgataccgcgagtgatggtggacattaaaatccagatcagccagGCGCa                                                                                                                                            <  2:572870‑M2/5‑1 (MQ=255)
    acagcccgctacccgccagagcggtcgccagaccaccaccttccgtattgggcggcggcgtcatttcagtgaacagcgccagcgacataccgtcgataccgcgagtgatggtggacattaaaatccagatcagccagGCGCaa                                                                                                                                           <  1:643416‑M2/6‑1 (MQ=255)
        cccgctacccgccagagcgttcgccagaccaccaccttccgtactgggcggcggcgtcatttcagtgaacagcgccagcgacataccgtcgataccgcgagtgctgggggacatcgaaatccagaccagccagGGGCGagcac                                                                                                                                       >  1:575422‑M2/134‑141 (MQ=255)
          cgctacccgccagagcgttcgccagaccaccaccttccgtattgggcggcggcgtcatttcagtgaacagcgccagcgacataccgtcgataccgcgagtgatggtggacattaaaatccagatcagccagGCGCAAGCtcag                                                                                                                                     <  2:687720‑M2/12‑1 (MQ=255)
                     agagcgttcgccagaccaccaccttccgtattgggcggcggcgtcatttcagtgaacagcgccggcgacataccgtcgataccgcgagtgatggtggacattaaaatccagatcagccagGCGCAAGCTCAGTCAGGAACAgg                                                                                                                          >  2:9040‑M2/121‑143 (MQ=255)
                       ggggttcgccagacacccgcctgcgctattgggcggcggcgtgatgtcagtgaacagggccagcgacataccgtcgataccgcgagtgatggtggacattaaaatccagatcagccagGCGCAAGCTCAGTCAGGAACAGGGc                                                                                                                        <  2:523408‑M2/25‑1 (MQ=255)
                            tcgccagaccaccaccttccgtattgggcggcggcgtcatttcagtgaacagcgccagcgacataccgtcgataccgcgagtgatggtggacattaaaatccagatcagccagGCGCAAGCTCAGTCAGGAACAGGGCGATAc                                                                                                                   <  2:559093‑M2/30‑1 (MQ=255)
                               ccagacccccaccttccgtatggggcggcggcgtcatttcagtgaacagcgccagcgacataccgtcgataccgcgagtgatggtggacattaaaatccagatcagccagGCGCAAGCTCAGTCAGGAACAGGGCGATACCGa                                                                                                                <  2:466834‑M2/33‑1 (MQ=255)
                                 agaccaccaccttccggattgggcggcggcgtcatttcagtgaacagcgccagcgacataccgtcgataccgcgagtgatggtggacattaacatccagatcagccagGCGCAAGCGCAGTCAGgaaccgggcgcgacgggca                                                                                                              >  1:564427‑M2/109‑134 (MQ=255)
                                  gaccaccaccttccgtattgggcggcggcgtcatttcagtgaacagcgccagcgacataccgtcggtaccgcgagtgatggtggacattaaaatccagatcagccagGCGCAAGCTCAGTCAGGAACAGGGCGACACCGAAAc                                                                                                             >  1:357226‑M2/108‑143 (MQ=255)
                                   accaccaccttccgtattgggcggcggcgtcatttcagtgaacagcgccagcgacataccgtcgataccgcgcgtgatggtggacattaaaagccagctcagccagGCGCAGGCTCAGTCAGGACCAGGGCGataccgcgacg                                                                                                            >  2:195502‑M2/107‑138 (MQ=255)
                                                         ggcggcgtcatttcagtgaacagcgccagcgacataccgtcgataccgcgagtgatggtggacattaaaatccagatcagccagGCGCAAGCTCAGTCAGGAACAGGGCGATACCGAAACTCACCGGGACGGCGACCagcagc                                                                                      >  2:514962‑M2/85‑143 (MQ=255)
                                                                tcatttcagtgaacagcgccagcgacataccgtcgataccgcgagtgatggtggacattaaaatccagatcagccagGCGCAAGCTCAGTCAGGAACAGGGCGATACCGAAACTCACCGGGACGGCGGTCAGCCGCGCGATaa                                                                               >  1:73371‑M2/78‑143 (MQ=255)
                                                                              agcgccagcgccataccgtcgataccgcgagtgatggtggccattaaaatccagatcagccagGCGCCAGCGCAGTCAGGAACAGGGCGATACCGAAACTCACCGGGACGGCGATCAGCAGCGCGATAAACGAGGccacccac                                                                 >  2:246871‑M2/64‑143 (MQ=255)
                                                                              agcgccagcgacataccgtcgataccgcgagtgatggtggacattaaaatccagatcagccagGCGCAAGCTCAGTCAGGAACAGGGCGATACCGAAACTCACCGGGACGGCGATCAGCAGCGCGATAAACGAAGTCACCAAc                                                                 <  2:27196‑M2/80‑1 (MQ=255)
                                                                                 gccagcgacataccgtcgataccgcgagtgatggtggacattaaaatccagatcagccagGCGCAAGCTCAGTCAGGAACACGGCGATACCGAAACTCACCGGGACGGCGCTCAGCAGCGCGATAAACGAAGCCACCAccgcc                                                              >  1:183705‑M2/61‑141 (MQ=255)
                                                                                                       cgcgagtgatggtggacattaaaatccagatcagccagGCGCAAGCTCAGTCAGGAACAGGGCGATACCGAAACTCACCGGGACGGCGAGCAGCAGCGCGATAAACGAAGTCACCAACGTACCGTAGATCGGCACCAGCGccc                                        >  1:651714‑M2/39‑143 (MQ=255)
                                                                                                             tgatggtggacattaacatccagatcagccagGCGCAAGCTCAGTCAGGAACAGGGCGATACCGAAACTCACCGGGACGGCGATCAGCAGCGCGATAAACGAAGTCACCAACGTACCGTAGATCGGCACCAGCGCCCCGTAGa                                  <  2:193784‑M2/111‑1 (MQ=255)
                                                                                                             cgatggtgcacagtaacagccagatcagccagGCGCAAGCGCAGTCAGGAGCAGGGCGATACCGAAACTCACCGGGACGGCGATCAGCAGCGCGATAAACGAAGTCACCAACGTACCGTAGATCGGCACCAGCGCCCCGTAGa                                  <  1:862‑M2/111‑1 (MQ=255)
                                                                                                                  gtggacattaaaatccagatcagccagGCGCAAGCTCAGTCAGGCGCAGGGCGATACCGGCCCTCACCGGGACGGCGATCAGCAGCGCGATAAGCGAAGTCACCAACGGACCGTAGATCGGCACCAGCGCCCCGTAGATATCg                             <  2:272435‑M2/116‑1 (MQ=255)
                                                                                                                  gtggacattaaaagccagatcagccagGCGCAAGCTCAGTCAGGAACAGGGCGATACCGAAACTCACCGGGACGGCGATCAGCAGCGCGATAAACGAAGTCACCAACGTACCGTAGATCGGCACCAGCGCCCCGTAGATATCg                             <  2:656913‑M2/116‑1 (MQ=255)
                                                                                                                     gacattaaaatccagatcagccagGCGCAAGCTCAGTCAGGAACAGGGCGATACCGAAACTCACCGGGACGGCGATCAGCAGCGCGATAAACGAAGTCACCAACGTACCGTAGATCGGCACCAGCGCCCCGTAGATATCGTTc                          >  2:519480‑M2/25‑143 (MQ=255)
                                                                                                                       cattaaaatccagatcagccagGCGCAAGCTCAGTCAGGAACAGGGCGATACCGAAACTCACCGGGACGGCGATCAGCAGCGCGATAAACGAAGTCACCAACGTACCGTAGATCGGCACCAGCGCCCCGTAGATATCGTTCgg                        >  2:362207‑M2/23‑143 (MQ=255)
                                                                                                                        actaagagccagatccgccagGCGCAAGCTCAGTCAGGCACAGGGCGATACCGAAACTCACCGGGACGGCGATCAGCAGCGCGATAAACGAAGTCACCAACGTACCGTAGATCGGCACCAGCGCCCCGTAGATATCGTTCGGt                       <  1:72047‑M2/122‑1 (MQ=255)
                                                                                                                         ttaaaatccagatcagccagGCGCAAGCTCAGTCAGGAACAGGGCGATACCGAAACTCACCGGGACGGCGATCAGCAGCGCGATAAACGAAGTCACCAACGTACCGTAGATCGGCACCAGCGCCCCGTAGATATCGTTCGGTg                      <  2:694599‑M2/123‑1 (MQ=255)
                                                                                                                                ccagatcagccagGCGCAAGCTCAGTCAGGAACAGGGCGATACCGAAACTCACCGGGACGGCGATCAGCAGCGCGATAAACGAAGTCACCAACGTACCGTAGATCGGCACCAGCGCCCCGTAGATATCGTTCGGTGCATCCCa               <  2:663545‑M2/130‑1 (MQ=255)
                                                                                                                                             gcgcAAGCTCAGTCAGGAACAGGGCGATACCGAAACTCACCGGGACGGCGAACAGCAGCGCGATAAACGAAGTCACCCACGTACCGTAGATCGCCACCAGCGCCCCGAAGAAAACGTTCGGCGCCTCCCACCCGTTGGTCCAt  >  1:9398/1‑143 (MQ=255)
                                                                                                                                             |                                                                                                                                              
ACATGATATTGCCGACCGGCTCCTGAAAGTAAACGGCGAACAGCGGCGCAAAGATAAACAGGCCCCACATGCCGTAAACGATACTTGGAATGGCTGCCAGCAGCTCAATGGCGATACCCAGCGGGCGTTTCAGCCAGCCAGGCGCAAGCTCAGTCAGGAACAGGGCGATACCGAAACTCACCGGGACGGCGATCAGCAGCGCGATAAACGAAGTCACCAACGTACCGTAGATCGGCACCAGCGCCCCGTAGATATCGTTCGGTGCATCCCACTCTTTGGTCCAT  >  NC_000913/3909944‑3910227

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: