Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,176,721 G→A T38I (ACC→ATC)  cpdA ← 3',5' cAMP phosphodiesterase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,176,7210GA100.0% 62.5 / NA 21T38I (ACC→ATC) cpdA3',5' cAMP phosphodiesterase
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (11/10);  total (11/10)

CAGCGAAATGCTGATAGGCCGCAGAGGATTGATCCTGCGCTAAATCACCTGTCGCGACAATCAGGTCGAATTCGTGCTGGTGTGGCCGAATCGCCTCCAGCACCGCCTGGTAACTCTCCCAGGTGTTTACCCCTAACAGGGCTTCGTGCTTTTGTGCAAACAGGTGAGTGTCGGTAATTTGTAAAATCCTGACTCTGGCCTCACCAGCCAGAGGAAGGGTTAACAGGCTTTCCAAATGGTGTCCTTAGGTTTCACGA  >  NC_000913/3176599‑3176855
                                                                                                                          |                                                                                                                                      
cAGCGAAATGCTGATAGGCCGCAGAGGATTGATCCTGCGCTAAATCACCTGTCGCGACAATCAGGTCGAATTCGTGCTGGTGTGGCCGAATCGCCTCCAGCACCGCCTGGTAACTCTCCCAGATGTTTACCCCTAACAg                                                                                                                        <  1:611625/139‑1 (MQ=255)
    gAAATGCTGATAGGCCGCAGAGGATTGATCCTGCGCTAAATCACCTGTCGCGACAATCAGGTCGAATTCGTGCTGGTGTGGCCGAATCGCCTCCAACACCGCCTGGTAACTCTCCCAGATGTTTACCCCTAACAGGGCt                                                                                                                    <  1:668619/139‑1 (MQ=255)
                  ccGCAGAGGATTGATCCTGCGCTAAATCACCTGTCGCGACACTCAGGTCGAATTCGTGCTGGTGTGGCCGAATCGCCTCCAGCACCGCCTGGTAACTCTCCCAGATGttt                                                                                                                                   <  1:521113/110‑1 (MQ=255)
                  ccGCAGAGGATTGATCCTGCGCTAAATCACCTGTCGCGACAATCAGGTCGAATTCGTGCTGGTGTGGCCGAATCGCCTCCAGCACCGCCTGGTAACTCTCCCAGATGttt                                                                                                                                   >  2:521113/1‑110 (MQ=255)
                                    gcgcTAAATCACCTGTCGCGACAATCAGGTCGAATTCGTGCTGGTGTGGCCGAATCGCCTCCAGCACCGCCTGGTAACTCTCCCAGATGTTTAcc                                                                                                                                >  1:105102/1‑95 (MQ=255)
                                    gcgcTAAATCACCTGTCGCGACAATCAGGTCGAATTCGTGCTGGTGTGGCCGAATCGCCTCCAGCACCGCCTGGTAACTCTCCCAGATGTTTAcc                                                                                                                                <  2:105102/95‑1 (MQ=255)
                                    gcgcTAAATCACCTGTCGCGACAATCAGGTCGAATTCGTGCTGGTGTGGCCGAATCGCCTCCAGCACCGCCTGGTAACTCTCCCAGATGTTTACCCCTAACAGGGCTTCGTGCTTTTGTGCAAACAGGTGAGTGTCGGt                                                                                    >  1:599865/1‑139 (MQ=255)
                                     cgcTAAATCACCTGTCGCGACAATCAGGTCGAATTCGTGCTGGTGTGGCCGAATCGCCTCCAGCACCGCCTGGTAACTCTCCCAGATGTTTACCCCTAACAGGGCTTCGTGCTTTTGTGCAAACAGGTGAGTTTCGGTa                                                                                   >  2:420063/1‑139 (MQ=255)
                                                             cAGGTCGAATTCGTGCTGGTGTGGCCGAATCGCCTCCAGCACCGCCTGGTAACTCTCCCAGATGTTTACCCCTAACAGGGCTTCGTGCTTTTGTGCAAACAGGTg                                                                                             <  2:513855/105‑1 (MQ=255)
                                                             cAGGTCGAATTCGTGCTGGTGTGGCCGAATCGCCTCCAGCACCGCCTGGTAACTCTCCCAGATGTTTACCCCTAACAGGGCTTCGTGCTTTTGTGCAAACAGGTg                                                                                             >  1:513855/1‑105 (MQ=255)
                                                                          tGCTGGTGTGGCCGAATCGCCTCCAGCACCGCCTGGTAACTCTCCCAGATGTTTACCCCTAACAGGGCTTCGTGCTTTTGTGCAAACAGGTGAGTGTCGGTAATTTGTAAAATCCTGACTCTGGCCTCACCAGCCagag                                              >  1:167041/1‑139 (MQ=255)
                                                                                         aTCGCCTCCAGCACCGCCTGGTAACTCTCCCAGATGTTTACCCCTAACAGGGCTTCGTGCTTTTGTGCAAACAGGTGAGTGTCGGTAATTTGTAAAATCCTGACTCTGGCCTCACCAGCCAGAGGAAGGGTTAACAGGc                               >  2:225381/1‑139 (MQ=255)
                                                                                         aTCGCCTCCAGCACCGCCTGGTAACTCTCCCAGATGTTTACCCCTAACAGGGCTTCGTGCTTTTGTGCAAACAGGTGAGTGTCGGTAATTTGTAAAATCCTGACTCTGGCCTCACCAGCCAGAGGAAGGGTTAACAGGc                               >  1:318929/1‑139 (MQ=255)
                                                                                                  aGCACCGCCTGGTAACTCTCCCAGATGTTTACCCCTAACAGGGCTTCGTGCTTTTGTGCAAACAGGTGAGTGTCGGTAATTTGTAAAATCCTGACTCTGGCCTCACCAGCCAGAGGAAGGGTTAACAGGCTTTCCAAAt                      <  2:280416/139‑1 (MQ=255)
                                                                                                       cGCCTGGTAACTCTCCCAGATGTTTACCCCTAACAGGGCTTCGTGCTTTTGTGCAAACAGGTGAGTGTCGGTAATTTGTAAAATCCTGACTCTGGCCTCACCAGCCAGAGGAAGGGTTAACAGGCTTTCCAAATGgtgt                 >  2:292895/1‑139 (MQ=255)
                                                                                                         ccTGGTAACTCTCCCAGATGTTTACCCCTAACAGGGCTTCGTGCTTTTGTGCAAACAGGTGAGTGTCGGTAATTTGTAAAATCCTGACTCTGGCCTCACCAGCCAGAGGAAGGGTTAACAGGCTTTCCAAATgg                    <  1:643674/134‑1 (MQ=255)
                                                                                                         ccTGGTAACTCTCCCAGATGTTTACCCCTAACAGGGCTTCGTGCTTTTGTGCAAACAGGTGAGTGTCGGTAATTTGTAAAATCCTGACTCTGGCCTCACCAGCCAGAGGAAGGGTTAACAGGCTTTCCAAATgg                    >  2:643674/1‑134 (MQ=255)
                                                                                                           tGGTAACTCTCCCAGATGTTTACCCCTAACAGGGCTTCGTGCTTTTGTGCAAACAGGTGAGTGTCGGTAATTTGTAAAATCCTGACTCTGGCCTCACCAGCCAGAGGAAGGGTTAACAGGCTTTCCAAATGGTGTCCtt             <  1:292895/139‑1 (MQ=255)
                                                                                                                    tcCCAGATGTTTACCCCTAACAGGGCTTCGTGCTTTTGTGCAAACAGGt                                                                                              <  1:631093/49‑1 (MQ=255)
                                                                                                                    tcCCAGATGTTTACCCCTAACAGGGCTTCGTGCTTTTGTGCAAACAGGt                                                                                              >  2:631093/1‑49 (MQ=255)
                                                                                                                      ccAGATGTTTACCCCTAACAGGGCTTCGTGCTTTTGTGCAAACAGGTGAGTGTCGGTAATTTGTAAAATCCTGACTCTGGCCTCACCAGCCAGAGGAAGGGTTAACAGGCTTTCCAAATGGTGTCCTTAGGTTTCacga  <  2:318929/139‑1 (MQ=255)
                                                                                                                          |                                                                                                                                      
CAGCGAAATGCTGATAGGCCGCAGAGGATTGATCCTGCGCTAAATCACCTGTCGCGACAATCAGGTCGAATTCGTGCTGGTGTGGCCGAATCGCCTCCAGCACCGCCTGGTAACTCTCCCAGGTGTTTACCCCTAACAGGGCTTCGTGCTTTTGTGCAAACAGGTGAGTGTCGGTAATTTGTAAAATCCTGACTCTGGCCTCACCAGCCAGAGGAAGGGTTAACAGGCTTTCCAAATGGTGTCCTTAGGTTTCACGA  >  NC_000913/3176599‑3176855

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