Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,568,449 C→T 8.1% G311G (GGC→GGT)  pcaK → 4‑hydroxybenzoate transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,568,4490CT 8.1% 141.4 / 2.9 49G311G (GGC→GGT) pcaK4‑hydroxybenzoate transporter
Reads supporting (aligned to +/- strand):  ref base C (23/22);  new base T (2/2);  total (25/24)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.21e-01

GGGCACCGTGTTGCTGTGGCTGACCTATTTCATGGGCCTGGTGATCGTCTACCTCCTGACCAGCTGGCTGCCAACCCTGATGCGTGACAGCGGCGCGAGCATGGAGCAGGCGGCGTTCATCGGCGCGCTGTTCCAGTTTGGCGGCGTGCTAAGTGCAGTCGGTGTGGGTTGGGCCATGGACCGTTTCAATC‑CGCACAAGGTAATTGGCATTTTCTATCTGCTGGCCGGGGTGTTCGCCTACGCAGTGGGGCAGAGCCTGGGCAACATCACTCTGCTAGCCAC  >  NC_002947/1568305‑1568586
                                                                                                                                                |                                                                                                                                          
gggCACCGTGTTGCTGTGGCTGACCTATTTCATGGGCCTGGTGATCGTCTACCTCCTGACCAGCTGGCTGCCAACCCTGATGCGTGACAGCGGCGCGAGCATGGAGCAGGCGGCGTTCATCGGCGCGCTGTTCCAGTTTGGCGGCGt                                                                                                                                          <  8:33714/147‑1 (MQ=255)
 ggCACCGTGTTGCTGTGGCTGACCTATTTCATGGGCCTGGTGATCGTCTACCTCCTGACCAGCTGGCTGCCAACCCTGATGCGTGACAGCGGCGCGAGCATGGAGCAGGCGGCGTTCATCGGCGCGCTGTTCCAGTTTGGCGGCGTGCt                                                                                                                                       <  5:365959/149‑1 (MQ=255)
  gCACCGTGTTGCTGTGGCTGACCTATTTCATGGGCCTGGTGATCGTCTACCTCCTGACCAGCTGGCTGCCAACCCTGATGCGTGACAGCGGCGCGAGCATGGAGCAGGCGGCGTTCATCGGCGCGCTGTTCCAGTTTGGCGGCGTGCt                                                                                                                                       >  2:111974/1‑148 (MQ=255)
   cACCGTGTTGCTGTGGCTGACCTATTTCATGGGCCTGGTGATCGTCTACCTCCTGACCAGCTGGCTGCCAACCCTGATGCGTGACAGCGGCGCGAGCATGGAGCAGGCGGCGTTCATCGGCGCGCTGTTCCAGTTTGGCGGCGTGCt                                                                                                                                       >  3:48514/1‑147 (MQ=255)
     ccGTGTTGCTGTGGCTGACCTATTTCATGGGCCTGGTGATCGTCTACCTCCTGACCAGCTGGCTGCCAACCCTGATGCGTGACAGCGGCGCGAGCATGGAGCAGGCGGCGTTCATCCGCGCGCTGTTCCAGTTTGGCGGCGTGCTAAGt                                                                                                                                   <  3:14309/149‑1 (MQ=255)
            gCTGTGGCTGACCTATTTCATGGGCCTGGTGATCGTCTACCTCCTGACCAGCTGGCTGCCAACCCTGATGCGTGACAGCGGCGCGAGCATGGAGCAGGCGGCGTTCATCGGCGCGCTGTTCCAGTTTGGCGGTGTGCTAAGTGCAGTCg                                                                                                                            >  4:234407/1‑149 (MQ=255)
              tgtgGCTGACCTATTTCATGGGCCTGGTGATCGTCTACCTCCTGACCAGCTGGCTGCCAACCCTGATGCGTGACAGCGGCGCGAGCATGGAGCAGGCGGCGTTCATCGGCGCGCTGTTCCAGTTTGGCGGCGt                                                                                                                                          >  5:20057/1‑133 (MQ=255)
              tgtgGCTGACCTATTTCATGGGCCTGGTGATCGTCTACCTCCTGACCAGCTGGCTGCCAACCCTGATGCGTGACAGCGGCGCGAGCATGGAGCAGGCGGCGTTCATCGGCGCGCTGTTCCAGTTTGGCGGCGt                                                                                                                                          <  6:20057/133‑1 (MQ=255)
                        cTATTTCATGGGCCTGGTGATCGTCTACCTCCTGACCAGCTGGCTGCCAACCCTGATGCGTGACAGCGGCGCGAGCATGGAGCAGGCGGCGTTCATCGGCGCGCTGTTCCAGTTTGGCGGTGTGCTAAGTGCAGTCGGTGTGGGTTggg                                                                                                                >  3:144439/1‑149 (MQ=255)
                        cTATTTCATGGGCCTGGTGATCGTCTACCTCCTGACCAGCTGGCTGCCAACCCTGATGCGTGACAGCGGCGCGAGCATGGAGCAGGCGGCGTTCATCGGCGCGCTGTTCCAGTTTGGCGGCGTGCTAAGt                                                                                                                                   <  8:106863/130‑1 (MQ=255)
                        cTATTTCATGGGCCTGGTGATCGTCTACCTCCTGACCAGCTGGCTGCCAACCCTGATGCGTGACAGCGGCGCGAGCATGGAGCAGGCGGCGTTCATCGGCGCGCTGTTCCAGTTTGGCGGCGTGCTAAGt                                                                                                                                   >  7:106863/1‑130 (MQ=255)
                        cTATTTCATGGGCCTGGTGATCGTCTACCTCCTGACCAGCTGGCTGCCAACCCTGATGCGTGACAGCGGCGCGAGCATGGAGCAGGCGGCGTTCATCGGCGCGCTGTTCCAGTTTGGCGGCGTGCTAAGTGCAGTCGGTGTGGGTTggg                                                                                                                >  1:133453/1‑149 (MQ=255)
                                     cTGGTGATCGTCTACCTCCTGACCAGCTGGCTGCCAACCCTGATGCGTGACAGCGGCGCGAGCATGGAGCAGGCGGCGTTCATCGGCGCGCTGTTCCAGTTTGGCGGTGTGCTAAGTGCAGTCGGTGTGGGTTGGGCCATGGACCGttt                                                                                                   <  3:234407/149‑1 (MQ=255)
                                       ggTGATCGTCTACCTCCTGACCAGCTGGCTGCCAACCCTGATGCGTGACAGCGGCGCGAGCATGGAGCAGGCGGCGTTCATCGGCGCGCTGTTCCAGTTTGGCGGCGTGCTAAGTGCAGTCGGTGTGGGTTGGGCCATGGACCGTTTc                                                                                                  >  4:124053/1‑148 (MQ=255)
                                        gTGATCGTCTACCTCCTGACCAGCTGGCTGCCAACCCTGATGCGTGACAGCGGCGCGAGCATGGAGCAGGCGGCGTTCATCGGCGCGCTGTTCCAGTTTGGCGGCGTGCTAAGTGCAGTCGGTGTGGGTTGGGCCATGGACCGTTTCaa                                                                                                <  5:372222/149‑1 (MQ=255)
                                           aTCGTCTACCTCCTGACCAGCTGGCTGCCAACCCTGATGCGTGACAGCGGCGCGAGCATGGAGCAGGCGGCGTTCATCGGCGCGCTGTTCCAGTTTGGCGGCGTGCTAAGTGCAGTCGGTGTGGGTTGGGCCATGGACCGTTTCAATc‑c                                                                                            <  4:38702/149‑1 (MQ=255)
                                              gTCTACCTCCTGACCAGCTGGCTGCCAACCCTGATGCGTGACAGCGGCGCGAGCATGGAGCAGGCGGCGTTCATCGGCGCGCTGTTCCAGTTTGGCGGTGTGCTAAGTGCAGTCGGTGTGGGTTGGGCCATGGACCGTTTCAATC‑CGca                                                                                         <  4:144439/149‑1 (MQ=255)
                                              gTCTACCTCCTGACCAGCTGGCTGCCAACCCTGATGCGTGACAGCGGCGCGAGCATGGAGCAGGCGGCGTTCATCGGCGCGCTGTTCCAGTTTGGCGGCGTGCTAAGTGCAGTCGGTGTGGGTTGGGCCATGGACCGTTTCAATC‑CGca                                                                                         >  7:150114/1‑149 (MQ=255)
                                              gTCTACCTCCTGACCAGCTGGCTGCCAACCCTGATGCGTGACAGCGGCGCGAGCATGGAGCAGGCGGCGTTCATCGGCGCGCTGTTCCAGTTTGGCGGCGTGCTAAGTGCAGTCGGTGTGGGTTGGGCCATGGACCGTTTCAATC‑CGca                                                                                         >  1:68979/1‑149 (MQ=255)
                                                  aCCTCCTGACCAGCTGGCTGCCAACCCTGATGCGTGACAGCGGCGCGAGCATGGAGCAGGCGGCGTTCATCGGCGCGCTGTTCCAGTTTGGCGGCGTGCTAAGTGCAGTCGGTGTGGGTTGGGCCATGGACCGtt                                                                                                    >  3:305111/1‑135 (MQ=255)
                                                  aCCTCCTGACCAGCTGGCTGCCAACCCTGATGCGTGACAGCGGCGCGAGCATGGAGCAGGCGGCGTTCATCGGCGCGCTGTTCCAGTTTGGCGGCGTGCTAAGTGCAGTCGGTGTGGGTTGGGCCATGGACCGtt                                                                                                    <  4:305111/135‑1 (MQ=255)
                                                                    tgCCAACCCTGATGCGTGACAGCGGCGCGAGCATGGAGCAGGCGGCGTTCATCGGCGCGCTGTTCCAGTTTGGCGGCGTGCTAAGTGCAGTCGGTGTGGGTTGGGCCATGGACCGTTTCAATC‑CGCACGAGGTAATTGGCATTTTCTAt                                                                   >  6:248326/1‑149 (MQ=255)
                                                                      ccAACCCTGATGCGTGACAGCGGCGCGAGCATGGAGCAGGCGGCGTTCATCGGCGCGCTGTTCCAGTTTGGCGGCGTGCTAAGTGCAGTCGGTGTGGGTTGGGCCATGGACCGTTTCAATC‑CGCACAAGGTAATTGGCATTTTCTAt                                                                   <  2:68979/147‑1 (MQ=255)
                                                                        aaCCCTGATGCGTGACAGCGGCGCGAGCATGGAGCAGGCGGCGTTCATCGGCGCGCTGTTCCAGTTTGGCGGCGTGCTAAGTGCAGTCGGTGTGGGTTGGGCCATGGACCGTTTCAATC‑CGCACAAGGTAATTGGCATTTTCTATctgc                                                               <  3:124053/149‑1 (MQ=255)
                                                                         aCCCTGATGCGTGACAGCGGCGCGAGCATGGAGCAGGCGGCGTTCATCGGCGCGCTGTTCCAGTTTGGCGGCGTGCTAAGTGGAGTCGGTGTGGGTTGGGCCATGGa                                                                                                         <  6:98532/107‑1 (MQ=255)
                                                                         aCCCTGATGCGTGACAGCGGCGCGAGCATGGAGCAGGCGGCGTTCATCGGCGCGCTGTTCCAGTTTGGCGGCGTGCTAAGTGCAGTCGGTGTGGGTTGGGCCATGGa                                                                                                         >  5:98532/1‑107 (MQ=255)
                                                                          cccTGATGCGTGACAGCGGCGCGAGCATGGAGCAGGCGGCGTTCATCGGCGCGCTGTTCCAGTTTGGCGGCGTGCTAAGTGCAGTCGGTGTGGGTTGGGCCATGGACCGTTTCAATC‑CGcac                                                                                        >  7:91648/1‑122 (MQ=255)
                                                                          cccTGATGCGTGACAGCGGCGCGAGCATGGAGCAGGCGGCGTTCATCGGCGCGCTGTTCCAGTTTGGCGGCGTGCTAAGTGCAGTCGGTGTGGGTTGGGCCATGGACCGTTTCAATC‑CGcac                                                                                        <  8:91648/122‑1 (MQ=255)
                                                                           ccTGATGCGTGACAGCGGCGCGAGCATGGAGCAGGCGGCGTTCATCGGCGCGCTGTTCCAGTTTGGCGGCGTGCTAAGTg                                                                                                                                  >  1:221828/1‑80 (MQ=255)
                                                                           ccTGATGCGTGACAGCGGCGCGAGCATGGAGCAGGCGGCGTTCATCGGCGCGCTGTTCCAGTTTGGCGGCGTGCTAAGTg                                                                                                                                  <  2:221828/80‑1 (MQ=255)
                                                                           ccTGATGCGTGACAGCGGCGAGAGCATGGAGCAGGCGGCGTTCATCGGCGCGCTGTTCCAGTTTGGCGGCGTGCTAAGTGCAGTCGGTGTGGGTTGGGCCATGGACCGTTTCAATC‑CGCACAAGGTAATTGGCATTTTCTATCTGCTgg                                                            >  6:95551/1‑149 (MQ=255)
                                                                                tGCGTGACAGCGGCGCGAGCATGGAGCAGGCGGCGTTCATCGGCGCGCTGTTCCAGTTTGGCGGCGTGCTAAGTGCAGTCGGTGTGGGTTGGGCCATGGACCGTTTCAATC‑CGCACAAGGTAATTGGCATTTTCTATCTGCTGGCCggg                                                       <  5:95551/149‑1 (MQ=255)
                                                                                   gTGACAGCGGCGCGAGCATGGAGCAGGCGGCGTTCATCGACGCGCTGTTCCAGTTTGGCGGCGTGCTAAGTGCAGTCGGTGTGGGTTGGGCCATGGACCGTTTCAATC‑CGCACAAGGTAATTGGCATTTTCTATCTGCTGGCCGGGgtg                                                    >  3:190091/1‑149 (MQ=255)
                                                                                      aCAGCGGCGCGAGCATGGAGCAGGCGGCGTTCATCGGCGCGCTGTTCCAGTTTGGCGGCGTGCTAAGTGCAGTCGGTGTGGGTTGGGCCATGGACCGTTTCAATC‑CGCACAAGGTAATTGGCATTTTCTATCTGCTGGCCGGGGTGTTc                                                 <  2:196049/149‑1 (MQ=255)
                                                                                        aGCGGCGCGAGCATGGAGCAGGCGGCGTTCATCGGCGCGCTGTTCCAGTTTGGCGGCGTGCTAAGTGCAGTCGGTGTGGGTTGGGCCATGGACCGTTTCAATC‑CGCACAAGGTAATTGGCATTTTCTATCTGCTGGCCGGGGTGTTCGc                                               <  4:48514/149‑1 (MQ=255)
                                                                                             cgcgAGCATGGAGCAGGCGGCGTTCATCGGCGCGCTGTTCCAGTTTGGCGGCGTGCTAAGTGCAGTCGGTGTGGGTTGGGCCATGGACCGTTTCAATC‑CGCACAAGGTAATTGGCATTTTCTATCTGCTGGCCGGGGTGTTCGCCTACg                                          >  4:31401/1‑149 (MQ=255)
                                                                                                         gCAGGCGGCGTTCATCGGCGCGCTGTTCCAGTTTGGCGGCGTGCTAAGTGCAGTCGGTGTGGGTTGGGCCATGGACCGTTTCAATC‑CGCACAAGGTAATTGGCATTTTCTATCTGCTGGCCGGGgtg                                                    <  2:32285/127‑1 (MQ=255)
                                                                                                         gCAGGCGGCGTTCATCGGCGCGCTGTTCCAGTTTGGCGGCGTGCTAAGTGCAGTCGGTGTGGGTTGGGCCATGGACCGTTTCAATC‑CGCACAAGGTAATTGGCATTTTCTATCTGCTGGCCGGGgtg                                                    >  1:32285/1‑127 (MQ=255)
                                                                                                                       tCGGCGCGCTGTTCCAGTTTGGCGGCGTGCTAAGTGCAGTCGGTGTGGGTTGGGCCATGGACCGTTTCAATC‑CGCACAAGGTAATTGGCATTTTCTATCTGCTGGCCGGGGTGTTCGCCTACGCAGTGGGGCAGAGCCTgg                        >  3:163144/1‑141 (MQ=255)
                                                                                                                       tCGGCGCGCTGTTCCAGTTTGGCGGCGTGCTAAGTGCAGTCGGTGTGGGTTGGGCCATGGACCGTTTCAATC‑CGCACAAGGTAATTGGCATTTTCTATCTGCTGGCCGGGGTGTTCGCCTACGCAGTGGGGCAGAGCCTgg                        <  4:163144/141‑1 (MQ=255)
                                                                                                                        cGGCGCGCTGTTCCAGTTTGGCGGCGTGCTAAGTGCAGTCGGTGTGGGTTGGGCCATGGACCGTTTCAATC‑CGCACAAGGTAATTGGCATTTTCTATCTGCTGGCCGGGGTGTTCGCCTACGCAGTGGGGCAGAGCCTGGGCAACATc                >  6:123991/1‑148 (MQ=255)
                                                                                                                        cGGCGCGCTGTTCCAGTTTGGCGGCGTGCTAAGTGCAGTCGGTGTGGGTTGGGCCATGGACCGTTTCAATC‑CGCACAAGGTAATTGGCATTTTCTATCTGCTGGCCGGGGTGTTCGCCTACGCAGTGGGGCAGAGCCTGGGCAACATCa               <  1:111974/149‑1 (MQ=255)
                                                                                                                           cgcgcTGTTCCAGTTTGGCGGCGTGCTAAGTGCAGTCGGTGTGGGTTGGGCCATGGACCGTTTCAATC‑CGCACAAGGTAATTGGCATTTTCTATCTGCTGGCCGGGGTGTTCGCCTACGCAGTGGGGCAGAGCCTGGGCAACATCActc            <  5:123991/149‑1 (MQ=255)
                                                                                                                            gcgcTGTTCCTGTTTGGCGGCGTGCTAAGTGCAGTCGGTGTGGGTTGGGCCATGGACCGTTTCAATC‑CGCACAAGGTAATTGGCATTTTCTATCTGCTGGCCGGGGTGTTCGCCTACGCAGTGGGGCAGAGCCTGGGCAACATCActct           <  8:150114/149‑1 (MQ=255)
                                                                                                                                 gTTCCAGTTTGGCGGCGTGCTAAGTGCAGTCGGTGTGGGTTGGGCCATGGACCGTTTCAATC‑CGCACAAGGTAATTGGCATTTTCTATCTGCTGGCCGGGGTGTTCGCCTACGCAGTGGGGCAGAGCCTgg                        >  7:70954/1‑131 (MQ=255)
                                                                                                                                 gTTCCAGTTTGGCGGCGTGCTAAGTGCAGTCGGTGTGGGTTGGGCCATGGACCGTTTCAATC‑CGCACAAGGTAATTGGCATTTTCTATCTGCTGGCCGGGGTGTTCGCCTACGCAGTGGGGCAGAGCCTgg                        <  8:70954/131‑1 (MQ=255)
                                                                                                                                  ttCCAGTTTGGCGGCGTGCTAAGTGCAGTCGGTGTGGGTTGGGCCATGGACCGTTTCAATCTGGCACAAGGTAATTGGCATTTTCTATCTGATGGCCGGGGTGTTCGCCTACGAAGTGGGGAAGAGCCTGGGCAACATCAGTATGCTAg      >  1:298478/1‑149 (MQ=255)
                                                                                                                                  ttCCAGTTTGGCGGCGTGCTAAGTGCAGTCGGTGTGGGTTGGGCCATGGACCGTTTCAATC‑CGCACAAGGTAATTGGCATTTTCTATCTGCTGGCCGGGGTGTTCGCCTACGCAGTGGGGCAGAGCCTGGGCAACATCACTCTGCTAGc     >  8:214571/1‑149 (MQ=255)
                                                                                                                                     cAGTTTGGCGGCGTGCTAAGTGCAGTCGGTGTGGGTTGGGCCATGGACCGGTTCAATC‑CGCACAAGGTAATTGGCATTTTCTATCTGCTGGCCGGGGTGTTCGCCTACGCAGTGGGGCAGAGCCTGGGCAACATCACTCTGCTAGCCAc  >  6:341033/1‑149 (MQ=255)
                                                                                                                                                |                                                                                                                                          
GGGCACCGTGTTGCTGTGGCTGACCTATTTCATGGGCCTGGTGATCGTCTACCTCCTGACCAGCTGGCTGCCAACCCTGATGCGTGACAGCGGCGCGAGCATGGAGCAGGCGGCGTTCATCGGCGCGCTGTTCCAGTTTGGCGGCGTGCTAAGTGCAGTCGGTGTGGGTTGGGCCATGGACCGTTTCAATC‑CGCACAAGGTAATTGGCATTTTCTATCTGCTGGCCGGGGTGTTCGCCTACGCAGTGGGGCAGAGCCTGGGCAACATCACTCTGCTAGCCAC  >  NC_002947/1568305‑1568586

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