Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,166,099 G→C 26.3% G163A (GGC→GCC)  PP_2780 → 3‑oxoacyl‑ACP synthase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,166,0990GC26.3% 25.8 / 13.9 23G163A (GGC→GCC) PP_27803‑oxoacyl‑ACP synthase
Reads supporting (aligned to +/- strand):  ref base G (10/7);  new base C (3/3);  total (13/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 = 1.40e-01

GTTTGCAGGACATGCTGCTCAGCCTGGCGGATAACCCGGCGCAGGCGCCGTCCTGGGCTCAATTGATGGCGCGACAGGCGCAGGTCGATCCTGACAGTTCACTGCGTCAGAGCAACGACCGCCTTGCGCGGGTGATTGCCGATGGCATCGGCAGTGAAGGCCCGGTGATCAATGTCAGCAGCGCCTGCGCGGGAGCGTCCCAGGCCATCGGCAATGCATTCCAGATGATCCGTCGCGGCGAGGTCGAGCTGGCCATTGCCGGGGGCGCCGATTCGGTGCTTAGCCTGGACACCA  >  NC_002947/3165955‑3166248
                                                                                                                                                |                                                                                                                                                     
gTTTGCAGGACATGCTGCTCAGCCTGGCGGATAACCCGGCGCAGGCGCCGTCCTGGGCTCAATTGATGGCGCGACAGGCGCAGGTCGATCCTGACAGTTCACTGCGTCAGAGCAACGACCGCCTTGCGCGGGTGATTGCCGATGGCATCg                                                                                                                                                  >  2:571840/1‑150 (MQ=255)
    gCAGGACATGCTGTTCAGCCTGGCGGATAACCCGGCGCAGGCGCCGTCCTGGGCTCAATTGATGGCGCGACAGGCGCAGGTCGATCCTGACAGTTCACTGCGTCAGAGCAACGACCGCCTTGCGCGGGTGATTGCCGATGCCATCGGCaa                                                                                                                                              <  1:127995/150‑2 (MQ=255)
    gCAGGACATGCTGCTCAGCCTGGCGGATAACCCGGCGCAGGCGCCGTCCTGGGCTCAATTGATGGCGCGACAGGCGCAGGTCGATCCTGACAGTTCACTGCGTCAGAGCAACGACCGCCTTGCGCGGGTGATTGCCGATGCCATCGGCaa                                                                                                                                              <  1:530512/150‑2 (MQ=255)
       ggACATGCTGCTCAGCCTGGCGGATAACCCGGCGAAGGCGCCGTCCTGGGCTCAATTGATGGCGCGACAGGCGCCGGTCGATCCTGACAGTTCACTGCGTCAGAGCAACGACCGCCTTGCGCGGGTGATTGCCGATGGCATCGGCAGTGa                                                                                                                                           >  1:357776/1‑150 (MQ=255)
           aTGCTGCTCAGCCTGGCGGATAACCCGGCGCAGGCGCCGTCCTGGGCTCAATTGATGGCGCGACAGGCGCAGGTCGATCCTGACAGTTCACTGCGTCAGAGCAACGACCGCCTTGCGCGGGTGATTGCCGATGGCATCGGCAGTGAAGGc                                                                                                                                       >  2:76375/1‑150 (MQ=255)
                        tGGCGGATAACCCGGCGCAGGCGCCGTCCTGGGCTCAATTGATGGCGCGACAGGCGCAGGTCGATCCTGACAGTTCACTGCGTCAGAGCAACGACCGCCTTGCGCGGGTGATTGCCGATGGCATCGGCAGTGAAGGCCCGGTGATCAATg                                                                                                                          >  1:445869/1‑150 (MQ=255)
                         ggCGGATAACCCGGCGCAGGCGCCGTCCTGGGCTCAATTGATGGCGCGACAGGCGCAGGTCGATCCTGACAGTTCACTACGTCAGAGCAACGACGACCTTGCGCATTTGATTGCCCATGGCATCGGCACTGAAGGCCCGTTGATCAATGt                                                                                                                         <  2:86558/150‑1 (MQ=255)
                                   ccGGCGCAGGCGCCGTCCTGGGCTCAATTGATGGCGCGACAGGCGCAGGTCGATCCTGACAGTTCACTGCGTCAGAGCAACGACCGCCTTGCGCGGGTGATTGCCGATGGCATCGGCAGTGAAGGCCCGGTGATCAATGTCAGCAgcgcg                                                                                                               >  2:438223/1‑149 (MQ=255)
                                                  tCCTGGGCTCAATTGATGGCGCGACAGGCGCAGGTCGATCCTGACAGTTCACTGCGTCAGAGCAACGACCGCCTTGCGCGGGTGATTGCCGATGGCATCGGCAGTGAAGGCCCGGTGATCAATGTCAGCAGCGCCTGCGCGGGAGCGTcc                                                                                                >  1:593133/1‑150 (MQ=255)
                                                         cTCAATTGATGGCGCGACAGGCGCAGGTCGATCCTGACAGTTCACTGCGTCAGAGCAACGACCGCCTTGCGCGGGTGATTGCCGATGGCATCGGCAGTGAAGGCCCGGTGATCAATGTCAGCAGCGCCTGCGCGGGAGCGTCCCAgg                                                                                            >  1:313131/1‑147 (MQ=255)
                                                         cTCAATTGATGGCGCGACAGGCGCAGGTCGATCCTGACAGTTCACTGCGTCAGAGCAACGACCGCCTTGCGCGGGTGATTGCCGATGGCATCGGCAGTGAAGGCCCGGTGATCAATGTCAGCAGCGCCTGCGCGGGAGCGTCCCAgg                                                                                            <  2:313131/147‑1 (MQ=255)
                                                                gATGGCGCGACAGGCGCACGTCGATCCTGACAGTTCACTGCGTCAGGGCAACGACCGCCTTGCGCGGGTGATTGCCGATGGCATCGGCAGGGAAGGGCAGGTGATCAATGTGAGCAGGGCCTGCGCGGGAGCGTCCCCGGCCATCGGCaa                                                                                  <  2:540099/150‑1 (MQ=255)
                                                                      gcgACAGGCGCAGGTCGATCCTGACAGTTCCTTGCGTCAGAGCAACGACCGCCTTGCGCGGGTGATTGCCGATGGCATCGGCAGTGAAGGCCCGGTGATCAATGTCAGCAGCGCCTGCGCGGGAGCGTCCCAGGCCATCGGCAATGCCtt                                                                            >  2:378531/1‑150 (MQ=255)
                                                                                                     cTGCGTCAGAGCAACGACCGCCTTGCGCGGGTGATTGCCGATGGCATCGGCAGTGAAGGCCCGGTGATCAATGTCAGCATCGCCTGCGCGGGATTGTCCCAGGCCATCGGCAATGCATTCCAGATGATCCGTCGCGGCGAGGTCGAGCTg                                             <  1:76375/150‑1 (MQ=255)
                                                                                                                          cTTGCGCGGGTGATTGCCGATGGCATCGGCAGTGAAGGCCCGGTGATCAATGTCAGCAGCGCCTGCGCGGGAGCGTCCCAGGCCATCGGCAATGCATTCCAGATGATCCGTCGCGGCGAGGTCGAGCTGGCCATTGCCGGGGGCGCCGAt                        >  2:193362/1‑150 (MQ=255)
                                                                                                                            tGCGCGGGTGATTGCCGATGGCATCGGCAGTGAAGGCCCGGTGATCAATGTCAGCAGCGCCTGCGCGGGAGCGTCCCAGGCCATCGGTAATGCATTCCAGATGATCCGTCGCGGCGAGGTCGAGCTGGCCATTGCCGGGGGCGCCGATTc                      <  1:393184/150‑1 (MQ=255)
                                                                                                                                  cttcactGCCGATGCCATCGGCAGTGAAGGCCCGGTGATCAATGTCAGCAGCGCCTGCGCGGGAGCGTCCCAGGCCATCGGCAATGCATTCCAGATGATCCGTCGCGGCGAGGTCGAGCTGGCCATTGCCGGGGGCGCCGATTCGGTGCt                >  1:129550/7‑150 (MQ=255)
                                                                                                                                   ttcactGCCGATGCCATCGGCAGTGAAGGCCCGGTGATCAATGTCAGCAGCGCCTGCGCGGGAGCGTCCCAGGCCATCGGCAATGCATTCCAGATGATCCGTCGCGGCGAGGTCGAGCTGGCCATTGCCGGGGGCGCCGATTCGGTGCtt               >  2:70210/6‑150 (MQ=255)
                                                                                                                                   gTGATTGCCGATGGCATCGGCAGTGAAGGCCCGGTGATCAATGTCAGCAGCGCCTGCGCGGGAGCGTCCCAGGCCATCGGCAATGCATTCCAGATGATCCGTCGCGGCGAGGTCGAGCTGGCCATTGCCGGGGgc                              >  2:276825/1‑135 (MQ=255)
                                                                                                                                   gTGATTGCCGATGGCATCGGCAGTGAAGGCCCGGTGATCAATGTCAGCAGCGCCTGCGCGGGAGCGTCCCAGGCCATCGGCAATGCATTCCAGATGATCCGTCGCGGCGAGGTCGAGCTGGCCATTGCCGGGGgc                              <  1:276825/135‑1 (MQ=255)
                                                                                                                                        tGCCGATGGCATCGGCAGTGAAGGCCCGGTGATAAATGTCATCAGCGCCTGCGCGGGAGCGTCCCAGGCCATCGGCAATGCATTCCAGATGATCCGTCGCGGCGAGGTCGAGCTGGCCATTGCCGGGGGCGCCGATTCGGTGCTTAGCCt          <  2:357776/150‑1 (MQ=255)
                                                                                                                                          ccGATGCCATCGGCAGTGAAGGCCCGGTGATCAATGTCAGCAGCGCCTGCGCGGGAGCGTCCCAGGCCATCGGCAATGCATTCCAGATGATCCGTCGCGGCGAGGTCGAGCTGGCCATTGCCGGGGCGCCGATTCGGTGCTTAGCCTGGa       >  1:548312/1‑150 (MQ=255)
                                                                                                                                           cGATGCCATCGGCAGTGAAGGCCCGGTGATCAATGTCAGCAGCGCCTGCGCGGGAGCGTCCCAGGCCATCGGCAATGCATTCCAGATGATCCGTCGCGGCGAGGTCGAGCTGGCCATTGCCGGGGGCGCCGATTCGGTGCTTAGCCTGGa       <  2:129550/150‑1 (MQ=255)
                                                                                                                                                gCATCGGCAGTGAAGGCCCGCTGATCAATGTCAGCAGCGCCTGCGCGGGAGCGTCCCAGGCCATCCGCAATGCATTCCAGATGATCCGTCGCGGCGAGGTCGAGCTGGCCATTGCCGGGGGCGCCGATTCGGTGCTTAGCCTGGACACCa  <  2:217203/150‑1 (MQ=255)
                                                                                                                                                |                                                                                                                                                     
GTTTGCAGGACATGCTGCTCAGCCTGGCGGATAACCCGGCGCAGGCGCCGTCCTGGGCTCAATTGATGGCGCGACAGGCGCAGGTCGATCCTGACAGTTCACTGCGTCAGAGCAACGACCGCCTTGCGCGGGTGATTGCCGATGGCATCGGCAGTGAAGGCCCGGTGATCAATGTCAGCAGCGCCTGCGCGGGAGCGTCCCAGGCCATCGGCAATGCATTCCAGATGATCCGTCGCGGCGAGGTCGAGCTGGCCATTGCCGGGGGCGCCGATTCGGTGCTTAGCCTGGACACCA  >  NC_002947/3165955‑3166248

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