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] [1] 35 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.690) 31/474 0.5 100% coding (644/1416 nt) phoA alkaline phosphatase
?NC_000913 403114 = 0 (0.000)coding (1368/1416 nt) phoA alkaline phosphatase

GGCGCTGGAATTCGCTAAAAAGGAGGGTAACACGCTGGTCATAGTCACCGCTGATCACGCCCACGCCAGCCAGATTGTTGCGCCGGATACCAAAGCTCCGGGCCTCACCCAGGCGCTAAATACCAAAGATGGCGCAGTGATGGTGATGAGTTACGGGAACTCCGAAGAGGATTCACAAGAACATACCGGCAGTCAGTTGCGTATTGCGGCGTATGGCCCGCATGCCGCCAATGTTGTTGGACTGACCGACCAGACCGATCTCTTCTACACCATGAAAGCCGCTCTGGGGCTGAAATAAAACCGCGCCCGGCAGTGAATTTTCGCTGCCGGGTGGTTTTTTTGCTGTTAGCAACCAGACTTAATGGCAGATCACGGGCGCATACGCTCATGGTTAAAACATGAAGAGGGATGGTGCTATGAAAATAACATTACTGGTTACCTTGCTTTTCGGTCTGGTTTTTTTAACCACCGTCGGCGCTGCCGAGAGAACTTTAACC  >  NC_000913/402865‑403361
                                                                                                                                                                                                                                                         |                                                                                                                                                                                                                                                       
aaccggtgtcaaaacctataacggcgcgctgggcgtcgatattcacgaaaaagatcacccaacgattctggaaatggcaaaagccgcaggtctggcgaccggtaacgtttctaccgcagagttgcaggatgccacgcccgctgcgctggtggcacatgtgacctcgcgcaaatgctacggtccgagcgcgaccagtgaaaaatgtccgggtaacgctctggaaaaaggcggaaaaggatcgattaccgacc                                                                                                                                                                                                                                                        <  1:307945‑M2/2‑1 (MQ=255)
            aacctataacggcgcggtgggcgtcgatagtcaagaaaaagatccccccacgattctggaaattgcaaaagccgcaggtctggcgaccggtaacgtttctaccgcagaggtgcaggatgccacgcccgctgcgctggtggcacatgtgacctcgcgcaaatgctacggtccgagcgcgaccagtgaaaaatggccggttaacgctctggaaaaaggcggaaaaggatcgatttccgaCCAGACCGATctct                                                                                                                                                                                                                                            <  2:439922‑M2/14‑1 (MQ=255)
               ctataacggcgcgctgggcgtcgatattcacgaaaaagatcacccaacgattctggaaatggcaaaagccgcaggtctggcgaccggtaacgtttctaccgcagagttgcaggatgccacgcccgctgcgctggtggcacatgtgacctcgcgcaaatgctacggtccgagcgcgaccagtgaaaaatgtccgggtaacgctctggaaaaaggcggaaaaggatcgattaccgaCCAGACCGATCtcttct                                                                                                                                                                                                                                         >  2:227155‑M2/235‑251 (MQ=255)
                            ctgggcgtcgatattcacgaaaaagatcacccaacgattctggaaatggcaaaagccgcaggtctggcgaccggtaacgtttctaccgcagagttgcaggatgccacgcccgctgcgctggtggcacatgtgacctcgcgcaaatgctacggtccgagcgcgaccagtgaaaaatgtccgggtaacgctctggaaaaaggcggaaaaggatcgattaccgaCCAGACCGATCTCTTCTACACCATGAAAGc                                                                                                                                                                                                                            >  1:48269‑M2/222‑251 (MQ=255)
                                 cgtctatatgcacgaaaaagatcacccaacgattctggaaatggcaaaagccgcaggtctggcgaccggtaacgtttctaccgcagaggtgcaggatgccacgcccgctgcgctggtggcacatgtgacctcgcgcaaatgctacggtccgagcgcgaccagtgaaaaatgtccgggtaacgctctggaaaaaggcggaaaaggatcgattaccgaCCAGTCCGATCTCTTCTACACCACGAAAGCCGctc                                                                                                                                                                                                                       <  2:360862‑M2/35‑1 (MQ=255)
                                    cgatattcacgaaaaagatcacccaacgattctggaaatggcaaaagccgcaggtctggcgaccggtaacgtttctaccgcagagttgcaggatgccacgcccgctgcgctggtggcacatgtgacctcgcgcaaatgctacggtccgagcgcgaccagtgaaaaatgtccgggtaacgctctggaaaaaggcggaaaaggatcgattaccgaCCAGACCGATCTCTTCTACACCATGAAAGCCGCTCTg                                                                                                                                                                                                                     <  1:35061‑M2/37‑1 (MQ=255)
                                                aaaagatcacccaacgattctggaaatggcaaaagccgcaggtctggcgaccggtaacgtttctaccgcagagttgcaggatgccacgcccgctgcgctggtggcacatgtgacctcgcgcaaatgctacggtccgagcgcgaccagtgaaaaatgtccgggtaacgctctggaaaaaggcggaaaaggatcgattaccgaCCAGACCGATCTCTTCTACACCATGAAAGCCGCTCTGGGGCTGAAATaaa                                                                                                                                                                                                        <  1:227155‑M2/50‑1 (MQ=255)
                                                 aaagatcacccaacgattctggaaatggcaaaagccgcaggtctggcgaccggtaacgtttctaccgcagagttgcaggatgccacgcccgctgcgctggtggcacatgtgacctcgcgcaaatgctacggtccgagcgcgaccagtgaaaaatgtccgggtaacgctctggaaaaaggcggaaaaggatcgattaccgaCCAGACCGATCTCTTCTACACCATGAAAGCCGCTCTGGGGCTGAAATaaaa                                                                                                                                                                                                       <  1:270292‑M2/51‑1 (MQ=255)
                                                   agatcacccaacgattctggaaatggcaaaagccgcaggtctggcgaccggtaacgtttctaccgcagagttgcaggatgccacgcccgctgcgctggtggcacatgtgacctcgcgcaaatgctacggtccgagcgcgaccagtgaaaaatgtccgggtaacgctctggaaaaaggcggaaaaggatcgattaccgaCCAGACCGATCTCTTCTACACCATGAAAGCCGCTCTGGGGCTGAAATAAAAcc                                                                                                                                                                                                     <  1:40332‑M2/53‑1 (MQ=255)
                                                                   ctggaaatggcaaaagccgcaggtctggcgaccggtaacgtttctaccgcagagttgcaggatgccacgcccgctgcgctggtggcacatgtgacctcgcgcaaatgctacggtccgagcgcgaccagtgaaaaatgtccgggtaacgctctggaaaaaggcggaaaaggatcgattaccgaCCAGACCGATCTCTTCTACACCATGaa                                                                                                                                                                                                                               >  2:73145‑M2/183‑209 (MQ=255)
                                                                   ctggaaatggcaaaagccgcaggtctggcgaccggtaacgtttctaccgcagagttgcaggatgccacgcccgctgcgctggtggcacatgtgacctcgcgcaaatgctacggtccgagcgcgaccagtgaaaaatgtccgggtaacgctctggaaaaaggcggaaaaggatcgattaccgaCCAGACCGATCTCTTCTACACCATGaa                                                                                                                                                                                                                               <  1:73145‑M2/27‑1 (MQ=255)
                                                                   ctggaaatggcaaaagccgcaggtctggcgaccggtaacgtttctaccgcagagttgcaggatgccacgcccgctgcgctggtggcacatgtgacctcgcgcaaatgctacggtccgagcgcgaccagtgaaaaatgtccgggtaacgctctggaaaaaggcggaaaaggatcgattaccgaCCAGACCGATCTCTTCTACACCATGAAaaccgctctggggctgaaataacacagagcccagacgtgaat                                                                                                                                                                                     >  2:167729‑M2/183‑231 (MQ=255)
                                                                              aaaagccgcaggtctggcgaccggtaacgtttctaccgcagagttgcaggatgccacgcccgctgcgctggtggcacatgtgacctcgcgcaaatgctacggtccgagcgcgaccagtgaaaaatgtccgggtaacgctctggaaaaaggcggaaaaggatcgattacctaCCAGACCGATCTCTTCTACACCATGAAAGCCGCTCTGGGGCTGAAATAAAACcg                                                                                                                                                                                                    <  1:254712‑M2/54‑1 (MQ=255)
                                                                              aaaagccgcaggactggcgaccggaaacgtttctaccgcagcgttgcaggatgccacccccgatgcgctgggggcacatgcgacctcgcgcaaatgctacggtccgagcgcgaccagtgaaaaatttccgggtaacgctctggaaaaaggcggcaaaggatcgattaccgaCCAGGCCGGTCTCTTCTACACCATGAAagccgctctggggctgaaataaaaccgagatcggacgaccctcgtgtcgggt                                                                                                                                                                           >  2:254712‑M2/172‑225 (MQ=255)
                                                                               aaaggcgcatgtctggcgaccgggatcgtttctaccgcagaggtgcaggttgccacggccgctgtgctggtggcacttgtgacctcgcgcaaatgatacggtccgagcgcgaccagtgaaaaatgtccgggtaacgctctggaaaaaggcggaaaaggatggattaccgaCCAGTCCGATCTCTTCTACACCATGAAAGCCGCTCTGGGGCTGAAATAAAACCGCGCCCGGCAGTGAATTTTCGctgctgg                                                                                                                                                                         <  2:261032‑M2/81‑4 (MQ=255)
                                                                                                        acgtttctaccgcagagttgcaggatgccacgcccgctgcgctggtggcacatgtgacctcgcgcaaatgctacggtccgagcgcgaccagtgaaaaatgtccgggtaacgctctggaaaaaggcggaaaaggatcgattaccgaCCAGACCGCTCTCTTCTACACCATGAAAGCCGCTCTGGGGCTGAAATAAAACCGCGCCCGGCAGTGAATTTTCGCTGCCGGGTGGTTTTTTAGCTGTTAGCAAcc                                                                                                                                                 >  2:562541‑M2/146‑250 (MQ=255)
                                                                                                                     agagttgcaggatgccacgcccgctgcgctggtggcacatgtgacctcgcgcaaatgctacggtccgagcgcgaccagtgcaaaatgtccgggtaacgctctggaaaaaggcggaaaaggatcgattaccgaCCAGACCGCTCTCTTCTACACCATGAAAGCCGCTCTGGGGCTGAAATAAAACCGCGccc                                                                                                                                                                                               >  1:41081‑M2/133‑191 (MQ=255)
                                                                                                                     agagttgcaggatgacacgcccgctgcgctgtgtgcacatgtgacctcgcgcaaaggctacggtccgcgcgcgaccagtgaaaaatgcccgggtaacgctctggaaaaaggcggaaaaggatcgattaccgaCCAGACCGACCTCTGCTACGCCATGAAAGCCGCTCTGGGGCTGAAATAAACCCTCGccc                                                                                                                                                                                               <  2:41081‑M2/59‑1 (MQ=255)
                                                                                                                          tgcaggatgccacgcccgctgcgaaggtggcacatgtgacctcgcgcaaatgctacggtccgagcgcgaccagtgaaaaatgtccgggtaacgctctggaaaaaggcggaaaaggttcgattaccgaCCAGACCGATCTCTTCTACACCATGAAAGCCGCTCTGGGGCTGAAATACAACCGCGCCCGGCAGTGAATTTTCGCTGCCGGGTGGTTTTTTTGCTGTTAGCAACCAGACTTAATAGCAGATCAc                                                                                                                              >  2:224126‑M2/128‑251 (MQ=255)
                                                                                                                                           gccgatctgggggcacatgtgccctctcgctaaggctacggtccgagcgggacctgtgaaaaatgtccggggaacgctctggaaaaaggctgaaaaggatcgattaccgaCCAGACCGAGCTCGTCTACTCCATGAACGCCGCTCTGGGGCTGAAATAAAACCGCGCCCGGCAGTGAATTTTCGCTGCCGGGTGGTTTTGTTGCTGTTAGCATCCAGACTTAATGGCAGATCACGGGCTCATACGTTCATg                                                                                                             <  2:564452‑M2/141‑1 (MQ=255)
                                                                                                                                               cgctggtggcacatgtgacctcgcgcaaatgctacggtccgagcgcgaccagtgaaaaatgtccgggtaacgctctggaaaaaggcggaaaaggatcgattaccgaCCAGACCGATCTCTTCTACACCATGAAAGCCGCTCTGGGGCTGAAATAAAACCGCGCCCGGCAGTGAATTTTCGCTGCCGGGTGGTTTTTTTGCTGTTAGCAACCAGACTTAATGGCAGATCACGGGCGCATACGCTCATGGTTa                                                                                                         <  1:562541‑M2/145‑1 (MQ=255)
                                                                                                                                                           atgtgacctcgcgcaaatgctacggtccgagcgcgaccagtgaaaaatgtccgggtaacgctctggaaaaaggcggaaaaggatcgattaccgaCCAGACCGATCTCTTCTACACCATGAAAGCCGCTCTGGGGCTGAAATAAAACCGCGCCCGGCAGTGAATTTTCGCTGCCGGGTGGTTTTTTTGCTGTTAGCAACCAGACTTAATGGCAGATCACGGGCGCATACGCTCATGGTTAAAACATGAagag                                                                                             <  1:305031‑M2/157‑1 (MQ=255)
                                                                                                                                                                   ccgagcaactgctacgggccgagcgcgaccagtgaaaaatgtccgggtaacgctctggaaaaaggcggaaaaggatcgattacctaCCAGGCCGATCTCTTCTACACCATGAAAGCCGCTCTGGGGCTGAAATAAATCCGCGCCCGGCAGTGAATTTTCGCTGCCGGGTGGTTTTTTTGCTGTTAGCAACCAGACTGAATGGCAGATCACGGTCGCATACGCTCATGGTTAAAACATGAAGAGGGATGGTg                                                                                     <  2:96286‑M2/165‑1 (MQ=255)
                                                                                                                                                                     gcgcaaatgctacggtccgagcgcgaccagtgaaaaatgtccgggtaacgctctggaaaaaggcggaaaaggatcgattaccgaCCAGACCGATCTCTTCTACACCATGAAAGCCGCTCTGGGGCTGAAATAAAACCGCGCCCGGCAGTGAATTTTCGCTGCCGGGTGGTTTTTTTGCTGTTAGCAACCAGACTTAATGGCAGATCACGGGCGCATTCGCTCATGGTTAAAACATGAAAAGGGATGGTGCt                                                                                   >  2:254636‑M2/85‑251 (MQ=255)
                                                                                                                                                                                                      aaaaggttcgggttacgctctggaaaaaggtggaaatggatccattaccgaCCAGACCGATCTCTTCTACACTCTGAAAGTCGCTCTGGGGCTGGAATAAAACCGCGCCTGGCAGTGATTTTTCGCTGCCGGGTGGTTTTTTTGCTGTTAGCAACCAGACTGAATGGGAGATCACTGGCGCATACGCTTATGGTTAAAACATGAAGAGGGATGGTGCTATGAAAATAACATTTCTGGTTACCTTGCTTTTc                                                  <  2:62045‑M2/200‑1 (MQ=255)
                                                                                                                                                                                                       aaatgtccgggtaacgctctggaaaaaggcggaaaaggatcgattaccgaCCAGACCGATCTCTTCTACACCATGAAAGCCGCTCTGGGGCTGAAATAAAACCGCGCCCGGCAGTGAATTTTCGCTGCCGGGTGGTTTTTTTGCTGTTAGCAACCAGACTTAATGGCAGATCACGGGCGCATACGCTCATGGTTAAAACATGAAGAGGGATGGTGCTATGAAAATAACATTACTGGTTACCTTGCTTTTCg                                                 >  1:509037‑M2/51‑251 (MQ=255)
                                                                                                                                                                                                             ccgggtaacgctctggaaaaaggcggaaaaggatcgattaccgaCCAGACCGATCTCTTCTACACCATGAAAGCCGCTCTGGGGCTGAAATAAAACCGCGCCCGGCAGTGAATTTTCGCTGCCGGGTGGTTTTTTTGCTGTTAGCAACCAGACTTAATGGCAGATCACGGGCGCATACGCTCATGGTTAAAACATGAAGAGGGATGGTGCTATGAAAATAACATTACTGGTTACCTTGCTTTTCGGTCTgg                                           >  1:224982‑M2/45‑251 (MQ=255)
                                                                                                                                                                                                                ggtaacgctctggaagaaggcggaaaaggatcgattaccgaCCAGACCGATCTCTGCTACACCATGAAAGCCGCTCTGGGGTTGAAATAAAACCGCGCCCGGCAGTGAATTTTCGCTGCCGGGTGGTTTTTTTGCTGTTAGCAACCAGACTTAATGGCAGATCACGGGCGCATACGCTCATGGTTAAATCATGAAGAGGGATGGTGCTATGAAAATAACATTACTGGTTACCTTGCttt                                                    <  2:396422‑M2/198‑1 (MQ=255)
                                                                                                                                                                                                                ggtaacgctctggaaaaaggcggaaaaggatcgattaccgaCCAGACCGATCTCTTCTACACCATGAAAGCCGCTCTGGGGCTGAAATAAAACCGCGCCCGGCAGTGAATTTTCGCTGCCGGGTGGTTTTTTTGCTGTTAGCAACCAGACTTAATGGCAGATCACGGGCGCATACGCTCATGGTTAAAACATGAAGAGGGATGGTGCTATGAAAATAACATTACTGGTTACCTTGCttt                                                    >  1:396422‑M2/42‑239 (MQ=255)
                                                                                                                                                                                                                ggtaacgctctggaaaaaggcggaaaaggatcgattaccgaCCAGACCGATCTCTTCTACACCATGAAAGCCGCTCTGGGGCTGAAATAAAACCGCGCCCGGCAGTGAATTTTCGCTGCCGGGTGGTTTTTTTGCTGTTAGCAACCAGACTTAATGGCAGATCACGGGCGCATACGCTCATGGTTAAAACATGAAGAGGGATGGTGCTATGAAAATAACATTACTGGTTACCTTGCTTGTCGGTCTGGttt                                        <  1:138529‑M2/210‑1 (MQ=255)
                                                                                                                                                                                                                             aaaacggcggaaacggatcgattaacgaACAGACCGATATCTTCTACACCATGAAAGCCGCTCTGGGGCTGAACTAACACCGCGCCCGGCAGTGAATTTTCGCTGACGGGTGGTTATTTTGCTGTTAGAGACCAGGCTTAATGGCAGATAACGGGCGCATAAGCTCATGGTTAAAACCTGAAGTGGGATGGTGCTATGAAAATAACATTACTGGTTACCTTGCTTTTCGGTCTGGTTTTTTtacccacctt                           >  2:36495‑M2/29‑249 (MQ=255)
                                                                                                                                                                                                                              aaaaggcggaaaaggagcgattaccgcCCAGACCGATCTCTTCTACACCATGAAAGCCGCTCTGGGGCTGAAATAAACCCGCGCCCGGCAGTGAATTTTCGCTGCCGGGTGGTTTTTTTGCTGTTAGCAACCAGACTTAATGGCAGATCACGGGCGCATACGCTCATGGTTAAAACATGAAGAGGGATGGTGCTATGAAAATAACATTACTGGTTACCTTGCTTTTCGGTCTGGTTTTTTTAACCACCGTc                          <  2:479777‑M2/224‑1 (MQ=255)
                                                                                                                                                                                                                                     ggaaaaggagcgcttaccgaCCAGAACGATCTCTTCTACACCATGAAAGCCGCTCTGGGGCTGAAATAAAACCGCGCCCGGAAGTGAAGTTTCGCTGCCGGGTGGTTTTTTTGCTGTTAGCAACCAGTCTTAATGGAAGATCACGGGCGCATAAGCCCATGGTTAAAACATGAAGAGGGATGGAAATTTGAAAATAACATTCCTGGTTACCTTGCTTTTCGGTCTGTTTTTTTTAACCaccgtcgacggtg                   >  2:556202‑M2/21‑245 (MQ=255)
                                                                                                                                                                                                                                       aaaaggatcgattaccgaCCAGACCGATCTCTTCTACACCATGAAAGCCGCTCTGGGGCTGAAATAAAACCGCGCCCGGCAGTGAATTTTCGCTGCCGGGTGGTTTTTTTGCTGTTAGCAACCAGACTTAATGGCAGATCACGGGCGCATACGCTCATGGTTAAAACATGAAGAGGGATGGTGCTATGAAAATAACATTACTGGTTACCTTGCTTTTCGGTCTGGTTTTTTTAACCACCGTCGGCGCTGcc                 >  1:244312‑M2/19‑251 (MQ=255)
                                                                                                                                                                                                                                                      cgaccAGACCGATCTCTTCTACACCATGAAAGCCGCTCTGGGGCTGAAATAAAACCGCGCCCGGCAGTGAATTTTCGCTGCCGGGTGGTTTTTTTGCTGTTAGCAACCAGACTTAATGGCAGATCACGGGCGCATACGCTCATGGTTAAAACATGAAGAGGGATGGTGCTATGAAAATAACATTACTGGTTACCTTGCTTTTCGGTCTGGTTTTTTTAACCACCGTCGGCGCTGCCGAGAGAACTTTAAcc  >  1:561996/1‑251 (MQ=255)
                                                                                                                                                                                                                                                         |                                                                                                                                                                                                                                                       
GGCGCTGGAATTCGCTAAAAAGGAGGGTAACACGCTGGTCATAGTCACCGCTGATCACGCCCACGCCAGCCAGATTGTTGCGCCGGATACCAAAGCTCCGGGCCTCACCCAGGCGCTAAATACCAAAGATGGCGCAGTGATGGTGATGAGTTACGGGAACTCCGAAGAGGATTCACAAGAACATACCGGCAGTCAGTTGCGTATTGCGGCGTATGGCCCGCATGCCGCCAATGTTGTTGGACTGACCGACCAGACCGATCTCTTCTACACCATGAAAGCCGCTCTGGGGCTGAAATAAAACCGCGCCCGGCAGTGAATTTTCGCTGCCGGGTGGTTTTTTTGCTGTTAGCAACCAGACTTAATGGCAGATCACGGGCGCATACGCTCATGGTTAAAACATGAAGAGGGATGGTGCTATGAAAATAACATTACTGGTTACCTTGCTTTTCGGTCTGGTTTTTTTAACCACCGTCGGCGCTGCCGAGAGAACTTTAACC  >  NC_000913/402865‑403361

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 31 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: