Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,166,255 C→A 16.6% A215E (GCG→GAG)  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,2550CA16.6% 80.8 / 6.4 36A215E (GCG→GAG) PP_27803‑oxoacyl‑ACP synthase
Reads supporting (aligned to +/- strand):  ref base C (12/18);  new base A (3/3);  total (15/21)
Fisher's exact test for biased strand distribution p-value = 6.77e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.57e-01

GTGAAGGCCCGGTGATCAATGTCAGCAGCGCCTGCGCGGGAGCGTCCCAGGCCATCGGCAATGCATTCCAGATGATCCGTCGCGGCGAGGTCGAGCTGGCCATTGCCGGGGGCGCCGATTCGGTGCTTAGCCTGGACACCATGACCGCGCTGTATCTGCTCGGCGCCGCCAGCAGCGAACAGCGTTGGGGCACCGAACTGTGTCGGCCGTTCGATCGCCATCGCAGTGGCTTGATTGCGGGTGAGGGCGGTGGTTTTGTCGTGCTCGAAACCCTCGAACGGGCGCTCGCCCGAGGGG  >  NC_002947/3166108‑3166404
                                                                                                                                                   |                                                                                                                                                     
gTGAAGGCCCGGTGATCAATGTCAGCAGCGCCTGCGTGGGAGCGTGCCAGGCCATCGGCAATGCATTGCAGATGGCCCGTCGCGGCGACGTCGAGCACGCCATTGCCGGGGGCGCCGATTCCGTGCTTAGCCTCGACACCATGGCCgcgc                                                                                                                                                     <  2:42604/150‑1 (MQ=255)
gTGAAGGCCCGGTGATCAATGTCAGCAGCGCCTGCGCGGGAGCGTCCCAGGCCATCCGCAATGCATTCCAGATGATCCGTCGCGGCGAGGTCGAGCTGGCCATTGCCGGGGGCGCCGATTCGGTGCTTAGCCTGGACACCATGACCgcgc                                                                                                                                                     <  2:755036/150‑1 (MQ=255)
 tGAAGGCCCGGTGATCAATGTCAGCAGCGCCTGCGCGGGAGCGTCCCAGGCCATCGGCAATGCATTCCAGATGATCCGTCGCGGCGAGGTCGAGCTGGCCATTGCCGGGGGCGCCGATTCGGTGCTTAGCCTGGACACCATGACCGCGCt                                                                                                                                                    >  2:662363/1‑150 (MQ=255)
                        gcagcGCCTGCGCGGGAGCGTCCCAGGCCATCGGCAATGCATTCCAGATGATCCGTCGCGGCGAGGTCGAGCTGGCCATTGCCGGGGGCGCCGATTCGGTGCTTAGCCTGGACACCATGACCGCGCTGTATCTGCTCGGCGCCGCcagca                                                                                                                             <  2:685863/150‑1 (MQ=255)
                                   gcgGGAGCGTCCCAGGCCATCGGCAATGCATTCCAGATGATCCGTCGCGGCGAGGTCGAGCTGGCCATTGCCGGGGGCGCCGATTCGGTGCTTAGCCTGGACACCATGACCGCGCTGTATCTGCTCGGCGCCGCCAGCAGCGAACAGCGt                                                                                                                  <  2:714805/150‑1 (MQ=255)
                                     gggAGCGTCCCAGGCCATCGGCAATGCATTCCAGATGATCCGTCGCGGCGAGGTCGAGCTGGCCATTGCCGGGGGCGCCGATTCGGTGCTTAGCCTGGACACCATGACCGCGCTGTATCTGCTCGGCGCCGCCAGCAGCGAACAGCGTTg                                                                                                                <  1:548454/150‑1 (MQ=255)
                                      ggAGCGTCCCAGGCCATCGGCAATGCATTCCAGATGATCCGTCGCGGCGAGGTCGAGCTGGCCATTGCCGGGGGCGCCGATTCGGTGCTTAGCCTGGACACCATGACCGCGCTGTATCTGCTCGGCGCCGCCAGCAGCGAACAGCGTTgg                                                                                                               <  1:694257/150‑1 (MQ=255)
                                       gAGCGTCCCAGGCCATCGGCAATGCATTTCAGATGATCCGTCGCGGCGAGGTCGAGCTGGCCATTGCCGGGGGCGCCGATTCGGTCCTTAGCCTGGACACCATGACCGCGCTGTATCTGCTCGGCGCCGCCAGCAGCGAACAGCGTTggg                                                                                                              <  2:269997/150‑1 (MQ=255)
                                         gCGTCCCAGGCCATCGGCAATGCATTCCAGATGATCCGTCGCGGCGAGGTCGAGCTGGCCATTGCCGGGGACGCGGATTCGGTGCTCAGGCTGGACACCATGACCGCGCTGTATATGCTCGGCGCCGCCAGCAGCGAACAGCGTTGGGGc                                                                                                            >  2:339347/1‑150 (MQ=255)
                                               cAGGCCATCGGCAATGCATTCCAGATGATCCGTCGCGGCGAGGTCGAGCTGGCCATTGCCGGGGGCGCCGATTCGGTGCTTAGCCTGGACACCATGACCGAGCTGTATCTGCTCGGCGCCGCCAGCAGCGAACAGCGTTgg                                                                                                               >  2:237617/1‑141 (MQ=255)
                                               cAGGCCATCGGCAATGCATTCCAGATGATCCGTCGCGGCGAGGTCGAGCTGGCCATTGCCGGGGGCGCCGATTCGGTGCTTAGCCTGGACACCATGACCGAGCTGTATCTGCTCGGCGCCGCCAGCAGCGAACAGCGTTgg                                                                                                               >  2:237840/1‑141 (MQ=255)
                                               cAGGCCATCGGCAATGCATTCCAGATGATCCGTCGCGGCGAGGTCGAGCTGGCCATTGCCGGGGGCGCCGATTCGGTGCTTAGCCTGGACACCATGACCGAGCTGTATCTGCTCGGCGCCGCCAGCAGCGAACAGCGTTgg                                                                                                               >  2:237815/1‑141 (MQ=255)
                                               cAGGCCATCGGCAATGCATTCCAGATGATCCGTCGCGGCGAGGTCGAGCTGGCCATTGCCGGGGGCGCCGATTCGGTGCTTAGCCTGGACACCATGACCGAGCTGTATCTGCTCGGCGCCGCCAGCAGCGAACAGCGTTgg                                                                                                               <  1:237840/141‑1 (MQ=255)
                                               cAGGCCATCGGCAATGCATTCCAGATGATCCGTCGCGGCGAGGTCGAGCTGGCCATTGCCGGGGGCGCCGATTCGGTGCTTAGCCTGGACACCATGACCGAGCTGTATCTGCTCGGCGCCGCCAGCAGCGAACAGCGTTgg                                                                                                               <  1:237815/141‑1 (MQ=255)
                                               cAGGCCATCGGCAATGCATTCCAGATGATCCGTCGCGGCGAGGTCGAGCTGGCCATTGCCGGGGGCGCCGATTCGGTGCTTAGCCTGGACACCATGACCGAGCTGTATCTGCTCGGCGCCGCCAGCAGCGAACAGCGTTgg                                                                                                               <  1:237617/141‑1 (MQ=255)
                                                   ccATCGGCAATGCATTCCAGATGATCCGTCGCGGCGAGGTCGAGCTGGCCATTGCCGGGGGCGCCGATTCGGTGCTTAGCCTGGACACCATGACCGCGCTGTATCTGCTCGGCGCCGCCAGCAGCGAACAGCGTTGGGGCACCGAACtgt                                                                                                  <  1:308295/150‑1 (MQ=255)
                                                         gCAATGCATTCCAGATGATCCGTCGCGGCGAGGTCGAGCTGGCCATTGCCGGGGGCGCCGATTCGGTGCTTAGCCTGGACACCATGACCGCGCTGTATCTGCGCGGCGCCGCCAGCAGCGAACAGCGTTGGGGCACCGAACTGTGTCGGc                                                                                            <  1:662363/150‑1 (MQ=255)
                                                                aTTCCAGATGATCCGTCGCGGCGAGGTCGAGCTGGCCATTGCCGGGGGCGCCGATTCGGTGCTTAGCCTGGACACCATGACCGCGCTGTATCTGCTCGGCGCCGCCAGCAGCGAACAGCGtt                                                                                                                 >  2:371377/1‑122 (MQ=255)
                                                                aTTCCAGATGATCCGTCGCGGCGAGGTCGAGCTGGCCATTGCCGGGGGCGCCGATTCGGTGCTTAGCCTGGACACCATGACCGCGCTGTATCTGCTCGGCGCCGCCAGCAGCGAACAGCGtt                                                                                                                 <  1:371377/122‑1 (MQ=255)
                                                                  tCCAGATGATCCGTCGCGGCGAGGTCGAGCTGGCCATTGCCGGGGGCGCCGATTCGGTGCTTAGCCTGGACACCATGACCGCGCTGTATCTGCTCGGCGCCGCCAGCAGCGAACAGCGTTGGGGTACCGAACTGTGTCGGCGGTTCGATc                                                                                   <  2:291486/150‑1 (MQ=255)
                                                                     aGATGATCCGTCGCGGCGAGGTCGAGCTGGCCATTGCCGGGGGCGCCGATTCGGTGCTTAGCCTGGACACCATGACCGCGCTGTATCTGCTCGGCGCCGCCAGCAGCGAACAGCGTTGGGGCACCGAACTGTGTCGGCCGTTCGATCGcc                                                                                >  2:752404/1‑150 (MQ=255)
                                                                            ccGTCGCGGCGAGGTCGAGCTGGCCATTGCCGGGGGCGCCGATTCGGTGCTTAGCCTGGACACCATGACCGCGCTGTATCTGCTCGgc                                                                                                                                       >  1:478878/1‑88 (MQ=255)
                                                                            ccGTCGCGGCGAGGTCGAGCTGGCCATTGCCGGGGGCGCCGATTCGGTGCTTAGCCTGGACACCATGACCGCGCTGTATCTGCTCGgc                                                                                                                                       <  2:478878/88‑1 (MQ=255)
                                                                                cgcgGCGAGGTCGAGCTGGCCATTGCCGGGGGCGCCGATTCGGTGCTTAGCCTGGACACCATGACCGCGCTGTATCTGCTCGGCGCCGCCAGCAGCGAACAGCGTTGGGGCACCGAACTGTGTCGGCCGTTCGATCGCCATCGCAGTGGc                                                                     <  2:383187/150‑1 (MQ=255)
                                                                                     cgAGGGCGAGCTGGCCATTGCCGGGGGCGCCGATTCGGTGCTTAGCCTGGACACCATGACCGCGCTGTATCTGCTCGGCGCCGCCAGCAGCGAACAGCGTTGGGGCACCGAACTGTGTCGGCCGTTCGATCGCCATCGCAGTGGCTTGAt                                                                >  2:688031/1‑150 (MQ=255)
                                                                                      gAGGTCGAGCTGGCCATTGCCGGGGGCGCCGATTCGGTGCTTAGCCTGGACACCATGACCGCGCTGTATCTGCTCGGCGCCGCCAGCAGCGAACAGCGTTGGGGCACCGAACTGTGTCGGCCGTTCGATCGCCATCGCGGTGGCTTGAtt                                                               >  1:389882/1‑150 (MQ=255)
                                                                                          tCGAGCTGGCCATTGCCGGGGGCGCCGATTCGGTGCTTAGCCTGGACACCATGACCGCGCTGTATCTGCTCGGCGCCGCCAGCAGCGAACAGCGTTGGGGCACCGAACTGTGTCGGCCGTTCGATCGCCATCGCAGTGGCTTGATTGCgg                                                           <  1:194280/150‑1 (MQ=255)
                                                                                             aGCTGGCCATTGCCGGGGGCGCCGATTCGGTGCTTAGCCTGGACACCATGACCGCGCt                                                                                                                                                    >  1:154468/1‑58 (MQ=255)
                                                                                             aGCTGGCCATTGCCGGGGGCGCCGATTCGGTGCTTAGCCTGGACACCATGACCGCGCt                                                                                                                                                    <  2:154468/58‑1 (MQ=255)
                                                                                                         ccGGGGGCGCCGATTCGGTGCTTAGCCTGGACACCATGACCGCGCTGTATCTGCTCGGCGCCGCCAGCAGCGAACAGCGTTGGGGCACCGAACTGTGTCGGCCGTTCGATCGCCATCGCAGTGGCTTGATTGCGGGTGAGGGCGGTGGtt                                            >  1:186184/1‑150 (MQ=255)
                                                                                                                 gcCGATTCGGTGCTTAGCCTGGACACCATGACCGCGCTGTATCTGCTCGGCGCCGCCAGCAGCGAACAGCGTTGGGGCACCGAACTGTGTCGGCCGTTCGATCGCCATCGCAGTGGCTTGATTGCGGGTGAGCGCCGTGGTTTTGTCGTg                                    >  2:687700/1‑150 (MQ=255)
                                                                                                                             cTTAGCCTGGACACCATGACCGCGCTGTATCTGCTCGGCGCCGCCAGCAGCGAACAGCGTTGGGGCACCGAACTGTGTCGGCCGTTCGATCGCCATCGCAGTGGCTTGATTGCGGGTGAGGGCGGTGGTTTTGTCGTGCTCGAAACCCTc                        <  2:402908/150‑1 (MQ=255)
                                                                                                                                     ggACACCATGACCGCGCTGTATCTGCTCGGCGCCGCCAGCAGCGAACAGCGTTGGGGCACCGAACTGTGTCGGCCGTTCGATCGCCATCGCAGTGGCTTGATTGCGGGTGAGGGCGGTGGTTTTGTCGTGCTCGAAACCCTCGAACGGgc                <  2:483594/150‑1 (MQ=255)
                                                                                                                                     ggACACCATGACCGCGCTGTATCTGCTCGGCGCCGCCAGCAGCGAACAGCGTTGGGGCACCGAACTGTGTCGGCCGTTCGATCGCCATCGCAGTGGCTTGATTGCGGGTGAGGGCGGTGGTTTTGTCGTGCTCGAAACCCTCGAACGGgc                <  1:299621/150‑1 (MQ=255)
                                                                                                                                     ggACACCATGACCGCGCTGTATCTGCTCGGCGCCGCCAGCAGCGAACAGCGTTGGGGCACCGAACTGTGTCGGCCGTTCGATCGCCATCGCAGTGGCTTGATTGCGGGTGAGGGCGGTGGTTTTGTCGTGCTCGAAACCCTCGAACGGgc                <  1:688031/150‑1 (MQ=255)
                                                                                                                                         acCATGACCGCGCTGTATCTGCTCGGCGCCGCCAGCAGCCAACCGCGTTGGCGTACCGAACTGTGTCGGCCGGTCGATCGCCATCGCAGTGGCTTGATTGCGGGTGAGGGCGGTGGTTTTATTGTACTCGAAACCCTCGAACGGGCGCTc            >  1:787443/1‑150 (MQ=255)
                                                                                                                                           cATGACCGCGCTGTATCTGCTCGGCGCCGCCAGCAGCGAACAGCGTTGGGGCACCGAACTGTGTCGGCCGTTCGATCGCCATCGCAGTGGCtt                                                                   >  2:365777/1‑93 (MQ=255)
                                                                                                                                           cATGACCGCGCTGTATCTGCTCGGCGCCGCCAGCAGCGAACAGCGTTGGGGCACCGAACTGTGTCGGCCGTTCGATCGCCATCGCAGTGGCtt                                                                   <  1:365777/93‑1 (MQ=255)
                                                                                                                                                 tgcgcTGTATCTGCTCGGCGCCGCCAGCAGCGAACAGCGTTGGGGCACCGAACTGTGTCGGCCGTTCGATCGCCATCGCAGTGGCTTGATTGCGGGTGAGGGCGGTGGTTTTGTCGTGCTCGAAACCCTCGAACGGGCGCTCGCCCGAgg    >  2:356333/2‑150 (MQ=255)
                                                                                                                                                 cgcgCTGTATCTGCTCGGCGCCGCCAGCAGCGAACAGCGTTGGGGCACCGAACTGTGTCGGCCGTTCGATCGCCATCGCAGTGGCTTGATTGCGGGTGAGGGCGGTGGTTTTGTCGTGCTCGAAACCCTCGAACGGGCGCTCGCCCGAgg    <  2:186184/150‑1 (MQ=255)
                                                                                                                                                   cgcTGTATCTGCTCGGCGCCGCCAGCAGCGAACAGCGTTGGGGCACCGAACTGTGTCGGCCGTTCGATCGCCATCGCAGTGGCTTGATTTCGGGTGAGGGCGGTGGTTTTGTCGTGCTCGAAACCCTCGAACGGGCGCTCGCCCGAgggg  >  1:538619/1‑150 (MQ=255)
                                                                                                                                                   |                                                                                                                                                     
GTGAAGGCCCGGTGATCAATGTCAGCAGCGCCTGCGCGGGAGCGTCCCAGGCCATCGGCAATGCATTCCAGATGATCCGTCGCGGCGAGGTCGAGCTGGCCATTGCCGGGGGCGCCGATTCGGTGCTTAGCCTGGACACCATGACCGCGCTGTATCTGCTCGGCGCCGCCAGCAGCGAACAGCGTTGGGGCACCGAACTGTGTCGGCCGTTCGATCGCCATCGCAGTGGCTTGATTGCGGGTGAGGGCGGTGGTTTTGTCGTGCTCGAAACCCTCGAACGGGCGCTCGCCCGAGGGG  >  NC_002947/3166108‑3166404

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: ‑