Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 6,172,417 C→T 8.8% W222* (TGG→TAG) ‡ 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,4170CT 8.8% 121.0 / 2.7 46W222* (TGG→TAG) ‡PP_5410DeoR family transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base C (22/19);  new base T (2/2);  total (25/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 = 4.81e-01

GCTCGGCGACAGCACCAACGATTCCACACATGTTTTATTGCTCCTGGCTGATCGCGGCGCAGATCAGGTTGATGCCGCGGGCTTGAATTTGTTCGCGCGCTGCTGCTGGCAGGCGTTCATCTGTAATAAGGGTGTGCACACTGCCCCAGGGCAGCTCGAGGTTAGGGATCTTGCGCCCCACCTTGTCCGACTCGACCATCACGATCACTTCACGGGCAACTTCCGCCATGACCCGGCTCAGCCCAAGCAACTCGTTGAAGGTAGTGGTACCCCGGCCCAGGTCGATGCCGTCG  >  NC_002947/6172271‑6172563
                                                                                                                                                  |                                                                                                                                                  
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                                                                >  3:282334/1‑149 (MQ=255)
                                                                                  ttGAATTTGTTCGCGCGCTGCTGCTGGCAGGCGTTCATCTGTAATAAGGGTGTGCACACTGCCCCAGGGCAGCTCGAGGTTAGGGATCTTGCGCCCCACCTTGTCCGACTCGACCATCACGATCACTTCACGGGCAACTTCCGCCATGa                                                                >  7:241960/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                               <  5:186219/145‑1 (MQ=255)
                                                                                                                       tCTGTAATAAGGGTGTGCACACTGCCCCAGGGCAGCTCGAGGTTAGGGATCTTGCGCCCCACCTTGTCCGACTCGACCATCACGATCACTTCACGGGCAACTTCCGCCATGACCCGGCTCAGCCCAAGCAACTCGTTGAAGGTAg                               >  6:186219/1‑145 (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)
                                                                                                                                                  |                                                                                                                                                  
GCTCGGCGACAGCACCAACGATTCCACACATGTTTTATTGCTCCTGGCTGATCGCGGCGCAGATCAGGTTGATGCCGCGGGCTTGAATTTGTTCGCGCGCTGCTGCTGGCAGGCGTTCATCTGTAATAAGGGTGTGCACACTGCCCCAGGGCAGCTCGAGGTTAGGGATCTTGCGCCCCACCTTGTCCGACTCGACCATCACGATCACTTCACGGGCAACTTCCGCCATGACCCGGCTCAGCCCAAGCAACTCGTTGAAGGTAGTGGTACCCCGGCCCAGGTCGATGCCGTCG  >  NC_002947/6172271‑6172563

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: