Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,824,314 C→G 100% A279P (GCG→CCG)  kguT ← 2‑ketogluconate transporter, putative

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,824,3140CG87.5% 39.4 / ‑0.2 16A279P (GCG→CCG) kguT2‑ketogluconate transporter, putative
Reads supporting (aligned to +/- strand):  ref base C (2/0);  new base G (5/9);  total (7/9)
Fisher's exact test for biased strand distribution p-value = 1.75e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.51e-01
Rejected as polymorphism: E-value score below prediction cutoff.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

TGCTCGCTGCCCAAGGCGTAGGAACCATAGAACGCCAACGCGGCAATCAGCAGCGGCGGCCAGACGAAGCGCTTGCGCTTTTGCAGGCGGTCGGACGCCCAGGACACGCCCACCATTGCCAGCACCGCCGCCAGGTAAGGCACCGAGGCCAGCCAGCCGGCCTGCACGATATCGACGTTCGCCGCCTGCTTGAGGATCGACGGCAGCCACAGCACGAAGCCATACACGCCGATGCTCCAGCAGAAGTACTGCAGCGCCAGGAT  >  NC_002947/3824184‑3824446
                                                                                                                                  |                                                                                                                                    
tGCTCGCTGCCCAAGGCGTAGGAACCATAGAACGCCAACGCGGCAATCAGCAGCGGCGGCCAGACGAAGCGCTTGCGCTTTTGCAGGCGGTCGGACGCCCAGGACACGCCCACCATTGCCAGCACCGCCGGCAGGt                                                                                                                                 <  2:485801/136‑1 (MQ=255)
                    ggAACCATAGAACGCCAACGCGGCAATCAGCAGCGGCGGCCAGACGAAGCGCTTGCGCTTTTGCAGGCGGTCGGACGCCCAGGACACGCCCACCATTGCCAGCACCGCCGGCAGGTAAGGCACCGAGGccagccag                                                                                                             <  1:217196/136‑1 (MQ=255)
                      aaCCATAGAACGCCAACGCGGCAATCAGCAGCGGCGGCCAGACGAAGCGCTTGCGCTTTTGCAGGCGGTCGGACGCCCAGGACACGCCCACCATTGCCAGCACCGCCGGCAGGTAAGGCACCGAGGCCagccagcc                                                                                                           >  1:164503/1‑136 (MQ=255)
                      aaCCATAGAACGCCAACGCGGCAATCAGCAGCGGCGGCCAGACGAAGCGCTTGCGCTTTTGCAGGCGGTCGGACGCCCAGGACACGCCCACCATTGCCAGCACCGCCGGCAGGTAAGGCACCGAGGCCagccagcc                                                                                                           >  2:270219/1‑136 (MQ=255)
                      aaCCATAGAACGCCAACGCGGCAATCAGCAGCGGCGGCCAGACGAAGCGCTTGCGCTTATGCAGGCGGTCGGACGCCCAGGACACGCCCACCATTCCCAGCACCGCCGGCAGGTAAGGCACCGAGGCCagccagcc                                                                                                           >  1:169691/1‑136 (MQ=255)
                                          gCAATCAGCAGCGGCGGCCAGACGAAGCGCTTGCGCTTTTGCAGGCGGTCGGACGCCCAGGACACGCCCACCATTGCCAGCACCGCCGGCAGGTAAGGCACCGAGGCCAGCCAGCCGGCCTGCACGATATCGACGt                                                                                       <  1:153301/136‑1 (MQ=255)
                                                             aGACGAAGCGCTTGCGCTTTTGCAGGCGGTCGGACGCCCAGGACACGCCCACCATTGCCAGCACCGCCGGCAGGTAAGGCACCGAGGCCAGCCAGCCGGCCTGCACGATATCGACGTTCGCCGCCTGCTTGAGGAt                                                                    >  1:79470/1‑136 (MQ=255)
                                                                   tgcgcttgggcttTTGCTGGCGGTCGGACGCCCTGGTCACGCCCACCATTGCCAGCACGGCCGGCAGGTAAGGCACGGAGGCCAGCCTGCCGGCCTGCACGTTATCGACGTTCGCCGCCTGCTTGAGGATCGACgg                                                              <  1:58874/135‑1 (MQ=255)
                                                                      gcttgcgcttTTGCAGGCGGTCGGACGCCCAGGACACGCCCACCATTGCCAGCACCGCCGGCAGGTAAGGCACCGAGGCCAGCCAGCCGGCCTGCACGATATCGACGTTCGCCGCCTGCTTGAGGATCGACGGCAg                                                           <  2:248526/136‑1 (MQ=255)
                                                                        ttgcgcttTTGCAGGCGGTCGGACGCCCAGGACACGCCCACCATTGCCAGCACCGCCGGCAGGTAAGGCACCGAGGCCAGCCAGCCGGCCTGCACGATATCGACGTTCGCCGCCTGCTTGAGGATCGACGGCAGcc                                                         <  2:448396/136‑1 (MQ=255)
                                                                             cttTTGCAGGCGGTCGGACGCCCAGGACACGCCCACCATTGCCAGCACCGCCGGCAGGTAAGGCACCGAGGCCAGCCAGCCGGCCTGCACGATATCGACGTTCGCCGCCTGCTTGAGGATCGACGGCAGCCACAGc                                                    <  1:161/136‑1 (MQ=255)
                                                                                ttGCAGGCGGTCGGACGCCCAGGACACGCCCACCATTGCTAGCACCGCCGGCAGGTAAGGCACCGAGGCCAGCCAGCCGGCCTGCACGATATCGACGTTCGCCGCCTGCTTGAGGATCGACGGCAGCCACAGCACg                                                 <  2:117059/136‑1 (MQ=255)
                                                                                                   cAGCACACGCCCACCATTGCCAGCACCGCCGGCAGGTAAGGCACCGAGCCCAGCCAGCCGGCCTGCACGATATCGACGTTCGCCGCCTGCTTGAGGATCGACGGCAGCCACAGCACGAAGCCATACACGCCGATGc                              <  1:33471/136‑1 (MQ=255)
                                                                                                              caccaTTGCCAGCACCGCCGCCAGGTAAGGCACCGAGGCCAGCCAGCCGGCCTGCACGATATCGACGTTCGCCGCCTGCTTGAGGATCGACGGCAGCCACAGCACGAAGCCATACACGCCGGTGCTCCAGCAGAAg                   >  2:22258/1‑136 (MQ=255)
                                                                                                                   ttGCCAGCACCGCCGGCAGGTAAGGCACCGAGGCCAGCCAGCCGGCCTGCACGATATCGACGTTCGCCGCCTGCTTGAGGATCGACGGCAGCCACAGCACGAAGCCATACACGCCGATGCTCCAGCAGAAGTACTg              >  2:19166/1‑136 (MQ=255)
                                                                                                                      ccAGCACCGCCGCCAGGTAAGGCACCGAGGCCAGCCAGCCGGCCTGCACGATATCGACGTTCGCCGCCTGCTTGAGGATCGACGGCAGCCACAGCACGAAGCCATACACGCCGGTGCTCCAGCAGAAGTACTGCAg           >  1:474372/1‑136 (MQ=255)
                                                                                                                               ccggcAGGTAAGGCACCGAGGCCAGCCAGCCGGCCTGCACGATATCGACGTTCGCCGCCTGCTTGAGGATCGACGGCAGCCACAGCACGAAGCCATACACGCCGATGCTCCAGCAGAAGTACTGCAGCGCCAGgat  <  2:130708/136‑1 (MQ=255)
                                                                                                                                  |                                                                                                                                    
TGCTCGCTGCCCAAGGCGTAGGAACCATAGAACGCCAACGCGGCAATCAGCAGCGGCGGCCAGACGAAGCGCTTGCGCTTTTGCAGGCGGTCGGACGCCCAGGACACGCCCACCATTGCCAGCACCGCCGCCAGGTAAGGCACCGAGGCCAGCCAGCCGGCCTGCACGATATCGACGTTCGCCGCCTGCTTGAGGATCGACGGCAGCCACAGCACGAAGCCATACACGCCGATGCTCCAGCAGAAGTACTGCAGCGCCAGGAT  >  NC_002947/3824184‑3824446

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