Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,434,657 C→T G637R (GGG→AGG)  cpdB ← 2',3'‑cyclic‑nucleotide 2'‑phosphodiesterase/3'‑nucleotidase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,434,6570CT96.7% 80.8 / ‑6.7 30G637R (GGG→AGG) cpdB2',3'‑cyclic‑nucleotide 2'‑phosphodiesterase/3'‑nucleotidase
Reads supporting (aligned to +/- strand):  ref base C (1/0);  new base T (19/10);  total (20/10)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.98e-01

GAAAAGGTTGCTGCACTGGCGGACTGGATTGAGTTGAATTTGCCCGAGGTGTTGGCGGTGCGGGATGAACGTGCGGCATAACTGCCGGATGCGCTGAGCTTATCTGGTCTACGGTGCATGGCTTGTAGGGCGGATAAGATGCGCCAGCATCGCATCCGCCAATAAATGCCGGATGCGGTTGATGATTTATAGGCCGAAAAAGAAGTGTTTTACTTACTCAAATCCACCTGATAAATCGCAAACCCGATATCATCGGTCGCGACTTTATTCATCGGATACTGCCCTTTCTCTTTAATAAACGCTGCGGCTTTATCTGACGGAGAGGTTTCGAAACGGATATCCAGTTTCTTATCGCCAGCTATCGGTGCTAAACGCCAGTTGTTATCTGCCGCCGGGTGAATTTCCCCCGCACGTTTCGACTCATCAGCAATCCACGCTGCCAGCACCGAGCGGTTCTCATCCGGTGAAGC  >  NC_000913/4434413‑4434882
                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                 
gAAAAGGTTGCTGCACTGGCGGACTGGATTGAGTTGAATTTGCCCGAGGTGTTGGCGGTGCGGGATGAACGTGCGGCATAACTGCCGGATGCGCTGAGCTTATCTGGTCTACGGTGCATGGCTTGTAGGGCGGATAAGATGCGCCAGCATCGCATCCGCCAATAAATGCCGGATGCGGTTGATGATTTATAGGCCGAAAAAGAAGTGTTTTACTTACTCAAATCCACCTGATAAATCGCAAACCTGatat                                                                                                                                                                                                                              >  1:321933/1‑250 (MQ=255)
              aCTGGCGGACTGGATTGAGTTGAATTTGCCCGAGGTGTTGGCGGGGCGGGATGAACGTGCGGCATAACTGCCGGATGCGCTGGGCTTATCTGGTCTACGGGGCATGGCTTGGAGGGGGGATAAGATGCGCCAGCATCGCATCCGCCAATAAATGCCGGATGCGGTTGATGATTTATAGGCCGACAAAGACGTGTTTTACTTACTAAAATCCCCCTGCTAAATCGCCAAACTGCTATCATCGGTAGCGACtt                                                                                                                                                                                                               >  2:329997/1‑251 (MQ=255)
                   cGGACTGGATTGAGTTGCATTTGCCCGAGGTGTTGGCGGTGCGGGATGAACGTGCGGCATAACTGCCGGATGCGCTGAGCTTATATGGTCTACGGTGCATGGCATGTAGGGCGGATAAGAAGCGCCAGCATCGCATCCGCCCATAAATGCCGGATGCGGGTGATGCTTTATAGGCCGAAAAAGAAGGGTTTTACTTACTCAAATCCACATGATAAATCGCAAACCTTATATCATCGGTCGAGACTTTATTc                                                                                                                                                                                                          >  1:417019/1‑251 (MQ=255)
                                           ccGAGGTGTTGGCGGTGCGGGATGAACGTGCGGCATAACTGCCGGATGCGCTGAGCTTATCTGGTCTACGGTGCATGGCTTGTAGGGCGGATAAGATGCGCCAGCATCGCATCCGCCAATAAATGCCGGATGCGGTTGATGATTTATAGGCCGAAAAAGAAGTGTTTTACTTACTCAAATCCACCTGATAAATCGCAAACCTGATATCATCGGTCGCGACTTTATTCATCGGATACTGccc                                                                                                                                                                                            <  1:91333/241‑1 (MQ=255)
                                           ccGAGGTGTTGGCGGTGCGGGATGAACGTGCGGCATAACTGCCGGATGCGCTGAGCTTATCTGGTCTACGGTGCATGGCTTGTAGGGCGGATAAGATGCGCCAGCATCGCATCCGCCAATAAATGCCGGATGCGGTTGATGATTTATAGGCCGAAAAAGAAGTGTTTTACTTACTCAAATCCACCTGATAAATCGCAAACCTGATATCATCGGTCGCGACTTTATTCATCGGATAATGccc                                                                                                                                                                                            >  2:91333/1‑241 (MQ=255)
                                            cGAGGAGTTGGCGGTGCGGGATGAACGTGCGGCATAACTGCCGGATGCGCTGCGCTTATCTGGTCTACGGGGCATGGCTTGTAGGGCGGATAAGATGCGCCAGCATCGCCTCCGCCAATAAATGCCGGGTGCGGTTGATGATTTATAGGCCGAAAAAGAAGTGTTTTACTTACTCAAATCCCCCTGTTAAATCGCAACCCTGATATCATCGGTCGCGACTTTAGTCATAGGATACTGCCCTTTCTCtttt                                                                                                                                                                                  >  2:278101/1‑249 (MQ=255)
                                                  gtTGGCGGTGCGGGATGAACGTGCGGCATAACTGCCGGATGCGCTGAGCTTATCTGGTCTACGGTGCATGGCTTGTAGGGCGGATAAGATGCGCCAGCATCGCATCCGCCAATAAATGCCGGATGCGGTTGATGATTTATAGGCCGAAAAAGAAGTGTTTTACTTACTCAAATCCACCTGATAAATCGCAAACCTGATATCATCGGTCGCGACTTTATTCATCGGATACTGCCCTTTCTCTTTAATAAACg                                                                                                                                                                           >  1:116886/1‑251 (MQ=255)
                                                          tGCGGGATGAACGTGCGGCATAACTGCCGGATGCGCTGAGCTTATCTGGTCTACGGTGCATGGCTTGTAGGGCGGATAAGATGCGCCAGCATCGCATCCGCCAATAAATGCCGGATGCGGTTGATGATTTATAGGCCGAAAAAGAAGTGTTTTACTTACTCAAATCCACCTGCTAAATCGCAAACCCGATATCATCGGTCGCGACTTTATTCATCGGATACTGCCCTTTCTCTTTAATAAACGCTGCGGcc                                                                                                                                                                   >  2:117012/1‑250 (MQ=255)
                                                                      gTGCGGCATAACTGCCGGATGCGCTGAGCTTATCTGGTCTACGGTGCATGGCTTGTAGGGCGGATAAGATGCGCCAGCATCGCATCCGCCAATAAATGCCGGATGCGGTTGATGATTTATAGGCCGAAAAAGAAGTGTTTTACTTACTCAAATCCACCTGATAAATCGCAAACCTGATATCATCGGTCGCGACTTTATTCATCGGATACTGCCCTTTCTCTTTAATAAACGCTGCGGCTTTATCTTACggg                                                                                                                                                       >  2:340352/1‑250 (MQ=255)
                                                                                     cGGATGCGCTGAGCTTATCTGGTCTACGGTGCATGGCTTGTAGGGCGGATAAGATGCGCCAGCATCGCATCCGCCAATAAATGCCGGATGCGGTTGATGATTTATAGGCCGCAAAAGAAGTGTTTTACTTACTCAAATCCACCTGATAAATCGCAAACCTGATATCATCGGTCGCGACTTTATTCATCGGATACTGCCCTTTCTCTTTAATAAACGCTGCGGCTTTATCTGACGGAGAGGTTTCGAAACgg                                                                                                                                        >  2:251832/1‑251 (MQ=255)
                                                                                                  cTTATCTGGTCTACGGTGCATGGCTTGTAGGGCGGCTACGATGCGCCAGCATCGCATCCGCCAATAAATGCCGGATGCGGTTGATGATTTATAGGCCGAAAAAGAAGTGTTTTACTTACTCAAATCCACCTGATAAATCGCAAACCTGATATCATCGGTCGCGACTTTATTCATCGGATACTGCCCTTTCTCTTTAATAAACGCTTCGGCTTTATCTGACGGAGCGGTTTCGACACGGATATCCAGTTTCt                                                                                                                           >  2:172802/1‑251 (MQ=255)
                                                                                                                                    gATAAGAGGGGCCAGCATCGCTGCCGCTAATAAAGGCCGGATGCGGTTGATGATTTATAGGCCGAAAAAGAAGTGTTTTACGTACTCAAATCCACCTGATAAATCGCAAACCTGATATCATCGGTCGCGACTTTATTCATCGGATACTGCCCTTTCTCTTTAATAAACGCTGCGGCTTTATCTGACGGAGAGGTTTCGAAACGGATATCCAGTTTCTTATCGCCAGCTATCGGTGGTAAACGCCAgttgtt                                                                                         <  1:388107/251‑1 (MQ=255)
                                                                                                                                                     tcgcatcCGCCAATAAATGCCGGATGCGGTTGATGATTTATAGGCCGAAAAAGAAGTGTTTTACTTACTCAAATCCACCTGATAAATCGCAAACCTGATATCAACGGTCGCGACTTTATTCATCGGCTACTGCCCTTTATCTTTAATAAACGCTGCGGCTTTATCTGACGGAGAGGTTTCGAAGCGGATATCCAGTTACTTATCGCCAGCTCTCGGTGCTAAACGCCAGTTGTTATCttcacccggagga                                                                         >  2:119423/1‑238 (MQ=255)
                                                                                                                                                      cgcatcCGCCAATAAATGCCGGATGCGGTTGATGATTTATAGGCCGAAAAAGAAGTGTTTTACTTACTCAAATCCACCTGATAAATCGCAAACCTGATATCATCGGTCGCGACTTTATTCATCGGATACTGCCCTTTCTCTTTAATAAACGCTGCGGCTTTATCTGACGGAGAGGTTTCGAAACGGCTATCCAGTTTCTTATCGCCAGCTATCGGTGCTAAACGCCAGTTGTTATCTGCCGCCGGGGGAAt                                                                       >  1:396669/1‑251 (MQ=255)
                                                                                                                                                         atcCGCCAATCAATGCAGGATGAGGTTGCTGATTTATAGGCCGAAAAAGAAGTGTTTTACTTACTCAACTCCACCTGATAAATCGCAAACCTGATATCATCGGTCGCGACTTTATTCATCGGATACTGCCCTTTCTCTTTAATAAACGCTGCGGCTTTATCTGACGGAGAGGTTTCGAAACGGAGATCCCGTTTCTTATCGCCAGCTATCGGTGCTAAACGCCAGTTGTTATCTGCCACCGGGGGAATTTc                                                                    >  2:389046/1‑251 (MQ=255)
                                                                                                                                                                   aaaTGCCGGATGCGGTTGATGATTTATAGGCCGAAAAAGAAGTGTTTTACTTACTCAAATCCACCTGATAAATCGCAAACCTGATATCATCGGTCGCGACTTTATTCATCGGATACTGCCCTTTCTCTTTAATAAACGCTGCGGCTTTATCTGACGGAGAGGTTTCGAAACGGATATCCAGTTTCTTATCGCCAGCTATCGGTGCTAAACGCCAGTTGTTATCTGCCGCCGGGTGAATTTCCCCCGCACg                                                           <  1:340352/250‑1 (MQ=255)
                                                                                                                                                                      ggCCGGATGGGGTTGATGATTTATAGGCCGAAAAAGAAGTGTTTTACTTACTCAAAGCCACCTGATAAATCGCAAACCTGATATCATCGGTCGCGACTTTATTCATCGGATACTGCCCTTTCTCTTTAATAAACGCTGCGGCTTTATCTGACGGAGAGGTTTCGAACCGGATATCCAGTTTCTTATCGCCAGCTATCGGTGTTAAACGCCAGTTGTTATCTGCCGCCGGGTGAATTTCCCCCGCACGTTTc                                                       <  2:222589/250‑1 (MQ=255)
                                                                                                                                                                                      ggctttATAGGCCGAAAAAGAAGGGTTTTACTTACTAACATCCACCTGATAAATCGCAAACCTGATATCATCGGTCGCGACTTTATTAATCGGATACTGCCCTTTCTCTTTAATAACCGCTGCGGCTTTATCTGACGGAGAGGTTTCGAAACGGATATCCAGTTTCTTATCGCCACCTATAGGTGCTAAACGCCAGTTGTTATCTGCCGCCGGGTGAATTTCCGCCGCACGTTTCGCCtcatacgcaaacc                                       >  2:288048/4‑242 (MQ=255)
                                                                                                                                                                                       gatTTATAGGCCGAAAAAGAAGTGTTTTACTTACTCAAATCCACCTGATAAATCGCAAACCTGATATCATCGGTCGCGACTTTATTCATCGGATACTGCCCTTTCTCTTTAATAAAc                                                                                                                                                                            <  2:343062/117‑1 (MQ=255)
                                                                                                                                                                                       gatTTATAGGCCGAAAAAGAAGTGTTTTACTTACTCAAATCCACCTGATAAATCGCAAACCTGATATCATCGGTCGCGACTTTATTCATCGGATACTGCCCTTTCTCTTTAATAAAc                                                                                                                                                                            >  1:343062/1‑117 (MQ=255)
                                                                                                                                                                                            ataGGCCGAAAAAGAAGTGTTTTACTTACTCAAATCCACCTGATAAATCGCAAACCTGATATCATCGGTCGCGACTTTATTCATCGGATACTGCCCTTTCTCTTTAATAAACGCTGCGGCTTTATCTGACGGAGAGGTTTCGAAACGGATATCCAGTTTCTTATCGCCAGCTATCGGTGCTAAACGCCAGTTGTTATCTGCCGCCGGGTGAATTTCCCCCGCACGTTTCGACTCATCAGCAATCCACGCTg                                 <  1:292578/251‑1 (MQ=255)
                                                                                                                                                                                            ataGGCCGAAAAAGAAGTGTTTTACTTACTCAAATCCACCTGATAAATCGCAAACCTGATATCATCGGTCGCGACTTTATTCATCGGATACTGCCCTTTCTCTTTAATAAACGCTGCGGCTTTATCTGACGGAGAGGTTTCGAAACGGATATCCAGTTTCTTATCGCCAGCTATCGGTGCTAAACGCCAGTTGTTATCTGCCGCCGGGTGAATTTCCCCCGCACGTTTCGACTCATCAGCAATCCACGCTg                                 <  1:218981/251‑1 (MQ=255)
                                                                                                                                                                                               ggCCGAAAAAGAAGTGTTTTACTTACTCAAATCCGCCTGATAAATCGCAAACCTGATATCATCGGGCGCGACTTTATTCATCGGATACTGCCCTTTCTCTTTAATAAACGCTGCGGCTTTATCTGACGGAGAGGTTTCGAAACGGATAGCCAGTTTCTTATCGCCAGCTATCGGTGCTAAACGCCAGTTGTTATCTGCCGCCGGGTGAATTTCCCCCGCACGTTTCGACTCATCAGCAATCCACGCTGCCa                              <  2:305430/251‑1 (MQ=255)
                                                                                                                                                                                               ggCCGAAAAAGAAGTGTTTTACTTACTCAAATCCACCTGATAAATCGCAAACCTGATATCATCGGTCGCGACTTTATTCATCGGATACTGCCCTTTCTCTTTAATAAACGCTGCGGCTTTATCTGACGGAGAGGTTTCGAAACGGATATCCAGTTTCTTATCGCCAGCTATCGGTGCTAAACGCCAGTTGTTATCTGCCGCCGTGTGAATTTCCCCCGCACGTTTCGACTCATCAGCAATCCACGCTGcc                               >  1:176346/1‑250 (MQ=255)
                                                                                                                                                                                               ggCCGAAAAAGAAGTGTTTTACTTACTCAAATCCACCTGATAAATCGCAAACCTGATATCATCGGTCGCGACTTTATTCATCGGATACTGCCCTTTCTCTTTAATAAACGCTGCGGCTTTATCTGACGGAGAGGTTTCGAAACGGATATCCAGTTTCTTATCACCAGCTATCGGTGCTAAACGCCAGTTGTTATCTGCCGCCGGGTGAATTTCCCCCGCACGTTTCGACTCATCAGCAATCCACGGTGcc                               >  2:192102/1‑250 (MQ=255)
                                                                                                                                                                                                         gaagTGTTTTACTTACTCAAATCCACCTGATAAATCGCAAACCTGATATCATCGGTCGCGACTTTATTCATCGGATACTGCCCTTTCTCTTTAATAAACGCTGCGGCTTTATCTGACGGAGAGGTTTCGAAACGGATATCCAGTTTCTTATCGCCAGCTATCGGTGCTAAACGCCAGTTGTTATCTGCCGCCGGGTGAATTTCCCCCGCACGTTTCGACTCATCAGCAATCCACGCTGCCAGCACCGAGCg                    >  1:78155/1‑251 (MQ=255)
                                                                                                                                                                                                          aagTGTTTTACTTCCTCAAAACCACCTGATAAATCGCCACCCTGATATCATCGGTCGCGACTTTATTCATCGGATACGGCCCGTTCTCTTTAATAACCGCCGCGGCTTTATCTGACGGAGAGGTGTCGAAACGGATATCCAGTTTCTTATCGCCAGCTATCGGTGCTAACCGCCAGTTGTTAGCTGCCGCCGGGTGAATTTCCCCCGCCCGTTTCGTCTCCCCAGCAAGCCTCGATGCCAGCCCCGAGCgg                   <  2:64718/251‑1 (MQ=255)
                                                                                                                                                                                                            gtgtTTTACTTATTCAAAGCCCCCTGAATAATCGCAAACCTTATATCATCGGTCGCCACTTCTTTCACCGGATCCTGCCTTTTCTCTTTAATAAACGCTGCGGCTTTATCGGACGGAGAGGTGTCGAAACGGATATCCAGTTACTTATCGGCAGCTATCTGTGCTAAACGCCAGTGGTTATCTGCCGCCGGGTGAAGTTCCCCCGCACGTTTCGACTCATCAGCAATCCACGCTGCCAGTACCGAGCTGtt                 <  2:371994/251‑1 (MQ=255)
                                                                                                                                                                                                                         tCAAATCCACCTGATAAATCGCAAACCTGATATCATCGGTCGCGACTTTATTCATCGGATACTGCCCTTTCTCTTTAATAAACGCTGCGGCTTTATCTGACGGAGAGGTTTCGAAACGGATATCCAGTTTCTTATCGCCAGCTATCGGTGCTAAACGCCAGTTGTTATCTGCCGCCGGGTGAATTTCCCCCGCACGTTTCGACTCATCAGCAATCCACGCTGCCAGCACCGAGCGGTTCTCATCCGGTGa     >  1:333292/1‑250 (MQ=255)
                                                                                                                                                                                                                            aaTCCACCTGATAAATCGCAAACCTGATATCATCGGTCGCGACTTTATTCATCGGATACTGCCCTTTCTCTTTAATAAACGCTGCGGCTTTATCTGACGGAGAGGTTTCGAAACGGATATCCAGTTTCTTATCGCCAGCTATCGGTGCTAAACGCCAGTTGTTATCTGCCGCCGGGTGAATTTCCCCCGCACGTTTCGACTCATCAGCAATCCACGCTGCCAGCACCGAGCGGTTCTCATCCGGTGAAGc  >  1:273404/1‑250 (MQ=255)
                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                 
GAAAAGGTTGCTGCACTGGCGGACTGGATTGAGTTGAATTTGCCCGAGGTGTTGGCGGTGCGGGATGAACGTGCGGCATAACTGCCGGATGCGCTGAGCTTATCTGGTCTACGGTGCATGGCTTGTAGGGCGGATAAGATGCGCCAGCATCGCATCCGCCAATAAATGCCGGATGCGGTTGATGATTTATAGGCCGAAAAAGAAGTGTTTTACTTACTCAAATCCACCTGATAAATCGCAAACCCGATATCATCGGTCGCGACTTTATTCATCGGATACTGCCCTTTCTCTTTAATAAACGCTGCGGCTTTATCTGACGGAGAGGTTTCGAAACGGATATCCAGTTTCTTATCGCCAGCTATCGGTGCTAAACGCCAGTTGTTATCTGCCGCCGGGTGAATTTCCCCCGCACGTTTCGACTCATCAGCAATCCACGCTGCCAGCACCGAGCGGTTCTCATCCGGTGAAGC  >  NC_000913/4434413‑4434882

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