Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 6,172,418 A→G 8.1% W222R (TGG→CGG) ‡ PP_5410 ← DeoR family transcriptional regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029476,172,4180AG 8.1% 135.1 / 2.7 49W222R (TGG→CGG) ‡PP_5410DeoR family transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base A (26/19);  new base G (2/2);  total (28/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

GCTCGGCGACAGCACCAACGATTCCACACATGTTTTATTGCTCCTGGCTGATCGCGGCGCAGATCAGGTTGATGCCGCGGGCTTGAATTTGTTCGCGCGCTGCTGCTGGCAGGCGTTCATCTGTAATAAGGGTGTGCACACTGCCCCAGGGCAGCTCGAGGTTAGGGATCTTGCGCCCCACCTTGTCCGACTCGACCATCACGATCACTTCACGGGCAACTTCCGCCATGACCCGGCTCAGCCCAAGCAACTCGTTGAAGGTAGTGGTACCCCGGCCCAGGTCGATGCCGTCGGCG  >  NC_002947/6172271‑6172566
                                                                                                                                                   |                                                                                                                                                    
gCTCGGCGACAGCACCAACGATTCCACACATGTTTTATTGCTCCTGGCTGATCGCGGCGCAGATCAGGTTGATGCCGCGGGCTTGAATTTGTTCGCGCGCTGCTGCTGGCAGGCGTTCATCTGTAATAAGGGTGTGCACACTGCCCCAg                                                                                                                                                     >  3:136304/1‑149 (MQ=255)
            cACCAACGATTCCACACATGTTTTATTGCTCCTGGCTGATCGCGGCGCAGATCAGGTTGATGCCGCGGGCTTGAATTTGTTCGCGCGCTGCTGCTGGCAGGCGTTCATCTGTAATAAGGGTGTGCACACTGCCCCAGGGCAGCTCGAgg                                                                                                                                         <  1:247437/149‑1 (MQ=255)
               cAACGATTCCACACATGTTTTATTGCTCCTGGCTGATCGCGGCGCAGATCAGGTTGATGCCGCGGGCTTGAATTTGTTCGCGCGCTACTGCTGGCAGGCGTTCATCTGTAATAAGGGTGTGCACACTGCCCCAGGGCAGCTCGAGGTTa                                                                                                                                      <  4:184092/149‑1 (MQ=255)
               cAACGATTCCACACATGTTTTATTGCTCCTGGCTGATCGAGGCGCAGATCAGGTTGATGCCGCGGGCTTGAATTTGTTCGCGCGCTGCTGCTGGCAGGCGTTCATCTGTAATAAGGGTGTGCAAACTGCCCAAGGGCAGCTCGAGGt                                                                                                                                        >  3:128258/1‑147 (MQ=255)
                         acacaTGTTTTATTGCTCCTGGCTGATCGCGGCGCAGATCAGGTTGATGCCGCGGGCTTGAATTTGTTCGCGCGCTGCTGCTGGCAGGCGTTCATCTGTAATAAGGGTGTGCACACTGCCCTGGGGCAgtgt                                                                                                                                             <  4:234525/132‑4 (MQ=255)
                         acacaTGTTTTATTGCTCCTGGCTGATCGCGGCGCAGATCAGGTTGATGCCGCGGGCTTGAATTTGTTCGCGCGCTGCTGCTGGCAGGCGTTCATCTGTAATAAGGGTGTGCACACTGCCCTGGGGCAgtgt                                                                                                                                             >  3:234525/1‑129 (MQ=255)
                         acacaTGTTTTATTGCTCCTGGCTGATCGCGGCGCAGATCAGGTTGATGCCGCGGGCTTGAATTTGTTCGCGCGCTGCTGCTGGCAGGCGTTCATCTGTAATAAGGGTGTGCACACTGCCCCAGGGCAGCTCGAGGTTAGGGATCTTgc                                                                                                                            >  3:160660/1‑149 (MQ=255)
                              tGTTTTATTGCTCCTGGCTGATCGCGGCGCAGATCAGGTTGATGCCGCGGGCTTGAATTTGTTCGCGCGCTGCTGCTGGCAGGCGTTCATCTGTAATAAGGGTGTGCACACTGCCCCAGGGCAGCTCGAGGTTAGGGATCTTGCGcccc                                                                                                                       >  2:29388/1‑149 (MQ=255)
                                 tttATTGCTCCTGGCTGATCGCGGCGCAGATCAGGTTGATGCCGCGGGCTTGAATTTGTTCGCGCGCTGCTGCTGGCAGGCGTTCATCTGTAATAAGGGTGTGCACACTGCCCCAGGGCAGATCGAGGTTAGGGATCTTGCGCCCCAcc                                                                                                                    >  7:94074/1‑149 (MQ=255)
                                       gCTCCTGGCTGATCGCGGCGCAGATCAGGTTGATGCGGCGGGCTTGAATTTGTTCGCGCGCTGCTGCTGGCAGGCGTTCATCTGTAATAAGGGTGTGCACACTGCCCCAGGGCAGCTCGAGGTTAGGGATCTTGCGCCCCACCTTGTc                                                                                                               <  4:136304/148‑1 (MQ=255)
                                                  atggcggcgGAGATCAGGTTGATGCCGCGGGCTTGAATTTGTTCGCGCGCTGCTGCTGGCAGGCGTTCATCTGTAATAAGGGTGTGCACACTGCCCCAGGGCAGCTCGAGGTTAGGGATCTTGCGCCCCACCTTGTCCGACTCGACCa                                                                                                    <  6:183469/145‑1 (MQ=255)
                                                        gcgcAGATCAGGTTTATGCCGCGGGCTTGAATTTGTTCGCGCGCTGCTGCTGGCAGGCGTTCATCTGTAATAAGGGTGTGCACACTGCCCCAGGGCAGCTCGAGGTTAGGGATCTTGCGCCCCACCTTGTCCGACTCGACCATCACGAt                                                                                             <  4:233079/149‑1 (MQ=255)
                                                          gcAGATCAGGTTGATGCCGCGGGCTTGAATTTGTTCGCGCGCTGCTGCTGGCAGGCGTTCATCTGTAATAAGGGTGTGCACACTGCCCCAGGGCAGCTCGAGGTTAGGGATCTTGCGCCCCACCTTGTCCGACTCGACCATCACGATc                                                                                            >  4:33692/1‑148 (MQ=255)
                                                            aGATCAGGTTGATGCCGCGGGCTTGAATTTGTTCGCGCGCTGCTGCTGGCAGGCGTTCATCTGTAATAAGGGTGTGCACACTGCCCCAGGGCAGCTCGAGGTTAGGGATCTTGCGCCCCACCTTGTCCGACTCGACCATCACGATCACt                                                                                         >  2:301173/1‑149 (MQ=255)
                                                                  ggTTGATGCCGCGGGCTTGAATTTGTTCGCGCGCTGCTGCTGGCAGGCGTTCATCTGTAATAAGGGTGTGCACACTGCCCCAGGGCAGCTCGAGGTTAGGGATCTTGCGCCCCACCTTGTCCGACTCGACCATCACGATCACTTCACgg                                                                                   <  4:136200/149‑1 (MQ=255)
                                                                    ttGATGCCGCGGGCTTGAATTTGTTCGCGCGCTGCTGCTGGCAGGCGTTCATCTGTAATAAGGGTGTGCACACTGCCCCAGGGCAGCTCGAGGTTAGGGATCTTGCGCCCCACCTTGTCCGACTCGACCATCACGATCACTTc                                                                                       <  7:231533/143‑1 (MQ=255)
                                                                    ttGATGCCGCGGGCTTGAATTTGTTCGCGCGCTGCTGCTGGCAGGCGTTCATCTGTAATAAGGGTGTGCACACTGCCCCAGGGCAGCTCGAGGTTAGGGATCTTGCGCCCCACCTTGTCCGACTCGACCATCACGATCACTTc                                                                                       >  8:231533/1‑143 (MQ=255)
                                                                       aTGCCGCGGGCTTGAATTTGTTCGCGCGCTGCTGCTGGCAGGCGTTCATCTGTAATAAGGGTGTGCACACTGCCCCAGGGCAGCTCGAGGTTAGGGATCTTGCGCCCCACCTTGTCCGACTCGACCATCACGATCACTTCACGGGCAAc                                                                              >  2:214655/1‑149 (MQ=255)
                                                                            gcgGGCTTGAATTTGTTCGCGCGCTGCTGCTGGCAGGCGTTCATCTGTAATAAGGGTGTGCACACTGCCCCAGGGCAGCTCGAGGTTAGGGATCTTGCGCCCCACCTTGTCCGACTCGACCATCACGATCACTTCACGGGCAACTTCCg                                                                         >  7:226802/1‑149 (MQ=255)
                                                                                  ttGAATTTGTTCGCGCGCTGCTGCTGGCAGGCGTTCATCTGTAATAAGGGTGTGCACACTGCCCCAGGGCAGCTCGAGGTTAGGGATCTTGCGCCCCACCTTGTCCGACTCGACCATCACGATCACTTCACGGGCAACTTCCGCCATGa                                                                   >  7:241960/1‑149 (MQ=255)
                                                                                  ttGAATTTGTTCGCGCGCTGCTGCTGGCAGGCGTTCATCTGTAATAAGGGTGTGCACACTGCCCCAGGGCAGCTCGAGGTTAGGGATCTTGCGCCCCACCTTGTCCGACTCGACCATCACGATCACTTCACGGGCAACTTCCGCCATGa                                                                   >  3:282334/1‑149 (MQ=255)
                                                                                    gAATTTGTTCGCGCGCTGCTGCTGGCAGGCGTTCATCTGTAATAAGGGTGTGCACACTGCCCCAGGGCAGCTCGAGGTTAGGGATCTTGCGCCCCACCTTGTCCGACTCGACCATCACGATCACTTCACGGGCAACTTCCGCCATGAcc                                                                 >  2:115325/1‑149 (MQ=255)
                                                                                    gAATTTGTTCGCGCGCTGCTGCTGGCAGGCGTTCATCTGTAATAAGGGTGTGCACACTGCCCCAGGGCAGCTCGAGGTTAGGGATCTTGCGCCCCACCTTGTCCGACTCGACCATCACGATCACTTCACGGGCAACTTCCGCCATGAcc                                                                 >  5:276355/1‑149 (MQ=255)
                                                                                     aaTTTGTTCGCGCGCTGCTGCTGGCAGGCGTTCATCTGTAATAAGGGTGTGCACACTGCCCCAGGGCAGCTCGAGGTTAGGGATCTTGCGCCCCACCTTGTCCGACTCGACCATCACGATCACTTCACGGGCAACTTCCGCCATGAccc                                                                <  4:282334/149‑1 (MQ=255)
                                                                                      aTTTGTTCGCGCGCTGCTGCTGGCAGGCGTTCATCTGTAATAAGGGTGTGCACACTGCCCCAGGGCAGCTCGAGGTTAGGGATCTTGCGCCCCACCTTGTCCGACTCGACCATCACGATCACTTCACGGGCAACTTCCGCCATGACCCg                                                               <  2:142487/149‑1 (MQ=255)
                                                                                                 cgcTGCTGCTGGCAGGCGTTCATCTGTAATAAGGGTGTGCACACTGCCCCAGGGCAGCTCGAGGTTAGGGATCTTGCGCCCCACCTTGTCCGACTCGACCATCACGATCACTTCACGGGCAACTTCCGCCATGACCCGGCTCAGCCCaa                                                    <  8:241960/149‑1 (MQ=255)
                                                                                                  gctgctgctGGCATGCGTTCATCTGTATTAAGGGTGTGCACACTGCCCCAGGGCAGCTCGAGGTTAGGGATCTTGCGGCCCACCTTGTCGGACTCGACCATCAGGATCACTTCACGGGCAACTTCCGCCATGACCCGGCTCAGCCCAAg                                                   <  4:128258/149‑1 (MQ=255)
                                                                                                  gctgctgctGGCAGGCGTTCATCTGTAATAAGGGTGTGCACACTGCCCCAGGGCAGCTCGAGGTTAGGGATCTTGCGCCCCACCTTGTCCGACTCGACCATCACGATCACTTCACGGGCAACTTCCGCCATGACCCGGCTCAGCCCAAg                                                   <  8:226802/149‑1 (MQ=255)
                                                                                                    tgctgctgGCAGGCGTTCATCTGTAATAAGGGTGTGCACACTGCCCCAGGGCAGCTCGAGGTTAGGGATCTTGCGCCCCACCTTGTCCGACTCGACCATCACGATCACTTCACGGGCAACTTCCGCCATGACCCGGCTCAGCCCAAGc                                                  >  3:318750/1‑148 (MQ=255)
                                                                                                         ctgGCAGGCGTTCATCTGTAATAAGGGTGTGCACACTGCCCCAGGGCAGCTCGAGGTTAGGGATCTTGCGCCCCACCTTGTCCGACTCGACCATCACGATCACTTCACGGGCAACTTCCGCCATGACCCGGCTCAGCCCAAGCAACTc                                             <  3:33692/148‑1 (MQ=255)
                                                                                                               ggCGTTCATCTGTAATAAGGGTGTGCACACTGCCCCAGGGCAGCTCGAGGTTAGGGATCTTGCGCCCCACCTTGTCCGACTCGACCATCACGATCACTTCACGGGCAACTTCCGCCATGACCCGGCTCAGCCCAAGCAACTCGTTGAAg                                      <  4:160660/149‑1 (MQ=255)
                                                                                                                       tCTGTAATAAGGGTGTGCACACTGCCCCAGGGCAGCTCGAGGTTAGGGATCTTGCGCCCCACCTTGTCCGACTCGACCATCACGATCACTTCACGGGCAACTTCCGCCATGACCCGGCTCAGCCCAAGCAACTCGTTGAAGGTAg                                  >  6:186219/1‑145 (MQ=255)
                                                                                                                       tCTGTAATAAGGGTGTGCACACTGCCCCAGGGCAGCTCGAGGTTAGGGATCTTGCGCCCCACCTTGTCCGACTCGACCATCACGATCACTTCACGGGCAACTTCCGCCATGACCCGGCTCAGCCCAAGCAACTCGTTGAAGGTAg                                  <  5:186219/145‑1 (MQ=255)
                                                                                                                          gTAATAAGGGTGTGCACACTGCCCCAGGGCAGCTCGAGGTTAGGGATCTTGCGCCCCACCTTGTCCGACTCGACCATCACGATCACTTCACGGGCAACTTCCGCCATGACCCGGCTCAGCCCAAGCAACTCGTTGAAGGTAGTGGTAcc                           >  1:128959/1‑149 (MQ=255)
                                                                                                                               aaGGGTGTGCACACTGCCCCAGGGCAGCTCGAGGTTAGGGATCTTGCGCCCCACCTTGTCCGACTCGACCATCACGATCACTTCACGGGCAACTTCCGCCATGACCCGGCTCAGCCCAAGCAACTCGTTGAAGGTAGTg                                <  5:300858/139‑1 (MQ=255)
                                                                                                                               aaGGGTGTGCACACTGCCCCAGGGCAGCTCGAGGTTAGGGATCTTGCGCCCCACCTTGTCCGACTCGACCATCACGATCACTTCACGGGCAACTTCCGCCATGACCCGGCTCAGCCCAAGCAACTCGTTGAAGGTAGTg                                >  6:300858/1‑139 (MQ=255)
                                                                                                                                  ggTGTGCACACTGCCCCAGGGCAGCTCGAGGTAGGGGATCTTGCGCCCCACCTTGTCCGACTCGACCATCACGATCACGTCACGGGCAACTTCCGCCATGACCCGGCTCAGCCAAAGCAACTCGTTGAAGGTAGTGGTACCCCGGcc                     >  4:283261/1‑147 (MQ=255)
                                                                                                                                        cgagcTGCCCTGGGGCAGCTCGAGGTTAGGGATCTTGCGCCCCACCTTGTCCGACTCGACCATCACGATCACTTCACGGGCAACTTCCGCCATGACCCGGCTCAGCCCAAGCAACTCGTTGAAGGTAGTGGTACCCCGGCCCAGGTCGa             >  3:101215/5‑149 (MQ=255)
                                                                                                                                        cacaCTGCCCCAGGGCAGCTCGAGGTTAGGGATCTTGCGCCCCACCTTGTCCGACTCGACCATCACGATCACTTCACGTGCAACTTCCGCCAt                                                                     >  6:300984/1‑93 (MQ=255)
                                                                                                                                        cacaCTGCCCCAGGGCAGCTCGAGGTTAGGGATCTTGCGCCCCACCTTGTCCGACTCGACCATCACGATCACTTCACGGGCAACTTCCGCCAt                                                                     <  5:300984/93‑1 (MQ=255)
                                                                                                                                        cacaCTGCCCCAGGGCAGCTCGAGGTTAGGGATCTTGCGCCCCACCTTGTCCGACTCGACCATCACGATCACTTCACGGGCAACTTCCGCCATGACCCGGCTCAGCCCAAGCAACTCGTTGAAGGTAGTGGTACCCCGGCCCAGGTCGa             >  8:279109/1‑149 (MQ=255)
                                                                                                                                         gagcTGCCCTGGGGCAGCTCGAGGTTAGGGATCTTGCGCCCCACCTTGTCCGACTCGACCATCACGATCACTTCACGGGCAACTTCCGCCATGACCCGGCTCAGCCCAAGCAACTCGTTGAAGGTAGTGGTACCCCGGCCCAGGTCGAt            <  4:101215/146‑1 (MQ=255)
                                                                                                                                         acacTGCCCCAGGGCAGCTCGAGGTTAGGGATCTTGCGCCCCACCTTGTCCGACTCGACCATCACGATCACTTCACGGGCAACTTCCGCCATGACCCGGCTCAGCCCAAGCAACTCGTTGAAGGTAGTGGTACCCCGGCCCAGGTCGAt            <  1:301173/149‑1 (MQ=255)
                                                                                                                                         acacTGCCCCAGGGCAGCTCGAGGTTAGGGATCTTGCGCCCCACCTTGTCCGACTCGACCATCACGATCACTTCACGGGCAACTTCCGCCATGACCCGGCTCAGCCCAAGCAACTCGTTGAAGGTAGTGGTACCCCGGCCCAGGTCGAt            <  1:214655/149‑1 (MQ=255)
                                                                                                                                         acacTCCCCCAGGGCAGCTCGAGGTTAGGGATCTTGCGCCCCACCATGTCCGACTCCACCATCACGATCACTTCACGGGCAAATTCAGCCATGAACAGGCTCAGCCCAAGCAACTCGTTAAAGCTAGTGGTACCACCGCCCacgacta             >  4:26298/1‑142 (MQ=255)
                                                                                                                                            cTGCCCCAGGGCAGCTCGAGGTTAGGGATCTTGCGCCCCACCTTGTCCGACTCGACCATCACGATCACTTCACGGGCTACTTCCGCCATGACCCGGCTCAGCCCAAGCAACTCGTTGAAGGTAGTGGTACCCCGGCCCAGGTCGATGcc         >  1:104028/1‑149 (MQ=255)
                                                                                                                                               ccccAGGGCAGCTCGAGGTTAGGGATCTTGCGCCCCACCTTGTCCGACTCGACCATCACGATCACTTCACGGGCAACTTCCGCCATGACCCGGCTCAGCCCAAGCAACTCGTTGAAGGTAGTGGTACCCCGGCCCAGGTCGATGCCGTc      >  7:60170/1‑149 (MQ=255)
                                                                                                                                                cccAGGGCAGCTCGAGGTTAGGGATCTTGCGCCCCACCTTGTCCGACTCGACCATCACGATCACTTCACGGGCAACTTCCGCCATGACCCGGCTCAGCCCAAGCAACTCGTTGAAGGTAGTGGTACCCCGGCCCAGGTCGATGCCGTCg     >  5:91067/1‑149 (MQ=255)
                                                                                                                                                  cAGGGCAGCTCGAGGTTAGGGATCTTGCGAAACACCTTGTCCGACTAGACCATCACGATCACTTCACGGGCAACTACCGCAATGAACCGGCTCAGCCCAACCAACTCGTTGAAGGTAGTGGTAACCCGGCCCAGGTCGATGACGTCg     >  8:248832/1‑147 (MQ=255)
                                                                                                                                                   aGGGCAGCTCGAGGTTAGGGATCTTGCGCCCCACCTTGTCCGACTCGACCATCACGATCACTTCACGGGCAACTTCCGCCATGACCCGGCTCAGCCCAAGCAACTCGTTGAAGGTAGTGGTACCCCGGCCCAGGTCGATGCCGTCGgcg  >  2:193809/1‑149 (MQ=255)
                                                                                                                                                   |                                                                                                                                                    
GCTCGGCGACAGCACCAACGATTCCACACATGTTTTATTGCTCCTGGCTGATCGCGGCGCAGATCAGGTTGATGCCGCGGGCTTGAATTTGTTCGCGCGCTGCTGCTGGCAGGCGTTCATCTGTAATAAGGGTGTGCACACTGCCCCAGGGCAGCTCGAGGTTAGGGATCTTGCGCCCCACCTTGTCCGACTCGACCATCACGATCACTTCACGGGCAACTTCCGCCATGACCCGGCTCAGCCCAAGCAACTCGTTGAAGGTAGTGGTACCCCGGCCCAGGTCGATGCCGTCGGCG  >  NC_002947/6172271‑6172566

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