Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 6,079,777 C→A 12.1% G45G (GGC→GGA PP_5331 → acyl‑CoA thioesterase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029476,079,7770CA12.1% 88.3 / 3.0 33G45G (GGC→GGAPP_5331acyl‑CoA thioesterase
Reads supporting (aligned to +/- strand):  ref base C (14/15);  new base A (2/2);  total (16/17)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.86e-01

GAGGCCCCATCGCATGATTGAGCTCGAACAAGAAGATCCAATCCCGCAAGGCGACCTGGCCTTGCAAATCACCGCATTGCCGCGAGAAACCAACGGGTTTGGCGATATTTTCGGCGGTTGGCTGGTTGCGCAGATGGACCTGGCAGGCACGGCCATGGCCAGCCGTGTCGCCGGCGGCCGCGTAGCCACCGTCGCCATCGACCGCATGGCCTTCCTGGTACCCGTCGCGGTCGGCGCACAGCTGTCTTTCTATACCCAGACCCTGGAAATCGGCCGCAGCTCGATCCAGATGATGG  >  NC_002947/6079630‑6079925
                                                                                                                                                   |                                                                                                                                                    
gagGCCCCATCGCATGATTGAGCTCGAACAAGAAGATCCAATCCCGCAAGGCGACCTGGCCTTGCAAATCACCGCATTGCCGCGAGAAACCAACGGGTTTGGCGATATTTTCGGCGGTTGGCTGGTTGCGCAGATGGACCTGGCAGGCAc                                                                                                                                                    >  2:252431/1‑150 (MQ=255)
     cccATCGCATGATTGAGCTCGAACAAGAAGATCCAATCCCGCAAGGCGACCTGGCCTTGCAAATCACCGCATTGCCGCGAGAAACCAACGGGTTTGGCGATATTTTCGGCGGTTGGCTGGTTTCGCAGATGGACCTTGCAGGCACGGccg                                                                                                                                               >  1:257462/1‑149 (MQ=255)
      ccATCGCATGATTGAGCTCGAACAAGAAGATCCAATCCCGCAAGGCGACCTGGCCTTGCAAATCACCGCATTGCCGCGAGAAACCAACGGGTTTGGCGATATTTTCGGCGGTTGGCTGGTTGCGCAGATGGACCTGGCAGGCACGGCCAt                                                                                                                                              <  1:271280/150‑1 (MQ=255)
       cATCGCATGATTGAGCTCGAACAAGAAGATCCAATCCCGCAAGGCGACCTGGCCTTGCAAATCACCGCATTGCCGCGAGAAACCAACGGGTTTGGCGATATTTTCGGCGGTTGGCTGGTTGCGCAGATGGACCTGGCAGGCACGGCCATg                                                                                                                                             <  2:80799/150‑1 (MQ=255)
       cATCGCATGATTGAGCTCGAACAAGAAGATCCAATCCCGCAAGGCGACCTGGCCTTGCAAATCACCGCATTGCCGCGAGAAACCAACGGGTTTGGCGATATTTTCGGCGGTTGGCTGGTTGCGCAGATGGACCTGGCAGGCACGGCCATg                                                                                                                                             <  1:539808/150‑1 (MQ=255)
        aTCGCATGATTGAGCTCGAACAAGAAGATCCAATCCCGCAAGGCGACCTGGCCTTGCAAATCACCGCATTGCCGCGAGAAACCAACGGGTTTGGCGATATTTTCGGCGGTTGGCTGGTTGCGCAGATGGACCTGGCAGGCACGGCCATgg                                                                                                                                            <  2:418208/150‑1 (MQ=255)
                     gCTCGAACAAGAAGATCCAATCCCGCAAGGCGACCTGGCCTTGCAAATCACCGCATTGCCGCGAGAAACCAACGGGTTTGGCGATATTTTCGGCGGTTGGCTGGTTGCGCAGATGGACCTGGCAGGCACGGCCATGGCCAGCCGTGTCGc                                                                                                                               <  1:80696/150‑1 (MQ=255)
                         gAACAAGAAGATCCAATCACGCAAGTCGACCTGGCCTTGCAAATCACCGCATTGCCGCGAGAAACCAACGGGTTTGGCGATATTTTCGGCGGTTGGCTGGTTGCGCAGATGGACCTGGCAGGCACGGCCATGGCCAGCCGTGTCGCcggc                                                                                                                           >  2:333057/1‑150 (MQ=255)
                                              cAAGGCGACCTGGCCTTGCAAATCACCGCATTGCCGCGAGAAACCAACGGGTTTGGCGATCTTTTCGGCGGTTTGCGGTTTGCGCAGATGGACCTGGCAGGCACGGCCATGGCCAGCCGTGTCGCCGGCGGCCGCGTAGCCACCGTCGcc                                                                                                      <  1:351254/150‑1 (MQ=255)
                                                    gACCTGGCCTTGCAAATCACCGCATTGCCGCGAGAAACCAACGGGTTTGGCGATATTTTCGGCGGTTGGCTGGTTGCGCAGATGGACCTGGCAGGCACGGCCATGGCCAGCCGTGTCGCGTGCTTCCGCGTAGCCACCGTCGCCATCGAc                                                                                                <  1:298625/150‑1 (MQ=255)
                                                    gACCTGGCCTTGCAAATCACCGCATTGCCGCGAGAAACCAACGGGTTTGGCGATATTTTCGGCGGTTGGCTGGTTGCGCAGATGGACCTGGCAGGCACGGCCATGGCCAGCCGTGTCGCCGGCGGCCGCGTAGCCACCGTCGCCATCGAc                                                                                                <  1:413195/150‑1 (MQ=255)
                                                     aCCTGGCCTTGCAAATCACCGCATTGCCGCGAGAAACCAACGGGTTTGGCGATATTTTCGGCGGTTGGCTGGTTGCGCAGATGGACCTGGCAGGCACGGCCATGGCCAGCCGTGTCGCCGGCGGCCGCGTAGCCACCGTCGCCATCGAcc                                                                                               <  2:82435/150‑1 (MQ=255)
                                                     aCCTGGCCTTGCAAATCACCGCATTGCCGCGAGAAACCAACGGGTTTGGCGATATTTTCGGCGGTTGGCTGGTTGCGCAGATGGACCTGGCAGGCACGGCCATGGCCAGCCGTGTCGCCGGCGGCCGCGTAGCCACCGTCGCCATCGAcc                                                                                               >  1:464993/1‑150 (MQ=255)
                                                       cTGGCCTTGCAAATCACCGCATTGCCGCGAGAAACCAACGGGTTTGGCGATATTTTCGGCGGTTGGCTGGTTGCGCAGATGGACCTGGCAGGCACGGCCATGGCCAGCCGTGTCGCCGGCGGCCGCGTAGCCACCGTCGCCATCGACCGc                                                                                             >  1:557248/1‑150 (MQ=255)
                                                                   aTCACCGCATTGCCGCGAGAAACCAACGGGTTTGGCGATATTTTCGGCGGTTGGCTGGTTGCGCAGATGGACCTGGCAGGCACGGCCATGGCCAGCCGTGTCGCCGGCGGCCGCGTAGCCACCGTCGCCATCGACCGCATg                                                                                          >  1:242867/1‑141 (MQ=255)
                                                                   aTCACCGCATTGCCGCGAGAAACCAACGGGTTTGGCGATATTTTCGGCGGTTGGCTGGTTGCGCAGATGGACCTGGCAGGCACGGCCATGGCCAGCCGTGTCGCCGGCGGCCGCGTAGCCACCGTCGCCATCGACCGCATg                                                                                          <  2:242867/141‑1 (MQ=255)
                                                                    ccACCGCCTTGCTGCGAGAATCCAACGGGTTTGGCGCGCTTTTCGGCTGTTGGCTGGTTGCGCAGATGGTCCTGGCAGGCTCGGCCATGGCCAGCCGTGTCGCCGGCGGCCGCGTAGCCACCGTCGCCATCGACCGCATGGCCTTCCTgg                                                                                <  2:257462/149‑1 (MQ=255)
                                                                      aCCGCATTGCCGCGAGAAACCAACGGGTTTGGCGATATTTTCGGCGGTTGGCTGGTTGCGCAGATGGACCTGGCAGGCACGGCCATGGCCAGCCGTGTCGCCGGCGGCCGCGTAGCCACCGTCGCCATCGACCGCATGGCCTTCCTGGTa                                                                              >  1:471945/1‑150 (MQ=255)
                                                                          cATTGCCGCGAGAAACCAACGGGTTTGGCGATATTTTCGGCGGTTGGCTGGTTGCGCAGATGGACCTGGCAGGCACGGCCATGGCCAGCCGTGTCGCCGGCGGCCGCGTAGCCACCGTCGCCATCGACCGCATGGCCTTCCTGGTACCCg                                                                          >  2:153589/1‑150 (MQ=255)
                                                                                   gagaAACCAACGGGTTTGGCGATATTTTCGGCGGTTGGCTGGTTGCGCAGATGGACGTGGCAGGCACGGCCATGGCCAGCCGTGTCGCCGGCGGCCGCGTAGCCACCGTCGCCATCGACCGCATGGCCTTCCTGGTACCCGTCGCGGTCg                                                                 >  2:365157/1‑150 (MQ=255)
                                                                                           aaCGGGTTTGGCGATATTTTCGGCGGTTGGCTGGTTGCGCAGATGGACCTGGCAGGCACGGCCATGGCCAGCCGTGTCGCCGGCGGCCGCGTAGCCACCGTCGCCATCGACCGCATGGCCTTCCTGGTACCCGTCGCGGTCGGCGCACAg                                                         <  2:557248/150‑1 (MQ=255)
                                                                                                gTTTGGCGATATTTTCGGCGGTTGGCTGGTTGCGCAGATGGACCTGGCAGGAACGGCCATGGCCa                                                                                                                                         >  1:392558/1‑65 (MQ=255)
                                                                                                gTTTGGCGATATTTTCGGCGGTTGGCTGGTTGCGCAGATGGACCTGGCAGGAACGGCCATGGCCa                                                                                                                                         >  1:392415/1‑65 (MQ=255)
                                                                                                gTTTGGCGATATTTTCGGCGGTTGGCTGGTTGCGCAGATGGACCTGGCAGGAACGGCCATGGCCa                                                                                                                                         <  2:392415/65‑1 (MQ=255)
                                                                                                gTTTGGCGATATTTTCGGCGGTTGGCTGGTTGCGCAGATGGACCTGGCAGGAACGGCCATGGCCa                                                                                                                                         <  2:392558/65‑1 (MQ=255)
                                                                                                 tttGGCGATATTTTCGGCGGTTGGCTGGTTGCGCAGATGGACCTGGCAGGCACGGCCATGGCCAGCCGTGTCGCCGGCGGCCGCGTAGCCACCGTCGCCATCGACCGCATGGCCTTCCTGGTACCCGTCGCGGTCGGCGCACAGCTGTCt                                                   <  1:375246/150‑1 (MQ=255)
                                                                                                     gCGATATTTTCGGCGGTTGGCTGGTTGCGCAGATGGACCTGGCAGGCACGGCCATGGCCAGCCGTGTCGCCGGCGGCCGCGTAGCCACCGTCGCCATCGACCGCATGGCCTTCCTGGTACCCGTCGCGGTCGGCGCACAGCTGTCTTTCt                                               >  1:282723/1‑150 (MQ=255)
                                                                                                      cGATATTTTCGGCGGTTGGCTGGTTGCGCAGATGGACCTGGCAGGCACGGCCATGGCCAGCCGTGTCGCCGGCGGCCGCGTAGCCACCGTCGCCATCGACCGCATGGCCTTCCTGGTACCCGTCGCGGTCGGCGCACAGCTGTCTTTCta                                              <  2:464993/150‑1 (MQ=255)
                                                                                                            tttCGGCGGTTGGCTGGTTGCGCAGATGGACCTGGCAGGCACGGCCATGGCCAGCCGTGTCGCCGGCGGCCGCGTAGCCACCGTCGCCATCGACCGCATGGCCTTCCTGGTACCCGTCGCGGTCGGCGCACAGCTGTCTTTCTATACCCa                                        >  2:183611/1‑150 (MQ=255)
                                                                                                                                     aTGGACCTGGCAGGCACGGCCATGGCCAGCCGTGTCGCCGGCGGCCGCGTAGCCACCGTCGCCATCGACCGCATGGCCTTCCTGGTACCCGTCGCGGTCGGCGCACAGCTGTCTTTCTATACCCAGACCCTGGAAATCGGCCGCAGCTCg               >  1:498054/1‑150 (MQ=255)
                                                                                                                                       ggACCTGGCAGTCACGGCCATGGCCAGCAGTGTCGCCGGCGGCGGCGTAGCCACCGTCGCCATCGACCGCATGCTCTTTCTGGTACCCGTCGCGGTCGGCGCACAGCTGTCTTTCTATACCCAGACCCTGGAAATCGGCCGCAGCTCGAt             <  1:549135/150‑1 (MQ=255)
                                                                                                                                         aCCTGGCAGGCACGGCCATGGCCAGCCGTGTCGCCGGCGGCCGCGTAGCCACCGTCGCCATCGACCGCATGGCCTTCCTGGTACCCGTCGCGGTCGGCGCACAGCTGTCTTTCTATACCCAGACCCTGGAAATCGGCCGCAGCTCGATcc           >  1:324117/1‑150 (MQ=255)
                                                                                                                                           cTGGCAGGCACGGCCATGGCCAGCCGTGTCGCCGGCGGCCGCGTAGCCACCGTCGCCATCGACCGCATGGCCTTCCTGGTACCCGTCGCGGTCGGCGCACAGCTGTCTTTCTATACCCAGACCCTGGAAATCGGCCGCAGCTCGATCCAg         >  1:281458/1‑150 (MQ=255)
                                                                                                                                                aggcaCGGCCATGGCCAGCCGTGTCGCCGGCGGCCGCGTAGCCACCGTCGCCATCGACCGCATGGCCTTCCTGGTACCCGTCGCGGTCGGCGCACAGCTGTCTTTCTATACCCAGACCCTGGAAATCGGCCGCAGCTCGATCCAgatgat    >  1:528506/1‑150 (MQ=255)
                                                                                                                                                aggcaCGGCCATGGCCAGCCGTGTCGCCGGCGGCCGCGTAGCCACCGTCGCCATCGACCGCATGGCCTTCCTGGTACCCGTCGCGGTCGGCGCACAGCTGTCTTTCTATACCCAGACCCTGGAAATCGGCCGCAGCTCGATCCAgatgat    >  1:238346/1‑150 (MQ=255)
                                                                                                                                                  gcaCGGCCATGGCCAGCCGTGTCGCCGGCGGCCGCGTAGCCACCGTCGCCATCGACCGCATGGCCTTCCTGGTACCCGTCGCGGTCGGCGCACAGCTGTCTTTCTATACCCAGACCCTGGAAATCGGCCGCAGCTCGATCCAGATGATgg  >  2:506840/1‑150 (MQ=255)
                                                                                                                                                   |                                                                                                                                                    
GAGGCCCCATCGCATGATTGAGCTCGAACAAGAAGATCCAATCCCGCAAGGCGACCTGGCCTTGCAAATCACCGCATTGCCGCGAGAAACCAACGGGTTTGGCGATATTTTCGGCGGTTGGCTGGTTGCGCAGATGGACCTGGCAGGCACGGCCATGGCCAGCCGTGTCGCCGGCGGCCGCGTAGCCACCGTCGCCATCGACCGCATGGCCTTCCTGGTACCCGTCGCGGTCGGCGCACAGCTGTCTTTCTATACCCAGACCCTGGAAATCGGCCGCAGCTCGATCCAGATGATGG  >  NC_002947/6079630‑6079925

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