Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP000730 1,776,114 G→A 29.8% P137L (CCA→CTA)  clpX ← S14 family endopeptidase ClpX

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0007301,776,1140GA29.8% 50.5 / 30.7 45P137L (CCA→CTA) clpXS14 family endopeptidase ClpX
Reads supporting (aligned to +/- strand):  ref base G (15/15);  new base A (9/4);  total (24/21)
Fisher's exact test for biased strand distribution p-value = 3.24e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.46e-01

TTTCATCTACATAAATAATACCTTTTTCGGCTTTATCAATGTCAAAGTCAGCTGCTTGAATTAATCTCAACAAGATATTTTCAACATCATCGCCTACATAACCAGCTTCAGTTAAACTTGTCGCATCTGCAATTGCAAATGGTACATTCAACGTCTTGGCTAAGGTTTGAGCTAATAATGTTTTACCACTACCTGTTGGCCCAATTAATGCAATGTTACTTTTTTGTAATTCAACATCATCTTCTTTTGGTCCTAATTGTTGAATACGCTTATAGTGG  >  CP000730/1775974‑1776251
                                                                                                                                            |                                                                                                                                         
tttCATCTACATAAATAATACCTTTTTCGGCTTTATCAATGTCAAAGTCAGCTGCTTGAATTAATCTCAACAAGATATTTTCAACATCATCGCCTACCTAACCAGCTTCAGTTAAACTTGTCGCATCTGCCATTGCAACTg                                                                                                                                           >  2:465562/1‑141 (MQ=255)
      cTACATAAATAATACCTTTTTCGGCTTTATCAATGTCAAAGTCAGCTGCTTGAATTAATCTCAACAAGATATTTTCAACATCATCGCCTACATAACCAGCTTCAGTTAAACTTGTCGCATCTGCAATTGCAAATAGTACAt                                                                                                                                     <  1:76746/141‑1 (MQ=255)
        aCATAAATAATACCTTTTTCGGCTTTATCAATGTCAAAGTCAGCTGCTTGAATTAATCTCAACAAGATATTTTCAACATCATCGCCTACATAACCAGCTTCAGTTAAACTTGTCGCATCTGCAATTGCAAATAGTACATTc                                                                                                                                   <  1:560696/141‑1 (MQ=255)
        aCATAAATAATACCTTTTTCGGCTTTATCAATGTCAAAGTCAGCTGCTTGAATTAATCTCAACAAGATATTTTCAACATCATCGCCTACATAACCAGCTTCAGTTAAACTTGTCGCATCTGCAATTGCAAATAGTACATTc                                                                                                                                   >  2:399847/1‑141 (MQ=255)
                aatcCCTTTCTCGGCTCTACCAATGTCAAAGTCAGCTGCTTGAATTAATCCCAACAAGATATTTTCAACATCATCGCCTACATAACCATCTTCAGTTAAACTTGTCGCATCTGCAATTGCAAATGGTACATTCAACGTCtt                                                                                                                           <  2:337037/141‑1 (MQ=255)
                  taCCTTTTTCGGCTTTATCAATGTCAAAGTCAGCTGCTTGAATTAATCTCAACAAGATATTTTCAACATCATCGCCTACATAACCAGCTTCAGTTAAACTTGTCGCATCTGCAATTGCAAATGGTACATTCAACGTCTTgg                                                                                                                         <  1:29424/141‑1 (MQ=255)
                   aCCTTTTTCGGCTTTATCAATGTCAAAGTCAGCTGCTTGAATTAATCTCAACAAGATATTTTCAACATCATCGCCTACATAACCAGCATCAGTTAAACTTGTCGCATCTGCAATTGCCGATAGTACATTCAACGTCTTGGc                                                                                                                        >  2:231852/1‑141 (MQ=255)
                     cTTTTTCGGCTTTATCAATGTCAAAGTCAGCTGCTTGAATTAATCTCAACAAGATATTTTCAACATCATCGCCTACATAACCAGCTTCAGTTAAACTTGTCGCATCTGCAATTGCAAATGGTACATTCAACGTCTTGGCTa                                                                                                                      <  2:22041/141‑1 (MQ=255)
                       ttttCGGCTTTATCAATGTCAAAGTCAGCTGCTTGAATTAATCTCAACAAGATATTTTCAACATCATCGCCTACATAACCAGCTTCAGTTAAACTTGTCGCATCTGCAATTGCAAATGGTACATTCAACGTCTTGGCTAAg                                                                                                                    >  1:258416/1‑141 (MQ=255)
                       ttttCGGCTTTATCAATGTCAAAGTCAGCTGCTTGAATTAATCTCAACAAGATATTTTCAACATCATCGCCTACATAACCAGCTTCAGTTAAACTTGTCGCATCTGCAATTGCAAATGGTACATTCAACGTCTTGGCTAAg                                                                                                                    <  2:355790/141‑1 (MQ=255)
                          tCGGCTTTATCAATGTCAAAGTCAGCTGCTTGAATTAATCTCAACAAGATATTTTCAACATCATCGCCTACATAACCAGCTTCAGTTAAACTTGTCGCATCTGCAATTGCAAATGGTACATTCAACGTCTTGGCTAAGGtt                                                                                                                 <  1:65755/141‑1 (MQ=255)
                           cGGCTTTATCAATGTCAAAGTCAGCTGCTTGAATTAATCTCAACAAGATATTTTCAACATCATCGCCTACATAACCAGCTTCAGTTAAACTTGTCGCATCTGCAATTGCAAATGGTACATTCAACGTCTTGGCTAAGGttt                                                                                                                >  2:494901/1‑141 (MQ=255)
                            ttttttttttaTGTCAAAGTCAGCTGCTTGAATTAATCTCAACAAGATATTTTCAACATCATCGCCTACATAACCAGCTTCAGTTTAACTTGTCTCATCTGCAATTGCAAATGGTACATTCAACGTCTTGGCTTAGGtttt                                                                                                               >  2:196054/11‑140 (MQ=255)
                            ggCTTTATCAATGTCAAAGTCAGCTGCTTGAATTAATCTCAACAAGATATTTTCAACATCATCGCCTACATAACCAGCTTCAGTTAAACTTGTCGCATCTGCAATTGCAAATGGTACATTCAACGTCTTGGCTAAGGTTTg                                                                                                               >  2:436539/1‑141 (MQ=255)
                                    cAATGTCAAAGTCAGCTGCTTGAATTAATCTCAACAAGATATTTTCAACATCATCGCCTACATAACCAGCTTCAGTTAAACTTGTCGCATCTGCAATTGCAAATGGTACATTCAACGTCTTTGCTAAGGTTTGAGCtaata                                                                                                       <  1:320057/141‑1 (MQ=255)
                                    cAATGTCAAAGTCAGCTGCTTGAATTAATCTCAACAAGATATTTTCAACATCATCGCCTACATAACCAGCTTCAGTTAAACTTGTCGCATCTGCAATTGCAAATGGTACATTCAACGTCTTGGCTAAGGTTTGAGCtaata                                                                                                       <  1:319810/141‑1 (MQ=255)
                                        gTCAAAGTCAGCTGCTTGAATTAATCTCAACAAGATATTTTCAACATCATCGCCTACATAACCAGCTTCAGTTAAACTTGTCGCATCTGCAATTGCAAATGGTATATTCAACGTCTTGGCTAAGGTTTGAGCTAATAATGt                                                                                                   >  2:403466/1‑141 (MQ=255)
                                            aaGTCAGCTGCTTGAATTAATCTCAACAAGATATTTTCAACATCATCGCCTACATAACCAGCTTCAGTTAAACTTGTCGCATCTGCAATTGCAAATAGTACATTCAACGTCTTGGCTAAGGTTTGAGCTAATAATGTTTTa                                                                                               >  2:89800/1‑141 (MQ=255)
                                                cAGCTGCTTGAATTAATCTCAACAAGATATTTTCAACATCATCGCCTACATAACCAGCTTCAGTTAAACTTGTCGCATCTGCAATTGCAAATGGTACATTCAACGTCTTGGCTAAGGTTTGAGCTAATAATGTTTTACCAc                                                                                           >  2:89617/1‑141 (MQ=255)
                                                cAGCTGCTTGAATTAATCTCAACAAGATATTTTCAACATCATCGCCTACATAACCAGCTTCAGTTAAACTTGTCGCATCTGCAATTGCAAATGGTACATTCAACGTCTTGGCTAAGGTTTGAGCTAATAATGTTTTACCAc                                                                                           >  1:427046/1‑141 (MQ=255)
                                                 aGCTGCTTGAATTAATCTCAACAAGATATTTTCAACATCATCGCCTACATAACCAGCTTCAGTTAAACTTGTCGCATCTGCAATTGCAAATGGTACATTCAACGTCTTGGCTAAGGTTTGAGCTAATAATGTTTTACCACt                                                                                          >  2:75452/1‑141 (MQ=255)
                                                    tgctTGAATTAATCTCAACAAGATATTTTCAACATCATCGCCTACATAACCAGCTTCAGTTAAACTTGTCGCATCTGCAATTGCAAATGGTACATTCAACGTCTTGGCTAAGGTTTGAGCTAATAATGTTTTACCACTAcc                                                                                       >  2:524244/1‑141 (MQ=255)
                                                     gctTGAATTAATCTCAACAAGATATTTTCAACATCATCGCCTACATAACCAGCTTCAGTTAAACTTGTCGCATCTGCAATTGCAAATAGTACATTCAACGTCTTGGCTAAGGTTTGAGCTAATAATGTTTTACCACTACCt                                                                                      >  1:210728/1‑141 (MQ=255)
                                                        tGAATTAATCTCAACAAGATATTTTCAACATCATCGCCTACATAACCAGCTTCAGTTAAACTTGTCGCATCTGCAATTGCAAATGGTACATTCAACGTCTTGGCTAAGGTTTGAGCTAATAATGTTTTACCACTACCTGtt                                                                                   >  1:355889/1‑141 (MQ=255)
                                                                         gATATTTTCAACATCATCGCCTACATAACCAGCTTCAGTTAAACTTGTCGCATCTGCAATTGCAAATGGTACATTCAACGTCTTGGCTAAGGTTTGAGCTAATAATGTTTTACCACTACCTGTTGGCCCAATTAATGCAAt                                                                  >  1:58147/1‑141 (MQ=255)
                                                                         gATATTTTCAACATCATCGCCTACATAACCAGCTTCAGTTAAACTTGTCGCATCTGCAATTGCAAATGGTACATTCAACGTCTTGGCTAAGGTTTGAGCTAATAATGTTTTACCACTACCTGTTGGCCCAATTAATGCAAt                                                                  >  2:71302/1‑141 (MQ=255)
                                                                           tagtttCAACATCATCGCCTACATAACCAGCTTCAGTTAAACTTGGCGCATCTGCAATTGCAAATGGTACATTCAACGTCTTGGTTAAAGATTGAACTAATAACATTATCCCACCACCTGTTGGCTCATTTTGTGCAATGt                                                                <  2:51352/138‑1 (MQ=255)
                                                                                      tcatcGCCTACATAACCAGCTTCAGTTAAACTTGTCGCATCTGCAATTGCAAATAGTACATTCAACGTCTTGGCTAAGGTTTGAGCTAATAATGTTTTACCACTACCTGTTGGTCCAATTAATGCAATGTTACTTTTTTGt                                                     >  2:328856/1‑141 (MQ=255)
                                                                                       catcGCCTACATAACCAGCTTCAGTTAAACTTGTCGCATCTGCAATTGCAAATCGTATTATTAGCGTCTTGGCTAAGGTTTGAGCTAATAATGTTTTACCACTACCTGTTGGCCCAATTAATGCAATGTTACTTTTTTGTa                                                    <  1:415461/141‑1 (MQ=255)
                                                                                         tcGCCTACATAACCAGCTTCAGTTAAACTTGTCGCATCTGCAATTGCAAATGGTACATTCAACGTCTTGGCTAAGGTTTGAGCTAATAATGTTTTACCACTACCTGTTGGCCCAATTAATGCAATGTTACTTTTTTGTAAt                                                  <  2:408774/141‑1 (MQ=255)
                                                                                              tACATAACCAGCTTCAGTTAAACTTGTCGCATCTGCAATTGCAAATGGTACATTCAACGTCTTGGCTAAGGTTCGAGCTAATACTGTTTTACCACTACCTGTTGGCCCAATTAATGCAATGTTACTTTTTTCTAAttaaca                                             >  2:137261/1‑137 (MQ=255)
                                                                                                cATAACCAGCTTCAGTTAAACTTGTCGCATCTGCAATTGCAAATAGTACATTCAACGTCTTGGCTAAGGTTTGAGCTAATAATGTTTTACCACTACCTGTTGGCCCAATTAATGCAATGTTACTTTTTTGTAATTCAAcat                                           <  1:495633/141‑1 (MQ=255)
                                                                                                       aGCTTCAGTTAAACTTGTCGCATCTGCAATTGCAAATGGTACATTCAACGTCTTGGCTAAGGTTTGAGCTAATAATGTTTTACCACTACCTGTTGGCCCAATTAATGCAATGTTACTTTTTTGTAATTCAACATCAtcttc                                    <  1:494901/141‑1 (MQ=255)
                                                                                                        gCTTCATTTAACTTTGTCGCATCTGAAATTGCAAATGGTACATTCAACGTCTTGGCTAAGGTTTGAGCTAATAATGTTTTACCACTACCTGTTGGCCCAATTAATGCAATGTTACTTTTTTGTAATTCAACATCAtcttct                                   <  2:507199/141‑1 (MQ=255)
                                                                                                        gCTTCAGTTAAACTTGTCGCATCTGCAATTGCAAATGGTACATTCAACGTCTTGGCTAAGGTTTGAGCTAATAATGTTTTACCACTACCTGTTGGCCCAATTAATGCAATGTTACTTTTTTGTAATTCAACATCAtcttct                                   <  1:13969/141‑1 (MQ=255)
                                                                                                                 aaaCTTTTCGCATCTGCAATTGCAAATGGTACATTCAACGTCTGGGCTAAGGTTTGAGCTAATAATGTTTTACCACTACCTGTTGGCCCAATTAATGCAATGTTACTTTTTTGTAATTCAACATCATCTTCTTTTGGTCCt                          <  1:196054/141‑1 (MQ=255)
                                                                                                                 aaaCTTGTCGCATCTGCAATTGCAAATGGTACATTCAACGTCTTGGCTAAGGTTTGAGCTAATAATGTTTTACCACTACCTGTTGGCCCAATTAATGCAATGTTACTTTTTTGTAATTCAACATCATCTTCTTTTGGTCCt                          <  1:436539/141‑1 (MQ=255)
                                                                                                                 aaaCGTGTCGCATCTGCAATTGCAAATAGTACATTCAACGTCTTGGCTAAGGTTTGAGCTAATAATGTTTTACCACTACCTGTTGGCCCAATTAATGCAATGTTACTTTTTTGTAATTCAACATCATCTTCTTTTGGTCCt                          <  2:330427/141‑1 (MQ=255)
                                                                                                                       gTCGCATCTGCAATTGCAAATGGTACATTCAACGTCTTGGCTAAGGTTTGAGCTAATAATGTTTTACCACTACCTGTTGGCCCAATTAATGCAATGTTACTTTTTTGTAATTCAACATCATCTTCTTTTGGTCCTGAttgt                    >  2:199088/1‑141 (MQ=255)
                                                                                                                           cATCTGCAATTGCAAATAGTACATTCAACGTCTTGGCTAAGGTTTGAGCTAATAATGTTTTACCACTACCTGTTGGCCCAATTAATGCAATGTTACTTTTTTGTAATTCAACATCATCTTCTTTTGGTCCTAATTGTTGaa                >  1:197010/1‑141 (MQ=255)
                                                                                                                             tCTGCAATTGCAAATAGTACATTCAACGTCTTGGCTAAGGTTTGAGCTAATAATGTTTTACCACTACCTGTTGGCCCAATTAATGCAATGTTACTTTTTTGTAATTCAACATCATCTTCTTTTGGTCCTCATTGTTGAATa              >  2:63935/1‑141 (MQ=255)
                                                                                                                                 cAATTGCAAATGGTACATTCAACGTCTTGGCTAAGGTTTGACCTAATAATGTTTTACCACTACCTGTTGGCCCAATTACTGCAATGTTACTTTTTTGTAATTCAACATCATCTTCTTTTGGTCCTAATTGTTGAATACGCt          <  2:33715/141‑1 (MQ=255)
                                                                                                                                   aTTGCAAATAGTACATTCAACGTCTTGGCTAAGGTTTGAGCTAATAATGTTTTACCACTACCTGTTGGCCCAATTAATGCAATGTTACTTTTTTGTAATTCAACATCATCTTCTTTTGGTCCTAATTGTTGACTACGCTta        >  1:505373/1‑141 (MQ=255)
                                                                                                                                   aTTGCAAATAGTACATTCAACGTCTTGGCTAAGGTTTGAGCTAATAATGTTTTACCACTACCTGTTGGCCCAATTAATGCAATGTTACTTTTTTGTAATTCAACATCATCTTCTTTTGGTCCTAATTGTTGAATCCGCTta        >  2:403181/1‑141 (MQ=255)
                                                                                                                                       cAAATTGTACATTCAACGTCTTGGCTAAGGTTTGAGCTAATAATGTTTTACCACTACCTGTTGGCCCAATTAATGCAATGTTACTTTTTTGTAATTCAACATCATCTTCTTTTGGTCCTAATTGTTGAATACGCTTATAGt    <  1:359521/141‑1 (MQ=255)
                                                                                                                                         aaTGGTACATTCAACGTCTTGGCTAAGGTTTGAGCTAATAATGTTTTACCACTACCTGTTGGCCCAATTAATGCAATGTTACTTTTTTATCATTCAACCTCACCTTCTTTTGCTCCTAATTGTTGAATACGCTTATAGTga  >  2:14048/1‑140 (MQ=255)
                                                                                                                                            |                                                                                                                                         
TTTCATCTACATAAATAATACCTTTTTCGGCTTTATCAATGTCAAAGTCAGCTGCTTGAATTAATCTCAACAAGATATTTTCAACATCATCGCCTACATAACCAGCTTCAGTTAAACTTGTCGCATCTGCAATTGCAAATGGTACATTCAACGTCTTGGCTAAGGTTTGAGCTAATAATGTTTTACCACTACCTGTTGGCCCAATTAATGCAATGTTACTTTTTTGTAATTCAACATCATCTTCTTTTGGTCCTAATTGTTGAATACGCTTATAGTGG  >  CP000730/1775974‑1776251

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