Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,018,022 C→A 24.9% R634R (CGC→CGA bcsB → cellulose synthase and translocator subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,018,0220CA24.9% 25.8 / 3.1 16R634R (CGC→CGAbcsBcellulose synthase and translocator subunit
Reads supporting (aligned to +/- strand):  ref base C (6/6);  new base A (2/2);  total (8/8)
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.59e-01

CCTGCGCGACGCGCCCGACCTTGGCCTGTTGCTCAGCGCCCAGCGTGACTGGTTGCTGCAGGGCCGCAGCGCATCGTTGCCCGGTGGCCAGCGTTTCGATACTGAGCCGGTGGCTGCAAGCAGCCGCGTGGCGGTCAGCGCCCAGGCGCCGATTGCGGCCATCACCGGGCTGAAGTCGCCGTTCCACGAGCAGCGCAGTGTTGTGGCCTTGCTGGCCAACAGCGACAGCGATTACGCATTGCTGCGCGACATGCTTG  >  NC_002947/3017896‑3018152
                                                                                                                              |                                                                                                                                  
ccTGCGCGACGCGCCCGACCTTGGCCTGTTGCTCAGCGCCCAGCGTGACTGGTTGCTGCAGGGCCGCAGCGCATCGTTGCCCGGTGGCCAGCGTTTCGATACTGAGCCGGTGGCTGCAAGCAGCCGCGTGGCGGTCAGCGCCCAGGCGcc                                                                                                             >  2:37466/1‑150 (MQ=255)
    cgcgACGCGCCCGACCTTGGCCTGTTGCTCAGCGCCCAGCGTGACTGGTTGCTGCAGGGCCGCAGCGCATCGTTGCCCGGTGGCCAGCGTTTCGATACTGAGCCGGTGGCTGCAAGCAGCCGCGTGGCGGTCAGCGCCCAGGCGCCGAtt                                                                                                         <  1:79843/150‑1 (MQ=255)
      cgACGCGCCCGACCTTGGCCTGTTGCTCAGCGCCCAGCGTGACTGGTTGCTGCAGGGCCGCAGCGCATCGTTGCCCGGTGGCCAGCGTTTCGATACTGAGCCGGTGGCTGCAAGCAGCCGCGTGGCGGTCAGCGCCCAGGCGCCGATTGc                                                                                                       >  1:319/1‑150 (MQ=255)
                   cTTGGCCTGTTGCTCAGCGCCCAGCGTGACTGGTTGCTGCAGGGCCGCAGCGCATCGTTGCCCGGGGGCCAGCGTTTCGATACTGAGCCGGTGGCTGCAAGCAGCCGCGTGGCGGTCAGCGCCCAGGCGCCGATTGCGGCCATCACCggg                                                                                          >  2:7304/1‑150 (MQ=255)
                           gTTGCTCAGCGCCCAGCGTGACTGGTTGCTGCAGGGCCGCAGCGCATCGTTGCCCGGTGGCCAGCGTTTCGATACTGAGCCGGTGGCTGCAAGCAGCCGCGTGGCGGTCAGCGCCCAGGCGCCGATTGCGGCCATCACCGGGCTGAAGTc                                                                                  <  2:358101/150‑1 (MQ=255)
                                      cccAGCGTGACTGGTTGCTGCAGGGCCGCAGCGCATCGTTGCCCGGTGGCCAGCGTTTCGATACTGAGCCGGTGGCTGCAAGCAGCCGAGTGGCGGTCAGCGCCCAGGCGCCGATTGCGGCCATCACCGGGCTGAAGTCGCCGTTCCACg                                                                       >  1:217642/1‑150 (MQ=255)
                                       ccAGCGTGACTGGTTGCTGCAGGGCCGCAGCGCATCGTTGCCCGGTGGCCAGCGTTTCGATACTGAGCCGGTGGCTGCAAGCAGCCGCGTGGCGGTCAGCGCCCAGGCGCCGATTGCGGCCATCACCGGGCTGaa                                                                                     >  1:70012/1‑135 (MQ=255)
                                       ccAGCGTGACTGGTTGCTGCAGGGCCGCAGCGCATCGTTGCCCGGTGGCCAGCGTTTCGATACTGAGCCGGTGGCTGCAAGCAGCCGCGTGGCGGTCAGCGCCCAGGCGCCGATTGCGGCCATCACCGGGCTGaa                                                                                     <  2:70012/135‑1 (MQ=255)
                                                 tGGTTGCTGCAGGGCCGCAGCGCATCGTTGCCCGGTGGCCAGCGTTTCGATACTGAGCCGGTGGCTGCAAGCAGCCGAGTGGCGGTCAGCGCCCAGGCGCCGATTGCGGCCATCACCGGGCTGAAGTCGCCGTTCCACGAGCAGCGCAgt                                                            >  1:360111/1‑150 (MQ=255)
                                                          cAGGGCCGCAGCGCATCGTTGCCCGGTGGCCAGCGTTTCGATACTGAGCCGGTGGCTGCAAGCAGCCGCGTGGCGGTCAGCGCCCAGGCGCCGATTGCGGCCATCACCGGGCTGAAGTCGCCGTTCCACGAGCAGCGCAGTGTTGTGGcc                                                   <  2:319/150‑1 (MQ=255)
                                                                cGCAGCGCATCGTTGCCCGGTGGCCAGCGTTTCGATACTGAGCCGGTGGCTGCAAGCAGCCGCGTGGCGGTCAGCGCCCAGGCGCCGATTGCGGCCATCACCGGGCTGAAGTCGCCGTTCCACGAGCAGCGCAGTGTTGTGGCCTTGCTg                                             <  1:7304/150‑1 (MQ=255)
                                                                       cATCGTTGCCCGGTGGCCAGCGTTTCGATACTGAGCCGGTGGCTGCAAGCAGCCGAGTGGCGGTCAGCGCCCAGGCGCCGATTGCGGCCATCACCGGGCTGAAGTCGCCGTTCCACGAGCAGCGCAGTGTTGTGGCCTTGCTGGCCAaca                                      <  2:360111/150‑1 (MQ=255)
                                                                                              ttCGATACTGAGCCGGTGGCTGCAAGCAGCCGAGTGGCGGTCAGCGCCCAGGCGCCGATTGCGGCCATCACCGGGCTGAAGTCGCCGTTCCACGAGCAGCGCAGTGTTGTGGCCTTGCTGGCCAACAGCGACAGCGATTACGCATtgctg               <  2:217642/150‑1 (MQ=255)
                                                                                                          ccGGTGGCTGCAAGCAGCCGCGTGGCGGTCAGCGCCCAGGCGCCGATTGCGGCCATCACCGGGCTGAAGTCGCCGTTCCACGAGCAGCGCAGTGTTGTGGCCTTGCTGGCCAACAGCGACAGCGATTACGCATTGCTGCGCGACATGCtt   >  1:67199/1‑150 (MQ=255)
                                                                                                               ggCTGCAAGCAGCCGCGTGGCGGTCAGCGCCCAGGCGCCGATTGCGGCCATCACCGGGCTGAAGTCGCCGTTCCACGAGCAGCGCAGTGTTGTGGCCTTGCTGGCCAACAGCGACAGCGATTACGCATTGCTGCGCGACATGCTTg  <  1:368712/146‑1 (MQ=255)
                                                                                                               ggCTGCAAGCAGCCGCGTGGCGGTCAGCGCCCAGGCGCCGATTGCGGCCATCACCGGGCTGAAGTCGCCGTTCCACGAGCAGCGCAGTGTTGTGGCCTTGCTGGCCAACAGCGACAGCGATTACGCATTGCTGCGCGACATGCTTg  >  2:368712/1‑146 (MQ=255)
                                                                                                                              |                                                                                                                                  
CCTGCGCGACGCGCCCGACCTTGGCCTGTTGCTCAGCGCCCAGCGTGACTGGTTGCTGCAGGGCCGCAGCGCATCGTTGCCCGGTGGCCAGCGTTTCGATACTGAGCCGGTGGCTGCAAGCAGCCGCGTGGCGGTCAGCGCCCAGGCGCCGATTGCGGCCATCACCGGGCTGAAGTCGCCGTTCCACGAGCAGCGCAGTGTTGTGGCCTTGCTGGCCAACAGCGACAGCGATTACGCATTGCTGCGCGACATGCTTG  >  NC_002947/3017896‑3018152

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