Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,493,781 C→T 8.7% R286H (CGC→CAC)  PP_2187 ← universal stress protein family protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,493,7810CT 8.7% 126.5 / 2.4 46R286H (CGC→CAC) PP_2187universal stress protein family protein
Reads supporting (aligned to +/- strand):  ref base C (18/24);  new base T (2/2);  total (20/26)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AAGGTGCGGGAGGCGTTGGTCTAGACTCCCGGGCGGCAACGGAACCTTGGGGCGGGTTTGCCGGCGAATGTGGCGTTGAAATCAATGACGCATTCGCGGGTGAACCCGCTCCTACAAGGGGCCTGCATATTGCCTACTCAGCCGCGGGCCAGCTCCGCCAGGTGGGCGATGGCTTCCTCACTCTTGAGCACCAGCACATCACCCTGCAGCGCATCCAGTACCACTTCCGAGGTATTACCAATCAACGCCCCGGAAATGCCCGTGCGGCACACCGTGCCCATCACCGTC  >  NC_002947/2493637‑2493924
                                                                                                                                                |                                                                                                                                               
aaGGTGCGGGAGGCGTTGGTCTAGACTCCCGGGCGGCAACGGAACCTTGGGGCGGGTTTGCCGGCGAATGTGGCGTTGAAATCAATGACGCATTCGCGGGTGAACCCGCTCCTACAAGGGGCCTGCATATTGCCTACTCAGCCGCGGGc                                                                                                                                             <  4:120625/149‑1 (MQ=255)
 aGGTGCGGGAGGCGTTGGTCTAGACTCCCGGGCGGCAACGGAACTTTGGGGCGGGTTTGCCGGCGAATGTGGCGTTGAAATCAATGACGCATTCGCGGGTGAACCCGCTCCTACAAGGGGCCTGCATATTGCCTACTCAGCCGCGGGcc                                                                                                                                            <  2:77047/149‑1 (MQ=255)
            gCGTTGGTCTAGACTCCCGGGCGGCAACGGAACCTTGGGGCGGGTTTGCCGGCGAATGTGGCGTTGAAATCAATGACGCATTCGCGGGTGAACCCGCTCCTACAAGGGGCCTGCATATTGCCTACTCAGCCGCGGGCCAGCTCCGCCAg                                                                                                                                 <  3:131075/149‑1 (MQ=255)
             cGTTGGTCTAGACTCCCGGGCGGCAACGGAACCTTGGGGCGGGTTTGCCGGCGAATGTGGCGTTGAAATCAATGACGCATTCGCGGGTGAACCCGCTCCTACAAGGGGCCTGCATATTGCCTACTCAGCCGCGGGCCAGCTCCGCCAgg                                                                                                                                <  7:100243/149‑1 (MQ=255)
                     tAGACTCCCGGGCGGCAACGGAACCTTGGGGCGGGTTTGCCGGCGAATGTGGCGTTGAAATCAATGACGCATTCGCGGGTGAACCCGCTCCTACAAGGGGCCTGCATATTGCCTACTCAGCCGCGGGCCAGCTCCGCCAGGTGGGCGAt                                                                                                                        >  3:192167/1‑149 (MQ=255)
                      aGACTCCCGGGCGGCAACGGAACCTTGGGGCGGGTTTGCCGGCGAATGTGGCGTTGAAATCAATGACGCATTCGCGGGTGAACCCGCTCCTACAAGGGGCCTGCATATTGCCTACTCAGCCGCGGGCCAGCTCCGCCAGGTgg                                                                                                                             >  7:243571/1‑143 (MQ=255)
                      aGACTCCCGGGCGGCAACGGAACCTTGGGGCGGGTTTGCCGGCGAATGTGGCGTTGAAATCAATGACGCATTCGCGGGTGAACCCGCTCCTACAAGGGGCCTGCATATTGCCTACTCAGCCGCGGGCCAGCTCCGCCAGGTgg                                                                                                                             <  8:243571/143‑1 (MQ=255)
                       gACTCCCGGGCGGCAACGGAACCTTGGGGCGGGTTTGCCGGCGAATGTGGCGTTGAAATCAATGACGCATTCGCGGGTGAACCCGCTCCTACAAGGGGCCTGCATATTGCCTACTCAGCCGCGGGCCAGCTCCGCCAGGTGGGCGa                                                                                                                         <  5:114461/146‑1 (MQ=255)
                       gACTCCCGGGCGGCAACGGAACCTTGGGGCGGGTTTGCCGGCGAATGTGGCGTTGAAATCAATGACGCATTCGCGGGTGAACCCGCTCCTACAAGGGGCCTGCATATTGCCTACTCAGCCGCGGGCCAGCTCCGCCAGGTGGGCGATgg                                                                                                                      >  5:302544/1‑149 (MQ=255)
                       gACTCCCGGGCGGCAACGGAACCTTGGGGCGGGTTTGCCGGCGAATGTGGCGTTGAAATCAATGACGCATTCGCGGGTGAACCCGCTCCTACAAGGGGCCTGCATATTGCCTACTCAGCCGCGGGCCAGCTCCGCCAGGTGGGCGATg                                                                                                                       <  6:201178/148‑1 (MQ=255)
                        aCTCCCGGGCGGCAACGGAACCTTGGGGCGGGTTTGCCGGCGAATGTGGCGTTGAAATCAATGACGCATTCGCGGGTGAACCCGCTCCTACAAGGGGCCTGCATATTGCCTACTCAGCCGCGGGCCAGCTCCGCCAGGTGGGCGATGGc                                                                                                                     <  1:70918/149‑1 (MQ=255)
                                 cggcAACGGAACCTTGGGGCGGGTTTGCCGGCGAATGTGGCGTTGAAATCAATGACGCATTCGCGGGTGAACCCGCTCCTACAAGGGGCCTGCATATTGCCTACTCAGCCGCGGGCCAGCTCCGCCAGGTGGGCGATGGCTTCCTCAct                                                                                                            <  6:102269/149‑1 (MQ=255)
                                     aaCGGAACCTTGGGGCGGGTTTGCCGGCGAATGTGGCGTTGAAATCAATGACGCATTCGCGGGTGAACCCGCTCCTACAAGGGGCCTGCATATTGCCTACTCAGCCGCGGGCCAGCTCCGcc                                                                                                                                   >  8:98096/1‑122 (MQ=255)
                                     aaCGGAACCTTGGGGCGGGTTTGCCGGCGAATGTGGCGTTGAAATCAATGACGCATTCGCGGGTGAACCCGCTCCTACAAGGGGCCTGCATATTGCCTACTCAGCCGCGGGCCAGCTCCGcc                                                                                                                                   <  7:98096/122‑1 (MQ=255)
                                         gAACCTTGGGGCGGGTTTGCCGGCGAATGTGGCGTTGAAATCAATGACGCATTCGCGGGTGAACCCGCTCCTACAAGGGGCCTGCATATTGCCTACTCAGCCGCGGGCCAGCTCCGCCAGGTGGGCGATGGCTTCCTCACTCTTGAGCa                                                                                                    >  6:260753/1‑149 (MQ=255)
                                             cTTGGGGCGGGTTTGCCGGCGAATGTGGCGTTGAAATCAATGACGCATTCGCGGGTGAACCCGCTCCTACAAGGGGCCTGCATATTGCCTACTCAGCCGCGGGCCAGCTCCGCCAGGTGGGCGATGGCTTCCTCACTCTTGAGCACCAg                                                                                                >  2:314629/1‑149 (MQ=255)
                                                  ggCAGGTTTGCCGGCGAATGTGGCGTTGAAATCAATGACGCATTCGCGGGTGAACCCGCTCCTACAAGGGGCCTGCATATTGCCTACTCAGCCGCGGGCCAGCTCCGCCa                                                                                                                                  >  6:238340/1‑110 (MQ=255)
                                                  ggCAGGTTTGCCGGCGAATGTGGCGTTGAAATCAATGACGCATTCGCGGGTGAACCCGCTCCTACAAGGGGCCTGCATATTGCCTACTCAGCCGCGGGCCAGCTCCGCCa                                                                                                                                  <  5:238340/110‑1 (MQ=255)
                                                                 gAATGTGGCGTTGAAATCAATGACGCATTCGCGGGTGAACCCGCTCCTACAAGGGGCCTGCATATTGCCTACTCAGCCGCGGGCCAGCTCCGCCAGGTGGGCGATGGCTTCCTCACTCTTGAGCACCAGCACATCACCCTTCAGCGCAt                                                                            >  2:12021/1‑149 (MQ=255)
                                                                  aaTGTGGCGTTGAAATCAATGACGCATTCGCGGGTGAACCCGCTCCTACAAGGGGCCTGCATATTGCCTACTCAGCCGCGGGCCAGCTCCGCCAGGTGGGCGATGGCTTCCTCACTCTTGAGCACCAGCACATCACCCTGCAGCGCATc                                                                           <  1:314629/149‑1 (MQ=255)
                                                                          gTTGAAATCAATGACGCATTCGCGGGTGAACCCGCTCCTACAAGGGGCCTGCATATTGCCTACTCAGCCGCGGGCCAGCTCCGCCAGGTGGGCGATGGCTTCCTCACTCTTGAGCACCAGCACATCACCCTGCAGCGCATCCAGTACCa                                                                   <  3:129958/149‑1 (MQ=255)
                                                                              aaaTCAATGACGCATTCGCGGGTGAACCCGCTCCTACAAGGGGCCTGCATATTGCCTACTCAGCCGCGGGCCAGCTCCGCCAGGTGGGCGATGGCTTCCTCACTCTTGAGCACCAGCACATCACCCTGCAGCGCATCCAGTACCACTTc                                                               >  6:231150/1‑149 (MQ=255)
                                                                                 tCAATGACGCATTCGCGGGTGAACCCGCTCCTACAAGGGGCCTGCATATTGCCTACTCAGCCGTGGGCCAGCTCCGCCAGGTGGGCGATGGCTTCCTCACTCTTGAGCACCAGCACATCACACAGCAGCGCATCCAGTACCACTTcc                                                              >  8:270795/1‑147 (MQ=255)
                                                                                   aaTGACGCATTCGCGGGTGAACCCGCTCCTACAAGGGGCCTGCATATTGCCTACTCAGCCGCGGGCCAGCTCCGCCAGGTGGGCGATGGCTTCCTCACTCTTGAGCACCAGCACATCACCCTGCAGCGCATCCAGTACCACTTCCGAg                                                           <  2:268403/148‑1 (MQ=255)
                                                                                    aTGACGCATTCGCGGGTGAACCCGCTCCTACAAGGGGCCTGCATATTGCCTACTCAGCCGCGGGCCAGCTCCGCCAGGTGGGCGATGGCTTCCTCACTCTTGAGCACCAGCACATCACCCTGCAGCGCATCCAGTACCACTTCCGAGGt                                                         <  5:260753/149‑1 (MQ=255)
                                                                                     tGACGCATTCGCGGGTGAACCCGCTCCTACAAGGGGCCTGCATATTGCCTACTCAGCCGCGGGCCAGCTCCGCCAGGTGGGCGATGGCTTCCTCACTCTTGAGCACCAGCACATCACCCTGCAGCGCATCCAGTACCACTTCCGAGGTa                                                        >  6:137080/1‑149 (MQ=255)
                                                                                          taggagcgGGTTCAACCGCTCCTACAAGGGGCCTGCATATTGCCTACTCAGCCGTGGGCCAGCTCCGCCAGGTGGGCGATGGCTTCCTCACTCTTGAGCACCAGCACATCACCCTGCAGCGCATCCAGTACCACTTCCGAGGTATTAc                                                    >  8:153053/6‑148 (MQ=255)
                                                                                              agcgGGTTCACCCGCTCCTACAAGGGGCCTGCATATTGCCTACTCAGCCGTGGGCCAGCTCCGCCAGGTGGGCGATGGCTTCCTCACTCTTGAGCACCAGCACATCACCCTGCAGCGCATCCAGTACCACTTCCGAGGTATTACCAATc                                               <  7:153053/148‑1 (MQ=255)
                                                                                                cgGGTGAACCCGCTCCTACAAGGGGCCTGCATATTGCCTACTCAGCCGCGGGCCAGCTCCGCCAGGTGGGCGATGGCTTCCTCACTCTTGAGCACCAGCACATCACCCTGCAGCGCATCCAGTACCACTTCCGAGGTATTACCAATCaa                                             <  6:87042/149‑1 (MQ=255)
                                                                                                 gggTGAACCCGCTCCTACAAGGGGCCTGCATATTGCCTACTCAGCCGCGGGCCAGCTCCGCCAGGTGGGCGATGGCTTCCTCACTCTTGAGCACCAGCACATCACCCTGCAGCGCATCCAGTACCACTTCCGAGGTATTACCAATCAAc                                            >  8:114794/1‑149 (MQ=255)
                                                                                                            cTCCTACAAGGGGCCTGCATATTGCCTACTCAGCCGTGGGCCAGCTCCGCCAGGTGGGCGATGGCTTCCTCACTCTTGAGCACCAGCACATCACCCTGCAGCGCATCCAGTACCACTTCCGAGGTATTACCAATCAACGCCCCGGAAAt                                 <  7:270795/149‑1 (MQ=255)
                                                                                                                tACAAGGGGCCTGCATATTGCCTACTCAGCCGCGGGCCAGCTCCGCCAGGTGGGCGATGGCTTCCTCAct                                                                                                            <  4:153176/70‑1 (MQ=255)
                                                                                                                tACAACGGGCCTGCATATTGCCTACTCAGCCGCGGGCCAGCTCCGCCAGGTGGGCGATGGCTTCCTCAct                                                                                                            >  3:153176/1‑70 (MQ=255)
                                                                                                                  cAAGGGGCCCGCATATTGCCTACTCAGCCGCGGGCCAGCTCCGCCAGGTGGGCGATGGCTTCCTCACTCTTGAGCACCAGCACATCACCCTGCAGCGCATCCAGTACCACTTCCGAGGTATTACCAATCAACGCCCCGGAAATGCCCGt                           >  2:128490/1‑149 (MQ=255)
                                                                                                                  cAAGGGGCCCGCATATTGCCTACTCAGCCGCGGGCCAGCTCCGCCAGGTGGGCGATGGCTTCCTCACTCTTGAGCACCAGCACATCACCCTGCAGCGCATCCAGTACCACTTCCGAGGTATTACCAATCAACGCCCCGGAAATGCCCGt                           >  4:60180/1‑149 (MQ=255)
                                                                                                                     ggggCCTGCATATTGCCTACTCAGCCGCGGGCCAGCTCCGCCAGGTGGGCGATGGCTTCCTCACTCTTGAGCACCAGCACATCACCCTGCAGCGCATCCAGTACCACTTCCGAGGTATTACCAATCAACGc                                          <  5:292166/131‑1 (MQ=255)
                                                                                                                     ggggCCTGCATATTGCCTACTCAGCCGCGGGCCAGCTCCGCCAGGTGGGCGATGGCTTCCTCACTCTTGAGCACCAGCACATCACCCTGCAGCGCATCCAGTACCACTTCCGAGGTATTACCAATCAACGc                                          >  6:292166/1‑131 (MQ=255)
                                                                                                                     ggggCCTGCATATTGCCTACTCAGCCGCGGGCCAGCTCCGCCAGGTGGGCGATGGCTTCCTCACTCTTGAGCACCAGCACATCACCCTGCAGCGCATCCAGTACCACTTCCGAGGTATTACCAATCAACGCCCCGGAAATGACCGTGCg                        >  6:87154/1‑149 (MQ=255)
                                                                                                                         ccTGCATATTGCCTACTCAGCCGCGGGCCAGCTCCGCCAGGTGGGCGATGGCTTCCTCACTCTTGAGCACCAGCACATCACCCTGCAGCGCATCCAGTACCACTTTCGAGGTATTACCAATCAACGCCCCGGAAATGCCCGTg                          <  1:215873/143‑1 (MQ=255)
                                                                                                                         ccTGCATATTGCCTACTCAGCCGCGGGCCAGCTCCGCCAGGTGGGCGATGGCTTCCTCACTCTTGAGCACCAGCACATCACCCTGCAGCGCATCCAGTACCACTTCCGAGGTATTACCAATCAACGCCCCGGAAATGCCCGTg                          >  2:215873/1‑143 (MQ=255)
                                                                                                                               tatTGCCTACTCAGCCGCGGGCCAGCTCCGCCAGGTGGGCGATGGCTTCCTCACTCTTGAGCACCAGCACATCACCCTGCAGCGCATCCAGTACCACTTCCGAGGTATTACCAATCAACGCCCCGGAAATGCCCGTGCGGCACACCGTg              <  6:302544/149‑1 (MQ=255)
                                                                                                                                 ttGCCTACTCAGCCGCGGGCCAGCTCCGCCAGGTGGGCGATGGCTTCCTCACTCTTGAGCACCAGCACATCACCCTGCAGCGCATCCAGTACCACTTCCGAGGTATTACCAATCAACGCCCCGGAAATGCCCGTGCGGCACACCGTGc             <  1:128490/148‑1 (MQ=255)
                                                                                                                                 ttGCCTACTCAGCCGCGGGCCAGCTCCGCCAGGTGGGCGATGGCTTCCTCACTCTTGAGCACCAGCACATCACCCTGCAGCGCATCCAGTACCACTTCCGAGGTATTACCAATCAACGCCCCGGAAATGCCCGTGCGGCACACCGTGc             <  5:231150/148‑1 (MQ=255)
                                                                                                                                 ttGCCTACTCAGCCGCGGGCCAGCTCCGCCAGGTGGGCGATGGCTTCCTCACTCTTGAGCACCAGCACATCACCCTGCAGCGCATCCAGTACCACTTCCGAGGTATTACCAATCAACGCCCCGGAAATGCCCGTGCGGCACACCGTGc             <  3:60180/148‑1 (MQ=255)
                                                                                                                                        cTCAGCCGCGGGCCAGCTCCGCCAGGTGGGCGATGGCTTCCTCACTCTTGAGCACCAGCACATAACCCTGCAGCGCATCCAGTACCACTTCCGAGGTATTACCAATCAACGCCCCGGAAATGCCCGTGCGGCACACCGTGCCCATCAcc     >  6:21749/1‑149 (MQ=255)
                                                                                                                                           aGCCGCGGGCCAGCTCCGCCAGGTGGGCGATGGCTTCCTCACTCTTGAGCACCAGCACATCACCCTGCAGCGCATCCAGTACCACTTCCGAGGTATTACCAATCAACGCCCCGGAAATGCCCGTGCGGCACACCGTGCCCATCACCGTc  <  4:192167/149‑1 (MQ=255)
                                                                                                                                                |                                                                                                                                               
AAGGTGCGGGAGGCGTTGGTCTAGACTCCCGGGCGGCAACGGAACCTTGGGGCGGGTTTGCCGGCGAATGTGGCGTTGAAATCAATGACGCATTCGCGGGTGAACCCGCTCCTACAAGGGGCCTGCATATTGCCTACTCAGCCGCGGGCCAGCTCCGCCAGGTGGGCGATGGCTTCCTCACTCTTGAGCACCAGCACATCACCCTGCAGCGCATCCAGTACCACTTCCGAGGTATTACCAATCAACGCCCCGGAAATGCCCGTGCGGCACACCGTGCCCATCACCGTC  >  NC_002947/2493637‑2493924

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