Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 677,954 C→T 54.5% G356D (GGC→GAC)  PP_0581 ← 3‑oxoacyl‑ACP reductase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947677,9540CT54.5% 6.7 / 53.3 44G356D (GGC→GAC) PP_05813‑oxoacyl‑ACP reductase
Reads supporting (aligned to +/- strand):  ref base C (7/13);  new base T (12/12);  total (19/25)
Fisher's exact test for biased strand distribution p-value = 3.72e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.32e-01

TGGGCATGGCTGCAGTCATGTGGGTTTCGATGAAACCGGGAGCGATGGCGTTGATGCTGCCACCGCGCTCGGCGAGCCGTGGCGCCCAGGCCTGGGCCAAACCGATCAGCCCGGCCTTGCTCGCGGCATAGTTGGCCTGCCCGCGGTTGCCGGCAATGCCGCTGACCGAGGCCAGCAGGGTGATGCGCGCGTTTTCCCCCAGCGCACCGTTGTCATACAACGCCTGGGTGAGCACTTGCGGTGCCTTGAGGTTCACGGCCAGCACGGCGTCCCAGTATTCCGGCGT  >  NC_002947/677805‑678090
                                                                                                                                                     |                                                                                                                                        
tGGGCATGGCTGCAGTCATGTGGGTTTCGATGAAACCGGGAGCGATGGCGTTGATGCTGCCACCGCGCTCGGCGAGCCGTGGCGCCCAGGCCTGGGCCAAACCGATCAGCCCGGCCTTGCTCGCGGCATAGTTGGCCTGCCCGCGGTTGc                                                                                                                                          <  2:610615/150‑1 (MQ=255)
  ggCATGGCTGCCGTCATGTGGGTTCCGATGAAACCGGGAGCGATGGCGTTGATGCTGCCAACGCGCTCGGCGAGCCGTGGCGCCCAGGCCTGGGCCAAACCGATCAGCCCGGCCTTGCTCGCGGCATAGTTGGCCTGCCCGCGGTTGTCg                                                                                                                                        <  1:491931/150‑1 (MQ=255)
       ggCTGCAGTCATGTGGGTTTCGATGAAACCGGGAGCGATGGCGTTGATGCTGCCACCGCGCTCGGCGAGCCGTGGCGCCCAGGCCTGGGCCAAACCGATCAGCCCGGCCTTGCTCGCGGCATAGTTGGCCTGCCCGCGGTTGTCGGCAAt                                                                                                                                   <  2:212127/150‑1 (MQ=255)
         cTGCAGTCATGTGGGTTTCGATGAAACCGGGAGCGATGGCGTTGATGCTGCCACCGCGCTCGGCGAGCCGTGGCGCCCAGGCCTGGGCCAAACCGATCAGCCCGGCCTTGCTCGCGGCATAGTTGGCCTGCCCGCGGTTGTCGGCAATGc                                                                                                                                 >  1:407613/1‑150 (MQ=255)
               tCATGTGGGTTTCGATGAAACCGGGAGCGATGGCGTTGATGCTGCCACCGCGCTCGGCGAGCCGTGGCGCCCAGGCCTGGGCCAAACCGATCAGCCCGGCCTTGCTCGCGGCATAGTTGGCCTGCCCGCGGTTGTCGGCAATGCCGCTGa                                                                                                                           >  2:102243/1‑150 (MQ=255)
               tCATGTGGGTTTCGATGAAACCGGGAGCGATGGCGTTGATGCTGCCACCGCGCTCGGCGAGCCGTGGCGCCCAGGCCTGGGCCAAACCGATCAGCCCGGCCTTGCTCGAGACAGAGTTGGCCTGCCCGCGGTTGTCGGCAATGCCGCTGa                                                                                                                           >  1:209293/1‑150 (MQ=255)
                  tgtgGGTTTCGATGAAACCGGGAGCGATGGCGTTGATGCTGCCACAGCGCTCGGCGAGCCGTGGCGCCCAGGCCTGGGCCAAACCGATCAGCCCGGCCTTGCTCGCGGCATAGTTGGCCTGCCCGCGGTTGCCGGCAATGCCGCTGACCg                                                                                                                        >  2:310759/1‑150 (MQ=255)
                    tgGGTTTCGATGAAACCGGGAGCGATGGCGTTGATGCTGCCACCGCGCTCGGCGAGCCGTGGCGCCCAGGCCTGGGCCAAACCGATCAGCCCGGCCTTGCTCGCGGCATAGTTGGCCTGCCCGCGGTTGCCGGCAATGCCGCTGACCGAg                                                                                                                      <  2:565990/150‑1 (MQ=255)
                    tgGGTTTCGATGAAACCGGGAGCGATGGCATGGATGCTGCCACCGCGCGCGGCGAGCCGTGGCGCCCAGGCCTGGGCCAAACCGATCAGCCCGGCCTTGCTCGCGGCATAGTTGGCCTGCCCCCGGTTGTCTGCAATGGCGCTGACCGAg                                                                                                                      <  2:250622/150‑1 (MQ=255)
                    tgGGTTTCGACGAACCCGGGCGAGATGGCGTTGTTGCTGTCACCGCGCTCTTCGAGCCGTGGCGCCCAGGCCTGGGCCAAACCGATCAGCCCGGCCTTGCTCGCGGCATAGTTGGCCTGCCCGCGGTTGTCGGCAATGCCGCTGACCGAg                                                                                                                      <  2:209293/150‑1 (MQ=255)
                            gATGAAACCGGGAGCGATGGCGTTGATGCTGCCACCGCGCTCGGCGAGCCGTGGCGCCCAGGCCTGGGCCAAACCGATCAGCCCGGCCTTGCTCGCGGCATAGTTGGCCTGCCCGCGGTTGTCGGCAATGCCGCTGACCGAGGCcagcag                                                                                                              >  1:570181/1‑150 (MQ=255)
                            gATGAAACCGGGAGCGATGGCGTTGATGCTGCCACCGCGCTCGGCGAGCCGTGGCGCCCAGGCCTGGGCCAAACCGATCAGCCCGGCCTTGCTCGCGGCATAGTTGGCCTGCCCGCGGTTGCCGGCAATGCCGCTGACCGAGGCcagcag                                                                                                              >  1:428924/1‑150 (MQ=255)
                                      ggAGCGATGGCGTTGATGCTGCCACCGCGCTCGGCGAGCCGTGGCGCCCAGGCCTGGGCCAAACCGATCAGCCCGGCCTTGCTCGCGGCATTGTTGGCCTGCCCGCGGTTGCCGGCAATGCCGCTGACCGAGGCCAGCAGGGTGATgcgc                                                                                                    <  1:4168/150‑1 (MQ=255)
                                      ggAGCGATGGCGTTGATGCTGCCACCGCGCTCGGCGAGCCGTGGCGCCCAGGCCTGGGCCAAACCGATCAGCCCGGCCTTGCTCGCGGCATAGTTGGCCTGCCCGCGGTTGTCGGCAATGCCGCTGACCGAGGCCAGCAGGGTGATgcgc                                                                                                    <  1:618061/150‑1 (MQ=255)
                                           gATGGCGTTGATGCTGCCACCGCGCTCGGCGAGCCGTGGCGCCCAGGCCTGGGCCAAACCGATCAGCCCGGCCTTGCTCGCGGCATAGTTGGCCTGCCCGCGGTTGCCGGCAATGCCGCTGACCGAGGCCAGCAGGGTGATGCGCGCGtt                                                                                               >  1:627746/1‑150 (MQ=255)
                                            aTGGCGTTGATGCTGCCACCGCGCTCGGCGAGCCGTGGCGCCCAGGCCTGGGCCAAACCGATCAGCCCGGCCTTGCTCGCGGCATAGTTGGCCTGCCCGCGGTTGTCGGCAATGCCGCTGACCGAGGCCAGCAGGGTGATGCGCGCGttt                                                                                              <  2:796859/150‑1 (MQ=255)
                                                     aTGCTGCCACCGCGCTCGGCGAGCCGTGGCGCCCAGGCCTGGGCCAAACCGATCAGCCCGGCCTTGCTCGCGGCATAGTTGGCCTGCCCGCGGTTGCCGGCAATGCCGCTGACCGAGGCCAGCAGGGTGATGCGCGCGTTTTCCCCCAgc                                                                                     <  2:410874/150‑1 (MQ=255)
                                                            cACCGCGCTCGGCGAGCCGTGGCGCCCAGGCCTGGGCCAAACCGATCAGCCCGGCCTTGCTCGCGGCATAGTTGGCCTGCCCGCGGTTGCCGGCAATGCCGCTGACCGAGGCCAGCAGGGTGATGCGCGCGTTTTCCCCCAGCGCACCGt                                                                              >  2:649080/1‑150 (MQ=255)
                                                                  gcTCGGCGAGCCGTGGCGCCCAGGCCTGGGCCAAACCGATCAGCCCGGCCTTGCTCGCGGCATAGTTGGCCTGCCCGCGGTTGTCGGCAATGCCGCTGACCGAGGCCAGCAGGGTGATGCGCGCGTTTTCCCCCAGCGCACCGTTGTCAt                                                                        >  1:211059/1‑150 (MQ=255)
                                                                    tCGGCGAGCCGTGGCGCCCAGGCCTGGGCCAAACCGATCAGCCCGGCCTTGCTCGCGGCATAGTTGGCCTGCCCGCGGTTGCCGGCAATGCCGCTGACCGAGGCCAGCAGGGTGACGCGCGCGTTTTCCCCCAGCGCACCGTTGTCATAc                                                                      >  1:31385/1‑150 (MQ=255)
                                                                      ggCGAGCCGTGGCGCCCAGGCCTGGGCCAAACCGATCAGCCCGGCCTTGCTCGCGGCATAGTTGGCCTGCCCGCGGTTGTCGGCCATGCCACTGACCGAGGCCAGCAGGGTGATGCGCGCGTTTTCCCCCAGCGCACCGTTGTCATACaa                                                                    <  2:573497/150‑1 (MQ=255)
                                                                             cGTGGCGCCCAGGCCTGGGCCAAACCGATCAGCCCGGCCTTGCTCGCGGCATAGTTGGCCTGCCCGCGGTTGCCGGCAATGCCGCTGACCGAGGCCAGCAGGGTGACGCGCGCGTTTTCCCCCAGCGCACCGTTGTCATACAACGCCTgg                                                             <  2:31385/150‑1 (MQ=255)
                                                                                 gcgcCCAGGCCTGGGCCAAACCGATCAGCCCGGCCTTGCTCGCTGCATAGTTGGCCTGCCCGCGGTTGCCGGCAATGCCGCTGACCGAGGCCAGCAGGGTGATGCGCGCGTTTTCCCCCAGCGCACCGTGGTCATACAACGCCTGGGTGa                                                         <  2:771832/150‑1 (MQ=255)
                                                                                      cAGGCCTGGGCCAAACCGATCAGCCCGGCCTTGCTCGCGGCATAGTTGGCCTGCCCGCGGTTGTCGGCAATGCCGCTGACCGAGGCCAGCAGGGTGATGCGCGCGTTTTCCCCCAGCGCACCGTTGTCATACAACGCCTGGGTGAGCACt                                                    >  1:269997/1‑150 (MQ=255)
                                                                                         gCCTGGGCCAAACCGATCAGCCCGGCCTTGCTCGCGGCATAGTTGGCCTGCCCGCGGTTGTCGGCAATGCCGCTGACCGAGGCCAGCAGGGTGATGCGCGCGTTTTCCCCCAGCGCACCGTTGTCATACAACGCCTGGGTGAGCACTTGc                                                 <  2:211059/150‑1 (MQ=255)
                                                                                         gCCTGGGCCAAACCGATCAGCCCGGCCTTGCTCGCGGCATAGTTGGCCTGCCCGCGGTTGCCGGCAATGCCGCTGACCGAGGCCAGCAGGGTGATGCGCGCGTTTTCCCCCAGCGCACCGTTGTCATACAACGCCTGGGTGAGCACTTGc                                                 <  1:649080/150‑1 (MQ=255)
                                                                                          ccTGGGCCAAACCGATCAGCCCGGCCTTGCTCGCGGCATAGTTGGCCTGCCCGCGGTTGTCGGCAATGCCGCTGACCGAGGCCAGCAGGGTGATGCGCGCGTTTTCCCCCAGCGCACCGTTGTCATACAACGCCTGGGTGAGCACTTGCg                                                <  1:102243/150‑1 (MQ=255)
                                                                                           cTGGGCCAAACCGATCAGCCCGGCCTTGCTCGCGGCATAGTTGGCCTGCCCGCGGTTGTCGGCAATGCCGCTGACCGAGGCCAGCAGGGTGATGCGCGCGTTTTCCCCCAGCGCACCGTTGTCATACAACGCCTGGGTGAGCACTTGCgg                                               >  2:561147/1‑150 (MQ=255)
                                                                                            tGGGCCAAACCGATCAGCCCGGCCTTGCTCGCGGCATAGTTGGCCTGCCCGCGGTTGTCGGCAATGCCGCTGACCGAGGCCAGCAGGGTGCTGCGCGCGTTTTCCCCCAGCGCACCGTTGTCATACAACGCCTGGTTGAGCACTTGCGGt                                              <  2:407613/150‑1 (MQ=255)
                                                                                             gggCCAAACCGATCAGCCCGGCCTTGCTCGCGGCATAGTTGGCCTGCCCGCGGATGTTGGCAATGCCGCTGACCGAGGCCAGCAGGGTGATGCGCGCGTTTTCCCCCAGCGCACCGTTGTCATACAACGCCTGGGTGAGCACTTGCGGTg                                             >  1:496275/1‑150 (MQ=255)
                                                                                               gCCAAACCGATCAGCCCGGCCTTGCTCGCGGCATAGTTGGCCTGCCCGCGGTTGCCGGCAATGCCGCTGACCGAGGCCAGCAGGGTGATGCGCGCGTTTTCCCCCAGCGCACCGTTGTCATACAACGCCTGGGTGAGCACTTGCGGTGcc                                           <  1:515101/150‑1 (MQ=255)
                                                                                                    aCCGATCAGCCCGGCCTTGCTCGCGGCATAGTTGGCCTGCCCGCGGTTGCCGGCAATGCCGCTGACCGAGGCCAGCAGGGTGATGCGCGCGTTTTCCCCCAGCGCACCGTTGTCATACAACGCCTGGGTGAGCACTTGCGGTGCCTTGAg                                      >  1:412587/1‑150 (MQ=255)
                                                                                                    aCCGATCAGCCCGGCCTTGCTCGCGGCATAGTTGGCCTGCCCGCGGTTGCCGGCAATGCCGCTGACCGAGGCCAGCAGGGTGATGCGCGCGTTTTCCCCCAGCGCACCGTTGTCATACAACGCCTGGGTGAGCACTTGCGGTGCCTTGAg                                      <  2:412587/150‑1 (MQ=255)
                                                                                                    aCCGATCAGCCCGGCCTTGCTCGCGGCATAGTTGGACTGCCCGCGGTTGCCGGCCACGCCGCTGACCGAGGCCAGCAGGGTGATGCGCGCGTTTTACCCCAGCGCACCGTTGTCATACAACGCCTGGGTGAGCACTTGCGGTGCCTTGAg                                      <  1:407671/150‑1 (MQ=255)
                                                                                                    aCCGATCAGCCCGGCCTTGCTCGCGGCATAGTTAGCCTGCCCGCGGTTGCCGGCAATGCCGCTGACCGAGGCCAGCAGGGTGATGCGCGCGTTTTCCCCCAGCGCACCGTTCTCATACAACGCCTGGGTGAGCACTTGCGGTGCCTTGAg                                      <  2:461046/150‑1 (MQ=255)
                                                                                                            gCCCGGCCTTGCTCGCGGCATAGTTGGCCTGCCCGCGGTTGTCGGCAATGCCGCTGACCGAGGCCAGCAGGGTGATGCGCGCGTTTTCCCCCAGCGCACCGTTGTCATACAACGCCTGGGTGAGCACTTGCGGTGCCTTGAGGTTCACgg                              >  1:369181/1‑150 (MQ=255)
                                                                                                                       cTCGCGGCATAGTTGGCCTGCCCGCGGTTGTCGGCAATGCCGCTGACCGAGGCCAGCAGGGTGATGCGCGCGTTTTCCCCCAGCGCACCGTTGTCATACAACGCCTGGGTGAGCACTTGCGGTGCCTTGAg                                      <  1:139563/131‑1 (MQ=255)
                                                                                                                       cTCGCGGCATAGTTGGCCTGCCCGCGGTTGTCGGCAATGCCGCTGACCGAGGCCAGCAGGGTGATGCGCGCGTTTTCCCCCAGCGCACCGTTGTCATACAACGCCTGGGTGAGCACTTGCGGTGCCTTGAg                                      >  2:139563/1‑131 (MQ=255)
                                                                                                                          gcgGCATAGTTGGCCTGCCCGCGGTTGTCGGCAATGCCGCTGACCGAGGCCAGCAGGGTGATGCGCGCGTTTTCCCCCAGCGCACCGTTGTCATACAACGCCTGGGTCCTGACTTGCGGTGCCTTGAGGTTCACGGCCAGCACGGCGTcc                >  1:406769/1‑150 (MQ=255)
                                                                                                                           cgGCATAGTTGGCCTGCCCGCGGTTGCCGGCAATGCCGCTGACAGAGGGCGGCAGTGTGATGCGCGCGTTTTCCCCCAGCGCACCGTTGTCATACAACGCCTGGGTGAGCACTTGCGGTGCCTTGAGGTTCACGGCCAGCACGGCGTccc               <  2:222816/150‑1 (MQ=255)
                                                                                                                            ggCATAGTTGGCCTGCCCGCGGTTGCCGGCAATGCCGCTGACCGAGGCCAGCAGGGTGATGCGCGCGTTTTCCCCCAGCGCACCGTTGTCATACAACGCCTGGGTGAGCACTTGCGGTGCCTTGAGGTTCACg                               <  1:717515/133‑1 (MQ=255)
                                                                                                                            ggCATAGTTGGCCTGCCCGCGGTTGCCGGCAATGCCGCTGACCGAGGCCAGCAGGGTGATGCGCGCGTTTTCCCCCAGCGCACCGTTGTCATACAACGCCTGGGTGAGCACTTGCGGTGCCTTGAGGTTCACg                               >  2:717515/1‑133 (MQ=255)
                                                                                                                             gCATAGTTGGCCTGCCCGCGGTTGTCGGCAATGCCGCTGACCGAGGCCAGCAGGGTGATGCGCGCGTTTTCCCCCAGCGCACCGTTGTCATACAACGCCTGGGTGAGCACTTGCGGTGCCTTGAGGTTCACGGCCAGCACGGCGTCCCAg             <  2:570181/150‑1 (MQ=255)
                                                                                                                                      gCCTGCCCGCGGTTGCCGGCAATGCCGCTGACCGAGGCCAGCAGGGTGATGCGCGCGTTTTCCCCCAGCGCACCGTTGTCATACAACGCCTGGGTGAGCACTTGCGGTGCCTTGAGGTTCACGGCCAGCACGGCGTCCCAGTATTCCGGc    <  1:310759/150‑1 (MQ=255)
                                                                                                                                        cTGCCCGCGGTTGTCGGCAATGCCGCTGACCGAGGCCAGCAGGGTGATGCGCGCGTTTTCCCCCAGCGCACCGTTGTCATACAACGCCTGGGTGAGCACTTGCGGTGCCTTGAGGTTCACGGCCAGCACGGCGTCCCAGAATTCCGGCGt  >  1:22737/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                        
TGGGCATGGCTGCAGTCATGTGGGTTTCGATGAAACCGGGAGCGATGGCGTTGATGCTGCCACCGCGCTCGGCGAGCCGTGGCGCCCAGGCCTGGGCCAAACCGATCAGCCCGGCCTTGCTCGCGGCATAGTTGGCCTGCCCGCGGTTGCCGGCAATGCCGCTGACCGAGGCCAGCAGGGTGATGCGCGCGTTTTCCCCCAGCGCACCGTTGTCATACAACGCCTGGGTGAGCACTTGCGGTGCCTTGAGGTTCACGGCCAGCACGGCGTCCCAGTATTCCGGCGT  >  NC_002947/677805‑678090

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