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 25 [0] [0] 26 [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)25 (0.650) 22/270 0.9 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

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

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: