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 38 [0] [0] 37 [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)37 (0.750) 30/270 0.6 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

TGCCGACCGGCTCCTGAAAGTAAACGGCGAACAGCGGCGCAAAGATAAACAGGCCCCACATGCCGTAAACGATACTTGGAATGGCTGCCAGCAGCTCAATGGCGATACCCAGCGGGCGTTTCAGCCAGCCAGGCGCAAGCTCAGTCAGGAACAGGGCGATACCGAAACTCACCGGGACGGCGATCAGCAGCGCGATAAACGAAGTCACCAACGTACCGTAGATCGGCACCAGCGCCCCGTAGATATCGTTCGGTGC  >  NC_000913/3909953‑3910208
                                                                                                                                    |                                                                                                                           
tcgctacccgccggggcgttcgccagacctccaccctccgtatggggcggcggcgtcatttcagtgaacagcgccagcgacataccgtcgataccgcgagtgatggtggacattaaaatccagatcagccagGCGCAAGCtca                                                                                                                   <  2:18522‑M2/11‑1 (MQ=255)
    tacccgccagagcgttcgccagaccaccaccttccgtattgggcggcggcgtcatttcagtgaacagcgccagcgacataccgtcgataccgcgagtgatggtggacattaaaagccagatcagccagGCGCAAGctcagtcg                                                                                                               >  1:709390‑M2/129‑142 (MQ=255)
         gccagagcgttcgccagaccaccaccttccgtattgggcggcggcgtcatttcagtgaacagcgccagcgacataccgtcgataccgcgagtgatggtggacattaaaatccagatcagccagGCGCAAGCTCAGTCAGGAAc                                                                                                          >  1:136225‑M2/124‑143 (MQ=255)
          ccagggcgttcgccagcccaccaccttccgtatcgggcggcggcgtcatttcagtgaacagcgccagcgacataccgtcgataccgcgagtgatggtggacattaaaatccagatcagccagGCGCAAGCTCAGTCAGGAACa                                                                                                         <  2:571108‑M2/21‑1 (MQ=255)
              agcgttcgccagaccaccaccttccgtattgggcggcggcgtcatttcagtgaacagcgccagcgacataccgtagataccgcgagtgatggtggacattaaaatccagatcagccagGCGCAAGCTCAGTCAGGAACAGGGc                                                                                                     >  1:177989‑M2/119‑143 (MQ=255)
                cgttcgccagaccaccaccttccgtattgggcggcggcgtcatttcagtgaacagcgccagcgacataccgtcgataccgcgagtgatggtggacattaaaatccagatcagccagGCGCAAGCTCAGTCAGGAACAGGGCGa                                                                                                   <  1:749181‑M2/27‑1 (MQ=255)
                       cagcccaccatcctccgtattgggcggcggcgtcacggcagtgaacagcgccagcgacctaccgtcgataccgcgagtgatcgtggacattaaaatccagatcagccagGCGCAAGCTCAGTCAGGAACAGGACGATACCGaa                                                                                            <  2:448651‑M2/34‑1 (MQ=255)
                       cagaccaccaccttccgtattgggcggcggcgtcatttcagtgaacagcgccagcgacataccgtcgataccgcgagtgatggtggacattaaaatccagatcagccagGCGCAAGCTCAGTCAGGAACAggg                                                                                                      >  2:62032‑M2/110‑133 (MQ=255)
                       cagaccaccaccttccgtattgggcggcggcgtcatttcagtgaacagcgccagcgacataccgtcgataccgcgagtgatggtggacattaaaatccagatcagccagGCGCAAGCTCAGTCAGGAACAggg                                                                                                      <  1:62032‑M2/24‑1 (MQ=255)
                         cagcaccgccttccgtatggggcggcgccgccacttcagtgaacagcgccagcgacatcccgtggataccgcgagtgatggtggacattaaaatccatatcagccagGCGCAAGCTCAGTCAGGAACAGGGCGATACCGacac                                                                                          <  1:366892‑M2/36‑1 (MQ=255)
                                  cttccgtattgggcggcggcgtcatttcagtgaacagcgccagcgacataccgtcgataccgcgagtgatggtggacattaaaatccagatcagccagGCGCAAGCTCAGTCAGGAACAGGGCGATACCGAAACTCACCGGGa                                                                                 >  1:536161‑M2/99‑143 (MQ=255)
                                        tattgcgcggcggcgtcatttcagtgaacagcgccagcgacataccgtcgataccgcgagtgattgtgtacattataatccagatcagccagGCGCAAGCTCAGTCAGGAACAGGGCGATACCGAAACTCACCGGGACGGCGa                                                                           <  1:447853‑M2/51‑1 (MQ=255)
                                        gagtgggcggcgacggaagttcagtgaacggcgccagcgacataccgtcgacaccgcgagtgatggtggacattaaaatccagatcagccagGCGCAAGCTCAGTCAGGAACAGGGCGATACCGAAACTCACCGGGACGGCGa                                                                           <  1:878495‑M2/51‑1 (MQ=255)
                                               cggcggcgtcatttcagtgaacagcgccagcgacataccgtcgataccgcgagtgatggtggacattaaaatccagatcagccagGCGCAAGCTCAGTCAGGAACAGGGCGATACCGAAACTCACCGGGACGGCGATcagcag                                                                    >  1:297526‑M2/86‑143 (MQ=255)
                                                 gcggcgtcatttcagtgaacagcgccagcgacataccgtcgataccgcgagtgatggtggacattaaaatccagatcagccagGCGCAAGCTCAGTCAGGAACAGGGCGATACCGAAACTCACCGGGACGGCGATCAGCAgcg                                                                  <  1:721482‑M2/60‑1 (MQ=255)
                                                        catttggggggacagcgcccgcgacataccggcgataccgcgagtgatggtggacagtaaaatccagatcagccagGCGCAAGCTCAGTCAGGAACAGGGCGATACCGAAACTCACCGGGACGGCGATCAGCAGCGCGATaaa                                                           <  2:467026‑M2/67‑1 (MQ=255)
                                                            tcagtggactgcgccagagagatagggtcgataccgagagtgattgtggacattaaaatccagatcagccagGCGCAAGCTCAGTCAGGAACAGGGCGATACCGAAACTCACCGGGACGGCGATCAGCAGCGCGATAAACGaa                                                       <  2:318403‑M2/71‑1 (MQ=255)
                                                               gtgaacagcgccagcgacataccgtcgataccgcgagtgatggtggacattaaaatccagatcagccagGCGCAAGCTCAGTCAGGAACAGGGCGATACCGAAACTCACCGGGACGGCGATGAGCAGCGCGACAAACGAAGTc                                                    >  2:786944‑M2/70‑143 (MQ=255)
                                                                  aacagcgccagcgacatacggtcgataccgcgagtgatggtggacattaaaatccagatcagccagGCGCAAGCTCAGTCAGGAACAGGGCGATACCGAAACTCACCGGGACGGCGATCAGCAGCGCGATAAACGAAGTCAcc                                                 <  1:495887‑M2/77‑1 (MQ=255)
                                                                  aacagcgccagcgacataccgtcgataccgcgagtgatggtggacattaaactccagatcagccagGCGCAAGCTCAGTCAGGAACAGGGCGATACCGAAACTCACCGGGCCGGcgctccgcagcgcgaacaccgcaatcacc                                                 >  2:477952‑M2/67‑129 (MQ=255)
                                                                    gcacgccagcgccttagcgtcgctaccgggggtgatggtggacattaaaatccagatcagccagGCGCAAGGTCAGTCAGGAACAGGGCTATACCGAAACTCACCGGGACGGCGATCAGCAGCGCGATAAACGAAGTCACCaa                                               <  2:762033‑M2/79‑1 (MQ=255)
                                                                      gcgccagcgacagaccgtcgataccgcgagtgatggtggacattaaaatccagatcagccagGCGCAAGCTCAGTCAGGAACAGGGCGATACCGAAACTCACCGGGACGGCGATCAGCAGCGCGATAAACGAAGTCACCAACg                                             <  2:34100‑M2/81‑1 (MQ=255)
                                                                          cagcgacataccgtcgataccgcgagtgatggtggacattaaaatccagatcagccagGCGCAAGCTCAGTCAGGAACAGGGCGATACCGAAACTCACCGGGACGGCGATCAGCAGCGCGATAAACGAAGTCACCAACGTAcc                                         >  1:831780‑M2/59‑143 (MQ=255)
                                                                          aagcgacataccgtcgatacagcgagggatggtggacagtaaaatccagatcagccagGCGCAAGCTCAGTCAGGAACAGGGCGATACCGAAACTCACCGGGACGGCGATCAGCAGCGCGATAAACGAAGTCACCAACGTAcc                                         <  1:810648‑M2/85‑1 (MQ=255)
                                                                             cgacataccgtcgataccgcgagtgatggtggacattaaaatccagatcagccagGCGCAAGCTCAGTCAGGAACAGGGCGATACCGAAACTCACGGGGACGGCGATCAGCGGCGAGACAAAAGAAGACACCAACGTACGGTa                                      >  2:366751‑M2/56‑143 (MQ=255)
                                                                                 ataccgtcgataccgcgagtgatggtggacattaaaatccagatcagccagGCGCAAGCTCAGTCAGGAACCGGTCGGTACCGAAACTCACCGGGACGACGATCAGCAGAGGGCTCAACGACGTCACCAACGTACCGTagagc                                  >  1:563231‑M2/52‑141 (MQ=255)
                                                                                  gaccgtcgagaccgcgagtgatggtggacaagaaaatcgagatcagccagGCGCAAGCTCAGTCAGGAACAGGGCGATACCGAAACTCACCGGGACGGCGATCAGCAGCGCGATAAACGAAGTCACCAACGTACCGTAGATCg                                 <  1:294649‑M2/93‑1 (MQ=255)
                                                                                       tcgataccgcgagtgatggtggacattaaaatccagatcagccagGCGCAAGCTCAGTCAGGAACAGGGCGATACCGAAACTCACCGGGACGGCGATCAGCAGCGCGATAAACGAAGTCACCAACGTACCGTAGAGCGGCAcc                            >  1:542569‑M2/46‑143 (MQ=255)
                                                                                          ataccgcgagtgatggtggacattaaaatccagatcagccagGCGCAAGCTCAGTCAGGAACAGGGCGATACCGAAACTCACCGGGACGGCGATCAGCAGCGCGATAAACGAAGTCACCAACGTACCGTAGATCGGCACCAgc                         >  1:534311‑M2/43‑143 (MQ=255)
                                                                                          ataccgcgagtgatggtggacattaaaatccagatcagccagGCGCAAGCTCAGTCAGGAACAGGGCGATACCGAAACTCACCGGGACGGCGATCAACAGCGCGACAAACGAAGTCACCAACGTACGCTAGCTCGGCACCAgc                         >  2:416693‑M2/43‑143 (MQ=255)
                                                                                                  agggctgttgaacattaaaggccagagcagccagGCGCAAGCTCAGTCAGGATCAGGGCGATACCGAAACTCGCCGGGACGGCGATCAGCAGCGCGATAAACGAAGTCACCAACGTACCGTAGATCGGCACCAGCGCCCCGTa                 <  1:698106‑M2/109‑1 (MQ=255)
                                                                                                   gagatggtcgacactaaaagcccggcaggccagGCGCAAGCGGAGTCAGGAACAGGGCGAGACCGAAACTCACCGGGACGGCGATCAGCAGCGCGATAAACGAAGTCACCAACGTACCGTAGATCGGCACCAGCGCCCCGTAg                <  2:135688‑M2/110‑1 (MQ=255)
                                                                                                         gtggacattaaaatccagatcagccagGCGCAAGCTCAGTCAGGAACAGGGCGATACCGAAACTCACCGGGACGGCGATCAGCAGCGCGATAAACGAAGTCACCAACGTACCGTAGATCGGCACCAGCGCCCCGTAGATATCg          >  2:719445‑M2/28‑143 (MQ=255)
                                                                                                         gtggacattaaaatccagatcagccagGCGCAAGCTCAGTCAGGAACAGGGCGATACCGAAACTCACCGGGACGGCGATCAGCAGCGCGACAAACGCAGTCACCAACGTACCGTAGCTCGGCACCAGCGCCCCGGAGATATCg          >  2:719451‑M2/28‑143 (MQ=255)
                                                                                                               attaaactccagatcagccagGCGCAAGCTCAGTCAGGAACAGGGCGATACCGAAGCTCACCGAGCAGGCTGTGTTCAGCGCGAGAAACGAAGTCACCAACGTACCGTAGATCGGCACCAGCGCCCCGTAGATATCGTTCGGt    <  1:903679‑M2/122‑1 (MQ=255)
                                                                                                               atcatgacccacatcagccagGCGCAAGCTCAGTCAGGAACAGGGCGATACCGAAACTCACCGGGACGGCGATCAGCAGCGCGATAAACGAAGTCACCAACGTACCGTAGATCGGCACCAGCGCCCCGTAGATATCGTTCGGt    <  1:159888‑M2/122‑1 (MQ=255)
                                                                                                                 taaaatccagatcagccagGCGCAAGCTCAGTCAGGAACAGGGCGATACCGAAACTCACCGGGACGGCGATCAGCAGCGCGATAAACGAAGTCACCAACGTACCGTAGATCGACACCAGCGCCCCGTAGCTATCGTGCGgcgc  >  1:587969‑M2/20‑143 (MQ=255)
                                                                                                                                    |                                                                                                                           
TGCCGACCGGCTCCTGAAAGTAAACGGCGAACAGCGGCGCAAAGATAAACAGGCCCCACATGCCGTAAACGATACTTGGAATGGCTGCCAGCAGCTCAATGGCGATACCCAGCGGGCGTTTCAGCCAGCCAGGCGCAAGCTCAGTCAGGAACAGGGCGATACCGAAACTCACCGGGACGGCGATCAGCAGCGCGATAAACGAAGTCACCAACGTACCGTAGATCGGCACCAGCGCCCCGTAGATATCGTTCGGTGC  >  NC_000913/3909953‑3910208

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.