Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 875,327 G→A *783* (TAG→TAA yliE → putative membrane‑anchored cyclic‑di‑GMP phosphodiesterase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913875,3270GA100.0% 98.6 / NA 32*783* (TAG→TAAyliEputative membrane‑anchored cyclic‑di‑GMP phosphodiesterase
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (16/16);  total (16/16)

TTGCGCGCGAATTAATCAAATTAGGATGTGTTCAGGCTCAGGGGTATCTGTACCAGAAACCCATGCCATTCTCCGCCTGGGATAAAAGTGGAAAATTAGTAAAAGAGTAGTTTACGTATGTCCAGAATCAATAAGTTCGTACTTACAGTCAGTCTGCTGATTTTTATCATGATTTCAGCAGTTGCCTGCGGGATCTACACTCAAATGGTAAAGGAACGGGTGTATAGCCTGAAAC  >  NC_000913/875218‑875452
                                                                                                             |                                                                                                                             
ttGCGCGCGAATTAATCAAATTAGGATGTGTTCAGGCTCAGGGGTATCTGTACCAGAAACCCATGCCATTCTCCGCCTGGGATAAAAGTGGAAAATTAGTAAAAGAGTAATTTACGTATGTCCAGAATCAATAAGTTCg                                                                                                  <  2:500087/139‑1 (MQ=255)
                             gtTCAGGCTCAGGGGTATCTGTACCAGAAACCCATGCCATTCTCCGCCTGGGATAAAAGTGGAAAATTAGTAAAAGAGTAATTTACGTATGTCCAGAATCAATAAGTTCGTACTTACAGTCAGTCTGCTGATTTTTATc                                                                     >  2:47709/1‑139 (MQ=255)
                              ttCAGGCTCAGGGGTATCTGTACCAGAAACCCATGCCATTCTCCGCCTGGGATAAAAGTGGAAAATTAGTAAAAGAGTAAtt                                                                                                                             <  2:575543/82‑1 (MQ=255)
                              ttCAGGCTCAGGGGTATCTGTACCAGAAACCCATGCCATTCTCCGCCTGGGATAAAAGTGGAAAATTAGTAAAAGAGTAAtt                                                                                                                             >  1:575543/1‑82 (MQ=255)
                                   gCTCAGGGGTATCTGTACCAGAAACCCATGCCATTCTCCGCCTGGGATAAAAGTGGAAAATTAGTAAAAGAGTAATTTACGt                                                                                                                        <  1:147737/82‑1 (MQ=255)
                                   gCTCAGGGGTATCTGTACCAGAAACCCATGCCATTCTCCGCCTGGGATAAAAGTGGAAAATTAGTAAAAGAGTAATTTACGt                                                                                                                        >  2:147737/1‑82 (MQ=255)
                                                  tACCAGAAACCCATGCCATTCTCCGCCTGGGATAAAAGTGGAAAATTAGTAAAAGAGTAATTTACGTATGTCCAGAATCAATAAGTTCGTACTTACAGTCAGTCTGCTGATTTTTATCATGATTTCAGCAGTTGCCTGc                                                >  1:844390/1‑139 (MQ=255)
                                                   aCCAGAAACCCATGCCATTCTCCGCCTGGGATAAAAGTGGAAAATTAGTAAAAGAGTAATTTACGTATGTCCa                                                                                                                 >  2:580192/1‑73 (MQ=255)
                                                   aCCAGAAACCCATGCCATTCTCCGCCTGGGATAAAAGTGGAAAATTAGTAAAAGAGTAATTTACGTATGTCCa                                                                                                                 <  1:580192/73‑1 (MQ=255)
                                                   aCCAGAAACCCATGCCATTCTCCGCCTGGGATAAAAGTGGAAAATTAGTAAAAGAGTAATTTACGTATGTCCAGAATCAATAAGTTCGTACTTACAGTCAGTCTGCTGATTTTTATCATGATTTCAGCAGTTGCCTGCg                                               >  2:480043/1‑139 (MQ=255)
                                                          aCCCATGCCATTCTCCGCCTGGGATAAAAGTGGAAAATTAGTAAAAGAGTAATTTACGTATGTCCAGAATCAATAAGTTCGTACTTACAGTCAGTCTGCTGATTTTTATCATGATTTCAGCAGTTGCCTGCGGGATCTa                                        >  1:969191/1‑139 (MQ=255)
                                                            ccATGCCATTCTCCGCCTGGGATAAAAGTGGAAAATTAGTAAAAGAGTAATTTACGTATGTCCAGAATCAATAAGTTCGTACt                                                                                              <  2:666539/83‑1 (MQ=255)
                                                            ccATGCCATTCTCCGCCTGGGATAAAAGTGGAAAATTAGTAAAAGAGTAATTTACGTATGTCCAGAATCAATAAGTTCGTACt                                                                                              >  1:666539/1‑83 (MQ=255)
                                                                  cATTCTCCGCCTGGGATAAAAGTGGAAAATTAGTAAAAGAGTAATTTACGTATGTCCAGAATCAATAAGTTCGTACTTACAGTCAGTCTGCTGATTTTTATCATGATTTCAGCAGTTGCCTGCGGGATCTACACTCaaa                                <  1:47709/139‑1 (MQ=255)
                                                                    ttCTCCGCCTGGGATAAAAGTGGAAAATTAGTAAAAGAGTAATTTACGTATGTCCAGAATCAATAAGtt                                                                                                    <  1:487347/69‑1 (MQ=255)
                                                                    ttCTCCGCCTGGGATAAAAGTGGAAAATTAGTAAAAGAGTAATTTACGTATGTCCAGAATCAATAAGtt                                                                                                    >  2:487347/1‑69 (MQ=255)
                                                                    ttCTCCGCCTGGGATAAAAGTGGAAAATTAGTAAAAGAGTAATTTACGTATGTCCAGAATCAATAAGTTCGTACTTACAGTCAGTCTGCTGATTTTTATCATGATTTCAGCAGTTGCCTGCGGGATCTACACTCAAATg                              >  1:847485/1‑139 (MQ=255)
                                                                     tctcCGCCTGGGATAAAAGTGGAAAATTAGTAAAAGAGTAATTTACGTATGTCCAGAATCAATAAGTTCGTACTTACAGTCAGTCTGCTGATTTTTATCAt                                                                   >  2:150524/1‑101 (MQ=255)
                                                                     tctcCGCCTGGGATAAAAGTGGAAAATTAGTAAAAGAGTAATTTACGTATGTCCAGAATCAATAAGTTCGTACTTACAGTCAGTCTGCTGATTTTTATCAt                                                                   <  1:150524/101‑1 (MQ=255)
                                                                     tctcCGCCTGGGATAAAAGTGGAAAATTAGTAAAAGAGTAATTTACGTATGTCCAGAATCAATAAGTTCGTACTTACAGTCAGTCTGCTGATTTTTATCATGATTTCAGCAGTTGCCTGCGGGATCTACACTCAAATgg                             >  1:112652/1‑139 (MQ=255)
                                                                      ctcCGCCTGGGATAAAAGTGGAAAATTAGTAAAAGAGTAATTTACGTATGTCCAGAATCAATAAGTTCGTACTTACAGTCAGTCTGCTGATTTTTATCATGATTTCAGCAGTTGCCTGCGGGATCTACACTCAAATGGt                            <  1:480043/139‑1 (MQ=255)
                                                                              gggATAAAAGTGGAAAATTAGTAAAAGAGTAATTTACGTATGTCCAGAATCAATAAGTTCGTACTTACAGTCAGTCTGCTGATTTTTATCATGATTTCAGCAGTTGCCTGCGGGATCTACACTCAAATGGTAAAGGAAc                    <  2:969191/139‑1 (MQ=255)
                                                                               ggATAAAAGTGGAAAATTAGTAAAAGAGTAATTTACGTATGTCCAGAATCAATAAGTTCGTACTTACAGTCAGTCTGCTGATTTTTATCATGATTTCAGCAGTTGCCTGCGGGATCTACACTCAAATGGTAAAGGAACg                   <  2:112652/139‑1 (MQ=255)
                                                                                     aaGTGGAAAATTAGTAAAAGAGTAATTTACGTATGTCCAGAATCAATAAGTTCGTACTTACAGTCAGTCTGCTGATTTTTATCATGATTTCAGCAGTTGCCTGCGGGATCTACACTCAAATGGTAAAGGAACGGGTGta             <  2:732166/139‑1 (MQ=255)
                                                                                          gAAAATTAGTAAAAGAGTAATTTACGTATGTCCAGAATCAATAAGTTCg                                                                                                  <  1:707510/49‑1 (MQ=255)
                                                                                          gAAAATTAGTAAAAGAGTAATTTACGTATGTCCAGAATCAATAAGTTCg                                                                                                  >  2:707510/1‑49 (MQ=255)
                                                                                          gAAAATTAGTAAAAGAGTAATTTACGTATGTCCAGAATCAATAAGTTCGTACTTACAGTCAGTctgct                                                                               >  1:930424/1‑68 (MQ=255)
                                                                                          gAAAATTAGTAAAAGAGTAATTTACGTATGTCCAGAATCAATAAGTTCGTACTTACAGTCAGTctgct                                                                               <  2:930424/68‑1 (MQ=255)
                                                                                          gAAAATTAGTAAAAGAGTAATTTACGTATGTCCAGAATCAATAAGTTCGTACTTACAGTCAGTCTGCTGATTTTTATCATGATTTCAGCAGTTGCCTGCGGGATCTACACTCAAATGGTAAAGGAACGGGTGTATAGcc        >  2:134959/1‑139 (MQ=255)
                                                                                                tAGTAAAAGAGCAAGTTACGTATGTCCAAATTCTATAAGATCGTACTTGCTGTCAGTCTGCTTATTTTTATTATGATTGCATCAGTTGCCACCGGGATCTACACTCAAATGGTAAAGGAACGGGTGTATAGCCTGAAAc  <  1:470498/139‑1 (MQ=255)
                                                                                                  gTAAAAGAGTAATTTACGTATGTCCAGAATCAATAAGTTCGTACTTACAGTCAGTCTGCTGATTTTTATCATGATTTCAGCAGt                                                       >  2:569244/1‑84 (MQ=255)
                                                                                                  gTAAAAGAGTAATTTACGTATGTCCAGAATCAATAAGTTCGTACTTACAGTCAGTCTGCTGATTTTTATCATGATTTCAGCAGt                                                       <  1:569244/84‑1 (MQ=255)
                                                                                                             |                                                                                                                             
TTGCGCGCGAATTAATCAAATTAGGATGTGTTCAGGCTCAGGGGTATCTGTACCAGAAACCCATGCCATTCTCCGCCTGGGATAAAAGTGGAAAATTAGTAAAAGAGTAGTTTACGTATGTCCAGAATCAATAAGTTCGTACTTACAGTCAGTCTGCTGATTTTTATCATGATTTCAGCAGTTGCCTGCGGGATCTACACTCAAATGGTAAAGGAACGGGTGTATAGCCTGAAAC  >  NC_000913/875218‑875452

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 16 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: