Predicted mutation
evidence seq id position mutation freq annotation gene description
RA ACB143 3,387,164 C→T 87.3% A95T (GCC→ACC)  pcpR_4 ← PCP degradation transcriptional activationprotein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*ACB1433,387,1640CT87.3% 136.9 / 14.5 57A95T (GCC→ACC) pcpR_4PCP degradation transcriptional activationprotein
Reads supporting (aligned to +/- strand):  ref base C (3/4);  new base T (23/24);  total (26/31)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.37e-01

CGTCGATCCGCGCCGCCTCCAGGCGGATGCCCGGCGCATGGCTGCGCAGATGGTTGAGGATGCTGGGCAGCAAGGTGATGTGGCTGGCATCGCTGACACACAGGCGAAAACGCCGCTGCGCTGTGGCGGCGATGAACTGCGGCTCCCAGGCGGTCAGCCGGCGCAGCGATTCCAGCACCTCGCGGCACGGCCCTATCAGCTGATCGGCGCGCGGGGTAGGGGCCATGCCGCCTGGCGTGCGCACGAACAGCG  >  ACB143/3387038‑3387289
                                                                                                                              |                                                                                                                             
cGTCGATCCGCGCCGCCTCCAGGCGGATGCCCGGCGCATGGCTGCGCAGATGGTTGAGGATGCTGGGCAGCAAGGTGATGTGGCTGGCATCGCTGACACACAGGCGAAAACGCCGCTGCGCTGTGGTGGCGATGaa                                                                                                                      <  1:222897/136‑1 (MQ=255)
 gTCGATCCGCGCCGCCTCCAGGCGGATGCCCGGCGCATGGCTGCGCAGATGGTTGAGGATGCTGGGCAGCAAGGTGATGTGGCTGGCATCGCTGACACACAGGCGAAAACGCCGCTGCGCTGTGGTGGCGATGAAc                                                                                                                     >  5:122629/1‑136 (MQ=255)
     aTCCGCGCCGCCTCCAGGCGGATGCCCGGCGCATGGCTGCGCAGATGGTTGAGGATGCTGGGCAGCAAGGTGATGTGGCTGGCATCGCTGACACACAGGCGAAAACGCCGCTGCGCTGTGGTGGCGATGAACTGCg                                                                                                                 >  3:281731/1‑136 (MQ=255)
       ccGCGCCGCCTCCAGGCGGATGCCCGGCGCATGGCTGCGCAGATGGTTGAGGATGCTGGGCAGCAAGGTGATGTGGCTGGCATCGCTGACACACAGGCGAAAACGCCGCTGCGCTGTGGCGGCGATGAACTGCGGc                                                                                                               >  6:132069/1‑136 (MQ=255)
                      gCGGATGCCCGGCGCATGGCTGCGCAGATGGTTGAGGATGCTGGGCAGCAAGGTGATGTGGCTGGCATCGCTGACACACAGGCGAAAACGCCGCTGCGCTGTGGTGGCGATGAACTGCGGCTCCCAGGCGGTCAGc                                                                                                >  1:21175/1‑136 (MQ=255)
                             cccGGCGCATGGCTGCGCAGATGGTTGAGGATGCTGGGCAGCAAGGTGATGTGGCTGGCATCGCTGACACACAGGCGAAAACGCCGCTGCGCTGTGGCGGCGATGAACTGCGGCTCCCAGGCGGTCAGCCGgcgc                                                                                          <  5:107307/135‑1 (MQ=255)
                                ggCGCATGGCTGCGCAGATGGTTGAGGATGCTGGGCAGCAAGGTGATGTGGCTGGCATCGCTGACACACAGGCGAAAACGCCGCTGCGCTGTGGTGGCGATGAACTGCGGCTCCCAGGCGGTCAGCCGGCGCAGc                                                                                       >  2:187286/1‑135 (MQ=255)
                                               aGATGGTTGAGGATGCTGGGCAGCAAGGTGATGTGGCTGGCATCGCTGACACACAGGCGAAAACGCCGCTGCGCTGTGGTGGCGATGAACTGCGGCTCCCAGGCGGTCAGCCGGCGCAGCGATTCCAGCACCTcgc                                                                       <  3:32864/136‑1 (MQ=255)
                                               aGATGGTTGAGGATGCTGGGCAGCAAGGTGATGTGGCTGGCATCGCTGACACACAGGCGAAAACGCCGCTGCGCTGTGGTGGCGATGAACTGCGGCTCCCAGGCGGTCAGCCGGCGCAGCGATTCCAGCACCTcgc                                                                       <  3:173351/136‑1 (MQ=255)
                                               aGATGGTTGAGGATGCTGGGCAGCAAGGTGATGTGGCTGGCATCGCTGACACACAGGCGAAAACGCCGCTGCGCTGTGGTGGCGATGAACTGCGGCTCCCAGGCGGTCAGCCGGCGCAGCGATTCCAGCACCTcgc                                                                       >  3:105388/1‑136 (MQ=255)
                                               aGATGGTTGAGGATGCTGGGCAGCAAGGTGATGTGGCTGGCATCGCTGACACACAGGCGAAAACGCCGCTGCGCTGTGGTGGCGATGAACTGCGGCTCCCAGGCGGTCAGCCGGCGCAGCGATTCCAGCACCTcgc                                                                       >  7:118028/1‑136 (MQ=255)
                                               aGATGGTTGAGGAGGCTGGGCAGCAAGGTGATGTGGCTGGCATCGCTGACACACAGGCGAAAACGCCGCTGCGCTGTGGCGGCGATGAACTGCGGCTCCCAGGCGGTCAGCCGGCGCAGCGATTCCAGCACCTcgc                                                                       <  7:201667/136‑1 (MQ=255)
                                                  tggtggAGGATGCTGGGCAGCAAGGTGATGTGGCTGGCATCGCGGACACACAGGCGAAAACGCCGCTGCGCTGGGGTGGCGATGAACTGCGGCTCCCAGGCGGGCAGCCGGCGCAGCGATTCCAGCACCTCGCGGc                                                                    <  5:141027/136‑1 (MQ=255)
                                                  tGGTTGAGGATGCTGGGCAGCAAGGTGATGTGGCTGGCATCGCTGACACACAGGCGAAAACGCCGCTGCGCTGTGGTGGCGATGAACTGCGGCTCCCAGGCGGTCAGCCGGCGCAGCGATTCCAGCACCTCGCGGc                                                                    <  7:185067/136‑1 (MQ=255)
                                                  tGGTTGAGGATGCTGGGCAGCAAGGTGATGTGGCTGGCATCGCTGACACACAGGCGAAAACGCCGCTGCGCTGTGGTGGCGATGAACTGCGGCTCCCAGGCGGTCAGCCGGCGCAGCGATTCCAGCACCTCGCGGc                                                                    <  8:182632/136‑1 (MQ=255)
                                                  tGGTTGAGGATGCGGGGCAGCAAGGTGATGTGGCTGGCATCGCTGACACACAGGCGAAAACGCCGCTGCGCTGTGGTGGCGATGAACTGCGGCTCCCAGGCGGTCAGCCGGCGCAGCGATTCCAGCACCTCGCGGc                                                                    <  8:12760/136‑1 (MQ=255)
                                                   ggTTGAGGATGCTGGGCAGCAAGGTGATGTGGCTGGCATCGCTGACACACAGGCGAAAACGCCGCTGCGCTGTGGTGGCGATGAACTGCGGCTCCCAGGCGGTCAGCCGGCGCAGCGATTCCAGCACCTCGCGGc                                                                    <  2:230092/135‑1 (MQ=255)
                                                    gTTGAGGATGCTGGGCAGCAAGGTGATGTGGCTGGCATCGCTGACACACAGGCGAAAACGCCGCTGCGCTGTGGTGGCGATGAACTGCGGCTCCCAGGCGGTCAGCCGgcgc                                                                                          <  2:59303/112‑1 (MQ=255)
                                                    gTTGAGGATGCTGGGCAGCAAGGTGATGTGGCTGGCATCGCTGACACACAGGCGAAAACGCCGCTGCGCTGTGGTGGCGATGAACTGCGGCTCCCAGGCGGTCAGCCGgcgc                                                                                          >  1:59303/1‑112 (MQ=255)
                                                    gTTGAGGATGCTGGGCAGCAAGGTGATGTGGCTGGCATCGCTGACACACAGGCGAAAACGCCGCTGCGCTGTGGTGGCGATGAACTGCGGCTCCCAGGCGGTCAGCCGGCGCAGCGATTCCAGCACCTCGCGGCAc                                                                  >  8:28079/1‑136 (MQ=255)
                                                    gTTGAGGATGCTGGGCAGCAAGGTGATGTGGCTGGCATCGCTGACACACAGGCCAAAACGCCGCTGCGCTGTGGTGGCGATGAACTGCGGCTCCCAGGCGGTCAGCCGGCGCAGCGATTCCAGCACCTCGCGGCAc                                                                  <  2:189475/136‑1 (MQ=255)
                                                      tGAGGATGCTGGGCAGCAAGGTGATGTGGCTGGCATCGCTGACACACAGGCGAAAACGCCGCTGCGCTGTGGTGGCGATGAACTGCGGCTCCCAGGCGGTCAGCCGGCGCAGCGATTCCAGCACCTCGCGGCACgg                                                                >  8:6337/1‑136 (MQ=255)
                                                                gggCAGCAAGGTGATGTGGCTGGCATCGCTGACACACAGGCGAAAACGCCGCTGCGCTGTGGTGGCGATGAACTGCGGCTCCCAGGCGGTCAGCCGGCGCAGCGATTCCAGCACCTCGCGGCACGGCCCTATCAGc                                                      >  3:122560/1‑136 (MQ=255)
                                                                gggCAGCAAGGTGATGTGGCTGGCATCGCTGACACACAGGCGAAAACGCCGCTGCGCTGTGGTGGCGATGAACTGCGGCTCCCAGGCGGTCAGCCGGCGCAGCGATTCCAGCACCTCGCGGCACGGCCCTATCAGc                                                      >  6:37121/1‑136 (MQ=255)
                                                                  gcagcaAGGTGATGTGGCTGGCATCGCTGACACACAGGCGAAAACGCCGCTGCGCTGTGGTGGCGATGAACTGCGGCTCCCAGGCGGTCAGCCGGCGCAGCGATTCCAGCACCTCGCGGCACGGCCCTATCAGCt                                                     >  8:205572/1‑135 (MQ=255)
                                                                  gcagcaAGGTGATGTGGCTGGCATCGCTGACACACAGGCGAAAACGCCGCTGCGCTGTGGTGGCGATGAACTGCGGCTCCCAGGCGGTCAGCCGGCGCAGCGATTCCAGCACCTCGCGGCACGGCCCTATCAGCt                                                     >  7:238504/1‑135 (MQ=255)
                                                                  gcagcaAGGTGATGTGGCTGGCATCGCTGACACACAGGCGAAAACGCCGCTGCGCTGTGGTGGCGATGAACTGCGGCTCCCAGGCGGTCAGCCGGCGCAGCGATTCCAGCACCTCGCGGCACGGCCCTATCAGCt                                                     >  5:148376/1‑135 (MQ=255)
                                                                    agcaAGGTGATGTGGCTGGCATCGCTGACACACAGGCGAAAACGCCGCTGCGCTGTGGTGGCGATGAACTGCGGCTCCCAGGCGGTCAGCCGGCGCAGCGATTCCAGCACCTCGCGGCACGGCCCTATCAGCTGAt                                                  >  6:92295/1‑136 (MQ=255)
                                                                        tggTGAGGGGGCGGGCATCGCTGACCCACAGGCGAACACGCCGCTGCGCTGGGGTGGCGATGAACTGCGGCTCCCAGGCGGTCAGCCGGCGCAGCGATTCCAGCACCTCGCGGCACGGCCCTATCAGCTGATCGgc                                              <  4:128789/135‑1 (MQ=255)
                                                                        aGGTTATGTGTCTGGCATCGCTGACCCACAGCCGAAAACGCCGATGCGCTGTGGGGGCGATGAACTGCGGCTCCCAGGCGGTCAGCCGGCGCAGCGATTCCAGCACCTCGCGGCACGGCCCTATCAGCTGATCGgc                                              <  4:133214/136‑1 (MQ=255)
                                                                        aGGTGATGTGGCTGGCATCGCTGACACACAGGCGAAAACGCCGCTGCGCTGTGGTGGCGATGAACTGCGGCTCCCAGGCGGTCAGCCGGCGCAGCGATTCCAGCACCTCGCGGCACGGCCCTATCAGCTGATCGgc                                              <  7:77050/136‑1 (MQ=255)
                                                                        aGGTGATGTGGCTGGCATCGCTGACACACAGGCGAAAACGCCGCTGCGCTGTGGTGGCGATGAACTGCGGCTCCCAGGCGGTCAGCCGGCGCAGCGATTCCAGCACCTCGCGGCACGGCCCTATCAGCTGATCGgc                                              <  5:209058/136‑1 (MQ=255)
                                                                        aGGTGATGTGGCTGGCATCGCTGACACACAGGCGAAAACGCCGCTGCGCTGTGGTGGCGATGAACTGCGGCTCCCAGGCGGTCAGCCGGCGCAGCGATTCCAGCACCTCGCGGCACGGCCCTATCAGCTGATCGgc                                              <  2:59736/136‑1 (MQ=255)
                                                                        aGGTGATGTGGCTGGCATCGCTGACACACAGGCGAAAACGCCGCTGCGCTGTGGCGGCGATGAACTGCGGCTCCCAGGCGGTCAGCCGGCGCAGCGATTCCAGCACCTCGCGGCACGGCCCTATCAGCTGATCGgc                                              <  6:63524/136‑1 (MQ=255)
                                                                        aGGGGATGTGGCTGGCATCGCTGACACACAGGCGAAAACGCCGCTGCGCTGTGGTGGCGATGAACTGCGGCTCCCAGGCGGTCAGCCGGCGCAGCGATTCCAGCACCTCGCGGCACGGCCCTATCAGCTGATCGgc                                              <  4:240960/136‑1 (MQ=255)
                                                                        aGGGGATGTGGCTGGCATCGCTGACACACAGGCGAAAACGCCGCTGCGCTGTGGTGGCGATGAACTGCGGCTCCCAGGCGGTCAGCCGGCGCAGCGATTCCAGCACCTCGCGGCACGGCCCTATCAGCTGATCGgc                                              <  6:18520/136‑1 (MQ=255)
                                                                        aGGGGATGTGGCTGGCATCGCGGACCCACAGGCCAAAACGCCGCGGCGCTGTGGTGGCGATGAACTGCGGCTCCCAGGCGGTCAGCCGGCGCAGCGATTCCAGCACCTCGCGGCACGGCCCTATCAGCTGATCGgc                                              <  2:144486/136‑1 (MQ=255)
                                                                        aGGGGATGGGGCGGGCAGCGCTGACACACAGGCGAACACGCCGCTGCGCTGTGGTGGCGATGAACTGCGGCTCCCAGGCGGTCAGCCGGCGCAGCGATTCCAGCACCTCGCGGCACGGCCCTATCAGCTGATCGgc                                              <  6:24839/136‑1 (MQ=255)
                                                                         ggTGCTGTGGCTGGCATCGCTGACACACAGGCGAAAACGCCGCTGCGCTGTGGTGGCGATGAACTGCGGCTCCCAGGCGGTCAGCCGGCGCAGCGATTCCAGCACCTCGCGGCACGGCCCTATCAGCTGATCGgcg                                             <  6:222152/136‑1 (MQ=255)
                                                                          gTGATGTGGCTGGCATCGCTGACACACAGGCGAAAACGCCGCTGCGCTGTGGCGGCGATGAACTGCGGCTCCCAGGCGGTCAGCCGGCGCAGCGATTCCAGCACCTCGCGGCACGGCCCTATCAGCTGATCGgcg                                             <  6:241064/135‑1 (MQ=255)
                                                                                 ggctggcATCGCTGACACACAGGCGAAAACGCCGCTGCGCTGTGGTGGCGATGAACTGCGGCTCCCAGGCGGTCAGCCGGCGCAGCGATTCCAGCACCTCGCGGCACGGCCCTATCAGCTGATCTGCGCGCggcg                                      >  6:41126/1‑133 (MQ=255)
                                                                                     ggcATCGCTGACACACAGGCGAAAACGCCGCTGCGCTGTGGCGGCGATGAACTGCGGCTCCCAGGCGGTCAGCCGGCGCAGCGATTCCAGCACCTCGCGGCACGGCCCTATCAGCTGATCGGCGCGCGGGGTAggg                                 >  2:65608/1‑136 (MQ=255)
                                                                                       cATCGCTGACACACAGGCGAAAACGCCGCTGCGCTGTGGTGGCGATGAACTGCGGCTCCCAGGCGGTCAGCCGGCGCAGCGTTTCCAGCACCTCGCGGCACGGCCCTATCAGCTGATCGGCGCGCGGGGTAGGGGc                               <  8:86413/136‑1 (MQ=255)
                                                                                       cATCGATGACACACAGCCGAAAAGGCCGCTGCGCTGTGGTGGCGATGAACTGCGGCTCCCAGGCGGTCAGCCGGCGCAGCGATTCCAGCACCTCGCGGCACGGCCCTATCAGCTGATCGGCGCGCGGGGTAGGGGc                               <  5:37659/136‑1 (MQ=255)
                                                                                        aTCGCTGACACACAGGCGAAAACGCCGCTGCGCTGTGGTGGCGATGAACTGCGGCTCCCAGGCGGTCAGCCGGCGCAGCGATTCCAGCACCTCGCGGCACGGCCCTATCAGCTGATCGGCGCGCGGGGTAGGGGc                               >  7:135334/1‑135 (MQ=255)
                                                                                        aTCGCTGACACACAGGCGAAAACGCCGCTGCGCTGTGGTGGCGATGAACTGCGGCTCCCAGGCGGTCAGCCGGCGCAGCGATTCCAGCACCTCGCGGCACGGCCCTATCAGCTGATCGGCGCGCGGGGTAGGGGc                               >  7:9425/1‑135 (MQ=255)
                                                                                               acacacAGGCGAAAACGCCGCTGCGCTGTGGTGGCGATGAACTGCGGCTCCCAGGCGGTCAGCCGGCGCAGCGATTCCAGCACCTCGCGGCACGGCCCTATCAGCTGATCGGCGCGCGGGGTAGGGGCCATgccgc                       >  7:19196/1‑136 (MQ=255)
                                                                                               acacacAGGCGAAAACGCCGCTGCGCTGTGGTGGCGATGAACTGCGGCTCCCAGGCGGTCAGCCGGCGCAGCGATTCCAGCACCTCGCGGCACGGCCCTATCAGCTGATCGGCGCGCGGGGTAGGGGCCATgccgc                       >  4:56616/1‑136 (MQ=255)
                                                                                               acacacAGGCGAAAACGCCGCTGCGCTGTGGTGGCGATGAACTGCGGCTCCCAGGCGGTCAGCCGGCGCAGCGATTCCAGCACCTCGCGGCACGGCCCTATCAGCTGATCGGCGCGCGGGGTAGGGGCCATgccgc                       >  3:14131/1‑136 (MQ=255)
                                                                                                       gCGAAAACGCCGCTGCGCTGTGGTGGCGATGAACTGCGGCTCCCAGGCGGTCAGCCGGCGCAGCGATTCCAGCACCTCGCGGCACGGCCCTATCAGCTGATCGGCGCGCGGGGTAGGGGCCATGCCGCCTGGCGTg               <  1:183429/136‑1 (MQ=255)
                                                                                                       gCGAAAACGCCGCTGCGCTGTGGGGGCGATGAACTGCGGCTCCCAGGCGGTCAGCCGGCGCAGCGATTCCAGCACCTCGCGGCACGGCCCTATCAGCTGATCGGCGCGCGGGGTAGGGGCCATGCCGCCTGGCGt                <  6:188912/135‑1 (MQ=255)
                                                                                                        cGAAAACGCCGCTGCGCTGTGGTGGCGATGAACTGCGGCTCCCAGGCGGTCAGCCGGCGCAGCGATTCCAGCACCTCGCGGCACGGCCCTATCAGCTGATCGGCGCGCGGGGTAGGGGCCATGCCGCCTGGCGTgc              >  6:152618/1‑136 (MQ=255)
                                                                                                        cGAAAACGCCGCTGCGCTGTGGCGGCGATGAACTGCGGCTCCCAGGCGGTCAGCCGGCGCAGCGATTCCAGCACCTCGCGGCACGGCCCTATCAGCTGATCGGCGCGCGGGGTAGGGGCCATGCCGCCTGGCGTgc              >  1:131381/1‑136 (MQ=255)
                                                                                                                   ctgcgctgTGGTGGCGATGAACTGCGGCTCCCAGGCGGTCAGCCGGCGCAGCGATTCCAGCACCTCGCGGCACGGCCCTATCAGCTGATCGGCGCGCGGGGTAGGGGCCATGCCGCCTGGCGTGCGCACGAACAGc   <  6:215959/136‑1 (MQ=255)
                                                                                                                   ctgcgctgTGGTGGCGATGAACTGCGGCTCCCAGGCGGTCAGCCGGCGCAGCGATTCCAGCACCTCGCGGCACGGCCCTATCAGCTGATCGGCGCGCGGGGTAGGGGCCATGCCGCCTGGCGTGCGCACGAACAGc   <  1:230364/136‑1 (MQ=255)
                                                                                                                   ctgcgctgTGGGGGCGATGAACTGCGGCTCCCACGCGGGCAGCCGGCGCAGCGATTCCAGCACCTCGCGGCACGGCCCTATCAGCTGATCGGCGCGCGGGGTAGGGGCCATGCCGCCTGGCGTGCGCACGAACAGc   <  3:171113/136‑1 (MQ=255)
                                                                                                                    tgcgctgTGGTGGCGATGAACTGCGGCTCCCAGGCGGTCAGCCGGCGCAGCGATTCCAGCACCTCGCGGCACGGCCCTATCAGCTGATCGGCGCGCGGGGTAGGGGCCATGCCGCCTGGCGTGCGCACGAACAGCg  >  1:62360/1‑136 (MQ=255)
                                                                                                                              |                                                                                                                             
CGTCGATCCGCGCCGCCTCCAGGCGGATGCCCGGCGCATGGCTGCGCAGATGGTTGAGGATGCTGGGCAGCAAGGTGATGTGGCTGGCATCGCTGACACACAGGCGAAAACGCCGCTGCGCTGTGGCGGCGATGAACTGCGGCTCCCAGGCGGTCAGCCGGCGCAGCGATTCCAGCACCTCGCGGCACGGCCCTATCAGCTGATCGGCGCGCGGGGTAGGGGCCATGCCGCCTGGCGTGCGCACGAACAGCG  >  ACB143/3387038‑3387289

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