Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 680,295 C→G 9.2% N394K (AAC→AAG) ‡ PP_0582 → thiolase family protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947680,2950CG 9.2% 151.5 / 5.8 54N394K (AAC→AAG) ‡PP_0582thiolase family protein
Reads supporting (aligned to +/- strand):  ref base C (21/28);  new base G (3/2);  total (24/30)
Fisher's exact test for biased strand distribution p-value = 6.46e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.22e-01

GGCCTGGGAAGATGCCGATTACTGCCGTGAACGGTTGGGCCTGGAGGCGCCGCTTGGCGCGATCGACCGCAGCAAACTCAATGTCAGGGGCAGTTCGCTGGCGGCCGGGCACCCGTTTGCGGCAACCGGTGGGCGGATATTGGCCAACCTGGCCAAGCTGTTGGCAACGGCGGGTAAGGGGCGCGGGCTGATTTCGATCTGTGCGGCAGGCGGGCAGGGGGTGACCGTTATCGTCGAGCGCTGAAGTACCTCCAGGGTAACGGCAGAGGCGGCGAGTGCCTCGATCATGAT  >  NC_002947/680148‑680438
                                                                                                                                                   |                                                                                                                                                
ggCCTGGGAAGATGCCGATTACTGCCGTGAACGGTTGGGCCTGGAGGCGCCGCTTGGCGCGATCGACCGCAGCAAACTCAATGTCAGGGGCAGTTCGCTGGCGGCCGGGCACCCGTTTGCGGCAACCGGTGGGCGGATATTGGCCAAcc                                                                                                                                                 >  6:265130/1‑149 (MQ=255)
   cTGGGAAGATGCCGATTACTGCCGTGAACGGTTGGGCCTGGAGGCGCCGCTTGGCGCGATCGACCGCAGCAAACTCAATGTCAGGGGCAGTTCGCTGGCGGCCGGGCACCCGTTTGCGGCAACCGGTGGGCGGATATTGGCCAACCTgg                                                                                                                                              >  2:3643/1‑149 (MQ=255)
      ggAAGATGCCGATTACTGCCGTGAACGGTTGGGCCTGGAGGCGCCGCTTGGCGCGATCGACCGCAGCAAACTCAATGTCAGGGGCAGTTCGCTGGCGGCCGGGCACCCGTTTGCGGCAACCGGTGGGCGGATATTGGCCAACCTGGc                                                                                                                                             <  1:284546/147‑1 (MQ=255)
       gAAGATGCCGATTACTGCCGTGAACGGTTGGGCCTGGAGGCGCCGCTTGGCGCGATCGACCGCAGCAAACTCAATGTCAGGGGCAGTTCGCTGGCGGCCGGGCACCCGTTTGCGGCAACCGGTGGGCGGATATTGGCCAACCTGGCCaa                                                                                                                                          <  7:156951/149‑1 (MQ=255)
       gAAGATGCCGATTACTGCCGTGAACGGTTGGGCCTGGAGGCGCCGCTTGGCGCGATCGACCGCAGCAAACTCAATGTCAGGGGCAGTTCGCTGGCGGCCGGGCACCCGTTTGCGGCAACCGGTGGGCGGATATTGGCCAACCTGGCCaa                                                                                                                                          <  7:11805/149‑1 (MQ=255)
        aaGATGCCGATTACTGCCGTGAACGGTTGGGCCTGGAGGCGCCGCTTGGCGCGATCGACCGCAGCAAACTCAATGTCAGGGGCAGTTCGCTGGCGGCCGGGCACCCGTTTGCGGCAACCGGTGGGCGGATATTGGCCAACCTGGCCAAg                                                                                                                                         >  4:325072/1‑149 (MQ=255)
            tGCCGATTACTGCCGTGAACGGTTGGGCCTGGAGGCGCCGCTTGGCGCGATCGACCGCAGCAAACTCAATGTCAGGGGCAGTTCGCTGGCGGCCGGGCACCCGTTTGCGGCAACCGGTGGGCGGATATTGGCCAACCTGGCCAAGCTGt                                                                                                                                     <  5:265130/149‑1 (MQ=255)
             gCCGATTACTGCCGTGAACGGTTGGGCCTGGAGGCGCCGCTTGGCGCGATCGACCGCAGCAAACTCAATGTCAGGGGCAGTTCGCTGGCGGCCGGGCACCCGTTTGCGGCAACCGGTGGGCGGATATTGGCCAACCTGGCCAAGCTGtt                                                                                                                                    >  5:188172/1‑149 (MQ=255)
                      tGCCGTGAACGGTTGGGCCTGGAGGCGCCGCTTGGCGCGATCGACCGCAGCAAACTCAATGTCAGGGGCAGTTCGCTGGCGGCCGGGCACCCGTTTGCGGCAACCGGTGGGCGGATATTGGCCAACCTGGCCAAGCTGTTGGCAAcggc                                                                                                                           >  5:207510/1‑149 (MQ=255)
                        ccGTGAACGGTTGGGCCTGGAGGCGCCGCTTGGCGCGATCGACCGCAGCAAACTCAATGTCAGGGGCAGTTCGCTGGCGGCCGGGCACCCGTTTGCGGCAACCGGTGGGCGGATATTGGCCAACCTGGCCAAGCTGTTGGCAAcggcgg                                                                                                                         <  3:3906/149‑1 (MQ=255)
                        ccGTGAACGGTTGGGCCTGGAGGCGCCGCTTGGCGCGATCGACCGCAGCAAACTCAATGTCAGGGGCAGTTCGCTGGCGGCCGGGCACCCGTTTGCGGCAACCGGTGGGCGGATATTGGCCAACCTGGCCAAGCTGTTGGCAAcggcg                                                                                                                          <  3:331444/148‑1 (MQ=255)
                          gTGAACGGTTGGGCCTGGAGGCGCCGCTTGGCGCGATCGACCGCAGCAAACTCAATGTCAGGGGCAGTTCGCTGGCGGCCGGGCACCCGTTTGCGGCAACCGGTGGGCGGATATTGGCCAACCTGGCCAAGCTGTTGGCAACGGCggg                                                                                                                        >  5:276990/1‑148 (MQ=255)
                             aaCGGTTGGGCCTGGAGGCGCCGCTTGGCGCGATCGACCGCAGCAAACTCAATGTCAGGGGCAGTTCGCTGGCGGCCGGGCACCCGTTTGCGGCAACCGGTGGGCGGATATTGGCCAACCTGGCCAAg                                                                                                                                         <  4:63869/128‑1 (MQ=255)
                             aaCGGTTGGGCCTGGAGGCGCCGCTTGGCGCGATCGACCGCAGCAAACTCAATGTCAGGGGCAGTTCGCTGGCGGCCGGGCACCCGTTTGCGGCAACCGGTGGGCGGATATTGGCCAACCTGGCCAAg                                                                                                                                         >  3:63869/1‑128 (MQ=255)
                             aaCGGTTGGGCCTGGAGGCGCCGCTTGGCGCGATCGACCGCAGCAAACTCAATGTCAGGGGCAGTTCGCTGGCGGCCGGGCACCCGTTTGCGGCAACCGGTGGGCGGATATTGGCCAACCTGGCCAAGCTGTTGGCAACGGCGGGTAAg                                                                                                                   <  8:20057/149‑1 (MQ=255)
                               cGGTTGGGCCTGGAGGCGCCGCTTGGCGCGATCGACCGCAGCAAACTCAATGTCAGGGGCAGTTCGCTGGCGGCCGGGCACCCGTTTGCGGCAACCGGTGGGCGGATATTGGCCAACCTGGCCAAGCTGTTGGCAACGGCGGGTAAggg                                                                                                                 <  7:93521/149‑1 (MQ=255)
                               cGGTGGGGCCGGGAGGCGCCGCTTGGCGCGATCGACCGCAGCAAACTCAATGTCAGGGGCAGTTCGCGGGCGGCCGGGCACCCGTTTGCGGCAACCGGTGGGCGGATATTGGCCAACCTGGCCAAGCTGTTGGCAACGGCGGGTAAggg                                                                                                                 <  3:155940/149‑1 (MQ=255)
                                     ggCCTGGAGGCGCCGCTTGGCGCGATCGACCGCAGCAAACTCAATGTCAGGGGCAGTTCGCTGGCGGCCGGGCACCCGTTTGCGGCAACCGGTGGGCGGATATTGGCCAACCTGGCCAAGCTGTTGGCAACGGCGGGTAAGGGGCGCgg                                                                                                           <  3:325072/149‑1 (MQ=255)
                                          ggAGGCGCCGCTTGGCGCGATCGACCGCAGCAAACTCAATGTCAGGGGCAGTTCGCTGGCGGCCGGGCACCCGTTTGCGGCAACCGGTGGGCGGATATTGGCCAACCTGGCCAAGCTGTTGGCAACGGCGGGTAAGGGGCGCGGGCTGa                                                                                                      >  2:172465/1‑149 (MQ=255)
                                              gcgcCGCTTGGCGCGATCGACCGCAGCAAACTCAATGTCAGGGGCAGTTCGCTGGCGGCCGGGCACCCGTTTGCGGCAACCGGTGGGCGGATATTGGCCAACCTGGCCAAGCTGTTGGCAACGGCGGGTAAGGGGCGCGGGCTGAttt                                                                                                   <  6:188172/148‑1 (MQ=255)
                                                gccgcTTGGCGCGATCGACCGCAGCAAACTCAATGTCAGGGGCAGTTCGCTGGCGGCCGGGCACCCGTTTGCGGCAACCGGTGGGCGGATATTGGCCAACCTGGCCAAGCTGTTGGCAACGGCGGGTAAGGGGCGCGGGCTGATTTCGa                                                                                                >  1:178575/1‑149 (MQ=255)
                                                 ccgcTTGGCGCGATCGACCGCAGCAAACTCAATGTCAGGGGCAGTTCGCTGGCGGCCGGGCACCCGTTTGCGGCAACCGGTGGGCGGATATTGGCCAACCTGGCCAAGCTGTTGGCAACGGCGGGTAAGGGGCGCGGGCTGATTTCGAt                                                                                               <  3:276831/149‑1 (MQ=255)
                                                  cgcTTGGCGCGATCGACCGCAGCAAACTCAATGTCAGGGGCAGTTCGCTGGCGGCCGGGCACCCGTTTGCGGCAACCGGTGGGCGGATATTGGCCAACCTGGCCAAGCTGTTGGCAACGGCGGGTAAGGGGCGCGGGCTGATTTCGATc                                                                                              <  6:207510/149‑1 (MQ=255)
                                                     ttGGCGCGATCGACCGCAGCAAACTCAATGTCAGGGGCAGTTCGCTGGCGGCCGGGCACCCGTTTGCGGCAACCGGTGGGCGGATATTGGCCAACCTGGCCAAGCTGTTGGCAACGGCGGGTAAgg                                                                                                                  >  6:101272/1‑126 (MQ=255)
                                                     ttGGCGCGATCGACCGCAGCAAACTCAATGTCAGGGGCAGTTCGCTGGCGGCCGGGCACCCGTTTGCGGCAACCGGTGGGCGGATATTGGCCAACCTGGCCAAGCTGTTGGCAACGGCGGGTAAgg                                                                                                                  <  5:101272/126‑1 (MQ=255)
                                                       ggCGCGATCGACCGCAGCAAACTCAATGTCAGGGGCAGTTCGCTGGCGGCCGGGCACCCGTTTGCGGCAACCGGTGGGCGGATATTGGCCAACCTGGCCAAGCTGTTGGCAACGGCGGGTACAGGGGCGCGGGCTGATTTCGATCtgt                                                                                           <  1:3643/148‑1 (MQ=255)
                                                              tCGACCGCAGCAAACTCAATGTCAGGGGCAGTTCGCTGGCGGCCGGGCACCCGTTTGCGGCAACCGGTGGGCGGATATTGGCCAACCTGGCCAAGCTGTTGGCAACGGCGGGTAAGGGGCGCGGGCTGATTTCGATCTGTGCGGCAGGc                                                                                  >  3:229569/1‑149 (MQ=255)
                                                                 aCCGCAGCAAACTCAATGTCAGGGGCAGTTCGCTGGCGGCCGGGCACCCGTTTGCGGCAACCGGTGGGCGGATATTGGCCAACCTGGCCAAGCTGTTGGCAACGGCGGGTAAGGGGCGCGGGCTGATTTCGATCTGTGCGGCAGGCggg                                                                               >  6:248377/1‑149 (MQ=255)
                                                                    gcagcaAACTCAATGTCAGGGGCAGTTCGCTGGCGGCCGGGCACCCGTTTGCGGCAACCGGTGGGCGGATATTGGCCAACCTGGCCAAGCTGTTGGCAACGGCGGGTAAGGGGCGCGGGCTGATTTCGATCTGTGCGGCAGGCGGGCAg                                                                            <  4:229569/149‑1 (MQ=255)
                                                                       gcaAACTCAATGTCAGGGTCAGTTCGCTGGCGGCCGGGGACCCGTTTGCGGCAACCGGTGGGCGGATATTGGCCAACCTGGCCTAGCTGTTGGCAACGGCGGGTAAGGGGCGCGGGCTGATTTCGATCTGTGCGGCAGGCGGGCAgggg                                                                         <  7:320899/149‑1 (MQ=255)
                                                                          aaCTCAATGTCAGGGGCAGTTCGCTGGCGGCCGGGCACCCGTTTGCGGCAACCGGTGGGCGGATATTGGCCAACCTGGCCAAGCTGTTGGCAACGGCGGGTAAGGGGCGCGGGCTGATTTCGATCTGTGCGGCAGGCGGGCAGGGGGTg                                                                      <  1:172465/149‑1 (MQ=255)
                                                                           aCTCAATGTCAGGGGCAGTTCGCTGGCGGCCGGGCACCCGTTTGCGGCAACCGGTGGGCGGATATTTGCAAACCTGGCCAAGCTGTTGGCAACGGCGGGTAAGGGGCGCGGGCTGATTTCGATCTGTGCGGCAGGCGGGCAGGGGGTGa                                                                     <  2:86077/149‑1 (MQ=255)
                                                                               aaTGTCAGGGGCAGTTCGCTGGCGGCCGGGCACCAGTTTGCGGCAACCGGTGGGCGGATATTGGCCAACCTGGCCAAGCTGTTGGCAACGGCGGGTAAGGGGCGCGGGCTGATTTCGATCTGTGCGGCAGGCGGGCAGGGGGAGACCGt                                                                 >  8:232978/1‑149 (MQ=255)
                                                                                      ggggCAGTTCGCTGGCGGCCGGGCACCCGTTTGCGGCAACCGGTGGGCGGATATTGGCCAACCTGGCCAAGCTGTTGGCAACGGCGGGTAAGGGGCGCGGGCTGATTTCGATCTGTGCGGCAGGCGGGCAGGGGGTGACCGTTAtcgtc                                                          <  5:248377/149‑1 (MQ=255)
                                                                                       gggCAGTTCGCTGGCGGCCGGGCACCCGTTTGCGGCAACCGGTGGGCGGATATTGGCCAACCTGGCCAAGCTGTTGGCAACGGCGGGTAAGGGGCGCGGGCTGATTTCGATCTGTGCGGCAGGCGGGCAGGGGGTGACCGTTAtcgtc                                                          >  2:5673/1‑148 (MQ=255)
                                                                                                gCTGGCGGCCGGGCACCCGTTTGCGGCAACCGGTGGGCGGATATTGGCCAACCTGGCCAAGCTGTTGGCAACGGCGGGTAAGGGGCGCGGGCTGATTTCGATCTGTGCGGCAGGCGGGCAGGGGGTGACCGTTATCGTCGAGCGCTg                                                  >  2:4032/1‑147 (MQ=255)
                                                                                                       gcCGGGCACCCGTTTGCGGCAACCGGTGGGCGGATATTGGCCAACCTGGCCAAGCTGTTGGCAACGGCGGGTAAGGGGCGCGGGCTGATTTCGATCTGTGCGGCAGGCGGGCAGGGGGTGACCGTTATCGTCGAGCGCTGAAGTACCTc                                         >  2:219365/1‑149 (MQ=255)
                                                                                                              aCCCGTTTGCGGCAACCGGTGGGCGGATATTGGCCAACCTGGCCAAGCTGTTGGCAACGGCGGGTAAGGGGCGCGGGCTGATTTCGATCTGTGCGGCAGGCGGGCAGGGGGTGACCGTTATCGTCGAGCGCTGAAGTACCTCCAGGGt                                   >  5:277768/1‑148 (MQ=255)
                                                                                                               cccGTTTGCGGCAACCGGTGGGCGGATATTGGCCAACCTGGCCAAGCTGTTGGCAACGGCGGGTAAGGGGCGCGGGCTGATTTCGATCTTTGAGGCAGGCGGGCAGGGGGTGACCGTTAACGTCGAGCGCTGAAGTACCTCCAGGGTaa                                 >  7:195756/1‑149 (MQ=255)
                                                                                                                       cGGCAACCGGTGGGCGGATATTGGCCAACCTGGCCAAGCTGTTGGCAACGGCGGGTAAGGGGCGCGGGCTGATTTCGATCTGTGCGGCAGGCGGGCAGGGGGTGACCGTTATCGTCGAGCGCTGAAGTACCTCCAGGGTAACGGCagag                         >  8:89250/1‑149 (MQ=255)
                                                                                                                          cAACCGGTGGGCGGATATTGGCCAACCTGGCCAAGCTGTTGGCAACGGCGGGTAAGGGGCGCGGGCTGATTTCGATCTGTGCGGCAGGCGGGCAGGGGGTGACCGTTATCGTCGAGCGCTGAAGTACCTCCAGGGTAACGGCAGAggc                       >  6:276228/1‑148 (MQ=255)
                                                                                                                          cAACCGGTGGGCGGATATTGGCCAACCTGGCCAAGCTGTTGGCAACGGCGGGTAAGGGGCGCGGGCTGATTTCGATCTGTGCGGCAGGCGGGCAGGGGGTGACCGTTATCGTCGAGCGCTGAAGTACCTCCAGGGTAACGGCAGAggc                       <  5:276228/148‑1 (MQ=255)
                                                                                                                            aCCGGTGGGCGGATATTGGCCAACCTGGCCAAGCTGTTGGCAACGGCGGGTAAGGGGCGCGGGCTGATTTCGATCTGTGCGGCAGGCGGGCAGGGGGTGACCGTTATCGTCGAGCGCTGAAGTACCTCCAGGGTAACGGCAGAggcggc                    <  7:89250/149‑1 (MQ=255)
                                                                                                                              cGGTGGGCGGATATTGGCCAACCTGGCCAAGCTGTTGGCAACGGCGGGTAAGGGGCGCGGGCTGATTTCGATCTGTGCGGCAGGCGGGCAGGGGGTGACCGTTATCGTCGAGCGCTGAAGTACCTCCAGGGTAACGGCAGAGGCGGCGa                  <  6:276990/149‑1 (MQ=255)
                                                                                                                                gTGGGCGGATATTGGCCAACCTGGCCAAGCTGTTGGCAACGGCGGGTAAGGGGCGCTGGCTGATTTCGATCTGTGCTGCAGTCGGGCAGGGGGTGACCGTTATCGTCGAGCGCTGAAGTACCTCCAGGGTAACTGCAGAGGCGGCGAGt                <  2:178575/149‑1 (MQ=255)
                                                                                                                                gTGGGCGGATATTGGCCAACCTGGCCAAGCTGTTGGCAACGGCGGGTAAGGGGCGCGGGCTGATTTCGATCTGTGCGGCAGGCGGGCAGGGGGTGACCGTTATCGTCGAGCGCTGAAGTACCTCCAGGGTAACGGCAGAGGCGGCGAGt                <  1:5673/149‑1 (MQ=255)
                                                                                                                                    caacagcttGGCCAGGTTGGCCAAGCTGTTGGCAACGGCGGGTAAGGGGCGCGGGCTGATTTCGATCTGTGCGGCAGGCGGGCAGGGGGTGACCGTTATCGTCGAGCGCTGAAGTACCTCCAGGGTAACGGCAGAGGCGGCGAGTGCCt            >  2:126011/8‑149 (MQ=255)
                                                                                                                                      ggATATTGGCCAACCTGGCCAAGCTGTTGGCAACGGCGGGTAAGGGGCGCGGGCTGATTTCGATCTGTGCGGCAGGCGGGCAGGGGGTGACCGTTATCGTCGAGCGCTGAAGTACCTCCAGGGTAACGGCAGAGGCGGCGAGTGCCTc           >  7:42725/1‑148 (MQ=255)
                                                                                                                                          cttGGCCAGGTTGGCCAAGCTGTTGGCAACGGCGGGTAAGGGGCGTGGGCTGATTTCGATCTGTGCGGCAGGCGGGCAGGGGGTGACCGTTATCGTCGAGCGCTGAAGTACCTCCAGGGTAACGGCAGAGGCGGCGAGTGCCTCGATCa      >  7:32547/2‑149 (MQ=255)
                                                                                                                                          cttGGCCAGGTTGGCCAAGCTGTTGGCAACGGCGGGTAAGGGGCGTGGGCTGATTTCGATCTGTGCGGCAGGCGGGCAGGGGGTGACCGTTATCGTCGAGCGCTGAAGTACCTCCAGGGTAACGGCAGAGGCGGCGAGTGCCTCGATCa      <  8:32547/148‑1 (MQ=255)
                                                                                                                                          cttGGCCAGGTTGGCCAAGCTGTTGGCAACGGCGGGTAAGGGGCGCGGGCTGATTTCGATCTGTGCGGCAGGCGGGCAGGGGGTGACCGTTATCGTCGAGCGCTGAAGTACCTCCAGGGTAACGGCAGAGGCGGCGAg                 >  2:241628/2‑138 (MQ=255)
                                                                                                                                          cttGGCCAGGTTGGCCAAGCTGTTGGCAACGGCGGGTAAGGGGCGCGGGCTGATTTCGATCTGTGCGGCAGGCGGGCAGGGGGTGACCGTTAACGTCGAGCGCTGAAGTACCTCCAGGGTAACGGCAGAGGCGGCGAg                 <  1:241628/137‑1 (MQ=255)
                                                                                                                                           ttGGCCAACCTGGCCAAGCTGTTGGCAACGGCGGGTAAGGGGCGCGGGCTGATTTCGATCTGTGCGGCAGGCGGGCAGGGGGTGACCGTTATCGTCGAGCGCTGAAGTACCTCCAGGGTAACGGCAGAGGCGGCGAGTGCCTCGATCa      <  8:42725/148‑1 (MQ=255)
                                                                                                                                              gCCAACCTGGCCAAGCTGTTGGCAACGGCGGGTAAGGGGCGCGGGCTGATTTCGATCTGTGCGGCAGGCGGGCAGGGGGTGACCGTTATCGTCGAGCGCTGAAGTACCTCCAGGGTAACGGCAGAGGCGGCGAGTGCCTCGCTCATGAt  >  8:111052/1‑149 (MQ=255)
                                                                                                                                              gCCAACCTGGCCAAGCTGTTGGCAACGGCGGGTAAGGGGCGCGGGCTGATTTCGATCTGTGCGGCAGGCGGGCAGGGGGTGACCGTTATCGTCGAGCGCTGAAGTACCTCCAGGGTAACGGCAGAGGCGGCGAGTGCCTCGATCATGa   <  1:4032/148‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                
GGCCTGGGAAGATGCCGATTACTGCCGTGAACGGTTGGGCCTGGAGGCGCCGCTTGGCGCGATCGACCGCAGCAAACTCAATGTCAGGGGCAGTTCGCTGGCGGCCGGGCACCCGTTTGCGGCAACCGGTGGGCGGATATTGGCCAACCTGGCCAAGCTGTTGGCAACGGCGGGTAAGGGGCGCGGGCTGATTTCGATCTGTGCGGCAGGCGGGCAGGGGGTGACCGTTATCGTCGAGCGCTGAAGTACCTCCAGGGTAACGGCAGAGGCGGCGAGTGCCTCGATCATGAT  >  NC_002947/680148‑680438

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