Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,981,060 G→C 8.5% *577S (TGA→TCA)  rpsA → 30S ribosomal protein S1

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,981,0600GC 8.5% 137.7 / 3.3 47*577S (TGA→TCA) rpsA30S ribosomal protein S1
Reads supporting (aligned to +/- strand):  ref base G (24/19);  new base C (2/2);  total (26/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 = 8.12e-01

CCGCGTTATCAGCCTGTCGATCAAATCGAAAGACGACGCTGAAGAGCGCGAAGC‑TATCCAGAGCCTGAAAAACGCTCCGGAAGCTGCTGCCGACACCACCATGGCTGCGCTGCTGCGCGAAGCAATGGCCAAGCAGAACTGAGTTCTGTTTGATCGGTAAAAAGGGCGACCCTCGGGTCGCCCTTTTTTTTGCTCCGCGACCAGGTGGTGCGGCGGATGGCACTGCTGGAACACGCCAGCTTCAGGCTGGCGAAGCCCGTTGTAGGAGCGGGTTCACCCGC  >  NC_002947/1980920‑1981200
                                                                                                                                             |                                                                                                                                            
gttctttATCAGCCTGTCTATCATATCGAATGACGTCGCTGAAGAGCGCGAAGC‑TATCCAGTGCCTGAAAAACGCTCCGGAAGCTGCTGCCGAGACCACCATGGCTGCGCTGCTGCGCGAAGCAATGGCCAAGCAGAACTGAGTTCTGt                                                                                                                                      <  8:164711/144‑1 (MQ=255)
 cgcgTTATCAGCCTGTCGATCAAATCGAAAGACGACGCTGAAGAGCGCGAAGC‑TATCCAGAGCCTGAAAAACGCTCCGGAAGCTGCTGCCGACACCACCATGGCTGCGCTGCTGCGCGAAGCAATGGCCAAGCAGAACTGAGTTCTGtt                                                                                                                                     <  4:149830/149‑1 (MQ=255)
 cgcgTTATCAGCCTGTCGATCAAATCGAAAGACGACGCTGAAGAGCGCGAAGC‑TATCCAGAGCCTGAAAAACGCTCCGGAAGCTGCTGCCGACACCACCATGGCTGCGCTGCTGCGCGAAGCAATGGCCAAACAGAACTCAGTTCTgct                                                                                                                                     >  7:104257/1‑147 (MQ=255)
   cgTTATCAGCCTGTCGATCAAATCGAAAGACGACGCTGAAGAGCGCGAAGC‑TATCCAGAGCCTGAAAAACGCTCCGGAAGCTGCTGCCGACACCACCATGGCTGCGCTGCTGCGCGAAGCAATGGCCAAGCAGAACTCAGTTCTGCTTg                                                                                                                                   <  7:85564/149‑1 (MQ=255)
   cgTTATCAGCCTGTAGATCAAATCGAAAGACGACGCTGAAGAGCGCGAAGC‑TATCCAGAGCCTGAAAAACGCTCCGGAAGCTGCTGCCGACACCACCATGGCTGCGCTGCTGCGCGAAGCAATGGCCAAGCAGAACTGAGTTCTGTTTg                                                                                                                                   <  2:23073/149‑1 (MQ=255)
      tATCAGCCTGTCGATCAAATCGAAAGACGACGCTGAAGAGCGCGAAGC‑TATCCAGAGCCTGAAAAACGCACCGGAACCTCCTGCCGACACCACCATGGCAGCGCTGCTGCGCGAAGCAATGGCCAAGCAGAACTGAGTTCTGt                                                                                                                                      >  7:164711/1‑143 (MQ=255)
        tCAGCCTGTCGATCAAATCGAAAGACGACGCTGAAGAGCGCGAAGC‑TATCCAGAGCCTGAAAAACGCTCCGGAAGCTGCTGCCGACACCACCATGGCTGCGCTGCTGCGCGAAGCAATGGCCAAACAGAACTCAGTTCTGCTTggcca                                                                                                                               <  8:104257/148‑5 (MQ=255)
           gCCTGTCGATCAAATCGAAAGACGACGCTGAAGAGCGCGAAGCATATCCAGAGCCTGAAAAACGCTCCGGAAGCTGCTGCCGACACCACCATGGCTGCGCTGCTGCGCGAAGCAATGGCCAAGCAGAACTGAGTTCTGTTTGATCGGt                                                                                                                             >  1:67164/1‑148 (MQ=255)
               gTCGATCAAATCGAAAGACGACGCTGAAGAGCGCGAAGC‑TATCCAGAGCCTGAAAAACGCTCCGGAAGCTGCTGCCGACACCACCATGGCTGCGCTGCTGCGCGAAGCAATGGCCAAGCAGAACTGAGTTCTGTTTGATCGGTAAAAAg                                                                                                                       >  3:282561/1‑149 (MQ=255)
                    tCAAATCGAAAGACGACGCTGAAGAGCGCGAAGC‑TATCCAGAGCCTGAAAAACGCTCCGGAAGCTGCTGCCGACACCACCATGGCTGCGCTGCTGCGCGAAGCAATGGCCAAGCAGAACTGAGTTCTGTTTGATCGGTAAAAAGGGCGa                                                                                                                  >  2:30128/1‑149 (MQ=255)
                    tCAAATCGAAAGACGACGCTGAAGAGCGCGAAGC‑TATCCAGAGCCTGAAAAACGCTCCCGAAGCTGCTGCCGACACCAACATGGCTGAGCTGCTGCGCGAAGCAATGGCCAAGCAGAACTGAGTTCTGTATGATCGGTAAAAAGGCCGa                                                                                                                  >  1:172590/1‑149 (MQ=255)
                         tCGAAAGACGACGCTGAAGAGCGCGAAGC‑TATCCAGAGCCTGAAAAACGCTCCGGAAGCTGCTGCCGACACCACCATGGCTGCGCTGCTGCGCGAAGCAATGGCCAAGCAGAACTGAGTTCTGTTTGATCGGTAAAAAGGGCGACCCTc                                                                                                             >  4:167162/1‑149 (MQ=255)
                                         aaGAGCGCGAAGC‑TATCCAGAGCCTGAAAAACGCTCCGGAAGCTGCTGCCGACACCACCATGGCTGCGCTGCTGCGCGAAGCAATGGCCAAGCAGAACTGAGTTCTGTTTGATCGGTAAAAAGGGCGACCCTCGGGTCGCCCttttttt                                                                                             <  8:109416/149‑1 (MQ=255)
                                             gcgcGAAGC‑TATCCAGAGCCTGAAAAACGCTCCGGAAGCTGCTGCCGACACCACCATGGCTGCGCTGCTGCGCGAAGCAATGGCCAAGCAGAACTGAGTTCTGTTTGATCGGTAAAAAGGGCGACCCTCGGGTCGCCCTTTTTTTTGCt                                                                                         >  6:318479/1‑149 (MQ=255)
                                              cgcgAAGC‑TATCCAGAGCCTGAAAAACGCTCCGGAAGCTGCTGCCGACACCACCATGGCTGCGCTGCTGCGCGAAGCAATGGCCAAGCAGAACTGAGTTCTGTTTGATCGGTAAAAAGGGCGACCCTCGGGTCGCCCTTTTTTTTGCTc                                                                                        <  1:30128/149‑1 (MQ=255)
                                                             gagCCTGAAAAACGCTCCGGAAGCTGCTGCCGACACCACCATGGCTGCGCTGCTGCGCGAAGCAATGGCCAAGCAGAACTGAGTTCTGtt                                                                                                                                     >  6:4063/1‑90 (MQ=255)
                                                             gagCCTGAAAAACGCTCCGGAAGCTGCTGCCGACACCACCATGGCTGCGCTGCTGCGCGAAGCAATGGCCAAGCAGAACTGAGTTCTGtt                                                                                                                                     <  5:4063/90‑1 (MQ=255)
                                                                  tGAAAAACGCTCCGGAAGCTGCTGCCGACACCACCATGGCTGCGCTGCTGCGCGAAGCAATGGCCAAGCAGAACTGAGTTCTGTTTGATCGGTAAAAAGGGCGACCCTCGGGTCGCCCTTTTTTTTGCTCcgc                                                                                     <  6:96164/133‑1 (MQ=255)
                                                                  tGAAAAACGCTCCGGAAGCTGCTGCCGACACCACCATGGCTGCGCTGCTGCGCGAAGCAATGGCCAAGCAGAACTGAGTTCTGTTTGATCGGTAAAAAGGGCGACCCTCGGGTCGCCCTTTTTTTTGCTCcgc                                                                                     >  5:96164/1‑133 (MQ=255)
                                                                  tGAAAAACGCTCCGGAAGCTGCTGCCGACACCACCATGGCTGCGCTGCTGCGCGAAGCAATGGCCAAGCAGAACTGAGTTCTGTTTGATCGGTAAAAAGGGCGACCCTCGGGTCGCCCTTTTTTTTGCTCCGCGACCAGGTGGTgcggc                                                                     >  3:155485/1‑149 (MQ=255)
                                                                   gAAAAACGCTCCGGAAGCTGCTGCCGACACCACCATGGCTGCGCTGCTGCGCGAAGCAATGGCCAAGCAGAACTGAGTTCTGTTTGATCGGTAAAAAGGGCGACCCTCGGGTCGc                                                                                                      >  3:188753/1‑115 (MQ=255)
                                                                   gAAAAACGCTCCGGAAGCTGCTGCCGACACCACCATGGCTGCGCTGCTGCGCGAAGCAATGGCCAAGCAGAACTGAGTTCTGTTTGATCGGTAAAAAGGGCGACCCTCGGGTCGc                                                                                                      <  4:188753/115‑1 (MQ=255)
                                                                   gAAAAACGCTCCGGAAGCTGCTGCCGACACCACCATGGCTGCGCTGCTGCGCGAAGCAATGGCCAAGCAGAACTGAGTTCTGTTTGATCGGTAAAAAGGGCGACCCTCGGGTCGCCCTTTTTTTTGCTCCGCGACCAGGTGGTgcggcg                                                                    >  1:56280/1‑149 (MQ=255)
                                                                       aaCGCTCCGGAAGCTGCTGCCGACACCACCATGGCTGCGCTGCTGCGCGAAGCAATGGCCAAGCAGAACTGAGTTCTGTTTGATCGGTAAAAAGGGCGACCCTCGGGTCGCCCTTTTTTTTGCTCCGCGACCAGGTGGTGCGGCGGATg                                                                <  3:167162/149‑1 (MQ=255)
                                                                          gCTCCGGAAGCTGCTGCCGACACCACCATGGCTGCGCTGCTGCGCGAAGCAATGGCCAAGCAGAACTGAGTTCTGTTTGATCGGTAAAAAGGGCGACCCTCGGGt                                                                                                         <  5:154545/105‑1 (MQ=255)
                                                                          gCTCCGGAAGCTGCTGCCGACACCACCATGGCTGCGCTGCTGCGCGAAGCAATGGCCAAGCAGAACTGAGTTCTGTTTGATCGGTAAAAAGGGCGACCCTCGGGt                                                                                                         >  6:154545/1‑105 (MQ=255)
                                                                          gCTCCGGAAGCTGCTGCCGACACCACCATGGCTGCGCTGCTGCGCGAAGCAATGGCCAAGCAGAACTGAGTTCTGTTTGATCGGTAAAAAGGGCGACCCTCGGGTCGCCCTTTTTTTTGCTCCGCGACCAGGTGGTGCGGCGGATGGCa                                                             >  6:286231/1‑149 (MQ=255)
                                                                                       ctgcCGACACCACCATGGCTGCGCTGCTGCGCGAAGCAATGGCCAAGCAGAACTGAGTTCTGTTTGATCGGTAAAAAGGGCGACCCTCGGGTCGCCCTTTTTTTTGCTCCGCGACCAGGTGGTGCGGCGGATGGCACTGCTGGAACACg                                                >  6:87656/1‑149 (MQ=255)
                                                                                         gcCGACACCACCATGGCTGCGCTGCTGCGCGAAGCAATGGCCAAGCAGAACTGAGTTCTGTTTGATCGGTAAAAAGGGCGACCCGAGGGTCGCCCTTTTTTTTGCTCCGCGACCAGGTGGTGCGGCGGATGGCACTGCTGGAACACGc                                               >  7:58732/1‑148 (MQ=255)
                                                                                             acacCACCATGGCTGCGCTGCTGCGCGAAGCAATGGCCAAGCAGAACTGAGTTCTGTTTGATCGGTAAAAAGGGCGACCCTCGGGTCGCCCTTTTTTTTGCTCCGCGACCAGGTGGTGCGGCGGATGGCACTGCTGGAACACGCCAGcc                                          <  4:282561/149‑2 (MQ=255)
                                                                                                  accaTGGCTGCGCTGCTGCGCGAAGCAATGGCCAAGCAGAACTGAGTTCTGTTTGATCGGTAAAAAGGGCGACCCTCGGGTCGCCCTTTTTTTTGCTCCGCGACCAGGTGGTGCGGCGGATGGCACTGCTGGAACACGCCAGCTTCAgg                                     >  7:254610/1‑149 (MQ=255)
                                                                                                   ccaTGGCTGCGCTGCTGCGCGAAGCAATGGCCAAGCAGAACTGAGTTCTGTTTGATCGGTAAAAAGGGCGACCCTCGGGTCGCCCTTTTTTTTGCTCCGCGACCAGGTGGTGCGGCGGATGGCACTGCTGGAACACGCCAGCTTCAGGc                                    >  6:153070/1‑149 (MQ=255)
                                                                                                         ctgcgctgcTGCGCGAAGCAATGGCCAAGCAGAACTGAGTTCTGTTTGATCGGTAAAAAGGGCGACCCTCGGGTCGCCCTTTTTTTTGCTCCGCGACCAGGTGGTGCGGCGGATGGCACTGCTGGAACACGCCAGCCTGAAGCTGGCgt                              <  2:56280/149‑2 (MQ=255)
                                                                                                            cgctgcTGCGCGAAGCAATGGCCAAGCAGAACTGAGTTCTGTTTGATCGGTAAAAAGGGCGACCCTCGGGTCGCCCTTTTTTTTGCTCCGCGACCAGGTGGTGCGGCGGATGGCACTGCTGGAACACGCCAGCTTCAGGCTGGCGAAGc                           <  5:318479/149‑1 (MQ=255)
                                                                                                            cgctgcTGCGCGAAGCAATGGCCAAGCAGAACTGAGTTCTGTTTGATCGGTAAAAAGGGCGACCCTCGGGTCGCCCTTTTTTTTGCTCCGCGACCAGGTGGTGCGGCGGATGGCACTGCTGGAACACGCCAGCCTGAAGCTGGCgtgtt                           <  6:101057/149‑5 (MQ=255)
                                                                                                                     gcgAAGCAATGGCCAAGCAGAACTGAGTTCTGTTTGATCGGTAAAAAGGGCGACCCTCGGGTCGCCCTTTTTTTTGCTCCGCGACCAGGTGGTGCGGCGGATGGCACTGCTGGAACACGCCAGCTTCAGGCTGGCGAAGCCCGTTGTAg                  <  1:209045/149‑1 (MQ=255)
                                                                                                                            aaTGGCCAAGCAGAACTGAGTTCTGTTTGATCGGTAAAAAGGGCGACCCTCGGGTCGCCCTTTTTTTTGCTCCGCGACCAGGTGGTGCGGCGGATGGCACTGCTGGAACACGCCAGCTTCAGGCTGGCGAAGCCCGTTGTAGGAGCggg           >  3:77823/1‑149 (MQ=255)
                                                                                                                              tGGCCAAGCAGAACTGAGTTCTGTTTGATCGGTAAAAAGGGCGACCCTCGGGTCGCCCTTTTTTTTGCTCCGCGACCAGGTGGTGCGGCGGATGGCACTGCTGGAACACGCCAGCTTCAGGCTGGCGAAGCCCGTTGTAGGAGCGGGtt         <  8:254610/149‑1 (MQ=255)
                                                                                                                                gCCAAGCAGAACTGAGTTCTGTTTGATCGGTAAAAAGGGCGACCCTCGGGTCGCCCTTTTTTTTGCTCCGCGACCAGGTGGTGCGGCGGATGGCACTGCTGGAACACg                                                >  4:92423/1‑108 (MQ=255)
                                                                                                                                gCCAAGCAGAACTGAGTTCTGTTTGATCGGTAAAAAGGGCGACCCTCGGGTCGCCCTTTTTTTTGCTCCGCGACCAGGTGGTGCGGCGGATGGCACTGCTGGAACACg                                                <  3:92423/108‑1 (MQ=255)
                                                                                                                                 ccAAGCAGAACTGAGTTCTGTTTGATCGGTAAAAAGGGCGACCCTCGGGTCGCCCTTTTTTTTGCTCCGCGACCAGGTGGTGCGGCGGATGGCACTGCTGGAACACGCCAGCTTCAGGCTGGCGAAGCCCGTTGTAGGAGCGGGTTCAc      >  2:254445/1‑149 (MQ=255)
                                                                                                                                  cAAACAGAACTCAGTTCTGTTTGATCGGTAAAAAGGGCGACCCTCGGGTCGCCCTTTTTTTTGCTCCGCGACCAGGTGGTGCGGCGGATGGCACTGCTGGAACACGCCAGCTTCAGGCTGGCGAAGCCCGTTGTAGGAGCGGGTTCAcc     >  5:90882/1‑149 (MQ=255)
                                                                                                                                     gCAGAACTGAGTTCTGTTTGATCGGTAAAAAGGGCGACCCTCGGGTCGCCCTTTTTTTTGCTCCGCGACCAGGTGGTGCGGCGGATGGCACTGCTGGAACACGCCAGCTTCAGGCTGGCGAAGCCCGTTGTAGGAGCGGGTTCACCcgc  >  2:139321/1‑149 (MQ=255)
                                                                                                                                      cAGAACTGAGTTCTGTTTGATCGGTAAAAAGGGCGACCCTCGGGTCGCCCTTTTTTTTGCTCCGCGACCAGGTGGTGCGGCGGATGGCACTGCTGGAACACGCCAGCTTCAGGCTGGCGAAGCCCGTTGTAg                  <  7:284955/132‑1 (MQ=255)
                                                                                                                                      cAGAACTGAGTTCTGTTTGATCGGTAAAAAGGGCGACCCTCGGGTCGCCCTTTTTTTTGCTCCGCGACCAGGTGGTGCGGCGGATGGCACTGCTGGAACACGCCAGCTTCAGGCTGGCGAAGCCCGTTGTAg                  >  8:284955/1‑132 (MQ=255)
                                                                                                                                         aaCTGAGTTCTGTTTGATCGGTAAAAAGGGCGACCCTCGGGTCGCCCTTTTTTTTGCTCCGCGACCAGGTGGTGCGGCGGATGGCACTGCTGGAACACGCCAGCTTCAGGCTGGCGAAGCCCGTTGTAg                  <  7:319931/129‑1 (MQ=255)
                                                                                                                                         aaCTGAGTTCTGTTTGATCGGTAAAAAGGGCGACCCTCGGGTCGCCCTTTTTTTTGCTCCGCGACCAGGTGGTGCGGCGGATGGCACTGCTGGAACACGCCAGCTTCAGGCTGGCGAAGCCCGTTGTAg                  >  8:319931/1‑129 (MQ=255)
                                                                                                                                             |                                                                                                                                            
CCGCGTTATCAGCCTGTCGATCAAATCGAAAGACGACGCTGAAGAGCGCGAAGC‑TATCCAGAGCCTGAAAAACGCTCCGGAAGCTGCTGCCGACACCACCATGGCTGCGCTGCTGCGCGAAGCAATGGCCAAGCAGAACTGAGTTCTGTTTGATCGGTAAAAAGGGCGACCCTCGGGTCGCCCTTTTTTTTGCTCCGCGACCAGGTGGTGCGGCGGATGGCACTGCTGGAACACGCCAGCTTCAGGCTGGCGAAGCCCGTTGTAGGAGCGGGTTCACCCGC  >  NC_002947/1980920‑1981200

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