Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 221,321 G→A 36.4% E31E (GAG→GAA PP_0169 → TauD/TfdA family dioxygenase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947221,3210GA36.4% 5.9 / 6.2 11E31E (GAG→GAAPP_0169TauD/TfdA family dioxygenase
Reads supporting (aligned to +/- strand):  ref base G (3/4);  new base A (2/2);  total (5/6)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GTTGCGCGCCATTGAGGTCGCGCATTGCGTGGGAGTGAAAAATGTCAGCCGCCTTGAACGCCTTGTCGTCCATCGACAGCGCCCAGCCGCAAGTTTTTGAAGTTCGCCCCTTCTCGGGTGCCGTCGGCGCCGAGATCTTAGGCCTGGACCTGGCCAAACCGGTCAGCGCCGAGGATTTCACCCGTATCCATCGCGCCCACCTCGACC  >  NC_002947/221188‑221394
                                                                                                                                     |                                                                         
gTTGCGCGCCATTGAGGTCGCGCATTGCGTGGGAGTGAAAAATGTCAGCCGCCTTGAACGCCTTGTCGTCCATCGACAGCGCCCAGCCGCAAGTTTTTGAAGTTCGCCCCTTCTCGGGTGCCGTCGGCGCCGAGAt                                                                         >  2:378463/1‑136 (MQ=255)
     gcgcCATTGAGGTCGCGCATTGCGTGGGAGTGAAAAATGTCAGCCGCCTTGAACGCCTTGTCGTCCATCGACAGCGCCCAGCCGCAAGTTTTTGAAGTTCGCCCCTTCTCGGGTGCCGTCGGCGCCGAGATCTTAg                                                                    <  1:108399/136‑1 (MQ=255)
     gcgcCATTGAGGTCGCGCATTGCGTGGGAGTGAAAAATGTCAGCCGCCTTGAACGCCTTGTCGTCCATCGACAGCGCCCAGCCGCAAGTTTTTGAAGTTCGCCCCTTCTCGGGTGCCGTCGGCGCCGAGATCTTAg                                                                    <  1:486536/136‑1 (MQ=255)
                      cATTGCGTGGGAGTGAAAAATGTCAGCCGCCTTGAACGCCTTGTCGTCCATCGACAGCGCCCAGCCGCAAGTTTTTGAAGTTCGCCCCTTCTCGGGTGCCGTCGGCGCCGAAATCTTAGGCCTg                                                               >  1:110177/1‑124 (MQ=255)
                      cATTGCGTGGGAGTGAAAAATGTCAGCCGCCTTGAACGCCTTGTCGTCCATCGACAGCGCCCAGCCGCAAGTTTTTGAAGTTCGCCCCTTCTCGGGTGCCGTCGGCGCCGAAATCTTAGGCCTg                                                               >  1:112859/1‑124 (MQ=255)
                      cATTGCGTGGGAGTGAAAAATGTCAGCCGCCTTGAACGCCTTGTCGTCCATCGACAGCGCCCAGCCGCAAGTTTTTGAAGTTCGCCCCTTCTCGGGTGCCGTCGGCGCCGAAATCTTAGGCCTg                                                               <  2:110177/124‑1 (MQ=255)
                      cATTGCGTGGGAGTGAAAAATGTCAGCCGCCTTGAACGCCTTGTCGTCCATCGACAGCGCCCAGCCGCAAGTTTTTGAAGTTCGCCCCTTCTCGGGTGCCGTCGGCGCCGAAATCTTAGGCCTg                                                               <  2:112859/124‑1 (MQ=255)
                            gTGGGAGTGAAAAATGTCAGCCGCCTTGAACGCCTTGTCGTCCATCGACAGCGCCCAGCCGCAAGTTTTTGAAGTTCGCCCCTTCTCGGGTGCCGTCGGCGCCGAGATCTTAGGCCTGGACCTGGCCAAACCGGTc                                             >  2:126674/1‑136 (MQ=255)
                                                  gccTTGAACGCCTTGTCGTCCATCGACAGCGCCCAGCCGCAAGTTTTTGAAGTTCGCCCCTTCTCGGGTGCCGTCGGCGCCGAGATCTTAGGCCTGGACCTGGCCAAACCGGTCAGCGCCGAGGATTTCCCCCGtt                       >  2:531166/1‑135 (MQ=255)
                                                                  cgtcCATCGACAGCGCCCAGCCGCAAGTTTTTGAAGTTCGCCCCTTCTCGGGTGCCGTCGGCGCCGAGATCTTAGGCCTGGACCTGGCCAAACCGGTCAGCGCCGAGGATTTCACCCGTATCCATCGCGCCCACCt       <  1:33011/136‑1 (MQ=255)
                                                                       aTCGACAGCGCCCAGCCGCAAGTTTTTGAAGTTCGCCCCTTCTCGGGTGCCGTCGGCGCCGAGATCTTAGGCCTGGACCTGGCCAAACCGGTCAGCGCCGAGGATTTCACCCGTATCCATCGCGCCCACCTCGacc  <  1:378463/136‑1 (MQ=255)
                                                                                                                                     |                                                                         
GTTGCGCGCCATTGAGGTCGCGCATTGCGTGGGAGTGAAAAATGTCAGCCGCCTTGAACGCCTTGTCGTCCATCGACAGCGCCCAGCCGCAAGTTTTTGAAGTTCGCCCCTTCTCGGGTGCCGTCGGCGCCGAGATCTTAGGCCTGGACCTGGCCAAACCGGTCAGCGCCGAGGATTTCACCCGTATCCATCGCGCCCACCTCGACC  >  NC_002947/221188‑221394

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