Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 188,837 C→T 100% R316C (CGC→TGC)  PP_0165 → GGDEF domain‑containing protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947188,8370CT100.0% 59.3 / NA 21R316C (CGC→TGC) PP_0165GGDEF domain‑containing protein
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (4/17);  total (4/17)
Rejected as polymorphism: Frequency below/above cutoff threshold.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

TGGCTGGCCTGAACGCCCGCCTTGGCGGCCAGCGTACCGACCAGTTGTTGCAAGCGGTAGGTGAACAGCTGCGTCGTACTTGTGCCAGCTACCCGGAAACCAACGACCTCATTTCCCGCAGCCGCGGCGGCGAGTTTGCAGTGCTGGCACCGGGCATGGTGCATGAGGAGGCAGTACACCTTGCCCAGGCTCTGGAGGCAACGTTGCAGAGCCTGCACGAAACGGGGGCCAGTGATATCGACCCGGT  >  NC_002947/188721‑188967
                                                                                                                    |                                                                                                                                  
tggctggcCTGAACGCCCGCCTTGGCGGCCAGCGTACCGACCAGTTGTTGCAAGCGGTAGGTGAACAGCTGCGTCGTACTTGTGCCAGCTACCCGGAAACCAACGACCTCATTTCCTGCAGCCGCGGCGGCGAGtt                                                                                                                 <  2:477873/136‑1 (MQ=255)
   ctggcCTGAACGCCCGCCTTGGCGGCCAGCGTACCGACCAGTTGTTGCAAGCGGTAGGTGAACAGCTGCGTCGTACTTGTGCCAGCTACCCGGAAACCAACTACCTCATTTCCTGCAGCCGCGGCGGCGAGTTTGc                                                                                                              >  2:166452/1‑136 (MQ=255)
   ctggcCTGAACGCCCGCCTTGGCGGCCAGCGTACCGACCAGTTGTTGCAAGCGGTAGGTGAACAGCTGCGTCGTACTTGTGCCAGCTACCCGGAAACCAACGACCTCATTTCCTGCAGCCGCGGCGGCGAGTTTGc                                                                                                              >  2:342932/1‑136 (MQ=255)
      gcCTGAACGCCCGCCTTGGCGGCCAGCGTACCGACCAGTTGTTGCAAGCGGTAGGTGAACAGCTGCGTCGTACTTGTGCCAGCTACCCGGAAACCAACGACCTCATTTCCTGCAGCCGCGGCGGCGAGTTTGCAGt                                                                                                           <  1:56533/136‑1 (MQ=255)
                 cgccTTGGCGGCCAGCGTACCGACCAGTTGTTGCAAGCGGTAGGTGAACAGCTGCGTCGTACTTGTGCCAGCTACCCGGAAACCAACGACCTCATTTCCTGCAGCCGCGGCGGCGAGTTTGCAGTGCTGGCACCgg                                                                                                <  2:404351/136‑1 (MQ=255)
                 cgccTTGGCGGCCAGCGTACCGACCAGTTGTTGCAAGCGGTAGGTGAACAGCTGCGTCGTACTTGTGCCAGCTACCCGGAAACCAACGACCTCATTTCCTGCAGCCGCGGCGGCGAGTTTGCAGTGCTGGCACCgg                                                                                                <  1:280726/136‑1 (MQ=255)
                 cgccTTGGCGGCCAGCGTACCGACCAGTTGTTGCAAGCGGTAGGTGAACAGCTGCGTCGTACTTGTGCCAGCTACCCGGAAACCAACGACCTCATTTCCTGCAGCCGCGGCGGCGAGTTTGCAGTGCTGGCACCgg                                                                                                <  1:353455/136‑1 (MQ=255)
                        gcggcCAGCGTACCGACCAGTTGTTGCAAGCGGTAGGTGAACAGCTGCGTCGTACTTGTGCCAGCTACCGGGAAACCAAGGACCTCATTTCCTGCAGCCGCGGCGGCGAGTTTGCAGTGCTGGCGCCGGGccgggg                                                                                         >  2:200235/1‑131 (MQ=255)
                             cAGCGTACCGACCAGTTGTTGCAAGCGGTAGGTGAACAGCTGCGTCGTACTTGTGCCAGCTACCCGGAAACCAACGACCTCATTTCCTGCAGCCGCGGCGGCGAGTTTGCAGTGCTGGCACCGGGCATGGTGCATg                                                                                    <  2:332243/136‑1 (MQ=255)
                                    ccGACCAGTTGTTGCAAGCGGTAGGTGAACAGCTGCGTCGTACTTGTGCAAGCTACCCGGAAACCAACGACCTCATTTCCTGCAGCCGCGGCGGCGAGTTTGCAGTGCTGGCACCGGGCATGGTGCATGAGGAGGc                                                                             <  1:497114/136‑1 (MQ=255)
                                             tgttgCAAGCGGTAGGTGAACAGCTGCGTCGTACTTGTGCCAGCTACCCGGAAACCAACGACCTCATTTCCTGCAGCCGCGGCGGCGAGTTTGCAGTGCTGGCACCGGGCATGGTGCATGAGGAGGCAGTACACCt                                                                    <  1:342932/136‑1 (MQ=255)
                                              gttgCAAGCGGTAGGTGAACAGCTGCGTCGTACTTGTGCCAGCTACCCGGAAACCAACGTCCTCATTTCCTGCAGCCGCGGCGGCGAGTTTGCAGTGCTGGCACCGTGCATGGTTCATGAGGAGGCAGTTCACCtt                                                                   <  1:375349/136‑1 (MQ=255)
                                              gttgCAAGCGGTAGGTGAACAGCTGCGTCGTACTTGTGCCAGCTACCCGGAAACCAACGACCTCATTTCCTGCAGCCGCGGCGGCGAGTTTGCAGTGCTGGCACCGGGCATGGTGCATGAGGAGGCAGTACACCtt                                                                   <  2:380494/136‑1 (MQ=255)
                                                             tGAACAGCTGCGTCGTACTTGTGCCAGCTACCCGGAAACCAACGACCTCATTTCCTGCAGCCGCGGCGGCGAGTTTGCAGTGCTGGCACCGGGCATGGTGCATGAGGAGGCAGTACACCTTGCCCAGGCTCTGGAg                                                    <  1:229267/136‑1 (MQ=255)
                                                             tGAACAGCTGCGTCGTACTTGTGCCAGCTACCCGGAAACCAACGACCTCATTTCCTGCAGCCGCGGCGGCGAGTTTGCAGTGCTGGCACCGGGCATGGTGCATGAGGAGGCAGTACACCTTGCCCAGGCTCTGGAg                                                    <  1:216112/136‑1 (MQ=255)
                                                                     tGCGTCGTACTTGTGCCAGCTACCCGGAAACCAACTACCTCATTTCCTGCAGCCGCGGCGGCGAGTTTGCAGTGCTGGCACCGGGCATGGTGCATGAGGAGGCAGTACACCTTGCCCAGGCTCTGGAGGCAACGtt                                            <  1:166452/136‑1 (MQ=255)
                                                                      gCGTCGTACTTGTGCCAGCTACCCGGAAACCAACGACCTCATTTCCTGCAGCCGCGGCGGCGAGTTTGCAGTGCTGGCACCGGGCATGGTGCATGAGGAGGCAGTACACCTTGCCCAGGCTCTGGAGGCAACGTTg                                           >  2:476310/1‑136 (MQ=255)
                                                                                             cGGAAACCAACGACCTCATTTCCTGCAGCCGCGGCGGCGAGTTTGCAGTGCTGGCACCGGGCATGGTGCATGAGGAGGCAGTACACCTTGCCCAGGCTCTGGAGGCAACGTTGCAGAGCCTGCACGAAACGGGGGc                    <  1:164188/136‑1 (MQ=255)
                                                                                                          ccTCATTTCCTGCAGCCGCGGCGGCGAGTTTGCAGTGCTGGCACCGGGCATGGTGCATGAGGAGGCAGTACACCTTGCCCAGGCTCTGGAGGCAACGTTGCAGAGCCTGCACGAAACGGGGGCCAGTGATATCGAc       <  1:67527/136‑1 (MQ=255)
                                                                                                               tttCCTGCCGCCGCGGCGGCGAGTTTGCAGTGCTGGCACCGGGCATGGTGCATGAGGAGGCAGTACACCTTGCCCTGGCTCTGGAGGCAACGTTGCAGAGCCTGCACGAAACGGGGGCCAGTGATATCGACCCGGt  <  1:200235/136‑1 (MQ=255)
                                                                                                               tttCCTGCAGCCGCGGCGGCGAGTTTGCAGTGCTGGCACCGGGCATGGTGAATGAGGAGGCAGTACACCTTGCCCAGGCTCTGGAGGCAACGTTGCAGAGCCTGCACGAAACGGGGGCCAGTGATATCGACCCGGt  <  1:270283/136‑1 (MQ=255)
                                                                                                                    |                                                                                                                                  
TGGCTGGCCTGAACGCCCGCCTTGGCGGCCAGCGTACCGACCAGTTGTTGCAAGCGGTAGGTGAACAGCTGCGTCGTACTTGTGCCAGCTACCCGGAAACCAACGACCTCATTTCCCGCAGCCGCGGCGGCGAGTTTGCAGTGCTGGCACCGGGCATGGTGCATGAGGAGGCAGTACACCTTGCCCAGGCTCTGGAGGCAACGTTGCAGAGCCTGCACGAAACGGGGGCCAGTGATATCGACCCGGT  >  NC_002947/188721‑188967

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