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

AAAACCTATAACGGCGCGCTGGGCGTCGATATTCACGAAAAAGATCACCCAACGATTCTGGAAATGGCAAAAGCCGCAGGTCTGGCGACCGGTAACGTTTCTACCGCAGAGTTGCAGGATGCCACGCCCGCTGCGCTGGTGGCACATGTGACCTCGCGCAAATGCTACGGTCCGAGCGCGACCAGTGAAAAATGTCCGGGTAACGCTCTGGAAAAAGGCGGAAAAGGATCGATTACCGAACAGCTGCTTAACGCTCGTGCCGACGTTACGCTTGGCGGCGGCGCAAAAACCTTTGCTGAAACGGCAACCGCTGGTGAATGGCAGGGAAAAACGCTGCGTGAACAGGCACAGGCGCGTGGTTATCAGTTGGTGAGCGATGCTGCCTCACTGAATTCGGTGACGGAAGCGAATCAGCAAAAACCCCTGCTTGGCCTGTTTGCTGACGGCAATATGCCA  >  NC_000913/402152‑402607
                                                                                                                                                                                                                                              |                                                                                                                                                                                                                         
aaaaCCTATAACGGCGCGCTGGGCGTCGATATTCACGAAAACGATCACCCCACGATTCTGGAAATGGCAAAAGCCGCAGGTCTGGCGACCGGTAACGTTTCTACCGCAGAGTTGCAGGATGCCACGCCCGCTGCGCTGGTGGCACATGTGACCTCGCGCAAATGCTACGGTCCGAGCGCGACCAGTGAAAAATGTCCGGGTAACGCTCTGGAAAAAGGCGGAAAAGGATCGATTACCGAccagaccgatct                                                                                                                                                                                                               >  1:177318‑M1/1‑239 (MQ=255)
                                        aaGATCACCCAACTATTCTGGAAAGGGCAAAAGCCGCAGGTCTGGCGACCGGTAACGTTTCTACCGCAGAGTTGCAGGATGCCACGCCCGCTGCGCTGGTTGCACATGGGACCTCGCGCAAATGCGTCGGTCCGAGCGCGACCAGTGAAAAATGTCCTGGTATCTCTCTGTATAAAGGCGGAAAAGGATCGATCACCGAccagaccgatctcttctacaccatgaaagccgctctggggctgaaataaaac                                                                                                                                                                       <  1:423153‑M1/251‑53 (MQ=255)
                                                            naaatGGCAAAAGCCGCAGGTCTGGCGACCGGTAACGTTTCTACCGCAGAGTTGCAGGATGCCACGCCCGCAGCGCTGGTGGCACATGTGACCTCGCGCAAATGCTACGGTCCGAGCGCGACCAGTGAAAAATGTCCGGGTAACGCTCTGGAAAAAGGCGGAAAAGGATCGATTACCGAccagaccgatctcttctacaccatgaaagccgctctggggctgaaataaaaccgcgcccggcagtgaatttt                                                                                                                                                   >  1:281328‑M1/2‑179 (MQ=255)
                                                             aaaTGGCAAAAGCCGCAGGTCTGGCGACCGGTAACGTTTCTACCGCAGAGTTGCAGGATGCCACGCCCGCTGCGCTGGTGGCACATGTGACCTCGCGCAAATGCTACGGTCCGAGCGCGACCAGTGAAAAATGTCCGGGTAACGCTCTGGAAAAAGGCGGAAAAGGATCGATTACCGAccagaccgatctcttctacaccatgaaagccgctctggggctgaaataaaaccgcgcccggcagtgaattttc                                                                                                                                                  <  1:106852‑M1/251‑74 (MQ=255)
                                                                  gCAAAAGCCGCAGGTCTGGCGACCGGTAACGTTTCTACCGCAGAGTTGCAGGATGCCAAGCCCGCTGCGCTGGTGGCACATGTGACCTCGCGCAAATGCTACGGTCCGAGCGCGACCAGTGAAAAATGTCCGGGTAACGCTCTGGAAAAAGGCGGAAAAGGATCGATTACCGAccagaccgatctcttctacaccatgaaagccgctctggggctgaaataaaaccgcgcccggcagtgaattttcgctgc                                                                                                                                             >  1:304672‑M1/1‑173 (MQ=255)
                                                                       aGCCGCAGGTCTGGCGACCGGTAACGTTTATACCGCCGAGTTGCAGGATGCCACGCACGCTGCGCTGGTGGCACATGTGACCTCGCGCAAATGCTACGGTCCGAGCGCGACCAGTGAAAAATGTCCGGGTAACACTCTGGAAAACGGCGGAAAAGGATCGATTACCGCccagaccgatctcttctacacaatgaaagccgctctggggctgaaataaaacagcgcccggcagtgaattttcgctgccgggt                                                                                                                                        >  1:408869‑M1/1‑168 (MQ=255)
                                                                                      gACCGGTAACGTTTCTACCGCAGAGTTGCAGGATGCCACGCCCGCTGCGCTGGTGGAACATGTGACCTCGCGCAAATGCTCCGGTCCGAGCGCGACCAGTGAAAAATGTCCGGGTAACGCTCTGGAAAAAGGCGGAAAAGGATCGATTACCGAccagaccgatctcttctacaccatgaaagccgctctggggctgaaataacaccgcgcccggcagtgaattttcgctgccgggtggtttttttgctgtt                                                                                                                         >  1:77110‑M1/1‑153 (MQ=255)
                                                                                                    cTACCGCAGAGTTGCAGGATGCCACGCCCGCTGCGCTGGTGGCACATGTGACCTCGCGCAAATGCTACGGTACGAGCGCGACCAGTGAAAAATGTCCGGGTAACGCTCTGGAAAAAGGCGGAAAAGGATCGATTACCGAccagaccgatctcttctacaccatgaaagccgctctggggctgaaataaaaccgcgcccggcagtgaattttcgctgccgggtggtttttttgctgttagcaaccagactt                                                                                                            >  1:234782‑M1/1‑139 (MQ=255)
                                                                                                              gTTGCAGGATGCCACGCCCGCTGCGCTGGTGGCACATGTGACCTCGCGCAAATGCTACGGTCCGAGCGCGACCAGTGAAAAATGTCCGGGTAACGCTCTGGAAAAAGGCGGAAAAGGATCGATTACCGAccagaccgatctcttctacaccatg                                                                                                                                                                                                  <  1:285946‑M1/154‑26 (MQ=255)
                                                                                                                 nnnnnnncncanncccgngGCGCTGGTGGCACATGTGACCTCNCGCAAATGCTACGGTCCGAGCGCGACCAGTGNANAATGTCCGGGTAACGCTCTGGAAAAAGGCGGAAAAGGATCGATTACCGNccagaccgatntctnntacancntnaaagccgctctggggctgaaataaaaccgngcccggcagtgaattttcgctgcnnggnnnttnntttgc                                                                                                                             <  2:326822‑M1/211‑95 (MQ=255)
                                                                                                                     gATGCCACGCCCACTGCGCTGGTGGCACATGGGACCTCGCGCAAATGCTACGGTCCGAGCGCGACCCGTGCAAAATGTCCGGGTAACGCTCTGGAAAAAGGCGGAAAAGGATCGATTACCGAccagaccgatctcttctacaccatgaaagccgctctggggctgaaataaaaccgcg                                                                                                                                                                   >  1:195072‑M1/1‑122 (MQ=255)
                                                                                                                         ccACGCCCGCTGCGCTGGTGGCACATGTGACCTCGCGCACATGCTACGGTCCGAGCGCGACCAGTGAACAATGTCCGGGTAACGCTCTGGAAAAAGGCGGAAAAGGATCGATTACCGAccagaccgatctcttctaca                                                                                                                                                                                                       >  1:238682‑M1/1‑118 (MQ=255)
                                                                                                                         ccACGCCCGCTGCGCTGGTGGCACATGTGACCTCGCGCAAATGCTACGGTCCGAGCGCGACCAGTGAAAAATGTCCGGGTAACGCTCTGGAAAAAGGCGGAAAAGGATCGATTACCGAccagaccgatctcttctacaccatgaaagccgctctggggctgaaataaaaccgcgcccggcagtgaattttcgctgccgggtggtttttttgctgttagcaaccagacttaatggcagatcacgggcgcat                                                                                       >  1:120515‑M1/1‑118 (MQ=255)
                                                                                                                          cACGCCCGCTGCGCTGGTGGCACATGTGACCTCGCGCAAATGCTACGGTCCGAGCGCGACCAGTGAAAAATGTCCGGGTAACGCTCTGGAAAAAGGCGGAAAAGGATCGATTACCGAccagaccgatctcttctacaccatgaaagccgctctggggctgaaataaaaccgcgcccggcagtgaattttcgctgccgggtggtttttttgc                                                                                                                             >  1:326822‑M1/1‑117 (MQ=255)
                                                                                                                             ncccgCTGCGCTGGTGGCACATGTGACCTCGCGCAAATGCTACGGTCCGAGCGCGACCAGTGAAAAATGTCCGGGTAACGCTCTGGAAAAAGGCGGAAAAGGATCGATTACAGAccagaccgatctcttctacaccatgaaagccgctctggggctgaaataaaaccgcgcccggcagtgaattttcgctgccgggtggtttttttgctgttagcaaccagacttaatggcagatcacgggcgcatacgct                                                                                  >  1:49747‑M1/2‑114 (MQ=255)
                                                                                                                               cnGCTGCGCTGGTGGCACATGTGACCTCGCGCAAATGCTACGGTCCGAGCGCGACCAGTGAAAAATGTCCGGGTAACGCTCTGGAAAAAGGCGGAAAAGGATCGATTACCGAccagaccgatctcttctacaccatgaaagccgctctggggctgaaataaaac                                                                                                                                                                       >  1:349306‑M1/1‑112 (MQ=255)
                                                                                                                                         ggtggCACATGTGACCTCGCGCAAATGCTACGGTCCGAGCGCGACCAGTGAAAAATGTACGGGTAACGCTCGGGAAAAAGGCGGAAAAGGATCGATTACCGAccagaccgagctcttctacaccatgaaagccgctctggggctgaaataaaaccgcgcccggcagggaattttcgcttccgggtgggttttttgctgttagcaaccagacttaatggcagatcacgggcgcatacgctcatggttaaaac                                                                      <  1:234166‑M1/251‑150 (MQ=255)
                                                                                                                                            ggCACATGTGACCTCGCGCAAATGCTACGGTCCGAGCGCGACCAGTGAAAAATGTCCGGGAAACGCTCTGGAACAAGGCGGCAAAGGATCGATTACCGAACAGCTGCTTAACGCTCGTGCCGACGTTACGCTTGGCGGCGGCGCAAAAACCTTTGCTGAAACGGc                                                                                                                                                         >  1:214289/1‑165 (MQ=255)
                                                                                                                                              cacaTGTGACCTCGCGCAAATGCTACGGTCCGAGCGCGACCAGTGAAAAATGTCCGGGTAACGCTCTGGAAAAAGGCGGAAAAGGATCGATTACCGAccagaccgatctcttctac                                                                                                                                                                                                        >  1:313937‑M1/1‑97 (MQ=255)
                                                                                                                                                     gACCTCGCGAAAATGCTACGGTCCGANCGCGACNNGTNAAAAATGTCCGGGTAACGCNCTNNNANNNGNCGGAANAGNNTCGCTCACCNAccagaccgatctcttctacaccatgaaagccgctctggggctgacataaaccngcgaccggcnntgaattttcgctgccgggtggntttnntnntntttntan                                                                                                                    >  2:103909‑M1/1‑90 (MQ=255)
                                                                                                                                                                           ccGATCGCGACCAGTGAAAAATGTCCGTGTAACGCTCTTGAAAAAGGCGGAAAAGGATCGATTACCGAccagaccgatctcttctccaccatgaaagc                                                                                                                                                                                             <  1:321937‑M1/98‑31 (MQ=255)
                                                                                                                                                                           ccGAGCGCGACCAGTGAAAAATGTCCGGGTAACGCTCTGGAAAAAGGCGGAAAAGGATCGATTACCGAccagaccgatctcttctacaccatgaaagccgctctggggctgaaataaaaccgcgcccggcagtgaattttcgctgccgggtggtttttttgctgttagcaaccagacttaatggcagatcacgggcgcatacgctc                                                                                 >  1:307195‑M1/1‑68 (MQ=255)
                                                                                                                                                                                                   ccGGGTAACGCTCTGGAAAAAGGCGGAAAAGGATCGATTACCGAccagaccgatctcttctacaccatgaaagccgctctggggctgaaataaaaccgcgcccggcagtgaattttcgctgccgggtggtttttttgctgttagcaaccagacttaatggcagatcacgggcgcatacgctcatggttaaaacatgaagagggatgatgctatgaaaataacattactggttgccttgcttttcggtctgg            <  1:356215‑M1/251‑208 (MQ=255)
                                                                                                                                                                                                     gggTAACGCTCGGGAAAAAGGCGGAAAATGATCGATGACCGAccagaccgatctcttctacaccatgaaagccgctctggggctgaaataaaaccgcgcccggcagtgaattttcgctgccgggtggtttttttgctgttag                                                                                                                       <  1:188694‑M1/142‑101 (MQ=255)
                                                                                                                                                                                                               ctGGAAAAAGGCGGAAAAGGATCGATTACCGAccagaccgatctcttctacaccaggaaagccgctctggggctgaaataaaaccgcgcccggcagtgaattttcgctgccgggtggtttttttgctgttagcaaccagacttaatggcagatcacgggcgcatacgctcatggttaaaacatgaagagggatggtgctatgaaaataacattactggttaccgggcttttcggtctggtgtttttaac  <  1:103909‑M1/249‑218 (MQ=255)
                                                                                                                                                                                                                                              |                                                                                                                                                                                                                         
AAAACCTATAACGGCGCGCTGGGCGTCGATATTCACGAAAAAGATCACCCAACGATTCTGGAAATGGCAAAAGCCGCAGGTCTGGCGACCGGTAACGTTTCTACCGCAGAGTTGCAGGATGCCACGCCCGCTGCGCTGGTGGCACATGTGACCTCGCGCAAATGCTACGGTCCGAGCGCGACCAGTGAAAAATGTCCGGGTAACGCTCTGGAAAAAGGCGGAAAAGGATCGATTACCGAACAGCTGCTTAACGCTCGTGCCGACGTTACGCTTGGCGGCGGCGCAAAAACCTTTGCTGAAACGGCAACCGCTGGTGAATGGCAGGGAAAAACGCTGCGTGAACAGGCACAGGCGCGTGGTTATCAGTTGGTGAGCGATGCTGCCTCACTGAATTCGGTGACGGAAGCGAATCAGCAAAAACCCCTGCTTGGCCTGTTTGCTGACGGCAATATGCCA  >  NC_000913/402152‑402607

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.