Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,382,317 A→T 16.1% E124D (GAA→GAT smpB → SsrA‑binding protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,382,3170AT16.1% 100.6 / 13.8 44E124D (GAA→GATsmpBSsrA‑binding protein
Reads supporting (aligned to +/- strand):  ref base A (18/18);  new base T (3/4);  total (22/22)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.22e-01

CCCACGTCATCGCCGACCCTACACGCACCCGCAAGCTGCTGCTGAACAAGCGCGAGCTCGAGCGCGTGGAAGCGGCCGTGGCGCAAAAAGGCTACACCTGTGTGGCTTTGGCGCTGTACTGGAGCAAGCACCTGATCAAGTGCGAAATCGCACTGGGCAAGGGCAAGAAGGAATTCGACAAGCGCGACACCATGCGCGAACGCGATTCCAACCGCGAGCTGCAGCGGGCTGTGCGGAACAAGGGCAAGGAAGAGTAACTACTCTTCAACCTGTTCCGGCCTCTTCGCGG  >  NC_002947/5382172‑5382460
                                                                                                                                                 |                                                                                                                                               
cccACGTCATCGCCGACCCTACACGCACCCGCAAGCTGCTGCTGAACAAGCGCGAGCTCGAGCGCGTGGAAGCGGCCGTGGCGCAAAAAGGCTACACCTGTGTGGCTTTGGCGCTGTACTGGAGCAAGCACCTGATCAAGTGCGAAATCg                                                                                                                                             >  2:710792/1‑150 (MQ=255)
 ccACGTCATCGCCGACCCTACACGCACCCGCAAGCTGCTGCTGAACAAGCGCGAGCTCGAGCGCGTGGAAGCGGCCGTGGCACAAAAAGGCTACACCTGTGTGGCTTTGGCGCTGTACTGGAGCAAGCACCTGATCAAGTGCGAAATCGc                                                                                                                                            >  2:200111/1‑150 (MQ=255)
   aCGTCATCGCCGACCCTACACGCACCCGCAAGCTGCTGCTGAACAAGCGCGAGCTCGAGCGCGTGGAAGCGGCCGTGGCGCAAAAAGGCTACACCTGTGTGGCTTTGGCGCTGTACTGGAGCAAGCACCTGATCAAGTGCGATTTCGCAc                                                                                                                                          <  1:491095/150‑1 (MQ=255)
     gTCATCGCCGACCCTACACGCACCCGCAAGCTGCTGCTGAACAAGCGCGAGCTCGAGCGCGTGGAAGCGGCCGTGGCGCAAAAAGGCTACACCTGTGTGGCTTTGGCGCTGTACTGGAGCAAGCACCTGATCAAGTGCGATTTCGCACtt                                                                                                                                        <  1:508401/150‑2 (MQ=255)
                        gCACCCGCAAGCTGCTGCTGAACAAGCGCGAGCTTGAGCGCGTGGAAGCGGCCGTGGCGCAAAAAGGCTACACCTGTGTGGCTTTGGCGCTGTACTGGAGCAAGCACCTGAGCAAGTCCGAAATCGCACTGGGCAAGGGCAAGAAGGAAt                                                                                                                     >  1:516392/1‑150 (MQ=255)
                        gCACCCGCAAGCTGCTGCTGAACAAGCGCGAGCTCGAGCGCGTGGAAGCGGCCGTGGCGCAAAAAGGCTACACCTGTGTGGCTTTGGCGCTGTACTGGAGCAAGCACCAGATAAGGTGAGAAATCGCACTGGGCAAGGGCAAGAAGGAAt                                                                                                                     >  2:649194/1‑150 (MQ=255)
                               cAAGCTGCTGCTGAACAAGCGCGAGCTCGAGCGCGTGGAAGCGGCCGTGGCGCAAAAAGGCTACACCTGTGTGGCTTTGGCGCTGTACTGGAGCAAGCACCTGATCAAGTGCGAAATCGCACTGGGCAAGGGCAAGAAGGAATTCGACaa                                                                                                              <  1:710792/150‑1 (MQ=255)
                                    tgctgctgAACAAGCGCGAGCTCGAGCGGATGGAAGCGGCCGTGACGCAAAAAGGCTACAGCTGTGTGGGTTTGGCGCTGTCCTGGAGCAAGCACCTGATCAAGTGCGACATCGCACTGGGCAAGGGCAAGAGGGAATTCGACAAGcgcg                                                                                                         >  1:417386/1‑150 (MQ=255)
                                          tgAACAAGCGCGAGCTCGAGCGCGTGGAAGCGGCCGTGGAGCAAAAAGGCTACACCTGTGTGGCTTTGGCGCTGTACTGGAGCAAGCACCTGATCAAGTGCGAAATCGCACTGGGCAAGGGCAAGAAGGAATTCGACAAGCGCGACACCa                                                                                                   <  1:649194/150‑1 (MQ=255)
                                           gAACAAGCGCGAGCTCGAGCGCGTGGAAGAGGCCGTGGCGCAAAAAGGCTACACCTGAGTGGCTTTGGCGCTGTACTGGAGCAAGCACCTGATCAAGTGCGAAATCGCACTGGGCAAGGGCAAGAAGGAATTCGACAAGCGCGACACCAt                                                                                                  >  1:534178/1‑150 (MQ=255)
                                           gAACAAGCGCGAGCTCGAGCGCGTGCCCGGTGCCGTGGCGCAAAAAGGCTACACCTGTGTGGCTTTGGCGCTGTACTGGAGCAAGCACCTGATCAAGTGCGAAATCGCACTGGGCAAGGGCAAGAAGGAATTCGACAAGCGCGACACCAt                                                                                                  >  1:450684/1‑150 (MQ=255)
                                                 gcgcGAGCTCGAGCGCGTGGAAGCGGCCGTGGCGCAAAAAGGCTACACCTGTGTGGCTTTGGCGCTGTACTGGAGCAAGCACCTGATCAAGTGCGAAATCGCACTGGGCAAGGGCAAGAAGGAATTCGACAAGCGCGACACCATGCGCGa                                                                                            >  1:668218/1‑150 (MQ=255)
                                                         tCGAGCGCGTGGAAGCGGCCGTGGCGCAAAAAGGCTACACCTGTGTGGCTTTGGCGCTGTACTGGAGCAAGCACCTGATCAAGTGCGAAATCGCACTGGGCAAGGGCAAGAAGGAATTCGACAAGCGCGACACCATGCGCGAACGCGAtt                                                                                    <  2:326481/150‑1 (MQ=255)
                                                         tCGAGCGCGTGGAAGCGGCCGTGGCGCAAAAAGGCTACACCTGTGTGGCTTTGGCGCTGTACTGGAGCAAGCACCTGATCAAGTGCGAAATCGCACTGGGCAAGGGCAAGAAGGAATTCGACAAGCGCGACACCATGCGCGAACGCGAtt                                                                                    >  2:468525/1‑150 (MQ=255)
                                                          cGAGCGCGTGCAAGCGGCCGTGGCGCAAAAAGGCTACACCTGTGTGGCTTTGGCGCTGTACTGGAGCAAGCACCTGATCAAGTGCGAAATCGCACTGGGCAAGGGCAAGAAGGAATTCGACAAGCGCGACACCATGCGCGAACGCGATTc                                                                                   <  1:210314/150‑1 (MQ=255)
                                                             gcgcGTGGAAGCGGCCGTGGCGCAAAAAGGCTACACCTGTGTGGCTTTGGCGCTGTACTGGAGCAAGCACCTGATCAAGTGCGAAATCGCACTGGGCAAGGGCAAGAAGGAATTCGACAAGCGCGACACCATGCGCGAACGCGATTCCaa                                                                                >  1:229399/1‑150 (MQ=255)
                                                                    gAAGCGGCCGTGGCGCAAAAAGGCTACACCTGTGTGGCTTTGGCGCTGTACTGGAGCAAGCACCTGATCAAGTGCGAAATCGCACTGGGCAAGGGCAAGAAGGAATTCGACAAGCGCGACACCATGCGCGAACGCGATTCCAACCGCGAg                                                                         <  2:575091/150‑1 (MQ=255)
                                                                      aGCGGCCGTGGCGCAAAAAGGCTACACCTGTGTGGCTTTGGCGCTGTACTGGAGCAAGCACCTGATCAAGTGCGAAATCGCACTGGGCAAGGGCaagaag                                                                                                                         >  2:44985/1‑100 (MQ=255)
                                                                      aGCGGCCGTGGCGCAAAAAGGCTACACCTGTGTGGCTTTGGCGCTGTACTGGAGCAAGCACCTGATCAAGTGCGAAATCGCACTGGGCAAGGGCaagaag                                                                                                                         <  1:44985/100‑1 (MQ=255)
                                                                        cGGCCGTGGCGCAAAAAGGCTACACCTGTGTGGCTTTGGCGCTGTACTGGAGCAAGCACCTGATCAAGTGCGAAATCGCACTGGGCAAGGGCaagaag                                                                                                                         <  1:7238/98‑1 (MQ=255)
                                                                        cGGCCGTGGCGCAAAAAGGCTACACCTGTGTGGCTTTGGCGCTGTACTGGAGCAAGCACCTGATCAAGTGCGAAATCGCACTGGGCAAGGGCaagaag                                                                                                                         >  2:7238/1‑98 (MQ=255)
                                                                          gCCGTGGCGCAAAAAGGCTACACCTGTGTGGCTTTGGCGCTGTACTGGAGCAAGCACCTGATCAAGTGCGAAATCGCACTgggc                                                                                                                                     <  1:387490/84‑1 (MQ=255)
                                                                          gCCGTGGCGCAAAAAGGCTACACCTGTGTGGCTTTGGCGCTGTACTGGAGCAAGCACCTGATCAAGTGCGAAATCGCACTgggc                                                                                                                                     >  2:387490/1‑84 (MQ=255)
                                                                                gcgcAAAAAGGCTACACCTGTGTGGCTTTGGCGCTGTACTGGAGCAAGCACCTGATCAAGTGCGAAATCGCACTGGGCAAGGGCAAGAAGGAATTCGACAAGCGCGACACCATGCGCGCACGAGATTCCAACCGCGAGCTGCAGCGGGCt                                                             >  1:641587/1‑150 (MQ=255)
                                                                                           cTACACCTGTGTGGCTTTGGCGCTGTACTGGAGCAAGCACCTGATCAAGTGCGAAATCGCACTGGGCAAGGGCAAGAAGGAATTCGACAAGCGCGACACCATGCGCGAACGCGATTCCAACCGCGAGCTGCAGCGGGCTGTGCGGAACaa                                                  >  2:444541/1‑150 (MQ=255)
                                                                                                 cTGTGTGGCTTTGGCGCTGTACTGGAGCAAGCACCTGATCAAGTGCGAAATCGCACTGGGCAAGGGCAAGAAGGa                                                                                                                       >  1:607178/1‑75 (MQ=255)
                                                                                                 cTGTGTGGCTTTGGCGCTGTACTGGAGCAAGCACCTGATCAAGTGCGAAATCGCACTGGGCAAGGGCAAGAAGGa                                                                                                                       <  2:607178/75‑1 (MQ=255)
                                                                                                      tgGCTTTGGCGCTGTACTGGAGCAAGCACCTGATCAAGTGCGAAATCGCACTGGGCAAGGGCAAGAAGGAATTCGACAAGCGCGACACCATGCGCGAACGCGATTCCAACCGCGAGCTGCAGCGGGc                                                              >  2:108239/1‑127 (MQ=255)
                                                                                                      tgGCTTTGGCGCTGTACTGGAGCAAGCACCTGATCAAGTGCGAAATCGCACTGGGCAAGGGCAAGAAGGAATTCGACAAGCGCGACACCATGCGCGAACGCGATTCCAACCGCGAGCTGCAGCGGGc                                                              <  1:108239/127‑1 (MQ=255)
                                                                                                       ggCTTTGGCGCTGTACTGGAGCAAGCACCTGATCAAGTGCGAAATCGCACTGGGCAAGGGCAAGAAGGAATTCGACAAGCGCGACACCATGCGCGAACGCGATTCCAACCGCGAGCTGCAGCGGGCTGTGCGGAACAAGGGCAAGGAaga                                      <  1:468525/150‑1 (MQ=255)
                                                                                                       ggCTTTGGCGCTGTACTGGAGCAAGCACCTGATCAAGTGCGAAATCGCACTGGGCAAGGGCAAGAAGGAATTCGACAAGCGCGACACCATGCGCGAACGCGATTCCAACCGCGAGCTGCAGCGGGCTGTGCGGAACAAGGGCAAGGAaga                                      <  2:668218/150‑1 (MQ=255)
                                                                                                         cTTTGGCGCTGTACTGGAGCAAGCACCTGATCAAGTGCGAAATCGCACTGGGCAAGGGCAAGAAGGAATTCGACAAGCGCGACACCATGCGCGAACGCGATTCCAACCGCGAGCTGCAGCGGGCTGTGCGGAACAAGGGCAAGGAAGAGt                                    <  2:740230/150‑1 (MQ=255)
                                                                                                                 cTGTAATGGAGCAAGCACCTGATCAAGTGCGAAATCGCACTGGGCAAGGGCAAGAAGGAATTCGACAAGCGCGACACCATGCGCGAACGCGATTCCAACCGCGAGCTGCAGCGGGCTGTGCGGAACAAGGGCAAGGAAGAGTAACTActc                            <  1:200111/150‑1 (MQ=255)
                                                                                                                              cggcccTGATCAAGTGCCAAACCGAACTGGGTAAGGGCAAGCAGGAATTCGACAAGCGCGACACCATGCGTCAACGAGATTCCAACCGCGATGTGCAGCGGGCTGTGCGGAAAAAGGGCAAGGAAGAGTAACTTCTCTTCAACCTGtgac               <  2:516392/146‑4 (MQ=255)
                                                                                                                                 aCCTGATCAAGTGCGAAATCGCACTGGGCAAGGGCAAGAAGGAATTCGACAAGCGCGACACCATGCGCGAACGCGATTCCaa                                                                                >  2:325159/1‑82 (MQ=255)
                                                                                                                                 aCCTGATCAAGTGCGAAATCGCACTGGGCAAGGGCAAGAAGGAATTCGACAAGCGCGACACCATGCGCGAACGCGATTCCaa                                                                                <  1:325159/82‑1 (MQ=255)
                                                                                                                                   cTGATCAAGTGCGAAATCGCACTGGGCAAGGGCAAGAAGGAATTCGACAAGCGCGACACCATGCGCGAACGCGATTCCAACCGCGAGCTGCAGCGGGCTGTGCGGAACAAGGGCAAGGaa                                        <  2:623565/120‑1 (MQ=255)
                                                                                                                                   cTGATCAAGTGCGAAATCGCACTGGGCAAGGGCAAGAAGGAATTCGACAAGCGCGACACCATGCGCGAACGCGATTCCAACCGCGAGCTGCAGCGGGCTGTGCGGAACAAGGGCAAGGaa                                        >  1:623565/1‑120 (MQ=255)
                                                                                                                                       tCAAGTGCGAAATCGCACTGGGCAAGGGCAAGAAGGAATTCGAAGAGCCGTGGACCAGGCGCGAACGCGATTCCAACCGCGAGCTGCAGCGGGCTGTGCGGAACAAGGGCAAGGAAGAGTAACTACTCTTCAACCTGTTCCGGCCTCTTc      <  2:450684/150‑1 (MQ=255)
                                                                                                                                        ccaGTGCGATTTCGCACTTGGCAAGGGCAAGAAGGAATTCGACAAGCGCGACACCATGCGCGAACGCGATTCCAACCGCGAGCTGCAGCGGGCTGTGCGGAACAAGGGCaa                                            <  2:198508/109‑1 (MQ=255)
                                                                                                                                        ccaGTGCGATTTCGCACTTGGCAAGGGCAAGAAGGAATTCGACAAGCGCGACACCATGCGCGAACGCGATTCCAACCGCGAGCTGCAGCGGGCTGTGCGGAACAAGGGCaa                                            >  1:198508/3‑111 (MQ=255)
                                                                                                                                        ccaGTGCGATTTCGCACTGGGCAAGGGCAAGAAGGAATTCGACAAGCGCGACACCATGCGCGAACGCGATTCCAACCGCGAGCTGCAGCGGGCTGTGCGGAACAAGGGCAAGGAAGAGTAACTACTCTTCAACCTGTTCCGGCCTCTTcg     >  2:235343/3‑150 (MQ=255)
                                                                                                                                           gTGCGATTTCGCACTGGGCAAGGGCAAGAAGGAATTCGACAAGCGCGACACCATGCGCGAACGCGATTCCAACCGCGAGCTGCAGCGGGCTGTGCGGAACAAGGGCAAGGAAGAGTAACTACTCTTCAACCTGTTCCGGCCTCTTCGCgg  >  2:153099/1‑150 (MQ=255)
                                                                                                                                           gTGCGATTTCGCACTGGGCAAGGGCAAGAAGGAATTCGACAAGCGCGACACCATGCGCGAACGCGATTCCAACCGCGAGCTGCAGCGGGCTGTGCGGAACAAGGGCAAGGAAGAGTAACTACTCTTCAACCTGTTCCGGCCTCTTCGCgg  <  1:153099/150‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                               
CCCACGTCATCGCCGACCCTACACGCACCCGCAAGCTGCTGCTGAACAAGCGCGAGCTCGAGCGCGTGGAAGCGGCCGTGGCGCAAAAAGGCTACACCTGTGTGGCTTTGGCGCTGTACTGGAGCAAGCACCTGATCAAGTGCGAAATCGCACTGGGCAAGGGCAAGAAGGAATTCGACAAGCGCGACACCATGCGCGAACGCGATTCCAACCGCGAGCTGCAGCGGGCTGTGCGGAACAAGGGCAAGGAAGAGTAACTACTCTTCAACCTGTTCCGGCCTCTTCGCGG  >  NC_002947/5382172‑5382460

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