Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,785,539 G→T 9.8% G98V (GGC→GTC)  rpsB → 30S ribosomal protein S2

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,785,5390GT 9.8% 153.5 / 4.9 50G98V (GGC→GTC) rpsB30S ribosomal protein S2
Reads supporting (aligned to +/- strand):  ref base G (25/20);  new base T (2/3);  total (27/23)
Fisher's exact test for biased strand distribution p-value = 6.51e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.61e-01

TGTTCAACGACGCTCTGTCGTTCGTAGAGCGCCTGGCCCAGGGCAAGAACAAGATCCTGTTCGTCGGCACCAAGCGTTCCGCCGGCAAGATCGTCGCCGAGCAAGCAGCTCGTTGCGGTTCGCCGTACGTTGACCACCGTTGGTTGGGCGGCATGCTGACCAACTACAAGACCATCCGCGCTTCGATCAAGCGTCTGCGCGACCTGGAAACCCAGGCCGAAGACGGCACTTTCGCCAAGCTGACCAAGAAAGAAGCCCTGATGCGCTCCCGCGACCTGGAAAAACTGGATCGCAGC  >  NC_002947/1785392‑1785687
                                                                                                                                                   |                                                                                                                                                    
tGTTCAACGACGCTCTGTCGTTCGTAGAGCGCCTGGCCCAGGGCAAGAACAAGATCCTGTTCGTCGGCACCAAGCGTTCCGCCGGCAAGATCGTCGCCGAGCAAGCAGCTCGTTGCGGTTCGCCGTACGTTGACCACCGTTGGTTGggc                                                                                                                                                     <  2:226556/149‑1 (MQ=255)
         acgCTCTGTCGTTCGTAGAGCGCCTGGCCCAGGGCAAGAACAAGATCCTGTTCGTCGGCACCAAGCGTTCCGCCGGCAAGATCGTCGCCGAGCAAGCAGCTCGTTGCGGTTCGCCGTACGTTGACCACCGTTGGTTGGGCGGCATGCTg                                                                                                                                            <  6:273491/149‑1 (MQ=255)
           gCTCTGTCGTTCGTAGAGCGCCTGGCCCAGGGCAAGAACAAGATCCTGTTCGTCGGCACCAAGCGTTCCGCCGGCAAGATCGTCGCCGAGCAAGCAGCTCGTTGCGGTTCGCCGTACGTTGACCACCGTTGGTTGGGCGGCATGCTGAc                                                                                                                                          <  6:126768/149‑1 (MQ=255)
              ctGTCGTTCGTAGAGCGCCTGGCCCAGGGCAAGAACAAGATCCTGTTCGTCGGCACCAAGCGTTCCGCCGGCAAGATCGTCGCCGAGCAAGCAGCTCGTTGCGGTTCGCCGTACGTTGACCACCGTTGGTTGGGCGGCATGCTGACCaa                                                                                                                                       >  4:123191/1‑149 (MQ=255)
                 tcgttcgtAGAGCGCCTGGCCCAGGGCAAGAACAAGATCCTGTTCGTCGGCACCAAGCGTTCCGCCGGCAAGATCGTCGCCGAGCAAGCAGCTCGTTGCGGTTCGCCGTACGTTGACCACCGTTGGTTGGTCAGCATGCCGCCCaa                                                                                                                                       <  1:136145/146‑1 (MQ=255)
                   gttcgtAGAGCGCCTGGCCCAGGGCAAGAACAAGATCCTGTTCGTCGGCACCAAGCGTTCCGCCGGCAAGATCGTCGCCGAGCAAGCAGCTCGTTGCGGTTCGCCGTACGTTGACCACCGTTGGTTGGTCAGCATGCCGCCCAAccaa                                                                                                                                   <  4:141791/148‑4 (MQ=255)
                   gttcgtAGAGCGCCTGGCCCAGGGCAAGAACAAGATCCTGTTCGTCGGCACCAAGCGTTCCGCCGGCAAGATCGTCGCCGAGCAAGCAGCTCGTTGCGGTTCGCCGTACGTTGACCACCGTTGGTTGGGCGGCATGCTGACCAACTAc                                                                                                                                   >  4:73210/1‑148 (MQ=255)
                      cgtAGAGCGCCTGGCCCAGGGCAAGAACAAGATCCTGTTCGTCGGCACCAAGCGTTCCGCCGGCAAGATCGTCGCCGAGCAAGCAGCTCGTTGCGGTTCGCCGTACGTTGACCACCGTTGGTTGGGCGGCATGCTGACCAACTACAAg                                                                                                                                <  3:123191/148‑1 (MQ=255)
                         agagCGCCTGGCCCAGGGCAAGAACAAGATCCTGTTCGTCGGCACAAAGCGTTCCGCCGGCAAGATCGTCGCCGAGCAAGCAGCTCGTTGCGGTTCGCCGTACGTTGACCACCGTTGGTTCGTCGTCATGCTGACCAACTACAATAcc                                                                                                                             >  7:104171/1‑148 (MQ=255)
                           agCGCCTGGCCCAGGGCAAGAACAAGATCCTGTTCGTCGGCACCAAGCGTTCCGCCGGCAAGATCGTCGCCGAGCAAGCAGCTCGTTGCGGTTCGCCGTACGTTGACCACCGTTGGTTGGGCGGCATGCTGACCAACTACAAGACCAta                                                                                                                          >  1:245905/1‑148 (MQ=255)
                           agCGCCTGGCCCAGGGCAAGAACAAGATCCTGTTCGTCGGCACCAAGCGTTCCGCCGGCAAGATCGTCGCCGAGCAAGCAGCTCGTTGCGGTTCGCCGTACGTTGACCACCGTTGGTTGGGCGGCATGCTGACCAACTACAAGACCATc                                                                                                                          >  7:12534/1‑149 (MQ=255)
                           agCGCCTGGCCCAGGGCAAGAACAAGATCCTGTTCGTCGGCACCAAGCGTTCCGCCGGCAAGATCGTCGCCGAGCAAGCAGCTCGTTGCGGTTCGCCGTACGTTGACCACCGTTGGTTGGGCGGCATGCTGACCAACTACAAGACCATc                                                                                                                          >  6:217070/1‑149 (MQ=255)
                              gcCTGGCCCAGGGCAAGAACAAGATCCTGTTCGTCGGCACCAAGCGTTCCGCCGGCAAGATCATCGCCGAGCAAGCAGCTCGTTGCGGTTCGCCGTACGTTGACCACCGTTGGTTGGGCGGCATGCTGACCAACTACAAGACCATCcgc                                                                                                                       >  1:47701/1‑149 (MQ=255)
                                     ccAGGGCAAGAACAAGATCCTGTTCGTCGGCACCAAGCGTTCCGCCGGCAAGATCGTCGCCGAGCAAGCAGCTCGTTGCGGTTCGCCGTACGTTGACCACCGTTGGTTGGGCGGCATGCTGACCAACTACAAGACCATCCGCGCTTCGa                                                                                                                <  8:12534/149‑1 (MQ=255)
                                      cAGGGCAAGAACAAGATCCTGTTCGTCGGCACCAAGCGTTCCGCCGGCAAGATCGTCGCCGAGCAAGCAGCTCGTTGCGGTTCGCCGTACGTTGACCACCGTTGGTTGGTCAGCATGCCGCCCaa                                                                                                                                       >  8:168484/1‑125 (MQ=255)
                                      cAGGGCAAGAACAAGATCCTGTTCGTCGGCACCAAGCGTTCCGCCGGCAAGATCGTCGCCGAGCAAGCAGCTCGTTGCGGTTCGCCGTACGTTGACCACCGTTGGTTGGTCAGCATGCCGCCCaa                                                                                                                                       <  7:168484/125‑1 (MQ=255)
                                         ggCAAGAACAAGATCCTGTTCGTCGGCACCAAGCGTTCCGCCGGCAAGATCGTCGCCGAGCAAGCAGCTCGTTGCGGTTCGCCGTACGTTGACCACCGTTGGTTGGGCGGCATGCTGACCAACTAc                                                                                                                                   <  7:286408/126‑1 (MQ=255)
                                         ggCAAGAACAAGATCCTGTTCGTCGGCACCAAGCGTTCCGCCGGCAAGATCGTCGCCGAGCAAGCAGCTCGTTGCGGTTCGCCGTACGTTGACCACCGTTGGTTGGGCGGCATGCTGACCAACTAc                                                                                                                                   >  8:286408/1‑126 (MQ=255)
                                         ggCAAGAACAAGATCCTGTTCGTCGGCACCAAGCGTTCCGCCGGCAAGATCGTCGCCGAGCAAGCAGCTCGTTGCGGTTCGCCGTACGTTGACCACCGTTGGTTGGGCGGCATGCTGACCAACTACAAGACCATCCGCGCTTCGATCaa                                                                                                            <  1:68932/149‑1 (MQ=255)
                                                 cAAGATCCTGTTCGTCGGCACCAAGCGTTCCGCCGGCAAGATCGTCGCCGAGCAAGCAGCTCGTTGCGGTTCGCCGTACGTTGACCACCGTTGGTTGGGCGGCATGCTGACCAACTACAAGACCATCCGCGCTTCGATCAAGCGTCTgc                                                                                                    >  7:299731/1‑149 (MQ=255)
                                                  aaGATCCTGTTCGTCGGCACCAAGCGTTCCGCCGGCAAGATCGTCGCCGAGCAAGCAGCTCGTTGCGGTTCGCCGTACGTTGACCACCGTTGGTTGGGCGGCATGCTGACCAACTACAAGACCATCCGCGCTTCGATCAAGCGTCTgcg                                                                                                   <  8:299731/149‑1 (MQ=255)
                                                       ccTGTTCGTCGGCACCAAGCGTTCCGCCGGCAAGATCGTCGCCGAGCAAGCAGCTCGTTGCGGTTCGCCGTACGTTGACCACCGTTGGTTGGGCGGCATGCTGACCAACTACAAGACCATCCGCGCTTCGATCAAGCGTCTGCGCGAcc                                                                                              >  6:327414/1‑149 (MQ=255)
                                                         tGTTCGTCGGCACCAAGCGTTCCGCCGGCAAGATCGTCGCCGAGCAAGCAGCTCGTTGCGGTTCGCCGTACGTTGACCACCGTTGGTTGGGCGGCATGCTGACCAACTACAAGACCATCCGCGCTTCGATCAAGCGTCTGCGCGACCTg                                                                                            >  7:177050/1‑149 (MQ=255)
                                                           ttCGTCGGCACCAAGCGTTCCGCCGGCAAGATCGTCGCCGAGCAAGCAGCTCGTTGCGGTTCGCCGTACGTTGACCACCGTTGGTTGGGCGGCATGCTGACCAACTACAAGACCATCCGCGCTTCGATCAAGCGTCTGCGCGACCTg                                                                                            <  6:253397/147‑1 (MQ=255)
                                                              gtcgGCACCAAGCGTTCCGCCGGCAAGATCGTCGCCGAGCAAGCAGCTCGTTGCGGTTCGCCGTACGTTGACCACCGTTGGTTGGGCGGCATGCTGACCAACTACAAGACCATCCGCGCTTCGATCAAGCGTCTg                                                                                                     >  3:184991/1‑135 (MQ=255)
                                                              gtcgGCACCAAGCGTTCCGCCGGCAAGATCGTCGCCGAGCAAGCAGCTCGTTGCGGTTCGCCGTACGTTGACCACCGTTGGTTGGGCGGCATGCTGACCAACTACAAGACCATCCGCGCTTCGATCAAGCGTCTg                                                                                                     <  4:184991/135‑1 (MQ=255)
                                                                 ggCACCAAGCGTTCCGCCGGCAAGATCGTCGCCGAGCAAGCAGCTCGTTGCGGTTCGCCGTACGTTGACCACCGTTGGTTGGGCGGCATGCTGACCAACTACAAGACCATCCGCGCTTCGATCAAGCGTCTGCGCGACCTGGAAACCCa                                                                                    <  3:73210/149‑1 (MQ=255)
                                                                  gCACCAAGCGTTCCGCCGGCAAGATCGTCGCCGAGCAAGCAGCTCGTTGCGGTTCGCCGTACGTTGACCACCGTTGGTTGGGCGGCATGCTGACCAACTACAAGACCATCCGCGCTTCGATCAAGCGTCTGCGCGACCTGGAAACCCAg                                                                                   >  2:111489/1‑149 (MQ=255)
                                                                          cGTTCCGCCGGCAAGATCGTCGCCGAGCAAGCAGCTCGTTGCGGTTCGCCGTACGTTGACCACCGTTGGTTGGGCGGCATGCTGACCAACTACAAGACCATCCGCGCTTCGATCAAGAGTCTGAGCGAACTGGAAACCCAGGCCGAAGa                                                                           >  8:233042/1‑149 (MQ=255)
                                                                               cgccgGCAAGATCGTCGCCGAGCAAGCAGCTCGTTGCGGTTCGCCGTACGTTGACCACCGTTGGTTGGGCGGCATGCTGACCAACTACAAGACCATCCGCGCTTCGATCAAGCGTCTGCGCGACCTGGAAACCCAGGCCGAAGACGGCa                                                                      >  5:53905/1‑149 (MQ=255)
                                                                                 ccgGCAAGATCGTCGCCGAGCAAGCAGCTCGTTGCGGTTCGCCGTACGTTGACCACCGTTGGTTGGGCGGCATGCTGACCAACTACAAGACCATCCGCGCTTCGATCAAGCGTCTGCGCGACCTGGAAACCCAGGCCGAAGACGGCACt                                                                    <  8:177050/149‑1 (MQ=255)
                                                                                   ggCAAGATCGTCGCCGAGCAAGCAGCTCGTTGCGGTTCGCCGTACGTTGACCACCGTTGGTTGGGCGGCATGCTGACCAACTACAAGa                                                                                                                               >  1:111422/1‑88 (MQ=255)
                                                                                   ggCAAGATCGTCGCCGAGCAAGCAGCTCGTTGCGGTTCGCCGTACGTTGACCACCGTTGGTTGGGCGGCATGCTGACCAACTACAAGa                                                                                                                               <  2:111422/88‑1 (MQ=255)
                                                                                   ggCAAGATCGTCGCCGAGCAAGCAGCTCGTTGCGGTTCGCCGTACGTTGACCACCGTTGGTTGGGCGGCATGCTGACCAACTACAAGACCATCCGCGCTTCGATCAAGCGTCTGCGCGACCTGGAAACCCAGGCCGAAGACGGCACtt                                                                   >  6:139487/1‑148 (MQ=255)
                                                                                     cAAGATCGTCGCCGAGCAAGCAGCTCGTTGCGGTTCGCCGTACGTTGACCACCGTTGGTTGGGCGGCATGCTGACCAACTACAAGACCATCCGCGCTTCGATCAAGCGTCTGCGCGACCTGGAAACCCAGGCGGAAGACGGCACTTTCg                                                                <  2:245905/149‑1 (MQ=255)
                                                                                               gCCGAGCAAGCAGCTCGTTGCGGTTCGCCGTACGTTGACCACCGTTGGTTGGGCGGCATGCTGACCAACTACAAGACCATCCGCGCTTCGATCAAGCGTCTGCGCGACCTGGAAACCCAGGCCGAAGACGGCACTTTCGCCAAGCTGAc                                                      <  5:139487/149‑1 (MQ=255)
                                                                                                  gAGCAAGCAGCTCGTTGCGGTTCGCCGTACGTTGACCACCGTTGGTTGGGCGGCATGCTGACCAACTACAAGACCATCCGCGCTTCGATCAAGCGTCTGCGCGACCTGGAAACCCAGGCCGAAGACGGCACTTTCGCCAAg                                                           <  2:277625/141‑1 (MQ=255)
                                                                                                  gAGCAAGCAGCTCGTTGCGGTTCGCCGTACGTTGACCACCGTTGGTTGGGCGGCATGCTGACCAACTACAAGACCATCCGCGCTTCGATCAAGCGTCTGCGCGACCTGGAAACCCAGGCCGAAGACGGCACTTTCGCCAAg                                                           >  1:277625/1‑141 (MQ=255)
                                                                                                    gcaagcaGCTCGTTGCGGTTCGCCGTACGTTGACCACCGTTGGTTGGGCGGCATGCTGACCAACTACAAGACCATCCGCGCTTCGATCAAGCGTCTGCGCGACCTGGAAACCCAGGCCGAAGACGGCACTTTCGCCAAGCTGACCaaga                                                 <  2:47701/149‑1 (MQ=255)
                                                                                                     caagcaGCTCGTTGCGGTTCGCCGTACGTTGACCACCGTTGGTTGGGCGGCATGCTGACCAACTACAAGACCATCCGCGCTTCGATCAAGCGTCTGCGCGACCTGGAAACCCAGGCCGAAGACGGCACTTTCGCCAAGCTGACCaaga                                                 >  6:205099/1‑148 (MQ=255)
                                                                                                      aagcaGCTCGTTGCGGTTCGCCGTACGTTGACCACCGTTGGTTGGGCGGCATGCTGACCAACTACAAGACCATCCGCGCTTCGATCAAGCGTCTGCGCGACCTGGAAACCCAGGCCGAAGACGGCACTTTCGCCAAGCTGACCaaga                                                 >  2:116975/1‑147 (MQ=255)
                                                                                                         cagcTCGTTGCGGTTCGCCGTACGTTGACCACCGTTGGTTGGGCGGCATGCTGACCAACTACAAGACCATCCGCGCTTCGATCAAGCGTCTGCGCGACCTGGAAACCCAGGCCGAAGACGGCACTTTCGCCAAGCTGACCaagaaag                                              >  6:27786/1‑147 (MQ=255)
                                                                                                         cagcTCGTTGCGGTTCGCCGTACGTTGACCACCGTTGGTTGGGCGGCATGCTGACCAACTACAAGACCATCCGCGCTTCGATCAAGCGTCTGCGCGACCTGGAAACCCAGGCCGAAGACGGCACTTTCGCCAAGCTGACCAagaaagaa                                            <  5:27786/149‑1 (MQ=255)
                                                                                                                  gCGGTTCGCCGTACGTTGACCACCGTTGGTTGTGCGGCATGCTGACCAACTACAAGACCATCCGCGCTTCGATCAAGCGTCTGCGCGACCTGGAAACCCAGGCCGAAGACGGCACTTTCGCCAATCTGACCAAGAAAGAAGCCCTGAt                                    <  7:233042/148‑1 (MQ=255)
                                                                                                                  gCGGTTCGCCGTACGTTGACCACCGTTGGTTGGGCGGCATGCTGACCAACTACAAGACCATCCGCGCTTCGATCAAGCGTCTGCGCGACCTGGAAACCCAGGCCGAAGACGGCACTTTCGCCAAGCTGACCAAGAAAGAAGCCCTGAt                                    <  1:116975/148‑1 (MQ=255)
                                                                                                                   cGGTTCGCCCTACGTTGACCAAAGTTGCTTCCGCGGCATGCTGACCAACTACAAGACCAACCGCGCTTCGCTCAAGCGTCTGCGCGACCTGGAAACCCAGGCCGAAGACGCCACTTTCGACAAGCGGACCAAGAAAGAAGCCCTCATgc                                  >  4:43450/1‑149 (MQ=255)
                                                                                                                    ggTTCGCCGTACGTTGACCACCGTTGGTTGGGCGGCATGCTGACCAACTACAAGACCATCCGCGCTTCGATCAAGCGTCTGCGCGACCTGGAAACCCAGGCCGAAGACGGCACTTTCGCCAAGCTGACCAAGAAAGAAGCCCTGATgcg                                 >  8:303543/1‑149 (MQ=255)
                                                                                                                                 ttGACCACCGTTGGTTGGGCGGCATGCTGACCAACTACAAGACCATCCGCGCTTCGATCAAGCGTCTGCGCGACCTGGAAACCCAGGCCGAAGACGGCACTTTCGCCAAGCTGACCAAGAAAGAAGCCCTGATGCGCTCCCGCGACCTg                    >  4:105798/1‑149 (MQ=255)
                                                                                                                                 ttGACCACCGTTGGTTGGGCGGCATGCTGACCAACTACAAGACCATCCGCGCTTCGATCAAGCGTCTGCGCGACCTGGAAACCCAGGCCGAAGACGGCACTTTCGCCAAGCTGACCAAGAAAGAAGCCCTGATGCGCTCCCGCGACCTg                    >  1:319364/1‑149 (MQ=255)
                                                                                                                                  tGACCACCGTTGGTTGGGCGGCATGCTGACCAACTACAAGACCATCCGCGCTTCGATCAAGCGTCTGCGCGACCTGGAAACCCAGGCCGAAGACGGCACTTTCGCCAAGCTGACCAAGAAAGAAGCCCTGAt                                    <  6:187525/132‑1 (MQ=255)
                                                                                                                                  tGACCACCGTTGGTTGGGCGGCATGCTGACCAACTACAAGACCATCCGCGCTTCGATCAAGCGTCTGCGCGACCTGGAAACCCAGGCCGAAGACGGCACTTTCGCCAAGCTGACCAAGAAAGAAGCCCTGAt                                    >  5:187525/1‑132 (MQ=255)
                                                                                                                                                   gcggcATGCTGACCAACTACAAGACCATCCGCGCTTCGATCAAGCGTCTGCGCGACCTGGAAACCCAGGCCGAAGACGGCACTTTCGCCAAGCTGACCAAGAAAGAAGCCCTGATGCGCTCCCGCGACCTGGAAAAACTGGATCGCAGc  >  1:117023/1‑149 (MQ=255)
                                                                                                                                                   |                                                                                                                                                    
TGTTCAACGACGCTCTGTCGTTCGTAGAGCGCCTGGCCCAGGGCAAGAACAAGATCCTGTTCGTCGGCACCAAGCGTTCCGCCGGCAAGATCGTCGCCGAGCAAGCAGCTCGTTGCGGTTCGCCGTACGTTGACCACCGTTGGTTGGGCGGCATGCTGACCAACTACAAGACCATCCGCGCTTCGATCAAGCGTCTGCGCGACCTGGAAACCCAGGCCGAAGACGGCACTTTCGCCAAGCTGACCAAGAAAGAAGCCCTGATGCGCTCCCGCGACCTGGAAAAACTGGATCGCAGC  >  NC_002947/1785392‑1785687

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