Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,514,591 A→G 8.1% N1087D (AAC→GAC)  PP_3108 → rhs family protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,514,5910AG 8.1% 137.2 / 2.8 49N1087D (AAC→GAC) PP_3108rhs family protein
Reads supporting (aligned to +/- strand):  ref base A (23/22);  new base G (2/2);  total (25/24)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.86e-01

TTCAATGCGCGGCAATTCGACAAGGTCAAGGACAACCGGATCAAGCAGTGGCGGGACCAGGAGTACCGCTACGATCCGTGGGGCAACCTGATCGAGAAGCGCTCTGGGCACAGCAAGCTGCAGCACTTCAGCTACGACTGCGAGAACCGGCTGGTGCGGGCGGAGACGTTGGTCAATGGCAAGCTGGAAAGCCAAGGCGAGTATCGTTACGACAGCTTGGGGAGGCGAATAGCCAAGCAGGCTGAAATCAATGGCGAGGTCGAGCAGAAGCGCTTCCTCTGGCAAGGGCT  >  NC_002947/3514447‑3514736
                                                                                                                                                |                                                                                                                                                 
ttCAATGCGCGGCAATTCGACAAGGTCAAGGACAACCGGATCAAGCAGTGGCGGGACCAGGAGTACCGCTACGATCCGTGGGGCAACCTGATCGAGAAGCGCTCTGGGCACAGCAAGCTGCAGCACTTCAGCTACGACTGCGAGAACCg                                                                                                                                               >  6:276519/1‑149 (MQ=255)
   aaTGCGCGGCAATTCGACAAGGTCAAGGACAACCGGATCAAGCAGTGGCGGGACCAGGAGTACCGCTACGATCCGTGGGGCAACCTGATCGAGAAGCGCTCTGGGCACAGCAAGCTGCAGCACTTCAGCTACGACTGCGAGAACCGGc                                                                                                                                             <  2:9890/148‑1 (MQ=255)
         cgGCAATTCGACAAGGTCAAGGACAACCGGATCAAGCAGTGGCGGGACCAGGAGTACCGCTACGATCCGTGGGGCAACCTGATCGAGAAGCGCTCTGGGCACAGCAAGCTGCAGCACTTCAGCTACGACTGCGAGAACCGGCTGGTgcg                                                                                                                                      >  1:60022/1‑149 (MQ=255)
         cgGCAATTCGACAAGGTCAAGGACAACCGGATCAAGCAGTGGCGGGACCAGGAGTACCGCTACGATCCGTGGGGCAACCTGATCGAGAAGCGCTCTGGGCACAGCAAGCTGCAGCACTTCAGCTACGACTGCGAGAACCGGCTGGTgcg                                                                                                                                      <  8:209064/149‑1 (MQ=255)
          ggCAATTCGACAAGGTCAAGGACAACCGGATCAAGCAGTGGCGGGACCAGGAGTACCGCTACGATCCGTGGGGCAACCTGATCGAGAAGCGCTCTGGGCACAGCAAGCTGCAGCACTTCAGCTACGACTGCGAGAACCGGCTGGTgcgg                                                                                                                                     >  3:2363/1‑149 (MQ=255)
            cAATTCGACAAGGTCAAGGACAACCGGATCAAGCAGTGGCGGGACCAGGAGTACCGCTACGATCCGTGGGGCAACCTGATCGAGAAGCGCTCTGGGCACAGCAAGCTGCAGCACTTCAGCTACGACTGCGAGAACCGGCTGGTgcgggc                                                                                                                                   <  4:2363/149‑1 (MQ=255)
            cAATTCGACAAGGTCAAGGACAACCGGATCAAGCAGTGGCGGGACCAGGAGTACCGCTACGATCCGTGGGGCAACCTGATCGAGAAGCGCTCTGGGCACAGCAAGCTGCAGCACTTCAGCTACGACTGCGAGAACCGGCTGGTgcgggc                                                                                                                                   >  1:10589/1‑149 (MQ=255)
                  gACAAGGTCAAGGACAACCGGATCAAGCAGTGGCGGGACCAGGAGTACCGCTACGATCCGTGGGGCAACCTGATCGAGAAGCGCTCTGGGCACAGCAATCTGCAGCACTTCAGCTACGACTGCGAGAACCGGCTGGTGCGGGCGGAGAc                                                                                                                             >  8:50488/1‑149 (MQ=255)
                      aGGTCAAGGACAACCGGATCAAGCAGTGGCGGGACCAGGAGTACCGCTACGATCCGTGGGGCAACCTGATCGAGAAGCGCTCTGGGCACAGCAAGCTGCAGCACTTCAGCTACGACTGCGAGAACCGGCTGGTGCGGGCGGAGACGTTg                                                                                                                         <  2:64710/149‑1 (MQ=255)
                             ggACAACCGGATCAAGCAGTGGCGGGACCAGGAGTACCGCTACGATCCGTGGGGCAACCTGATCGAGAAGCGCTCTGGGCACAGCAAGCTGCAGCACTTCAGCTACGACTGCGAGAACCGGCTGGTgcggg                                                                                                                                    >  1:64662/1‑131 (MQ=255)
                             ggACAACCGGATCAAGCAGTGGCGGGACCAGGAGTACCGCTACGATCCGTGGGGCAACCTGATCGAGAAGCGCTCTGGGCACAGCAAGCTGCAGCACTTCAGCTACGACTGCGAGAACCGGCTGGTgcggg                                                                                                                                    <  2:64662/131‑1 (MQ=255)
                              gACAACCGGATCAAGCAGTGGCGGGACCAGGAGTACCGCTACGATCCGTGGGGCAACCTGATCGAGAAGCGCTCTGGGCACAGCAAGCTGCAGCACTTCAGCTACGACTGCGAGAACCGGCTGGTGCGGGCGGAGACGTTGGTCAATgg                                                                                                                 <  5:276519/149‑1 (MQ=255)
                                 aaCCGGATCAAGCAGTGGCGGGACCAGGAGTACCGCTACGATCCGTGGGGCAACCTGATCGAGAAGCGCTCTGGGCACAGCAAGCTGCAGCACTTCAGCTACGACTGCGAGAACCGGCTGGTGCGGGCGGAGACGTTGGTCAATGGc                                                                                                                >  7:128268/1‑147 (MQ=255)
                                                tGGCGGGACCAGGAGTACCGCTACGATCCGTGGGGCAACCTGATCGAGAAGCGCTCTGGGCACAGCAAGCTGCAGCACTTCAGCTACGACTGCGAGAACCGGCTGGTGCGGGCGGAGACGTTGGTCAATGGCAAGCTGGAAAGCCAAgg                                                                                               <  8:128268/149‑1 (MQ=255)
                                                  gCGGGACCAGGAGTACCGCTACGATCCGTGGGGCAACCTGATCGAGAAGCGCTCTGGGCACAGCAAGCTGCAGCACTTCAGCTACGACTGCGAGGACCGGCTGGTGCGGGCGGAGACGTTGGTCAATGGCAAGCTGGAAAGCCAAGGc                                                                                              >  6:260237/1‑148 (MQ=255)
                                                  gCGGGACCAGGAGTACCGCTACGATCCGTGGGGCAACCTGATCGAGAAGCGCTCTGGGCACAGCAAGCTGCAGCACTTCAGCTACGACTGCGAGGACCGGCTGGTGCGGGCGGAGACGTTGGTCAATGGCAAGCTGGAAAGCCAAGGc                                                                                              >  6:325703/1‑148 (MQ=255)
                                                    gggACCAGGAGTACCGCTACGATCCGTGGGGCAACCTGATCGAGAAGCGCTCTGGGCACAGCAAGCTGCAGCACTTCAGCTACGACTGCGAGAACCGGCTGGTGCGGGCGGAGACGTTGGTCAATGGCAAGCTGGAAAGCCAAGGCGAg                                                                                           >  8:84935/1‑149 (MQ=255)
                                                         cAGGAGTACCGCTACGATCCGTGGGGCAACCTGATCGAGAAGCGCTCTGGGCACAGCAAGCTGCAGCACTTCAGCTACGACTGCGAGAACCGGCTGGTGCGGGCGGAGACGTTGGTCAATGGCAAGCTGGAAAGCCAAGGCGAGTATCg                                                                                      <  2:60022/149‑1 (MQ=255)
                                                            gAGTACCGCTACGATCCGTGGGGCAACCTGATCGAGAAGCGCTCTGGGCACAGCAAGCTGCAGCACTTCAGCTACGACTGCGAGAACCGGCTGGTGCGGGCGGAGACGTTCGTCAATGGCAAGCTGGAAAGCCAAGGCGAGTATCGTTa                                                                                   >  5:73717/1‑149 (MQ=255)
                                                             aGTACCGCTACGATCCGTGGGGCAACCTGATCGAGAAGCGCTCTGGGCACAGCAAGCTGCAGCACTTCAGCTACGACTGCGAGAACCGGCTGGTGCGGGCGGAGACGTTGGTCAATGGCAAGCTGGAAAGCCAAGGCGAGTATCGTTAc                                                                                  <  7:84935/149‑1 (MQ=255)
                                                                 ccGCTACGATCCGTGGGGCAACCTGATCGAGAAGCGCTCTGGGCACAGCAAGCTGCAGCACTTCAGCTACGACTGCGAGAACCGGCTGGTGCGGGCGGAGACGTTg                                                                                                                         <  5:48795/106‑1 (MQ=255)
                                                                 ccGCTACGATCCGTGGGGCAACCTGATCGAGAAGCGCTCTGGGCACAGCAAGCTGCAGCACTTCAGCTACGACTGCGAGAACCGGCTGGTGCGGGCGGAGACGTTg                                                                                                                         >  6:48795/1‑106 (MQ=255)
                                                                        gATCCGTGGGGCAACCTGATCGAGAAGCGCTCTGGGCACAGCAAGCTGCAGCACTTCAGCTACGACTGCGAGAACCGGCTGGTGCGGGCGGAGACGTTGGTCAATGGCAAGCTGGAAAGCCAAGGCGAGTATCGTTACGACAGCTTggg                                                                       <  2:10589/149‑1 (MQ=255)
                                                                            cGTGGGGCAACCTGATCGAGAAGCGCTCTGGGCACAGCAAGCTGCAGCACTTCAGCTACGACTGCGAGGACCGGCTGGTGCGGGCGGAGACGTTGGTCAATGGCAAGCTGGAAAGCCAAGGCGAGTATCGTTACGACAGCTTGGGGAgg                                                                   <  5:260237/149‑1 (MQ=255)
                                                                            cGTGGGGCAACCTGATCGAGAAGCGCTCTGGGCACAGCAAGCTGCAGCACTTCAGCTACGACTGCGAGGACCGGCTGGTGCGGGCGGAGACGTTGGTCAATGGCAAGCTGGAAAGCCAAGGCGAGTATCGTTACGACAGCTTGGGGAgg                                                                   <  5:325703/149‑1 (MQ=255)
                                                                            cGTGGGGCAACCTGATCGAGAAGCGCTCTGGGCACAGCAAGCTGCAGCACTTCAGCTACGACTGCGAGAACCGGCTGGTGCGGGCGGAGACGTTGGTCAATGGCAAGCTGGAAAGCCAAGGCGAGTATCGTTACGACAGCTTGGGGAgg                                                                   <  6:73717/149‑1 (MQ=255)
                                                                             gTGGGGCAACCTGATCGAGAAGCGCTCTGGGCACAGCAAGCTGCAGCACTTCAGCTACGACTGCGAGAACCGGCTGGTGCGGGCGGAGACGTTGGTCAATGGCAAGCTGGAAAGCCAAGGCGAGTATCGTTACGACAGCTTGGGGAGGc                                                                  >  8:263770/1‑149 (MQ=255)
                                                                                 ggCAACCTGATCGAGAAGCGCTCTGGGCACAGCAAGCTGCAGCACTTCAGCTACGACTGCGAGAACCGGCTGGTGCGGGCGGAGACGTTGGTCAATGGCAAGCTGGAAAGCCAAGGCGAGTATCGTTACGACAGCTTGGGGAGGCGa                                                                <  5:171461/147‑1 (MQ=255)
                                                                                 ggCAACCTGATCGAGAAGCGCTCTGGGCACAGCAAGCTGCAGCACTTCAGCTACGACTGCGAGAACCGGCTGGTGCGGGCGGAGACGTTGGTCAATGGCAAGCTGGAAAGCCAAGGCGAGTATCGTTACGACAGCTTGGGGAGGCGa                                                                >  6:171461/1‑147 (MQ=255)
                                                                                  gCAACCTGATCGAGAAGCGCTCTGGGCACAGCAAGCTGCAGCACTTCAGCTACGACTGCGAGAACCGGCTGGTGCGGGCGGAGACGTTGGTCAATGGCAAGCTGGAAAGCCAAGGCGAGTATCGTTACGACAGCTTGGGGAGGCTAATa                                                             <  7:263770/149‑1 (MQ=255)
                                                                                        tGATCGAGAAGCGCTCTGGGCACAGCAAGCTGCAGCACTTCAGCTACGACTGCGAGAACCGGCTGGTGCGGGCGGAGACGTTGGTCAATgg                                                                                                                 <  6:57421/91‑1 (MQ=255)
                                                                                        tGATCGAGAAGCGCTCTGGGCACAGCAAGCTGCAGCACTTCAGCTACGACTGCGAGAACCGGCTGGTGCGGGCGGAGACGTTGGTCAATgg                                                                                                                 >  5:57421/1‑91 (MQ=255)
                                                                                        tGATCGAGAAGCGCTCTGGGCACAGCAAGCTGCAGCACTTCAGCTACGACTGCGAGAACCGACTGGTGCGGGCGGAGACGTTGGTCAATGGCAAGCTGGAAAGCCAAGGCGAGTATCGTTACGACAGCTTGGGGAGGCGAATAGCCAAg                                                       >  3:45545/1‑149 (MQ=255)
                                                                                         gATCGAGAAGCGCTCTGGGCACAGCAAGCTGCAGCACTTCAGCTACGACTGCGAGAACCGGCTGGTGCGGGCGGAGACGTTGGTCAATGGCAAGCTGGAAAGCCAAGGCGAGTATCGTTACGACAGCTTGGGGAGGCGAATAGCCAAg                                                       >  8:28451/1‑148 (MQ=255)
                                                                                          aTCGAGAAGCGCTCTGGGCACAGCAAGCTGCAGCACTTCAGCTACGACTGCGAGAACCGACTGGTGCGGGCGGAGACGTTGGTCAATGGCAAGCTGGAAAGCCAAGGCGAGTATCGTTACGACAGCTTGGGGAGGCGAATAGCCAAg                                                       <  4:45545/147‑1 (MQ=255)
                                                                                            cGAGAAGCGCTCTGGGCACAGCAAGCTGCAGCACTTCAGCTACGACTGCGAGAACCGGCTGGTGCGGGCGGAGACGTTGGTCAATGGCAAGCTGGAAAGCCAAGGCGAGTATCGTTACGACAGCTTGGGGAGGCGAATAGCCAAGCAgg                                                   >  4:157042/1‑149 (MQ=255)
                                                                                                aaGCGCTCTGGGCACAGCAAGCTGCAGCACTTCAGCTACGACTGCGAGAACCGGCTGGTGCGGGCGGAGACGTTGGTCAATGGCAAGCTGGAAAGCCAAGGCGAGTATCGTTACGACAGCTTGGGGAGGCGAATAGCCAAGCAGGCTg                                                >  8:303729/1‑148 (MQ=255)
                                                                                                  gcgcTCTGGGCACAGCAAGCTGCAGCACTTCAGCTACGACTGCGAGAACCGGCTGGTGCGGGCGGAGACGTTGGTCAATGGCAAGCTGGAAAGCCAAGGCGAGTATCGTTACGACAGCTTGGGGAGGCGAATAGCCAAGCAGGCTg                                                >  6:313437/1‑146 (MQ=255)
                                                                                                  gcgcTCTGGGCACAGCAAGCTGCAGCACTTCAGCTACGACTGCGAGAACCGGCTGGTGCGGGCGGAGACGTTGGTCAATGGCAAGCTGGAAAGCCAAGGCGAGTATCGTTACGACAGCTTGGGGAGGCGAATAGCCAAGCAGGCTg                                                <  5:313437/146‑1 (MQ=255)
                                                                                                     ctctGGGCACAGCAAGCTGCAGCACTTCAGCTACGACTGCGAGAACCGGCTGGTGCGGGCGGAGACGTTGGTCAATGGCAAGCTGGAAAGCCAAGGCGAGTATCGTTACGACAGCTTGGGGAGGCGAATAGCCAAGCAGGCTGAAATc                                           >  6:314619/1‑148 (MQ=255)
                                                                                                               aGCAAGCTGCAGCACTTCAGCTACGACTGCGAGAACCGGCTGGTGCGGGCGGAGACGTTGGTCAATGGCAAGCTGGAAAGCCAAGGCGAGTATCGTTACGACAGCTTGGGGAGGCGAATAGCCAAGCAGGCTGAAATCAATGGCGAGGt                                >  1:294041/1‑149 (MQ=255)
                                                                                                                    ggtGCATCACTTCAGCTACGACTGGGAGAACCGGCTGGTGCGGGCGGAGACGTTGGTCAATGGCAAGCTGGAATGCCAAGGCGAGTATCGTTACGACAGCTTGGGGAGGCGAATAGCCAAGCAGGCTGAAATCAATGGCGAGGTCGAGc                           <  7:50488/147‑1 (MQ=255)
                                                                                                                    gCTGCAGCACTTCAGCTACGACTGCGAGAACCGGCTGGTGCGGGCGGAGACGTTGGTCAATGGCAAGCTGGAAAGCCAAGGCGAGTATCGTTACGACAGCTTGGGGAGGCGAATAGCCAAGCAGGCTGAAATCAATGGCGAGGTCGAGc                           <  7:303729/149‑1 (MQ=255)
                                                                                                                      tGCAGCACTTCAGCTACGACTGCGAGAACCGGCTGGTGCGGGCGGAGACGTTGGTCAATGGCAAGCTGGAAAGCCAAgg                                                                                               <  4:248877/79‑1 (MQ=255)
                                                                                                                      tGCAGCACTTCAGCTACGACTGCGAGAACCGGCTGGTGCGGGCGGAGACGTTGGTCAATGGCAAGCTGGAAAGCCAAgg                                                                                               <  4:253153/79‑1 (MQ=255)
                                                                                                                      tGCAGCACTTCAGCTACGACTGCGAGAACCGGCTGGTGCGGGCGGAGACGTTGGTCAATGGCAAGCTGGAAAGCCAAgg                                                                                               >  3:253153/1‑79 (MQ=255)
                                                                                                                      tGCAGCACTTCAGCTACGACTGCGAGAACCGGCTGGTGCGGGCGGAGACGTTGGTCAATGGCAAGCTGGAAAGCCAAgg                                                                                               >  3:248877/1‑79 (MQ=255)
                                                                                                                                       gACTGCGAGAACCGGCTGGTGCGGGCGGAGACGTTGGTCAATGGCAAGCTGGAAAGCCAAGGCGAGTATCGTTACGACAGCTTGGGGAGGCGAATAGCCaa                                                        >  7:165188/1‑101 (MQ=255)
                                                                                                                                       gACTGCGAGAACCGGCTGGTGCGGGCGGAGACGTTGGTCAATGGCAAGCTGGAAAGCCAAGGCGAGTATCGTTACGACAGCTTGGGGAGGCGAATAGCCaa                                                        <  8:165188/101‑1 (MQ=255)
                                                                                                                                               gaACCGGCTGGTGCGGGCGGAGACGTTGGTCAATGGCAAGCTGGAAAGCCAAGGCGAGTATCGTTACGACAGCTTGGGGAGGCGAATAGCCAAGCAGGCTGAAATCAATGGCGAGGTCGAGCAGAAGCGCTTCCTCTGGCAAGGGCt  >  1:202711/1‑147 (MQ=255)
                                                                                                                                                |                                                                                                                                                 
TTCAATGCGCGGCAATTCGACAAGGTCAAGGACAACCGGATCAAGCAGTGGCGGGACCAGGAGTACCGCTACGATCCGTGGGGCAACCTGATCGAGAAGCGCTCTGGGCACAGCAAGCTGCAGCACTTCAGCTACGACTGCGAGAACCGGCTGGTGCGGGCGGAGACGTTGGTCAATGGCAAGCTGGAAAGCCAAGGCGAGTATCGTTACGACAGCTTGGGGAGGCGAATAGCCAAGCAGGCTGAAATCAATGGCGAGGTCGAGCAGAAGCGCTTCCTCTGGCAAGGGCT  >  NC_002947/3514447‑3514736

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