Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,382,318 A→T 11.9% I125F (ATC→TTC)  smpB → SsrA‑binding protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,382,3180AT11.9% 122.5 / 8.8 42I125F (ATC→TTC) smpBSsrA‑binding protein
Reads supporting (aligned to +/- strand):  ref base A (19/18);  new base T (2/3);  total (21/21)
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.93e-01

CACGTCATCGCCGACCCTACACGCACCCGCAAGCTGCTGCTGAACAAGCGCGAGCTCGAGCGCGTGGAAGCGGCCGTGGCGCAAAAAGGCTACACCTGTGTGGCTTTGGCGCTGTACTGGAGCAAGCACCTGATCAAGTGCGAAATCGCACTGGGCAAGGGCAAGAAGGAATTCGACAAGCGCGACACCATGCGCGAACGCGATTCCAACCGCGAGCTGCAGCGGGCTGTGCGGAACAAGGGCAAGGAAGAGTAACTACTCTTCAACCTGTTCCGGCCTCTTCGCGGGCAA  >  NC_002947/5382174‑5382464
                                                                                                                                                |                                                                                                                                                  
cACGTCATCGCCGACCCTACACGCACCCGCAAGCTGCTGCTGAACAAGCGCGAGCTCGAGCGCGTGGAAGCGGCCGTGGCGCAAAAAGGCTACACCTGTGTGGCTTTGGCGCTGTACTGGAGCAAGCACCTGATCAAGTGCGAAATCGCa                                                                                                                                               <  2:463863/150‑1 (MQ=255)
  cGTCATCGCCGACCCTACACGCACCCGCAAGCTGCTGCTGAACAAGCGCGAGCTCGAGCGCGTGGAAGCGGCCGTGGCGCAAAAAGGCTACACCTGTGTGGCTTTGGCGCTGTACTGGAGCAAGCACCTGATCAAGTGCGATTTCGCACt                                                                                                                                             <  2:248203/150‑1 (MQ=255)
    tCATCGCCGACCCTACACGCACCCGCAAGCTGCTGCTGAACAAGCGCGAGCTCGAGCGCGTGGAAGCGGCCGTGGCGCAAAAAGGCTACACCTGTGTGGCTTTGGCGCTGTACTGGAGCAAGCACCTGATCAAGTGCGAAATCGCACTgg                                                                                                                                           <  1:426507/150‑1 (MQ=255)
     taTCGCCGACCCTACACGCACCCGCAAGCTGCTGCTGAACAAGCGCGAGCTCGAGCGCGTGGAAGCGGCCGTGGCGCAAAAAGGCTACACCTGTGTGGCTTTGGCGCTGTACTGGAGCAAGCACCTGATCAAGTGCGATTTCGCACttga                                                                                                                                          <  2:257078/149‑4 (MQ=255)
      aTCGCCGACCCTACACGCACCCGCAAGCTGCTGCTGAACAAGCGCGAGCTCGAGCGCGTGGAAGCGGCCGTGGCGCAAAAAGGCTACACCTGTGTGGCTTTGGCGCTGTACTGGAGCAAGCACCTGATCAAGTGCGAAATCGCACTgggc                                                                                                                                         >  1:725935/1‑150 (MQ=255)
         gCCGACCCTACACGCACCCGCAAGCTGCTGCTGAACAAGCGCGAGCTCGAGCGCGTGGAAGCGGCCGTGGCGCAAAAAGGCTACACCTGTGTGGCTTTGGCGCTGTACTGGAGCAAGCACCTGATCAAGTGCGAAATCGCACTgggcaag                                                                                                                                      <  1:772973/150‑1 (MQ=255)
               ccTACACGCACCCGCAAGCTGCTGCTGAACAAGCGCGAGCTCGAGCGCGTGGAAGCGGCCGTGGCGCAAAAAGGCTACACCTGTGTGGCTTTGGCGCTGTACTGGAGCAAGCACCTGATCAAGTGCGAAATCGCACTGggcaagggcaag                                                                                                                                >  1:417565/1‑150 (MQ=255)
                  acacGCACCCGCAAGCTGCTGCTGAACAAGCGCGAGCTCGAGCGCGTGGAAGCGGCCGTGGCGCAAAAAGGCTACACCTGTGTGGCTTTGGCGCTGTACTGGAGCAAGCACCTGATCAAGTGCGAAATCGCACTGGGCAAGGGCaagaag                                                                                                                             <  1:29127/150‑1 (MQ=255)
                  acacGCACCCGCAAGCTGCTGCTGAACAAGCGCGAGCTCGAGCGCGTGGAAGCGGCCGTGGCGCAAAAAGGCTACACCTGTGTGGCTTTGGCGCTGTACTGGAGCAAGCACCTGATCAAGTGCGAAATCGCACTGGGCAAGGGCaagaag                                                                                                                             <  1:419065/150‑1 (MQ=255)
                        aCCCGCAAGCTGCTGCTGAACAAGCGCGAGCTCGAGCGCGTGGAAGCGGCCGTGGCGCAAAAAGGCTACACCTGTGTGGCTTTGGCGCTGTACTGGAGCAAGCACCTGATCAAGTGCGAAATCGCACTGGGCAAGGGCAAGAAGGAATTc                                                                                                                       <  2:991034/150‑1 (MQ=255)
                          ccGCAAGCTGCTGCTGTACAAGCGCGAGCTCGAGCGCGTGGAAGCGGCCGTGGCGCAAAAAGGCTACACCTGTGTGGCTTTGGCGCTGTACTGGAGCAAGCACCTGATCAAGTGCGAAATCGCACTGGGCAAGGGCAAGAAGGAATTCgc                                                                                                                     <  2:827594/150‑2 (MQ=255)
                               aGCTGCTGCTGAACAAGCGCGAGCTCGAGCGCGTGGAAGCGGCCGTGGCGCAAAAAGGCTACACCTGTGTGGCTTTGGCGCTGTACTGGAGCAAGCACCTGATCAAGTGCGAAATCGCACTGGGCAAGGGCAAGAAGGAATTCGACAAgc                                                                                                                <  2:108044/150‑1 (MQ=255)
                                gctgctgctGAACAAGCGCGAGCTCGAGCGCGTGGAAGCGGCCGTGGCGCAAAAAGGCTACACCTGTGTGGCTTTGGCGCTGTACTGGAGCAAGCACCTGATCAAGTGCGAAATCGCACTGGGCAAGGGCAAGAAGGAATTCGACAAgcg                                                                                                               >  1:136049/1‑150 (MQ=255)
                                   gcggctgAACAAGGGCGAGCTCGAGCGCGTGGAAGCGGCCGTGGCGCAAAAAGGCTACACCTGTGTGGCTTTGGCGCTGTACTGGAGCAAGCACCTGATCAAGTGCGAAATCGCACTGGGCAAGGGCAAGAAGGAATTCGACAAGCGCGa                                                                                                            <  2:66903/147‑1 (MQ=255)
                                      gctgAACAAGCGCGAGCTCGAGCGCGTGGAAGCGGCCGTGGCGCAAAAAGGCTACACCTGTGTGGCTTTGGCGCTGTACTGGAGCAAGCACCTGATCAAGTGCGAAATCGCACTGGGCAAGGGCAAGAAGGAATTCGACAAGCGCGacac                                                                                                         <  2:417565/150‑1 (MQ=255)
                                      gctgAACAAGCGCGAGCTCGAGCGCGTGGAAGCGGCCGTGGCGCAAAAAGGCTACACCTGTGTGGCTTTGGCGCTGTACTGGAGCAAGCACCTGATCAAGTGCGAAATCGCACTGGGCAAGGGCAAGAAGGAATTCGACAAGCGCGacac                                                                                                         >  2:1194011/1‑150 (MQ=255)
                                                                            tGGCGCAAAAAGGCTACACCTGTGTGGCTTTGGCGCTGTACTGGAGCAAGAACCTGATCAAGTGCGAAATCGCACTGGGCAAGGGCAAGAAGGAATTCGACAAGCGCGACACCATGCGCGAACGCGATTCCAACCGCGAGCTGCAGCggg                                                                   <  2:136049/150‑1 (MQ=255)
                                                                              gcgcAAAAAGGCTACACCTGTGTGGCTTTGGCGCTGTACTGGAGCAAGCACCTGATCAAGTGCGAAATCGCACTGGGCAAGGGCAAGAAGGAATTCGACAAGCGCGACACCATGCGCGAACGCGATTCCAACCGCGAGCTGCAGCGGGCt                                                                 >  2:114004/1‑150 (MQ=255)
                                                                                gcAAAAAGGCTACACCTGTGTGGCTTTGGCGCTGTACTGGAGCAAGCACCTGATCAAGTGCGAAATCGCACTGGGCAAGGGCAAGAAGGAATTCGACAAGCGCGACACCATGCGCGAACGCGATTCCAACCGCGAGCTGCAGCGGGCtgt                                                               >  1:771739/1‑150 (MQ=255)
                                                                                  aaaaaGGCTACACCTGTGTGGCTTTGGCGCTGTACTGGAGCAAGCACCTGATCAAGTGCGAAATCGCACTGGGCAAGGGCAAGAAGGAATTCGACAAGCGCGACACCATGCGCGAACGCGATTCCAACCGCGAGCTGCAGCGGGCTGTGc                                                             <  2:725935/150‑1 (MQ=255)
                                                                                    aaaGGCTACACCTGTGTGGCTTTGGCGCTGTACTGGAGCAAGCACCTGATCAAGTGCGAAATCGCACTGGGCAAGGGCAAGAAGGAATTCGACAAGCGCGACACCATGCGCGAACGCGATTCCAACCGCGAGCTGCAGCGGGCTGTGCgg                                                           <  1:876768/150‑1 (MQ=255)
                                                                                             acCTGTGTGGCTTTGGCGCTGTACTGGAGCAAGCACCTGATCAAGTGCGAAATCGCACTGGGCAAGGGCAAGAAGGAATTCGACAAGCGCGACACCATGCGCGAACGCGATTCCAACCGCGAGCTGCAGCGGGCTGTGCGGAACAAGGGc                                                  >  1:1122749/1‑150 (MQ=255)
                                                                                                    tgGCTTTGGCGCTGTACTGGAGCAAGCACCTGATCAAGTGCGAAATCGCACTGGGCAAGGGCAAGAAGGAATTCGACAAGCGCGACACCATGCGCGAACGCGATTCCAACCGCGAGCTg                                                                          <  1:259059/119‑1 (MQ=255)
                                                                                                    tgGCTTTGGCGCTGTACTGGAGCAAGCACCTGATCAAGTGCGAAATCGCACTGGGCAAGGGCAAGAAGGAATTCGACAAGCGCGACACCATGCGCGAACGCGATTCCAACCGCGAGCTg                                                                          >  2:259059/1‑119 (MQ=255)
                                                                                                    tgGCTTTGGCGCTGTACTGGAGCAAGCACCTGATCAAGTGCGAAATCGCACTGGGCAAGGGCAAGAAGGAATTCGACAAGCGCGACACCATGCGCGAACGCGATTCCAACCGCGAGCTGCAGCGGGCTGTGCGGAACAAGGGCAAGGAag                                           >  2:1141881/1‑150 (MQ=255)
                                                                                                      gCTTTGGCGCTGTACTGGAGCAAGCACCTGATCAAGTGCGAAATCGCACTGGGCAAGGGCAAGAAGGAATTCGACAAGCGCGACACCATGCGCGAACGCGATTCCAACCGCGAGCTGCAGCGGGCTGTGCGGAACAAGGGCAAGGAagag                                         >  1:60691/1‑150 (MQ=255)
                                                                                                       cTTTGGCGCTGTACTGGAGCAAGCACCTGATCAAGTGCGAAATCGCACTGGGCAAGGGCAAGAAGGAATTCGACAAGCGCGACACCATGCGCGAACGCGATTCCAACCGCGAGCTGCAGCGGGCTGTGCGGAACAAGGGCAAGGAAGAGt                                        >  1:176877/1‑150 (MQ=255)
                                                                                                        tttGGCGCTGTACTGGAGCAAGCACCTGATCAAGTGCGAAATCGCACTGGGCAAGGGCAAGAAGGAATTCGACAAGCGCGACACCATGCGCGAACGCGATTCCAACCGCGAGCTGCAGCGGGCTGTGCGGAACAAGGGCAAgg                                              <  2:514877/143‑1 (MQ=255)
                                                                                                        tttGGCGCTGTACTGGAGCAAGCACCTGATCAAGTGCGAAATCGCACTGGGCAAGGGCAAGAAGGAATTCGACAAGCGCGACACCATGCGCGAACGCGATTCCAACCGCGAGCTGCAGCGGGCTGTGCGGAACAAGGGCAAgg                                              >  1:514877/1‑143 (MQ=255)
                                                                                                             cgcTGTACTGGAGCAAGCACCTGATCAAGTGCGAAATCGCACTGGGCAAGGGCAAGAAGGAATTCGACAAGCGCGACACCATGCGCGAACGCGATTCCAACCGCGAGCTGCAGCGGGCTGTGCGGAACAAGGGCAAGGAAGAGTAactac                                  >  2:100330/1‑150 (MQ=255)
                                                                                                                tGTACTGGAGCAAGCACCTGATCAAGTGCGAAATCGCACTGGGCAAGGGCAAGAAGGAATTCGACAAGCGCGACACCAt                                                                                                      <  1:397696/79‑1 (MQ=255)
                                                                                                                tGTACTGGAGCAAGCACCTGATCAAGTGCGAAATCGCACTGGGCAAGGGCAAGAAGGAATTCGACAAGCGCGACACCAt                                                                                                      >  2:397696/1‑79 (MQ=255)
                                                                                                                  tACTGGAGCAAGCACCTGATCAAGTGCGAAATCGCACTGGGCAAGGGCAAGAAGGAATTCGACAAGCGCGACACCATGCGCGAACGCGATTCCAACCGCGAGCTGCAGCGGGCTGTGCGGAACAAGGGCaa                                                >  2:66498/1‑131 (MQ=255)
                                                                                                                  tACTGGAGCAAGCACCTGATCAAGTGCGAAATCGCACTGGGCAAGGGCAAGAAGGAATTCGACAAGCGCGACACCATGCGCGAACGCGATTCCAACCGCGAGCTGCAGCGGGCTGTGCGGAACAAGGGCaa                                                <  1:66498/131‑1 (MQ=255)
                                                                                                                            aGCACCTGATCAAGTGCGAAATCGCACTGGGCAAGGGCAAGAAGGAATTCGACAAGCGCGACACCATGCGCGAACGCGATTCCAACCGCGAGCTGCAGCGGGCTGTGCGGAACAAGGGCAAGGAAGAGTAACTACTCTTCAACCTGTTCC                   >  2:248332‑M1/1‑150 (MQ=255)
                                                                                                                                ccTGATCAAGTGCGAAATCGCACTGGGCAAGGGCAAGAAGGAATTCGACAAGCGCGACACCATGCGCGAACGCGATTCCAACCGCGAGCTGCAGCGGGCTGTGCGGAACAAGGGCAAGGAAGAGt                                        >  2:433185/1‑125 (MQ=255)
                                                                                                                                ccTGATCAAGTGCGAAATCGCACTGGGCAAGGGCAAGAAGGAATTCGACAAGCGCGACACCATGCGCGAACGCGATTCCAACCGCGAGCTGCAGCGGGCTGTGCGGAACAAGGGCAAGGAAGAGt                                        <  1:433185/125‑1 (MQ=255)
                                                                                                                                 cttgcccagtgcgatTTCGCACTGGGCAAGGGCAAGAAGGAATTCGACAAGCGCGACACCATGCGCGAACGCGATTCCAACCGCGAGCTGCAGCGGGCTGTGCGGAACAAGGGCAAGGAAGAGTAACTACTCTTCAACCTGTTCCGGCCT              >  1:539876‑M1/8‑150 (MQ=255)
                                                                                                                                      ccaGTGCGATTTCGCACTGGGCAAGGGCAAGAAGGAATTCGACAAGCGCGACACCATGCGCGAACGCGATTCCAACCGCGAGCTGCAGCGGGCTGTGCGGAACAAg                                                     >  1:840761/3‑106 (MQ=255)
                                                                                                                                      ccaGTGCGATTTCGCACTGGGCAAGGGCAAGAAGGAATTCGACAAGCGCGACACCATGCGCGAACGCGATTCCAACCGCGAGCTGCAGCGGGCTGTGCGGAACAAg                                                     <  2:840761/104‑1 (MQ=255)
                                                                                                                                       aaGTGCGAAATCGCACTGGGCAAGGGCAAGAATGAATTCGACAAGCGCGACACCATGCGCGAACGCGATTCCAACCGCGAGCTGCAGCGGGCTGTGCGGAACAAGGGCAAGGAAGAGTAACTACTCTTCAACCTGTTCCGGCCTCTTCGC        >  2:1117841‑M1/1‑150 (MQ=255)
                                                                                                                                             gAAATCGCACTGGGCAAGGGCAAGAAGGAATTCGACAAGCGCGACACCATGCGCGAACGCGATTCCAACAGCGAGCTGCAGCGGGCTGTGCGGAACAAGGGCAAGGAAGAGTAACTACTCTTCAACCTGTTCCGGCCTCTTCGCGGGCAA  >  1:857746‑M1/1‑150 (MQ=255)
                                                                                                                                                |                                                                                                                                                  
CACGTCATCGCCGACCCTACACGCACCCGCAAGCTGCTGCTGAACAAGCGCGAGCTCGAGCGCGTGGAAGCGGCCGTGGCGCAAAAAGGCTACACCTGTGTGGCTTTGGCGCTGTACTGGAGCAAGCACCTGATCAAGTGCGAAATCGCACTGGGCAAGGGCAAGAAGGAATTCGACAAGCGCGACACCATGCGCGAACGCGATTCCAACCGCGAGCTGCAGCGGGCTGTGCGGAACAAGGGCAAGGAAGAGTAACTACTCTTCAACCTGTTCCGGCCTCTTCGCGGGCAA  >  NC_002947/5382174‑5382464

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