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 24 [1] [0] 24 phoA alkaline phosphatase

New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913 = 4023901 (0.060)24 (1.410) 21/282 0.0 98.0% coding (644/1416 nt) phoA alkaline phosphatase
?NC_000913 403114 = 0 (0.000)coding (1368/1416 nt) phoA alkaline phosphatase

TTCGCTAAAAAGGAGGGTAACACGCTGGTCATAGTCACCGCTGATCACGCCCACGCCAGCCAGATTGTTGCGCCGGATACCAAAGCTCCGGGCCTCACCCAGGCGCTAAATACCAAAGATGGCGCAGTGATGGTGATGAGTTACGGGAACTCCGAAGAGGATTCACAAGAACATACCGGCAGTCAGTTGCGTATTGCGGCGTATGGCCCGCATGCCGCCAATGTTGTTGGACTGACCGACCAGACCGATCTCTTCTACACCATGAAAGCCGCTCTGGGGCTGAAATAAAACCGCGCCCGGCAGTGAATTTTCGCTGCCGGGTGGTTTTTTTGCTGTTAGCAACCAGACTTAATGGCAGATCACGGGCGCATACGCTCATGGTTAAAACATGAAGAGGGATGGTGCTATGAAAATAACATTACTGGTTACCTTGCTTTTCGGTCTGGTTTTTTTAAC  >  NC_000913/402875‑403330
                                                                                                                                                                                                                                               |                                                                                                                                                                                                                        
aaaacctataacggcgcgctgggcgtcgatattcacgaaaacgatcaccccacgattctggaaatggcaaaagccgcaggtctggcgaccggtaacgtttctaccgcagagttgcaggatgccacgcccgctgcgctggtggcacatgtgacctcgcgcaaatgctacggtccgagcgcgaccagtgaaaaatgtccgggtaacgctctggaaaaaggcggaaaaggatcgattaccgaCCAGACCGAtct                                                                                                                                                                                                               >  1:177318‑M2/240‑251 (MQ=255)
                                        aagatcacccaactattctggaaagggcaaaagccgcaggtctggcgaccggtaacgtttctaccgcagagttgcaggatgccacgcccgctgcgctggttgcacatgggacctcgcgcaaatgcgtcggtccgagcgcgaccagtgaaaaatgtcctggtatctctctgtataaaggcggaaaaggatcgatcaccgaCCAGACCGATCTCTTCTACACCATGAAAGCCGCTCTGGGGCTGAAATAAAAc                                                                                                                                                                       <  1:423153‑M2/52‑1 (MQ=255)
                                                            naaatggcaaaagccgcaggtctggcgaccggtaacgtttctaccgcagagttgcaggatgccacgcccgcagcgctggtggcacatgtgacctcgcgcaaatgctacggtccgagcgcgaccagtgaaaaatgtccgggtaacgctctggaaaaaggcggaaaaggatcgattaccgaCCAGACCGATCTCTTCTACACCATGAAAGCCGCTCTGGGGCTGAAATAAAACCGCGCCCGGCAGTGAAtttt                                                                                                                                                   >  1:281328‑M2/180‑251 (MQ=255)
                                                             aaatggcaaaagccgcaggtctggcgaccggtaacgtttctaccgcagagttgcaggatgccacgcccgctgcgctggtggcacatgtgacctcgcgcaaatgctacggtccgagcgcgaccagtgaaaaatgtccgggtaacgctctggaaaaaggcggaaaaggatcgattaccgaCCAGACCGATCTCTTCTACACCATGAAAGCCGCTCTGGGGCTGAAATAAAACCGCGCCCGGCAGTGAATTTTc                                                                                                                                                  <  1:106852‑M2/73‑1 (MQ=255)
                                                                  gcaaaagccgcaggtctggcgaccggtaacgtttctaccgcagagttgcaggatgccaagcccgctgcgctggtggcacatgtgacctcgcgcaaatgctacggtccgagcgcgaccagtgaaaaatgtccgggtaacgctctggaaaaaggcggaaaaggatcgattaccgaCCAGACCGATCTCTTCTACACCATGAAAGCCGCTCTGGGGCTGAAATAAAACCGCGCCCGGCAGTGAATTTTCGCTGc                                                                                                                                             >  1:304672‑M2/174‑251 (MQ=255)
                                                                       agccgcaggtctggcgaccggtaacgtttataccgccgagttgcaggatgccacgcacgctgcgctggtggcacatgtgacctcgcgcaaatgctacggtccgagcgcgaccagtgaaaaatgtccgggtaacactctggaaaacggcggaaaaggatcgattaccgcCCAGACCGATCTCTTCTACACAATGAAAGCCGCTCTGGGGCTGAAATAAAACAGCGCCCGGCAGTGAATTTTCGCTGCCGggt                                                                                                                                        >  1:408869‑M2/169‑251 (MQ=255)
                                                                                      gaccggtaacgtttctaccgcagagttgcaggatgccacgcccgctgcgctggtggaacatgtgacctcgcgcaaatgctccggtccgagcgcgaccagtgaaaaatgtccgggtaacgctctggaaaaaggcggaaaaggatcgattaccgaCCAGACCGATCTCTTCTACACCATGAAAGCCGCTCTGGGGCTGAAATAACACCGCGCCCGGCAGTGAATTTTCGCTGCCGGGTGGTTTTTTTGCTGtt                                                                                                                         >  1:77110‑M2/154‑251 (MQ=255)
                                                                                                    ctaccgcagagttgcaggatgccacgcccgctgcgctggtggcacatgtgacctcgcgcaaatgctacggtacgagcgcgaccagtgaaaaatgtccgggtaacgctctggaaaaaggcggaaaaggatcgattaccgaCCAGACCGATCTCTTCTACACCATGAAAGCCGCTCTGGGGCTGAAATAAAACCGCGCCCGGCAGTGAATTTTCGCTGCCGGGTGGTTTTTTTGCTGTTAGCAACCAGACtt                                                                                                            >  1:234782‑M2/140‑250 (MQ=255)
                                                                                                              gttgcaggatgccacgcccgctgcgctggtggcacatgtgacctcgcgcaaatgctacggtccgagcgcgaccagtgaaaaatgtccgggtaacgctctggaaaaaggcggaaaaggatcgattaccgaCCAGACCGATCTCTTCTACACCATg                                                                                                                                                                                                  <  1:285946‑M2/25‑1 (MQ=255)
                                                                                                                 nnnnnnncncanncccgnggcgctggtggcacatgtgacctcncgcaaatgctacggtccgagcgcgaccagtgnanaatgtccgggtaacgctctggaaaaaggcggaaaaggatcgattaccgnCCAGACCGATNTCTNNTACANCNTNAAAGCCGCTCTGGGGCTGAAATAAAACCGNGCCCGGCAGTGAATTTTCGCTGCNNGGNNNTTNNTTTGc                                                                                                                             <  2:326822‑M2/94‑1 (MQ=255)
                                                                                                                     gatgccacgcccactgcgctggtggcacatgggacctcgcgcaaatgctacggtccgagcgcgacccgtgcaaaatgtccgggtaacgctctggaaaaaggcggaaaaggatcgattaccgaCCAGACCGATCTCTTCTACACCATGAAAGCCGCTCTGGGGCTGAAATAAAACcgcg                                                                                                                                                                   >  1:195072‑M2/123‑178 (MQ=255)
                                                                                                                         ccacgcccgctgcgctggtggcacatgtgacctcgcgcacatgctacggtccgagcgcgaccagtgaacaatgtccgggtaacgctctggaaaaaggcggaaaaggatcgattaccgaCCAGACCGATCTCTTCTaca                                                                                                                                                                                                       >  1:238682‑M2/119‑138 (MQ=255)
                                                                                                                         ccacgcccgctgcgctggtggcacatgtgacctcgcgcaaatgctacggtccgagcgcgaccagtgaaaaatgtccgggtaacgctctggaaaaaggcggaaaaggatcgattaccgaCCAGACCGATCTCTTCTACACCATGAAAGCCGCTCTGGGGCTGAAATAAAACCGCGCCCGGCAGTGAATTTTCGCTGCCGGGTGGTTTTTTTGCTGTTAGCAACCAGACTTAATGGCAGATCACGGGCGCAt                                                                                       >  1:120515‑M2/119‑250 (MQ=255)
                                                                                                                          cacgcccgctgcgctggtggcacatgtgacctcgcgcaaatgctacggtccgagcgcgaccagtgaaaaatgtccgggtaacgctctggaaaaaggcggaaaaggatcgattaccgaCCAGACCGATCTCTTCTACACCATGAAAGCCGCTCTGGGGCTGAAATAAAACCGCGCCCGGCAGTGAATTTTCGCTGCCGGGTGGTTTTTTTGc                                                                                                                             >  1:326822‑M2/118‑211 (MQ=255)
                                                                                                                             ncccgctgcgctggtggcacatgtgacctcgcgcaaatgctacggtccgagcgcgaccagtgaaaaatgtccgggtaacgctctggaaaaaggcggaaaaggatcgattacagaCCAGACCGATCTCTTCTACACCATGAAAGCCGCTCTGGGGCTGAAATAAAACCGCGCCCGGCAGTGAATTTTCGCTGCCGGGTGGTTTTTTTGCTGTTAGCAACCAGACTTAATGGCAGATCACGGGCGCATACGCt                                                                                  >  1:49747‑M2/115‑251 (MQ=255)
                                                                                                                               cngctgcgctggtggcacatgtgacctcgcgcaaatgctacggtccgagcgcgaccagtgaaaaatgtccgggtaacgctctggaaaaaggcggaaaaggatcgattaccgaCCAGACCGATCTCTTCTACACCATGAAAGCCGCTCTGGGGCTGAAATAAAAc                                                                                                                                                                       >  1:349306‑M2/113‑164 (MQ=255)
                                                                                                                                         ggtggcacatgtgacctcgcgcaaatgctacggtccgagcgcgaccagtgaaaaatgtacgggtaacgctcgggaaaaaggcggaaaaggatcgattaccgaCCAGACCGAGCTCTTCTACACCATGAAAGCCGCTCTGGGGCTGAAATAAAACCGCGCCCGGCAGGGAATTTTCGCTTCCGGGTGGGTTTTTTGCTGTTAGCAACCAGACTTAATGGCAGATCACGGGCGCATACGCTCATGGTTAAAAc                                                                      <  1:234166‑M2/149‑1 (MQ=255)
                                                                                                                                              cacatgtgacctcgcgcaaatgctacggtccgagcgcgaccagtgaaaaatgtccgggtaacgctctggaaaaaggcggaaaaggatcgattaccgaCCAGACCGATCTCTTCTac                                                                                                                                                                                                        >  1:313937‑M2/98‑116 (MQ=255)
                                                                                                                                                     gacctcgcgaaaatgctacggtccgancgcgacnngtnaaaaatgtccgggtaacgcnctnnnannngncggaanagnntcgctcaccnaCCAGACCGATCTCTTCTACACCATGAAAGCCGCTCTGGGGCTGACATAAACCNGCGACCGGCNNTGAATTTTCGCTGCCGGGTGGNTTTNNtnntntttntan                                                                                                                    >  2:103909‑M2/91‑188 (MQ=255)
                                                                                                                                                                           ccgatcgcgaccagtgaaaaatgtccgtgtaacgctcttgaaaaaggcggaaaaggatcgattaccgaCCAGACCGATCTCTTCTCCACCATGAAAGc                                                                                                                                                                                             <  1:321937‑M2/30‑1 (MQ=255)
                                                                                                                                                                           ccgagcgcgaccagtgaaaaatgtccgggtaacgctctggaaaaaggcggaaaaggatcgattaccgaCCAGACCGATCTCTTCTACACCATGAAAGCCGCTCTGGGGCTGAAATAAAACCGCGCCCGGCAGTGAATTTTCGCTGCCGGGTGGTTTTTTTGCTGTTAGCAACCAGACTTAATGGCAGATCACGGGCGCATACGCTc                                                                                 >  1:307195‑M2/69‑206 (MQ=255)
                                                                                                                                                                                                   ccgggtaacgctctggaaaaaggcggaaaaggatcgattaccgaCCAGACCGATCTCTTCTACACCATGAAAGCCGCTCTGGGGCTGAAATAAAACCGCGCCCGGCAGTGAATTTTCGCTGCCGGGTGGTTTTTTTGCTGTTAGCAACCAGACTTAATGGCAGATCACGGGCGCATACGCTCATGGTTAAAACATGAAGAGGGATGATGCTATGAAAATAACATTACTGGTTGCCTTGCTTTTCGGTCTgg            <  1:356215‑M2/207‑1 (MQ=255)
                                                                                                                                                                                                     gggtaacgctcgggaaaaaggcggaaaatgatcgatgaccgaCCAGACCGATCTCTTCTACACCATGAAAGCCGCTCTGGGGCTGAAATAAAACCGCGCCCGGCAGTGAATTTTCGCTGCCGGGTGGTTTTTTTGCTGTTAg                                                                                                                       <  1:188694‑M2/100‑1 (MQ=255)
                                                                                                                                                                                                               ctggaaaaaggcggaaaaggatcgattaccgaCCAGACCGATCTCTTCTACACCAGGAAAGCCGCTCTGGGGCTGAAATAAAACCGCGCCCGGCAGTGAATTTTCGCTGCCGGGTGGTTTTTTTGCTGTTAGCAACCAGACTTAATGGCAGATCACGGGCGCATACGCTCATGGTTAAAACATGAAGAGGGATGGTGCTATGAAAATAACATTACTGGTTACCGGGCTTTTCGGTCTGGTGTTTTTAac  <  1:103909‑M2/217‑1 (MQ=255)
                                                                                                                                                                                                                                               |                                                                                                                                                                                                                        
TTCGCTAAAAAGGAGGGTAACACGCTGGTCATAGTCACCGCTGATCACGCCCACGCCAGCCAGATTGTTGCGCCGGATACCAAAGCTCCGGGCCTCACCCAGGCGCTAAATACCAAAGATGGCGCAGTGATGGTGATGAGTTACGGGAACTCCGAAGAGGATTCACAAGAACATACCGGCAGTCAGTTGCGTATTGCGGCGTATGGCCCGCATGCCGCCAATGTTGTTGGACTGACCGACCAGACCGATCTCTTCTACACCATGAAAGCCGCTCTGGGGCTGAAATAAAACCGCGCCCGGCAGTGAATTTTCGCTGCCGGGTGGTTTTTTTGCTGTTAGCAACCAGACTTAATGGCAGATCACGGGCGCATACGCTCATGGTTAAAACATGAAGAGGGATGGTGCTATGAAAATAACATTACTGGTTACCTTGCTTTTCGGTCTGGTTTTTTTAAC  >  NC_000913/402875‑403330

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 13 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 38 ≤ 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.