Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,644,438 G→T 15.3% E153D (GAG→GAT)  PP_3206 → hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,644,4380GT15.3% 81.8 / 2.0 26E153D (GAG→GAT) PP_3206hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (13/9);  new base T (2/2);  total (15/11)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.82e-01

TCTGTCCACCGGTGCTGGTGTACGCCAGCAGCGTCGCGGTGTACGGCGGCAAGTTGCCGACGCGCATGGACGAAGGCCAGCCAGCGTCGCCGCAGCTGTCCTATGCGGCGCACAAACGCATGGTGGAGATCGCGTTGCAGGATCTTGCCCGCCGTGGCGAGGTTGACGGCCGCGCCTTGCGTTTGCCCGGTATCGTTGCTCGCCCGCGCGAGCCGAATGGTCTGCGCTCGGCGTTCATGAGCGACCTGTT  >  NC_002947/3644311‑3644560
                                                                                                                               |                                                                                                                          
tctGTCCACCGGTGCTGGTGTACGCCAGCAGCGTCGCGGTGTACGGCGGCAAGTTGCCGACGCGCATGGACGAAGGCCAGCCAGCGTCGCCGCAGCTGTCCTATGCGGCGCACAAACGCATGGTGGAGATCGCGtt                                                                                                                    <  1:213283/136‑1 (MQ=255)
tctGTCCACCGGTGCTGGTGTACGCCAGCAGCGTCGCGGTGTACGGCGGCAAGTTGCCGACGCGCATGGACGAAGGCCAGCCAGCGTCGCCGCAGCTGTCCTATGCGGCGCACAAACGCATGGTGGAGATCGCGtt                                                                                                                    <  1:215421/136‑1 (MQ=255)
   gTCCACCGGTGCTGGTGTACGCCAGCAGCGTCGCGGTGTACGGCGGCAAGTTGCCGACGCGCATGGACGAAGGCCAGCCAGCGTCGCCGCAGCTGTCCTATGCGGCGCACAAACGCATGGTGGAGATCGCGTTGCa                                                                                                                 >  1:368901/1‑136 (MQ=255)
          ggTGCTGGTGTACGCCAGCAGCGTCGCGGTGTACGGCGGCAAGTTGCCGACGCGCATGGACGAAGGCCAGCCAGCGTCGCCGCAGCTGTCCTATGCGGCGCACAAACGCATGGTGGAGATCGCGTTGCAGGATCtt                                                                                                          >  2:7050/1‑136 (MQ=255)
             gCTGGTGTACGCCAGCAGCGTCGCGGTGTACGGCGGCAAGTTGCCGACGCGCATGGACGAAGGCCAGCCAGCGTCGCCGCAGCTGTCCTATGCGGCGCACAAACGCATGGTGGAGATCGCGTTGCAGGATCTTGcc                                                                                                       <  2:139880/136‑1 (MQ=255)
                                    cgGTGTACGGCGGCAAGTTGCCGACGCGCATGGACGAAGGCCAGCCAGCGTCGCCGCAGCTGTCCTATGCTGCGCACAAACGCATGGTGGATATCGCGTTGCAGGATCTTGCCCGCCGTGGCGAGGTTGACGGCcg                                                                                >  2:248160/1‑136 (MQ=255)
                                    cgGTGTACGGCGGCAAGTTGCCGACGCGCATGGACGAAGGCCAGCCAGCGTCGCCGCAGCTGTCCTATGCTGCGCACAAACGCATGGTGGATATCGCGTTGCAGGATCTTGCCCGCCGTGGCGAGGTTGACGGCcg                                                                                >  2:246005/1‑136 (MQ=255)
                                    cgGTGTACGGCGGCAAGTTGCCGACGCGCATGGACGAAGGCCAGCCAGCGTCGCCGCAGCTGTCCTATGCGGCGCACAAACGCATGGTGGAGATCGCGTTGCAGGATCTTGCCCGCCGTGGCGAGGTTGACGGCcg                                                                                >  1:285952/1‑136 (MQ=255)
                                                gcAAGTTGCCGACGCGCATGGACGAAGGCCAGCCAGCGTCGCCGCAGCTGTCCTATGCTGCGCACAAACGCATGGTGGATATCGCGTTGCAGGATCTTGCCCGCCGTGGCGAGGTTGACGGCCGCGCCTTGCGttt                                                                    <  1:246005/136‑1 (MQ=255)
                                                gcAAGTTGCCGACGCGCATGGACGAAGGCCAGCCAGCGTCGCCGCAGCTGTCCTATGCTGCGCACAAACGCATGGTGGATATCGCGTTGCAGGATCTTGCCCGCCGTGGCGAGGTTGACGGCCGCGCCTTGCGttt                                                                    <  1:248160/136‑1 (MQ=255)
                                                    gTTGCCGACGCGCATGGACGAAGGCCAGCCAGCGTCGCCGCAGCTGTCCTATGCGGCGCACAAACGCATGGTGGAGATCGCGTTGCAGGATCTTGCCCGCCGTGGCGAGGTTGACGGCCGCGCCTTGCGTTTGccc                                                                >  1:511377/1‑136 (MQ=255)
                                                    gTTGCCGACGCGCATGGACGAAGGCCAGACAGCGTCGACGCAGCTGTCCTATGCGGCGCACAATCGCATGCTGGAGATCGCGTAGCAGCATCTTGCCCGCCGTGGCGAGGTTGACGGCCGCGCCTTGCGTTTGcgc                                                                >  1:388935/1‑134 (MQ=255)
                                                      tGCCGACGCGCATGGACGAAGGCCAGCCAGCGTCGCCGCAGCTGTCCTATGCGGCGCACAAACGCATGGTGGAGATCGCGTTGCAGGATCTTGCCCGCCGTGGCGAGGTTGACGGCCGCGCCTTGCGTTTGCCCgg                                                              >  1:546277/1‑136 (MQ=255)
                                                           aCGCGCATGGACGAAGGCCAGCCAGCGTCGCCGCAGCTGTCCTATGCGGCGCACAAACGCATGGTGGAGATCGCGTTGCAGGATCTTGCCCGCCGTGGCGAGGTTGACGGCCGCGCCTTGCGTTTGCCCGGTATCg                                                         <  1:329/136‑1 (MQ=255)
                                                           aCGCGCATGGACGAAGGCCAGCCAGCGTCGCCGCAGCTGTCCTATGCGGCGCACAAACGCATGGTGGAGATCGCGTTGCAGGATCTTGCCCGCCGTGGCGAGGTTGACGGCCGCGCCTTGCGTTTGCCCGGTATCg                                                         <  1:232490/136‑1 (MQ=255)
                                                                                  agcGTCGCCGCAGCTGTCCTATGCGGCGCACAAACGCATGGTGGAGATCGCGTTGCAGGATCTTGCCCGCCGTGGCGAGGTTGACGGCCGCGCCTTGCGTTTGCCCGGTATCGTTGCTCGCCCGCGCGAGCCGAAt                                  >  1:462370/1‑136 (MQ=255)
                                                                                                          ggcgCACAAACGCATGGTGGAGATCGCGTTGCAGGATCTTGCCCGCCGTGGCGAGGTTGACGGCCGCGCCTTGCGTTTGCCCGGTATCGTTGCTCGCCCGCGCGAGCCGAATGGTCTGCGCTCGGCGTTCATGAGc          <  1:320921/136‑1 (MQ=255)
                                                                                                          ggcgCACAAACGCATGGTGGAGATCGCGTTGCAGGATCTTGCCCGCCGTGGCGAGGTTGACGGCCGCGCCTTGCGTTTGCCCGGTATCGTTGCTCGCCCGCGCGAGCCGAATGGTCTGCGCTCGGCGTTCATGAGc          >  2:320921/1‑136 (MQ=255)
                                                                                                          ggcgCACAAACGCATGGTGGAGATCGCGTTGCAGGATCTTGCCCGCCATGGCGAGGTTGACGGCCGCGCCTTGCGTTTGCCCGGTATCGTTGCTCGCCCGCGCGAGCCGAATGGTCTGCGCTCGGCGTTCATGAGc          <  1:286024/136‑1 (MQ=255)
                                                                                                               acaAACGCATGGTGGAGATCGCGTTGCAGGATCTTGCCCGCCGTGGCGAGGTTGACGGCCGCGCCTTGCGTTTGCCCGGTATCGTTGCTCGCCCGCGCGAGCCGAATGGTCTGCGCTCGGCGTTCATGAGCGACCt     >  1:296677/1‑136 (MQ=255)
                                                                                                               acaAACGCATGGTGGAGATCGCGTTGCAGGATCTTGCCCGCCGTGGCGAGGTTGACGGCCGCGCCTTGCGTTTGCCCGGTATCGTTGCTCGCCCGCGCGAGCCGAATGGTCTGCGCTCGGCGTTCATGAGCGACCt     >  1:260443/1‑136 (MQ=255)
                                                                                                               acaAACGCATGGTGGAGATCGCGTTGCAGGATCTTGCCCGCCGTGGCGAGGTTGACGGCCGCGCCTTGCGTTTGCCCGGTATCGTTGCTCGCCCGCGCGAGCCGAATGGTCTGCGCTCGGCGTTCATGAGCGACCt     >  1:16146/1‑136 (MQ=255)
                                                                                                               acaAACGCATGGTGGAGATCGCGTTGCAGGATCTTGCCCGCCGTGGCGAGGTTGACGGCCGCGCCTTGCGTTTGCCCGGTATCGTTGCTCGCCCGCGCGAGCCGAATGGACTGCGCTCGGCGTGCATGAGCGACCt     >  1:60372/1‑136 (MQ=255)
                                                                                                                caAACGCATGGTGGAGATCGCGTTGCAGGATCTTGCCCGCCGTGGCGAGGTTGACGGCCGCGCCTTGCGTTTGCCCGGTATCGTTGCTCGCCCGCGCGAGCCGAATGGTCTGCGCTCGGCGTTCATGAGCGACCTg    <  2:368901/136‑1 (MQ=255)
                                                                                                                 aaaCGCATGGTGGAGATCGCGTTGCAGGATCTTGCCCGCCGTGGCGAGGTTGACGGCCGCGCCTTGCGTTTGCCCGGTATCGTTGCTCGCCCGCGCGAGCCGAATGGTCTGCGCTCGGCGTTCATGAGCGACCTGt   >  1:227679/1‑136 (MQ=255)
                                                                                                                  aaCGCATGGTGGAGATCGCGTTGCAGGATCTTGCCCGCCGTGGCGAGGTGGACGGCCGCGCCTTGCGTTTGCCCGGTATCGTTGCTCGCCCGCGCGAGCCGAATGGTCTGCGCTCGGCGTTCATGAGCGACCTGtt  <  1:160737/136‑1 (MQ=255)
                                                                                                                               |                                                                                                                          
TCTGTCCACCGGTGCTGGTGTACGCCAGCAGCGTCGCGGTGTACGGCGGCAAGTTGCCGACGCGCATGGACGAAGGCCAGCCAGCGTCGCCGCAGCTGTCCTATGCGGCGCACAAACGCATGGTGGAGATCGCGTTGCAGGATCTTGCCCGCCGTGGCGAGGTTGACGGCCGCGCCTTGCGTTTGCCCGGTATCGTTGCTCGCCCGCGCGAGCCGAATGGTCTGCGCTCGGCGTTCATGAGCGACCTGTT  >  NC_002947/3644311‑3644560

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