Predicted mutation
evidence seq id position mutation annotation gene description
MC JC NC_000913 402,391 Δ723 bp coding (645‑1367/1416 nt) phoA → alkaline phosphatase

Missing coverage evidence...
   seq id start end size ←reads reads→ gene description
* * ÷ NC_000913 402391 403113 723 34 [0] [0] 34 phoA alkaline phosphatase

New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913 = 4023900 (0.000)34 (0.620) 31/482 0.7 100% coding (644/1416 nt) phoA alkaline phosphatase
?NC_000913 403114 = 0 (0.000)coding (1368/1416 nt) phoA alkaline phosphatase

TGGAATTCGCTAAAAAGGAGGGTAACACGCTGGTCATAGTCACCGCTGATCACGCCCACGCCAGCCAGATTGTTGCGCCGGATACCAAAGCTCCGGGCCTCACCCAGGCGCTAAATACCAAAGATGGCGCAGTGATGGTGATGAGTTACGGGAACTCCGAAGAGGATTCACAAGAACATACCGGCAGTCAGTTGCGTATTGCGGCGTATGGCCCGCATGCCGCCAATGTTGTTGGACTGACCGACCAGACCGATCTCTTCTACACCATGAAAGCCGCTCTGGGGCTGAAATAAAACCGCGCCCGGCAGTGAATTTTCGCTGCCGGGTGGTTTTTTTGCTGTTAGCAACCAGACTTAATGGCAGATCACGGGCGCATACGCTCATGGTTAAAACATGAAGAGGGATGGTGCTATGAAAATAACATTACTGGTTACCTTGCTTTTCGGTCTGGTTTTTTTAACCACCGTCGGCGCTGCCGAGAGAAC  >  NC_000913/402870‑403354
                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                                
gtgtcaaaacctataacggcgcgctgggcgtcgatattcacgaaaaagatcacccaacgattctggaaatggcaaaagccgcaggtctggcgaccggtaacgtttctaccgcagagttgcaggatgccacgcccgctgcgctggtggcacatgtgacctcgcgcaaatgctacggtccgagcgcgaccagtgaaaaatgtccgggtaacgctctggaaaaaggcggaaaaggatcgattaccgaCCAGAcc                                                                                                                                                                                                                                            <  2:152121‑M2/7‑1 (MQ=255)
    caaatcctataacggcgcgctggtcgtcgatattcacgaaaaagatctcccaacgattctggaaatggcaaaagccgcaggtctggcgaccggtaacgtttcttccgcagagttgcaggatgccacgcccgctgcgctggtggcacatgtgacctcgcgcaaatgctacggtccgagcgcgtccagtgaaaaatgtccgggtaacgctctggaaaaaggcggaaaaggatcgattaccgaCCAGACCGAtc                                                                                                                                                                                                                                        <  1:82920‑M2/11‑1 (MQ=255)
            ataacggcgcgctgggcgtcgatattcacgaaaaagatcacccaacgattctggaaatggcaaaagccgcaggtctggcgaccggtaacgtttctaccgcagagttgcaggatgccacgcccgctgcgctggtggcacatgtgacctcgcgcaaatgctacggtccgagcgcgaccagtgaaaaatgtccgggtaacgctctggaaaaaggcggaaaaggatcgattaccgaCCAGACCGATCTCTTCTac                                                                                                                                                                                                                                >  1:191841‑M2/233‑251 (MQ=255)
                         gggcgtcgatattcacgaaaaagatcacccaacgattctggaaatggcaaaagccgcaggtctggcgaccggtaacgtttctaccgcagagttgcaggatgccacgcccgctgcgctggtggcacatgtgacctcgcgcaaatgctacggtccgagcgcgaccagtgaaaaatgtccgggtaacgctctggaaaaaggcggaaaaggatcgattaccgaCCAGACCGATCTCTTCTACACCATGAAAGc                                                                                                                                                                                                                     <  1:375258‑M2/30‑1 (MQ=255)
                             gtcgatattcacgaaaaagatcacccaacgattctggaaatggcaaaagccgcaggtctggcgaccggtaacgtttctaccgcagagttgcaggatgccacgcccgctgcgctggtggcacatgtgacctcgcgcaaatgctacggtccgagcgcgaccagtgaaaaatgtccgggtaacgctctggaaaaaggcggaaaaggatcgattaccgaCCAGACCGATCTCTTCTACACCATGAAAGCCGctct                                                                                                                                                                                                               >  1:179416‑M2/216‑251 (MQ=255)
                                 atattcacgaaaaagatcacccaacgattctggaaatggcaaaagccgcaggtctggcgaccggtaacgtttctaccgcagagttgcaggatgccacgcccgctgcgctggtggcacatgtgacctcgcgcaaatgctacggtccgagcgcgaccagtgaaaaatgtccgggtaacgctctggaaaaaggcggaaaaggatcgattaccgaCCAGACCGATCTCTTCTACACCATGAAAGCCGCTCTgggg                                                                                                                                                                                                           >  1:12940‑M2/212‑251 (MQ=255)
                                 atattcacgaaaaagatcacccaacgattctggaaatggcaaaagccgcaggtctggcgaccggtaacgtttctaccgcagagttgcaggatgccacgcccgctgcgctggtggcacatgtgacctcgcgcaaatgctacggtccgagcgcgaccagtgaaaaatgtccgggtaacgctctggaaaaaggcggaaaaggatcgattaccgaCCAGACCGATCTCTTCTACACCATGAAAGCCGCTCTgggg                                                                                                                                                                                                           >  1:519468‑M2/212‑251 (MQ=255)
                                      cacgaaaaagatcccccaacgattctggaaatggcaaaagccgcaggtctggcgaccggtaacgtttctaccgcagagttgcaggatgccacgcccgctgcgctggtggcacatgtgacctcgcgcaaatgctacggtccgagcgcgaccagtgaaaaatgtccgggtaacgctctggaaaaaggcggaaaaggatcgattaccgaCCAGACCGATCTCTTCTACACCATGAAAGCCGCTCTGGGGCTGa                                                                                                                                                                                                       <  2:376745‑M2/44‑1 (MQ=255)
                                           aaaagatcacccaacgattctggaaatggcacaagccgcaggtctggcgaccggtaacgtttctaccgcagagttgcaggatgccacgcccgctgcgctggtggcacatgtgacctcgcgcaaatgctacggtccgagcgcgaccagtgaaaaatgtccgggtaacgctctggaaaaaggcggaaaaggatcgattaccgaCCAGCCCGATCTCTTCTACACCATGAAAGCCGCTCTGGGGCTGAAATaa                                                                                                                                                                                                  >  2:338477‑M2/202‑250 (MQ=255)
                                                         cgattctggaaatggcaaaagccgcaggtctggcgaccggtaacgtttctaccgcagagttgcaggatgccacgcccgctgcgctggtggcacatgtgacctcgcgcaaatgctacggtccgagcgcgaccagtgaaaaatgtccgggtaacgctctggaaaaaggcggaaaaggatcgattaccgaCCAGACCGATCTCTTCTACACCATGAAAGCCGCTCTGGGGCTGAAATAAAACCGCGCCCGGCAg                                                                                                                                                                                   <  2:314315‑M2/64‑1 (MQ=255)
                                                                         aaaagccgcaggtctggcgaccggtaacgtttctaccgcagagttgcaggatgccacgcccgctgcgctggtggcacatgtgacctcgcgcaaatgctacggtccgagcgcgaccagtgaaaaatgtccgggtaacgctctggaaaaaggcggaaaaggatcgattaccgaCCAGACCGATCTCTTCTACACCATGAAAGCCGCTCTGGGGCTGAAATAAAACCGCGCCCGGCAGTGAATTTTCGCTGCCg                                                                                                                                                                   >  1:486081‑M2/172‑251 (MQ=255)
                                                                                        ggcgaccggtaacgtttctaccgcagagttgcaggatgccacgcccgctgcgctggtggcacatgtgacctcgcgcaaatgctacggtccgagcgcgaccagtgaaaaatgtccgggtaacgctctggaaaaaggcggaaaaggatcgattaccgaCCAGACCGATCGCGTCTACACCATGAAAGCCGCTCTGGGGCTGAAATAAAACCGCGCCCGGCAGTGAATTTTCGCGGCCGGGGGGTTTTTTTg                                                                                                                                                      <  2:83442‑M2/93‑1 (MQ=255)
                                                                                                                      gcaggatgccacgcccgctgcgctggtggcacatgtgacctcgcgcaaatgctacggtccgagcgcgaccagtgaaaaatgtccgggtaacgctctggaaaaaggcggaaaaggatcgattaccgaCCAGACCGATCTCTTCTACACCATGAAAGCCGCTCTGGGGCTGAAATAAAACCGCGCCCGGCAGTGAATTTTCGCTGCCGGGTGGTTTTTTTGCTGTTAGCAACCAGACTTAATGGCAGATCACg                                                                                                                      <  1:555918‑M2/125‑1 (MQ=255)
                                                                                                                       caggatgccacgcccgctgcgctggtggcacatgtgacctcgcgcaaatgctacggtccgagcgcgaccagtgaaaaatgtccgggtaacgctctggaaaaaggcggaaaaggatcgattaccgaCCAGACCGATCTCTTCTACACCATGAAAGCCGCTCTGGGGCTGAAATAAAACCGCGCCCGGCAGTGAATTTTCGCTGCCGGGTGGTTTTTTTGCTGTTAGCAACCAGACTTAATGGCAGATCACgg                                                                                                                     >  1:191576‑M2/126‑251 (MQ=255)
                                                                                                                                 cgcccgctgcgctggtggcacatgtgacctcgcgcaaatgctacggtccgagcgcgaccagtgaaaaatgtccgggtaacgctctggaaaaaggcggaaaaggatcgattaccgaCCAGACCGATCTCTTCTACACCATGAAAGCCGCTCTGGGGCTGAAATAAAACCGCGCCCGGCAGTGAATTTTCGCTGCCGGGTGGTTTTTTTGCTGTTAGCAACCAGACTTAATGGCAGATCACGGGCGCATACGc                                                                                                           <  1:181545‑M2/136‑1 (MQ=255)
                                                                                                                                          cgctggtggcacatgtgaccgcgcgcaaatgctacggtccgagcgcgaccagtgaaaaatgtccgggtaacgctctggaaaaaggcggaaaaggatcgattaccgaCCAGACCGATCTCTTCTACACCATGAAAGCCGCTCTGGGGCTGAAATAAAACCGCGCCCGGCAGTGAATTTTCGCTGCCGGGTGGTTTTTTTGCTGTTAGCAACCAGACTTAATGGCAGATCACGGGCGCATACGCTCATGGTTa                                                                                                  <  2:114681‑M2/145‑1 (MQ=255)
                                                                                                                                               gtggcacatgtgacctcgcgcaaatgctacggtccgagcgcgaccagtgaaaaatgtccgggtaacgctctggaaaaaggcggaaaaggatcgattaccgaCCAGACCGATCTCTTCTACACCATGAAAGCCGCTCTGGGGCTGAAATAAAACCGCGCCCGGCAGTGAATTTTCGCTGCCGGGTGGTTTTTTTGCTGTTAGCAACCAGACTTAATGGCAGATCACGGGCGCATACGCTCATGGTTAAAAc                                                                                              >  1:106872‑M2/102‑250 (MQ=255)
                                                                                                                                                           acctcgcgcaaatgctacggtccgagcgcgaccagtgaaaaatgtccgggtaacgctctggaaaaaggcggaaaaggatcgattaccgaCCAGACCGATCTCTTCTACACCATGAAAGCCGCTCTGGGGCTGAAATAAAACCGCGCCAGGCAGTGAATTTTCGCTGCCGGGTGGTTTTTTTGCTGTTAGCAACCAGACTTAATGGCAGATCACGGGCGCATACGCTCATGGTTAAAACTTGAAGAGGGATg                                                                                 >  2:491037‑M2/90‑251 (MQ=255)
                                                                                                                                                               cgcgcaaatgctacggtccgagcgcgaccagtgaaaaatgtccgggtaacgctctggaaaaaggcggaaaaggatcgattaccgaCCAGACCGATCTCTTCTACACCATGAAAGCCGCTCTGGGGCTGAAATAAAACCGCGCCCGGCAGTGAATTTTCGCTGCCGGGTGGTTTTTTTGCTGTTAGCAACCAGACTTAATGGCAGATCACGGGCGCATACGCTCATGGTTAAAACATGAAGAGGGATGGTg                                                                              <  2:15388‑M2/165‑1 (MQ=255)
                                                                                                                                                                          tacggtccgagcgcgaccagtgaaaaatggccgggtaacgctctggaaaaaggcagaaaaggatcgattaccgaCCAGACCGATCTCTTCTACACCAGGAAAGCCGCTCTGGGGCTGAAATAAAACCGCGCCCGGCAGTGAATTTTCGCTGCCGGGTGGTTTTTTTGCTGTTAGCAACCAGACTTAATGGCAGATCACGGGCGCATACGCTCATGGTTAAAACATGAAGAGGGATGGTTCTATGAAAATa                                                                   <  2:59316‑M2/176‑1 (MQ=255)
                                                                                                                                                                                          ccagtgaaaaatgtccgggtaacgctctggaaaaaggcggaaaaggatcgattaccgaCCAGACCGATCTCTTCTACACCATGAAAGCCGCTCTGGGGCTGAAATAAAACCGCGCCCGGCAGTGAATTTTCGCTGCCGGGTGGTTTTTTTGCTGTTAGCAACCAGACTTAATGGCAGATCACGGGCGCATACGCTCATGGTTAAAACATGAAGAGGGATGGTGCTATGAAAATAACATTACTGGTTACCtt                                                  >  1:73417‑M2/59‑251 (MQ=255)
                                                                                                                                                                                             gggaaaaatggccgggtaacgctctggaaaaaggcggaaaaggatcgattaccgaCCAGACCGATCTCTTCTACACCATGAAAGCCGCTCTGGGGCTGAAATAAAACCGCGCCCGGCAGTGAATTTTCGCTGCCGGGTGGTTTTTTTGCTGTTAGCAACCAGACTTAATGGCAGATCACGGGCGCATACGCTCATGGTTAAAACAGGAAGAGGGATGGTGCTATGAAAATAACATTACTGGTTACCTTGCt                                               <  1:335791‑M2/196‑1 (MQ=255)
                                                                                                                                                                                               gaaaaatgtccgggtaacgctctggaaaaaggcggaaaaggatcgattaccgaCCAGACCGATCTCTTCTACACCATGAAAGCCGCTCTGGGGCTGAAATAAAACCGCGCCCGGCAGTGAATTTTCGCTGCCGGGTGGTTTTTTTGCTGTTAGCAACCAGACTTAATGGCAGATCACGGGCGCATACGCTCATGGTTAAAACATGAAGAGGGATGGTGCTATGAAAATAACATTACTGGTTACCTTGCtt                                              <  2:558091‑M2/197‑1 (MQ=255)
                                                                                                                                                                                                 aaaatgtccgggtcacgatctggacaaaggcggcaaaggatagattacagaCCAGACCGATCTCTTCTACACCATGAAAGCCGCTCTGGGGCTGAAATAAAACCGCGCCCGGCAGTGAATTTTCGCTGCCGGGTGGTTTTTTTGCTGTTAGCAACCAGACTTAATGGCAGATAACGGGCGCATAAGCTCATGGTTAAAACATGACGAGGGATGGTGCTAAGAAAATAACATTACTGGTTACCTTGcttttt                                           >  2:389642‑M2/52‑250 (MQ=255)
                                                                                                                                                                                                  aaatgtccgtgtaacgctctggaaaaaggcggaaaaggatcgattaccgaCCAGACCGATCTCTTCTACACCATGAAAGCCGCGCTGGGGCTGAAATAAAACCGCGCCCGGCAGTGAATGTTCGCTGCCGGGTGGTTTTTTGGCTGTTAGCAACCAGACTTAATGGCAGATCACGGGCGCAGACGCTCATGGTTAAAGCATTAAGAGGGATGGTGCTATGAAAATAACATTACTGGTTACCTTGCTTTTCg                                          <  1:82704‑M2/201‑1 (MQ=255)
                                                                                                                                                                                                        ccgggtaacgctcttgaaaaaggcggaaaaggatcgattaccgaCCAGACCGATCTCTTCTACACCATGAAAGCCGCTCTGGGGCTGAAATAAAACCGCGCCCGGCAGTGAATTTTCGCTGCCGGGTGGTTTTTTTGCTGTTAGCAACCAGACTTAATGGCAGATCACGGGCGCATACGCTCATGGTTAAAACATGAAGAGGGATGGTGCTATGAAAATAACATTACTGGTTACCTTGCTTTTCGGTCTgg                                    <  1:461624‑M2/207‑1 (MQ=255)
                                                                                                                                                                                                         cgggtaacgctctggaaaaaggcggaaaaggatcgattaccgaCCAGACCGATCTCTTCTACACCATGAAAGCCGCTCTGGGGCTGAAATAAAACCGCGCCCGGCAGTGAATTTTCGCTGCCGGGTGGTTTTTTTGCTGTTAGCAACCAGACTTAATGGCAGATCACGGgc                                                                                                                   <  1:498341‑M2/128‑1 (MQ=255)
                                                                                                                                                                                                         cgggtaacgctctggaaaaaggcggaaaaggatcgattaccgaCCAGACCGATCTCTTCTACACCATGAAAGCCGCTCTGGGGCTGAAATAAAACCGCGCCCGGCAGTGAATTTTCGCTGCCGGGTGGTTTTTTTGCTGTTAGCAACCAGACTTAATGGCAGATCACGGgc                                                                                                                   >  2:498341‑M2/44‑171 (MQ=255)
                                                                                                                                                                                                         cgggtaacgctctggaaaaaggcggaaaaggatcgattaccgaCCAGACCGATCTCTTCTACACCATGAAAGCCGCTCTGGGGCTGAAATAAAACCGCGCCCGGCAGTGAATTTTCGCTGCCGGGTGGTTTTTTTGCTGTTAGCAACCAGACTTAATGGCAGATCACGGGCGCATACGCTCATGGTTAAAACATGAAGAGGGATGGTGCTATGAAAATAACATTACTGGTTACCTTGCTTTTCGGTCTGGt                                   >  1:307919‑M2/44‑251 (MQ=255)
                                                                                                                                                                                                             taacgctctggaaaaaggcggaaaaggatcgattaccgaCCAGACCGATCTCTTCTACACCATGAAAGCCGCTCTGGGGCTGAAATAAAACCGCGCCCGGCAGTGAATTTTCGCTGCCGGGTGGTTTTTTTGCTGTTAGCAACCAGACTTAATGGCAGATCACGGGCGCATACGCTCATGGTTAAAACATGAAGAGGGATGGTGCTATGAAAATAACATTACTGGTTACCTTGCTTTTCGGTCTGGttttt                               >  1:37565‑M2/40‑251 (MQ=255)
                                                                                                                                                                                                             taacgctctggaaaaaggcggaaaaggatcgattaccgaCCAGACCGATCTCTTCTACACCATGAAAGCCGCTCTGGGGCTGAAATAAAACCGCGCCCGGCAGTGAATTTTCGCTGCCGGGTGGTTTTTTTGCTGTTAGCAACCAGACTTAATGGCAGATCACGGGCGCATACGCTCATGGTTAAAACATGAAGAGGGATGGTGCTATGAAAATAACATTACTGGTTACCTTGCTTTTCGGTCTGGttttt                               >  2:19152‑M2/40‑251 (MQ=255)
                                                                                                                                                                                                                                       gatcgattaccggCCAGACCGATCTCGTCTACGCCATGAAAGCCGCTCTGGGGCTGAAATAAACCCGCGCCCGGCAGTGAATTTTCGCTGCCGGGTGGTTTTTTTGATGTTAGCAACCAGACTTAATGGCAGATCACGGGCGCATACGCTCATGGTTAAAACATGAAGAGGGATGGTGCTATGAAAATAACATTACTGGTTACCTTGCTTTTCGGTCTGGTTTTTTTAACCACCGTCGGCGCTGCCgagag     <  1:96685‑M2/238‑1 (MQ=255)
                                                                                                                                                                                                                                        atcgattaccgaCCAGACCGATCTCTTCTACACCATGAAAGCCGCTCTGGGGCTGAAATAAAACCGCGCCCGGCAGTGAATTTTCGCTGCCGGGTGGTTTTTTTGCTGTTAGCAACCAGCCTTAATGGCAGATCACGGGCGCATACACTCATGGTTAAAACATGAAGAGGGATGGTGCTATGAAAATAACATTACTGGTTACCTTGCTTTTCGGTCTGGTATTTTTAACCACCGTCGGCGCTGCCgag       >  1:490881‑M2/13‑248 (MQ=255)
                                                                                                                                                                                                                                          cgattaccgaCCAGACCGATCTCTTCTACACCATGAAAGCCGCTCTGGGGCTGAAATAAAACCGCGCCCGGCAGTGAATTTTCGCTGCCGGGTGGTTTTTTTGCTGTTAGCAACCAGACTTAATGGCAGATCACGGGCGCATACGCTCATGGTTAAAACATGAAGAGGGATGGTGCTATGAAAATAACATTACTGGTTACCTTGCTTTTCGGTCTGGTTTTTTTAACCACCGTCGGCGCTGCCGAGAGAAc  <  1:428995‑M2/241‑1 (MQ=255)
                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                                
TGGAATTCGCTAAAAAGGAGGGTAACACGCTGGTCATAGTCACCGCTGATCACGCCCACGCCAGCCAGATTGTTGCGCCGGATACCAAAGCTCCGGGCCTCACCCAGGCGCTAAATACCAAAGATGGCGCAGTGATGGTGATGAGTTACGGGAACTCCGAAGAGGATTCACAAGAACATACCGGCAGTCAGTTGCGTATTGCGGCGTATGGCCCGCATGCCGCCAATGTTGTTGGACTGACCGACCAGACCGATCTCTTCTACACCATGAAAGCCGCTCTGGGGCTGAAATAAAACCGCGCCCGGCAGTGAATTTTCGCTGCCGGGTGGTTTTTTTGCTGTTAGCAACCAGACTTAATGGCAGATCACGGGCGCATACGCTCATGGTTAAAACATGAAGAGGGATGGTGCTATGAAAATAACATTACTGGTTACCTTGCTTTTCGGTCTGGTTTTTTTAACCACCGTCGGCGCTGCCGAGAGAAC  >  NC_000913/402870‑403354

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 24 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 39 ≤ 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.