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

ATTAACAGCCCGCTACCCGCCAGAGCGTTCGCCAGACCACCACCTTCCGTATTGGGCGGCGGCGTCATTTCAGTGAACAGCGCCAGCGACATACCGTCGATACCGCGAGTGATGGTGGACATTAAAATCCAGATCAGCCAGAACAGGCCGAAGGCCATCGTCGCCATCGAGAGCGTCAGCGCAATACGGTTTTTGAGGCGGCGACGCGCCTGCATTTTGCGGCGAGATTCAGCCAGCGCCGCAGTGGTTTGCATTTCAACCATAGCCATTA  >  NC_000913/3909170‑3909440
                                                                                                                                            |                                                                                                                                  
cttaACCGCCCGCGACCCGCCAGAGCGTTCGCCAGACCACCACCTTCCGTATTGGGCGGCGGCGTCATTTCAGTGAACCGCGCCAGCGACATACCGTCGATACCGCGAGTGATGGTGGACATTAAAATCCAGATCAGCCAGgc                                                                                                                                  <  2:349932‑M1/142‑3 (MQ=255)
   accagCCCGCCACCCGCCAGAGCGTTCGCCAGACCACCACCTTCCGTATTGGGCGGCGGCGTCATTTCAGTGAACAGCGCCAGCGACATACCGTCGATACCGCGAGTGATGGTGGACATTAAAATCCAGATCAGCCAGgcgca                                                                                                                               <  2:572870‑M1/141‑6 (MQ=255)
    aCAGCCCGCTACCCGCCAGAGCGGTCGCCAGACCACCACCTTCCGTATTGGGCGGCGGCGTCATTTCAGTGAACAGCGCCAGCGACATACCGTCGATACCGCGAGTGATGGTGGACATTAAAATCCAGATCAGCCAGgcgcaa                                                                                                                              <  1:643416‑M1/143‑7 (MQ=255)
        cccGCTACCCGCCAGAGCGTTCGCCAGACCACCACCTTCCGTACTGGGCGGCGGCGTCATTTCAGTGAACAGCGCCAGCGACATACCGTCGATACCGCGAGTGCTGGGGGACATCGAAATCCAGACCAGCCAGgggcgagcac                                                                                                                          >  1:575422‑M1/1‑133 (MQ=255)
          cGCTACCCGCCAGAGCGTTCGCCAGACCACCACCTTCCGTATTGGGCGGCGGCGTCATTTCAGTGAACAGCGCCAGCGACATACCGTCGATACCGCGAGTGATGGTGGACATTAAAATCCAGATCAGCCAGgcgcaagctcag                                                                                                                        <  2:687720‑M1/143‑13 (MQ=255)
                     agagCGTTCGCCAGACCACCACCTTCCGTATTGGGCGGCGGCGTCATTTCAGTGAACAGCGCCGGCGACATACCGTCGATACCGCGAGTGATGGTGGACATTAAAATCCAGATCAGCCAGgcgcaagctcagtcaggaacagg                                                                                                             >  2:9040‑M1/1‑120 (MQ=255)
                       ggggttcGCCAGACACCCGCCTGCGCTATTGGGCGGCGGCGTGATGTCAGTGAACAGGGCCAGCGACATACCGTCGATACCGCGAGTGATGGTGGACATTAAAATCCAGATCAGCCAGgcgcaagctcagtcaggaacagggc                                                                                                           <  2:523408‑M1/140‑26 (MQ=255)
                            tCGCCAGACCACCACCTTCCGTATTGGGCGGCGGCGTCATTTCAGTGAACAGCGCCAGCGACATACCGTCGATACCGCGAGTGATGGTGGACATTAAAATCCAGATCAGCCAGgcgcaagctcagtcaggaacagggcgatac                                                                                                      <  2:559093‑M1/143‑31 (MQ=255)
                               ccAGACCCCCACCTTCCGTATGGGGCGGCGGCGTCATTTCAGTGAACAGCGCCAGCGACATACCGTCGATACCGCGAGTGATGGTGGACATTAAAATCCAGATCAGCCAGgcgcaagctcagtcaggaacagggcgataccga                                                                                                   <  2:466834‑M1/143‑34 (MQ=255)
                                 aGACCACCACCTTCCGGATTGGGCGGCGGCGTCATTTCAGTGAACAGCGCCAGCGACATACCGTCGATACCGCGAGTGATGGTGGACATTAACATCCAGATCAGCCAGgcgcaagcgcagtcaggaaccgggcgcgacgggca                                                                                                 >  1:564427‑M1/1‑108 (MQ=255)
                                  gACCACCACCTTCCGTATTGGGCGGCGGCGTCATTTCAGTGAACAGCGCCAGCGACATACCGTCGGTACCGCGAGTGATGGTGGACATTAAAATCCAGATCAGCCAGgcgcaagctcagtcaggaacagggcgacaccgaaac                                                                                                >  1:357226‑M1/1‑107 (MQ=255)
                                   accaccaccTTCCGTATTGGGCGGCGGCGTCATTTCAGTGAACAGCGCCAGCGACATACCGTCGATACCGCGCGTGATGGTGGACATTAAAAGCCAGCTCAGCCAGgcgcaggctcagtcaggaccagggcgataccgcgacg                                                                                               >  2:195502‑M1/1‑106 (MQ=255)
                                                         ggcggcgTCATTTCAGTGAACAGCGCCAGCGACATACCGTCGATACCGCGAGTGATGGTGGACATTAAAATCCAGATCAGCCAGgcgcaagctcagtcaggaacagggcgataccgaaactcaccgggacggcgaccagcagc                                                                         >  2:514962‑M1/1‑84 (MQ=255)
                                                                tCATTTCAGTGAACAGCGCCAGCGACATACCGTCGATACCGCGAGTGATGGTGGACATTAAAATCCAGATCAGCCAGgcgcaagctcagtcaggaacagggcgataccgaaactcaccgggacggcggtcagccgcgcgataa                                                                  >  1:73371‑M1/1‑77 (MQ=255)
                                                                              agcgccagcgccATACCGTCGATACCGCGAGTGATGGTGGCCATTAAAATCCAGATCAGCCAGgcgccagcgcagtcaggaacagggcgataccgaaactcaccgggacggcgatcagcagcgcgataaacgaggccacccac                                                    >  2:246871‑M1/1‑63 (MQ=255)
                                                                              aGCGCCAGCGACATACCGTCGATACCGCGAGTGATGGTGGACATTAAAATCCAGATCAGCCAGgcgcaagctcagtcaggaacagggcgataccgaaactcaccgggacggcgatcagcagcgcgataaacgaagtcaccaac                                                    <  2:27196‑M1/143‑81 (MQ=255)
                                                                                 gcCAGCGACATACCGTCGATACCGCGAGTGATGGTGGACATTAAAATCCAGATCAGCCAGgcgcaagctcagtcaggaacacggcgataccgaaactcaccgggacggcgctcagcagcgcgataaacgaagccaccaccgcc                                                 >  1:183705‑M1/1‑60 (MQ=255)
                                                                                                       cgcgAGTGATGGTGGACATTAAAATCCAGATCAGCCAGgcgcaagctcagtcaggaacagggcgataccgaaactcaccgggacggcgagcagcagcgcgataaacgaagtcaccaacgtaccgtagatcggcaccagcgccc                           >  1:651714‑M1/1‑38 (MQ=255)
                                                                                                             tGATGGTGGACATTAACATCCAGATCAGCCAGgcgcaagctcagtcaggaacagggcgataccgaaactcaccgggacggcgatcagcagcgcgataaacgaagtcaccaacgtaccgtagatcggcaccagcgccccgtaga                     <  2:193784‑M1/143‑112 (MQ=255)
                                                                                                             cgaTGGTGCACAGTAACAGCCAGATCAGCCAGgcgcaagcgcagtcaggagcagggcgataccgaaactcaccgggacggcgatcagcagcgcgataaacgaagtcaccaacgtaccgtagatcggcaccagcgccccgtaga                     <  1:862‑M1/142‑112 (MQ=255)
                                                                                                                  gtggACATTAAAATCCAGATCAGCCAGgcgcaagctcagtcaggcgcagggcgataccggccctcaccgggacggcgatcagcagcgcgataagcgaagtcaccaacggaccgtagatcggcaccagcgccccgtagatatcg                <  2:272435‑M1/143‑117 (MQ=255)
                                                                                                                  gtggACATTAAAAGCCAGATCAGCCAGgcgcaagctcagtcaggaacagggcgataccgaaactcaccgggacggcgatcagcagcgcgataaacgaagtcaccaacgtaccgtagatcggcaccagcgccccgtagatatcg                <  2:656913‑M1/143‑117 (MQ=255)
                                                                                                                     gACATTAAAATCCAGATCAGCCAGgcgcaagctcagtcaggaacagggcgataccgaaactcaccgggacggcgatcagcagcgcgataaacgaagtcaccaacgtaccgtagatcggcaccagcgccccgtagatatcgttc             >  2:519480‑M1/1‑24 (MQ=255)
                                                                                                                       cATTAAAATCCAGATCAGCCAGgcgcaagctcagtcaggaacagggcgataccgaaactcaccgggacggcgatcagcagcgcgataaacgaagtcaccaacgtaccgtagatcggcaccagcgccccgtagatatcgttcgg           >  2:362207‑M1/1‑22 (MQ=255)
                                                                                                                        actaagagccagatccgcCAGgcgcaagctcagtcaggcacagggcgataccgaaactcaccgggacggcgatcagcagcgcgataaacgaagtcaccaacgtaccgtagatcggcaccagcgccccgtagatatcgttcggt          <  1:72047‑M1/135‑123 (MQ=255)
                                                                                                                         ttAAAATCCAGATCAGCCAGgcgcaagctcagtcaggaacagggcgataccgaaactcaccgggacggcgatcagcagcgcgataaacgaagtcaccaacgtaccgtagatcggcaccagcgccccgtagatatcgttcggtg         <  2:694599‑M1/143‑124 (MQ=255)
                                                                                                                                ccAGATCAGCCAGgcgcaagctcagtcaggaacagggcgataccgaaactcaccgggacggcgatcagcagcgcgataaacgaagtcaccaacgtaccgtagatcggcaccagcgccccgtagatatcgttcggtgcatccca  <  2:663545‑M1/143‑131 (MQ=255)
                                                                                                                                            |                                                                                                                                  
ATTAACAGCCCGCTACCCGCCAGAGCGTTCGCCAGACCACCACCTTCCGTATTGGGCGGCGGCGTCATTTCAGTGAACAGCGCCAGCGACATACCGTCGATACCGCGAGTGATGGTGGACATTAAAATCCAGATCAGCCAGAACAGGCCGAAGGCCATCGTCGCCATCGAGAGCGTCAGCGCAATACGGTTTTTGAGGCGGCGACGCGCCTGCATTTTGCGGCGAGATTCAGCCAGCGCCGCAGTGGTTTGCATTTCAACCATAGCCATTA  >  NC_000913/3909170‑3909440

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: 
Reads not counted as support for junction
read_name Not counted due to insufficient overlap past the breakpoint.
read_name Not counted due to not crossing MOB target site duplication.