Predicted mutation
evidence seq id position mutation freq annotation gene description
RA ACB143 4,509,674 G→A 9.5% S481S (TCC→TCT)  clpA ← ATP‑dependent Clp protease ATP‑binding subunitClpA

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*ACB1434,509,6740GA 9.5% 128.0 / 3.2 42S481S (TCC→TCT) clpAATP‑dependent Clp protease ATP‑binding subunitClpA
Reads supporting (aligned to +/- strand):  ref base G (16/22);  new base A (2/2);  total (18/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 = 1.00e+00

CATGTCGAAACGCACCAACTCCACACCCAGCGCCTTGGCCAATTGCCGCGCCGCTTCGGTCTTGCCGACGCCGGTAGGGCCGGCGAACAGGAACGAACCCACCGGCTTGTCTGGCGACTTGAGGCCGGCACGGGACAACTTGATCGCGGTGGCCAGCGAGTCGATTGCCTGGTCCTGGCCAAACACGGTCAGTTTCAGGTCGCGCTCAAGGTTGCGCAGCAGCTCCTTGTCGGAACTGGTGACGTGCTTCGGCGGAATTCGC  >  ACB143/4509542‑4509803
                                                                                                                                    |                                                                                                                                 
cATGTCGAAACGCACCAACTCCACACCCAGCGCCTTGGCCAATTGCCGCGCCGCTTCGGTCTTGCCGACGCCGGTAGGGCCGGCGAACAGGAACGAACCCACCGGCTTGTCTGGCGACTTGAGGCCGGCACGGGa                                                                                                                                 >  3:173446/1‑135 (MQ=255)
 aTGTCGAAACGCACCAACTCCACACCCAGCGCCTTGGCCAATTGCCGCGCCGCTTCGGTCTTGCCGACGCCGGTAGGGCCGGCGAACAGGAACGAACCCACCGGCTTGTCTGGCGACTTGAGGCCGGCACGGGAc                                                                                                                                >  1:119614/1‑135 (MQ=255)
     cGAAACGCACCAACTCCACACCCAGCGCCTTGGCCAATTGCCGCGCCGCTTCGGTCTTGCCGACGCCGGTAGGGCCCGCGAACAGGAACGAACCCACCGGCTTGTCTGGCGACTAGAGGCCGGCACGGGACAACtt                                                                                                                           >  2:121218/1‑136 (MQ=255)
      gAAACGCACCAACTCCACACCCAGCGCCTTGGCCAATTGCCGCGCCGCTTCGGTCTTGCCGACGCCGGTAGGGCCGGCGAACAGGAACGAACCCACCGGCTTGTCTGGCGACTTGAGGCCGGCACGGGACAACtt                                                                                                                           >  7:184485/1‑135 (MQ=255)
       aaaCGCACCAACTCCACACCCAGCGCCTTGGCCAATTGCCGCGCCGCTTCGGTCTTGCCGACGCCGGTAGGGCCGGCGAACAGGAACGAACCCACCGGCTTGTCTGGCGACTTGAGGCCGGCACGGGACAACTTGa                                                                                                                         >  2:5611/1‑136 (MQ=255)
                         cccAGCGCCTTGGCCAATTGCCGCGCCGCTTCGGTCTTGCCGACGCCGGTAGGGCCGGCGAACAGGAACGAACCCACCGGCTTGTCTGGCGACTTGAGGCCGGCACGGGACAACTTGATCGCGGTGGCCAGCGAGt                                                                                                       <  6:289539/136‑1 (MQ=255)
                         cccAGCGCCTTGGCCAATTGCCGCGCCGCTTCGGTCTTGCCGACGCCGGTAGGGCCGGCGAACAGGAACGAACCCACCGGCTTGTCTGGCGACTTGAGGCCGGCACGGGACAACTTGATCGCGGTGGCCAGCGAGt                                                                                                       <  4:134303/136‑1 (MQ=255)
                         cccAGCGCCTTGGCCAATTGCCGCGCCGCTTCGGTCTTGCCGACGCCGGTAGGGCCGGCGAACAGGAACGAACCCACCGGCTTGTCTGGCGACTTGAGGCCGGCACGGGACAACTTGATCGCGGTGGCCAGCGAGt                                                                                                       <  3:257212/136‑1 (MQ=255)
                              cgcCTTGGCCAATTGCCGCGCCGCTTCGGTCTTGCCGACGCCGGTAGGGCCGGCGAACAGGAACGAACCCACCGGCTTGTCTGGCGACTTGAGGCCGGCACGGGACAACTTGATCGCGGTGGCCAGCGAGTCGAtt                                                                                                  <  4:120255/136‑1 (MQ=255)
                              cgcCTTGGCCAATTGCCGCGCCGCTTCGGTCTTGCCGACGCCGGTAGGGCCGGCGAACAGGAACGAACCCACCGGCCTGTCTGGCGACTTGAGGCCGGCACGGGACAACTTGATCGCGGTGGCCAGCGAGTCGAtt                                                                                                  <  3:18448/136‑1 (MQ=255)
                               gcCTTGGCCAATTGCCGCGCCGCTTCGGTCTTGCCGACGCCGGTAGGGCCGGCGCACAGGAACGAACCCACCGGCTTGTCTGGCGACTTGAGGCCGGCACGGGACAACTTGATCGCGGTGGCCAGCGAGTCGATTg                                                                                                 >  8:122643/1‑136 (MQ=255)
                               gcCTTGGCCAATTGCCGCGCCGCTTCGGTCTTGCCGACGCCGGTAGGGCCGGCGAACAGGAACGAACCCACCGGCTTGTCTGGCGACTTGAGGCCGGCACGGGACAACTTGATCGCGGTGGCCAGCGAGTCGATTg                                                                                                 >  2:203355/1‑136 (MQ=255)
                                ccTTGGCCAATTGCCGCGCCGCTTCGGTCTTGCCGACGCCGGTAGGGCCGGCGAACAGGAACGAACCCACCGGCTTGTCTGGCGACTTGAGGCCGGCACGGGACAACTTGATCGCGGTGGCCAGCGAGTCGATTGc                                                                                                >  2:46066/1‑136 (MQ=255)
                                ccTTGGCCAATTGCCGCGCCGCTTCGGTCTTGCCGACGCCGGTAGGGCCGGCGAACAGGAACGAACCCACCGGCTTGTCTGGCGACTTGAGGCCGGCACGGGACAACTTGATCGCGGTGGCCAGCGAGTCGATTGc                                                                                                >  3:173002/1‑136 (MQ=255)
                                ccTTGGCCAATTGCCGCGCCGCTTCGGTCTTGCCGACGCCGGTAGGGCCGGCGAACAGGAACGAACCCACCGGCTTGTCTGGCGACTTGAGGCCGGCACGAGACAACTTGATCGCGGTGGCCAGCGAGTCGATTGc                                                                                                >  2:125858/1‑136 (MQ=255)
                                ccTTGGCCAATTGCCGCGCCGCTTCGGTCTTGCCGACGCCGGTAGGGCCGGCGAACAGGAACGAACCCACCGGCTTGTCTGGCGACTTGAGGCCGGCACGAGACAACTTGATCGCGGTGGCCAGCGAGTCGATTGc                                                                                                >  2:125861/1‑136 (MQ=255)
                                       cAATTGCCGCGCCGCTTCGGTCTTGCCGACGCCGGTAGGGCCGGCGCACAGGAACGAACCCACCGGCTTGTCTGGCGACTTGAGGCCGGCACGGGACAACGTGATCGCGGTGGCCAGCGAGTCGATTTCCTGGtgc                                                                                         >  3:123510/1‑134 (MQ=255)
                                       cAATTGCCGCGCCGCTTCGGTCTTGCCGACGCCGGTAGGGCCGGCGAACAGGAACGAACCCACCGGCTTGTCTGGCGACTTGAGGCCGGCACGGGACAACTTGATCGCGGTGGCCAGCGAGTCGATTGCCTGGTcc                                                                                         >  1:246559/1‑136 (MQ=255)
                                                  ccgcTTCGGTCTTGCCGACGCCGGTAGGGCCGGCGAACAGGAACGAACCCACCGGCTTGTCTGGCGACTTGAGGCCGGCACGGGACAACTTGATCGCGGTGGCCAGCGAGTCGATTGCCTGGTCCTGGCCCAacac                                                                              >  5:221405/1‑136 (MQ=255)
                                                  ccgcTTCGGTCTTGCCGACGCCGGTAGGGCCGGCGAACAGGAACGAACCCACCGGCTTGTCTGGCGACTTGAGGCCGGCACGGGACAACTTGATCGCGGTGGCCAGCGAGTCGATTGCCTGGTCCTGGCCAAacac                                                                              <  3:170025/136‑1 (MQ=255)
                                                  ccgcTTCGGTCTTGCCGACGCCGGTAGGGCCGGCGAACAGGAACGAACCCACCGGCTTGTCTGGCGACTTGAGGCCGGCACGGGACAACTTGATCGCGGTGGCCAGCGAGTCGATTGCCTGGTCCTGGCCAAacac                                                                              <  6:73915/136‑1 (MQ=255)
                                                  ccgcTTCGGTCTTGCCGACGCCGGTAGGGCCGGCGAACAGGAACGAACCCACCGGCTTGTCTGGCGACGTGAGGCCGGCACGGGACAACTTGATCGCGGTGGCCAGCGAGTCGATTGCCTGGTGCTGGCCaaaacc                                                                              >  2:200500/1‑133 (MQ=255)
                                                  ccgcTTCGGTCTTGCCGACGCCGGTAGGGCCGGCGAACAGGAACGAACCCACCGGCTTGGCTGGCGACTTGAGGCCGGCACGGGACAACTTGATCGCGGTGGCCAGCGAGTCGATTGCCTGGTCCTGGCCAAacac                                                                              <  6:230710/136‑1 (MQ=255)
                                                  ccgcGTCGGTCTTGCCGACGCCGGTAGGGCCGGCGAACAGGAACGAACCCACCGGCTTGTCTGGCGACTTGAGGCCGGCACGGGACAACTTGATCGCGGTGGCCAGCGAGTCGATTGCCTGGTCCTGGCCAAacac                                                                              <  2:258063/136‑1 (MQ=255)
                                                               gCCGACGCCGGTAGGGCCGGCGAACAGGAACGAACCCACCGGCTTGTCTGGCGACTTGAGGCCGGCACGGGACAACTTGATCGCGGTGGCCAGCGAGTCGATTGCCTGGTCCTGGCCAAACACGGTCAGTTTCa                                                                   <  5:3813/134‑1 (MQ=255)
                                                               gCCGACGCCGGTAGGGCCGGCGAACAGGAACGAACCCACCGGCTTGTCTGGCGACTTGAGGCCGGCACGGGACAACTTGATCGCGGTGGCCAGCGAGTCGATTGCCTGGTCCTGGCCAAACACGGTCAGTTTCa                                                                   <  5:69277/134‑1 (MQ=255)
                                                               gCCGACGCCGGTAGGGCCGGCGAACAGGAACGAACCCACCGGCTTGTCTGGCGACTTGAGGCCGGCACGGGACAACTTGATCGCGGTGGCCAGCGAGTCGATTGCCTGGTCCTGGCCAAACACGGTCAGTTTCAg                                                                  <  1:2831/135‑1 (MQ=255)
                                                                      ccGGTAGGGCCGGCGAACAGGAACGAACCCACCGGCTTGTCTGGCGACTTGAGGCCGGCACGGGACAACTTGATCGCGGTGGCCAGCGAGTCGATTGCCTGGTCCTGGCCAAACACGGTCAGTTTCAGGTCGCGCt                                                          <  4:173002/136‑1 (MQ=255)
                                                                            gggCCGGCGAACAGGAACGAACCCACCGGCTTGTCTGGCGACTTGAGGCCGGCACGGGACAACTTGATCGCGGTGGCCAGCGAGTCGATTGCCTGGTCCTGGCCAAACACGGTCAGTTTCAGGTCGCGCTCAAgg                                                     <  5:199593/135‑1 (MQ=255)
                                                                            gggCCGGCGAACAGGAACGAACCCACCGGCTTGTCTGGCGACTTGAGGCCGGCACGGGACAACTTGATCGCGGTGGCCAGCGAGTCGATTGCCTGGTCCTGGCCAAACACGGTCAGTTTCAGGTCGCGCTCAAgg                                                     <  3:201987/135‑1 (MQ=255)
                                                                                           aacgaacCGACCGGCTTGTCTGGCGACTTGAGGCCGGCACGGGACAACTTGATCGCGGTGGCCAGCGAGTCGATTGCCTGGTCCTGGCCAAACACGGTCAGTTTCAGGTCGCGCTCAAGGTTGCGCAGCAGCTCCt                                     <  5:202549/136‑1 (MQ=255)
                                                                                           aacgaacCCACCGGCTTGTCTGGCGACTTGAGGCCGGCACGGGACAACTTGATCGCGGTGGCCAGCGAGTCGATTGCCTGGTCCTGGCCAAACACGGTCAGTTTCAGGTCGCGCTCAAGGTTGCGCAGCAGCTCCt                                     <  7:207579/136‑1 (MQ=255)
                                                                                           aacgaacCCACCGGCTTGTCTGGCGACTTGAGGCCGGCACGGGACAACTTGATCGCGGTGGCCAGCGAGTCGATTGCCTGGTCCTGGCCAAACACGGTCAGTTTCAGGTCGCGCTCAAGGTTGCGCAGCAGCTCCt                                     <  4:13387/136‑1 (MQ=255)
                                                                                           aacgaacCCACCGGCTTGTCTGGCGACTTGAGGCCGGCACGGGACAACTTGATCGCGGTGGCCAGCGAGTCGATTGCCTGGTCCTGGCCAAACACGGTCAGTTTCAGGTCGCGCTCAAGGTTGCGCAGCAGCTCCt                                     >  5:206996/1‑136 (MQ=255)
                                                                                              gaacCCACCGGCTTGTCTGGCGACTTGAGGCCGGCACGGGACAACTTGATCGCGGTGGCCAGCGAGTCGATTGCCTGGTCCTGGCCAAACACGGTCAGTTTCAGGTCGCGCTCAAGGTTGCGCAGCAGCTCCTTGt                                  <  1:270772/136‑1 (MQ=255)
                                                                                              gaacCCACCGGCTTGTCTGGCGACTTGAGGCCGGCACGAGACAACTTGATCGCGGTGGCCAGCGAGTCGATTGCCTGGTCCTGGCCAAACACGGTCAGTTTCAGGTCGCGCTCAAGGTTGCGCAGCAGCTCCTTGt                                  <  1:125858/136‑1 (MQ=255)
                                                                                              gaacCCACCGGCTTGTCTGGCGACTTGAGGCCGGCACGAGACAACTTGATCGCGGTGGCCAGCGAGTCGATTGCCTGGTCCTGGCCAAACACGGTCAGTTTCAGGTCGCGCTCAAGGTTGCGCAGCAGCTCCTTGt                                  <  1:125861/136‑1 (MQ=255)
                                                                                                 cccACCGGCTTGTCTGGCGACTTGAGGCCGGCACGGGACAACTTGATCGCGGTGGCCAGCGAGTCGATTGCCTGGTCCTGGCCAAACACGGTCAGTTTCAGGTCGCGCTCAAGGTTGCGCAGCAGCTCCTTGTCgg                               <  6:4574/136‑1 (MQ=255)
                                                                                                       ggCTTGTCTGGCGACTTGAGGCCGGCACGGGACAACTTGATCGCGGTGGCCAGCGAGTCGATTGCCTGGTCCTGGCCAAACACGGTCAGTTTCAGGTCGCGCTCAAGGTTGCGCAGCAGCTCCTTGTCGGAACTgg                         >  7:133848/1‑136 (MQ=255)
                                                                                                         cTTGTCTGGCGACTTGAGGCCGGCACGGGACAACTTGATCGCGGTGGCCAGCGAGTCGATTGCCTGGTCCTGGCCAAACACGGTCAGTTTCAGGTCGCGCTCAAGGTTGCGCAGCAGCTCCTTGTCGGAACTGGTg                       <  3:238420/136‑1 (MQ=255)
                                                                                                                     cTTGAGGCCGGCACGGGACAACTTGATCGCGGTGGCCAGCGAGTCGATTGCCTGGTCCTGGCCAAACACGGTCAGTTTCAGGTCGCGCTCAAGGTTGCGCAGCAGCTCCTTGTCGGAACTGGTGACGTGCTTcggc           >  2:91991/1‑136 (MQ=255)
                                                                                                                              ggCACGGGACAACTTGATCGCGGTGGCCAGCGAGTCGATTGCCTGGTCCTGGCCAAACACGGTCAGTTTCAGGTCGCGCTCAAGGTTGCGCAGCAGCTCCTTGTCGGAACTGGTGACGTGCTTCGGCGGAATTcgc  <  7:60884/136‑1 (MQ=255)
                                                                                                                                    |                                                                                                                                 
CATGTCGAAACGCACCAACTCCACACCCAGCGCCTTGGCCAATTGCCGCGCCGCTTCGGTCTTGCCGACGCCGGTAGGGCCGGCGAACAGGAACGAACCCACCGGCTTGTCTGGCGACTTGAGGCCGGCACGGGACAACTTGATCGCGGTGGCCAGCGAGTCGATTGCCTGGTCCTGGCCAAACACGGTCAGTTTCAGGTCGCGCTCAAGGTTGCGCAGCAGCTCCTTGTCGGAACTGGTGACGTGCTTCGGCGGAATTCGC  >  ACB143/4509542‑4509803

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