Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,532,765 A→G 7.7% S4G (AGT→GGT)  PP_2224 → hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,532,7650AG 7.7% 150.7 / 2.8 52S4G (AGT→GGT) PP_2224hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (30/18);  new base G (2/2);  total (32/20)
Fisher's exact test for biased strand distribution p-value = 6.34e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.82e-01

GTGATGCTGGTATTCGCCGCGCAGCAGCAAGTGCAGGAAGGCCTGGGCAAGTTGCTTAAAGCAGTGACTGGCAAAGTCGCCAAATGGCGCAAGGGTGAACTGACCACCTTCAAGGCCAGTGAATGACCGACAGTCACGACTCACCAGCCTCGCGGGCCGAGTTTCATCGGCAGCACCAGGCAGATGCCGTGGCTGAGGCAGAACGTTTGCTGGCCGGGCGTGAGAAGCTGCAAGGCGCCTGGTTGACTTGGGTAGCGGCTGAGCTCTATCGCCTCAGC  >  NC_002947/2532634‑2532911
                                                                                                                                   |                                                                                                                                                  
gTGATGCTGGTATTCGCCGCGCAGCAGCAAGTGCAGGAAGGCCTGGGCAAGTTGCTTAAAGCAGTGACTGGCAAAGTCGCCAAATGGCGCAAGGGTGAACTGACCACCTTCAAGGCCAGTGAATGACCGACAGTCACGACTCACCAGcc                                                                                                                                   >  8:56692/1‑149 (MQ=255)
       tGGTATTCGCCGCGCAGCAGCAAGTGCAGGAAGGCCTGGGCAAGTTGCTTAAAGCAGTGACTGGCAAAGTCGCCAAATGGCGCAAGGGTGAACTGACCACCTTCAAGGCCAGTGAATGACCGACAGTCACGACTCACCAGCCTCGCggg                                                                                                                            >  1:307463/1‑149 (MQ=255)
       tGGTATTCGCCGCGCAGCAGCAAGTGCAGGAAGGCCTGGGCAAGTTGCTTAAAGCAGTGACTGGCAAAGTCGCCAAATGGCGCAAGGGTGAACTGACAACCTTCAAGGCCAGTGAATGACCGACAGTCACGACTCACCAGCCTCGCggg                                                                                                                            >  7:270073/1‑149 (MQ=255)
         gTATTCGCCGCGCAGCAGCAAGTGCAGGAAGGCCTGGGCAAGTTGCTTAAAGCAGTGACTGGCAAAGTCGCCAAATGGCGCAAGGGTGAACTGACCACCTTCAAGGCCAGTGAATGACCGACAGTCACGACTCACCAGCCTCGCGGGcc                                                                                                                          >  6:327310/1‑149 (MQ=255)
            ttCGCCGCGCAGCAGCAAGTGCAGGAAGGCCTGGGCAAGTTGCTTAAAGCAGTGACTGGCAAAGTCGCCAAATGGCGCAAGGGTGAACTGACCACCTTCAAGGCCAGTGAATGACCGACAGTCACGACTCACCAGCCTCGCGGGCCGAg                                                                                                                       <  4:142935/149‑1 (MQ=255)
               gccgcGCAGCAGCAAGTGCAGGAAGGCCTGGGCAAGTTGCTTAAAGCAGTGACTGGCAAAGTCGCCAAATGGCGCAAGGGTGAACTGACCACCTTCAAGGCCAGTGAATGACCGACAGTCACGACTCACCAGCCTCGCGGGCCGAGttt                                                                                                                    >  8:80415/1‑149 (MQ=255)
               gccgcGCAGCAGCAAGTGCAGGAAGGCCTGGGCAAGTTGCTTAAAGCAGTGACTGGCAAAGTCGCCAAATGGCGCAAGGGTGAACTGACCACCTTCAAGGCCAGTGAATGACCGACAGTCACGACTCACCAGCCTCGCGGGCCGAGttt                                                                                                                    <  2:307463/149‑1 (MQ=255)
                 cgcgCAGCAGCAAGTGCAGGAAGGCCTGGGCAAGTTGCTTAAAGCAGTGACTGGCAAAGTCGCCAAATGGCGCAAGGGTGAACTGACCACCTTCAAGGCCAGTGAATGACCGACAGTCACGACTCACCAGCCTCGCGGGCCGAGTTTCa                                                                                                                  >  6:9432/1‑149 (MQ=255)
                              gTGCAGGAAGGCCTGGGCAAGTTGCTTAAAGCAGTGACTGGCAAAGTCGCCAAATGGCGCAAGGGTGAACTGACCACCTTCAAGGCCAGTGAATGACCGACAGTCACGACTCACCAGCCTCGCGGGCCGAGTTTCATCGGCAGCACCAg                                                                                                     <  7:80415/149‑1 (MQ=255)
                               tGCAGGAAGGCCTGGGCAAGTTGCTTAAAGCAGTGACTGGCAAAGTCGCCAAATGGCGCAAGGGTGAACTGACCACCTTCAAGGCCAGTGAATGACCGACAGTCACGACTCACCAGCCTCGCGGGCCGAGTTTCATCGGCAGCACCAgg                                                                                                    >  5:205474/1‑149 (MQ=255)
                                 cAGGAAGGCCTGGGCAAGTTGCTTAAAGCAGTGACTGGCAAAGTCGCCAAATGGCGCAAGGGTGAACTGACCACCTTCAAGGCCAGTGAATGACCGACAGTCACGACTCACCAGCCTCGCGGGCCGAGTTTCATCGGCAGCACCAgg                                                                                                    >  8:298998/1‑147 (MQ=255)
                                         ccTGGGCAAGTTGCTTAAAGCAGTGACTGGCAAAGTCGCCAAATGGCGCAAGGGTGAACTGACCACCTTCAAGGCCAGTGAATGACCGACAGTCACGACTCACCAGCCTCGCGGGCCGAGTTTCATCGGCAGCACCAGGCAGATGCCGt                                                                                          <  7:56692/149‑1 (MQ=255)
                                         ccTGGGCAAGTTGCTTAAAGCAGTGACTGGCAAAGTCGCCAAATGGCGCAAGGGTGAACTGAACACCTTCAAGGCCAGTGAATGACCGACAGTCACGACTCACCAGCCTCGCGGGCCGAGTTTCATCGGCAGCACCAGGCAGATGCCGt                                                                                          >  6:254657/1‑149 (MQ=255)
                                                      cTTAAAGCAGTGACTGGCAAAGTCGCCAAATGGCGCAAGGGTGAACTGACCACCTTCAAGGCCAGTGAATGACCGACAGTCACGACTCACCAGCCTCGCGGGCCGAGTTTCATCGGCAGCACCAGGCAGATGCCGTGGCTGAGGCAGaa                                                                             >  7:73124/1‑149 (MQ=255)
                                                       ttAAAGCAGTGACTGGCAAAGTCGCCAAATGGCGCAAGGGTGAACTGACCACCTTCAAGGCCAGTGAATGACCGACAGTCACGACTCACCAGCCTCGCGGGCCGAGTTTCATCGGCAGCACCAGGCAGATGCCGTGGCTGAGGCAGAAc                                                                            >  7:177531/1‑149 (MQ=255)
                                                          aaGCAGTGACTGGCAAAGTCGCCAAATGGCGCAAGGGTGAACTGACCACCTTCAAGGCCAGTGAATGACCGACAGTCACGACTCACCAGCCTCGCGGGCCGAGTTTCATCGGCAGCACCAGGCAGATGCCGTGGCTGAGGCAGAACGtt                                                                         <  7:298998/149‑1 (MQ=255)
                                                                tGACTGGCAAAGTCGCCAAATGGCGCAAGGGTGAACTGCCCACCTTCAAGGCCAGTGAATGACCGACAGTCACGACTCACCAGCCTCGCGGGCCGAGTTTCATCGGCAGCACCAGGCAGATGCCGTGGCTGAGGCAGAACGTTTGCTgg                                                                   >  7:242557/1‑149 (MQ=255)
                                                                tGACTGGCAAAGTCGCCAAATGGCGCAAGGGTGAACTGACCACCTTCAAGGCCAGTGAATGACCGACAGTCACGACTCACCAGCCTCGCGGGCCGAGTTTCATCGGCAGCACCAGGCAGATGCCGTGGCTGAGGCAGAACGTTTGCTgg                                                                   >  3:71286/1‑149 (MQ=255)
                                                                tGACTGGCAAAGTCGCCAAATGGCGCAAGGGTGAACTGACCACCTTCAAGGCCAGTGAATGACCGACAGTCACGACTCACCAGCCTCGCGGGCCGAGTTTCATCGGCAGCACCAGGCAGATGCCGTGGCTGAGGCAGAACGTTTGCTgg                                                                   >  5:29479/1‑149 (MQ=255)
                                                                        ccttgcGCCATTTGGCGCAAGGGTGAACTGACCACCTTCAAGGCCAGTGAATGACCGACAGTCACGACTCACCAGCCTCGCGGGCCGAGTTTCATCGGCAGCACCAGGCAGATGCCGTGGCTGAGGCAGAACGTTTGCTGGCCGGGCGt                                                           >  2:137493/6‑149 (MQ=255)
                                                                                cAAATGGCGCAAGGGTGAACTGACCACCTTCAAGGCCAGTGAATGACCGACAGTCACGACTCACCAGCCTCGCGGGCCGAGTTTCATCGGCAGCACCAGGCAGATGCCGt                                                                                          <  7:1322/110‑1 (MQ=255)
                                                                                cAAATGGCGCAAGGGTGAACTGACCACCTTCAAGGCCAGTGAATGACCGACAGTAACGACTCACCAGCCTCGCGGGCCGAGTTTCATCGGCAGCACCAGGCAGATGCCGt                                                                                          >  8:1322/1‑110 (MQ=255)
                                                                                 aaaTGGCGCAAGGGTGAACTGACCACCTTCAAGGCCAGTGAATGACCGACAGTCACGACTc                                                                                                                                          >  3:274780/1‑61 (MQ=255)
                                                                                 aaaTGGCGCAAGGGTGAACTGACCACCTTCAAGGCCAGTGAATGACCGACAGTCACGACTc                                                                                                                                          <  4:274780/61‑1 (MQ=255)
                                                                                               tGAACTGACCACCTTCAAGGCCAGTGAATGACCGACAGTCACGACTCACCAGCCTCGCGGGCCGAGTTTCATCGGCAGCACCAGGCAGATGCCGTGGCTGAGGCAGAACGTTTGCTGGCCGGGCGTGAGAAGCTGCAAGGCGCCTGGtt                                    <  8:177531/149‑1 (MQ=255)
                                                                                                      accaccTTCAAGGCCAGTGAATGACCGACAGTCAAGACTCACCAGCCTCGCGGGCCGAGTTTCATCGGCAGCACCAGGCAGATGCCGTGGCTGAGGCAGAACGTTTGCTGGACGGGCGTGAGAAGCTGCAAGGCGCCTGGTTGACTTgg                             >  5:47917/1‑149 (MQ=255)
                                                                                                        caccTTCAAGGCCAGTGAATGACCGACAGTCACGACTCACCAGCCTCGCGGGCCGAGTTTCATCGGCAGCACCAGGCAGATGCCGTGGCTGAGGCAGAACGTTTGCTGGCCGGGCGTGAGAAGCTGCAAGGCGCCTGGTTGACTTggg                            >  8:175896/1‑148 (MQ=255)
                                                                                                         accTTCAAGGCCAGTGAATGACCGACAGTCACGACTCACCAGCCTCGCGGGCCGAGTTTCATCGGCAGCACCAGGCAGATGCCGTGGCTGAGGCAGAACGTTTGCTGGCCGGGCGTGAGAAGCTGCaa                                               >  5:311904/1‑128 (MQ=255)
                                                                                                         accTTCAAGGCCAGTGAATGACCGACAGTCACGACTCACCAGCCTCGCGGGCCGAGTTTCATCGGCAGCACCAGGCAGATGCCGTGGCTGAGGCAGAACGTTTGCTGGCCGGGCGTGAGAAGCTGCaa                                               <  6:311904/128‑1 (MQ=255)
                                                                                                          ccTTCTAGGCCAGTGAATGTGCGACAGTCACGACTCTCCAGGCTCGCGGGCCGTGTTTCATCGGCAGCACCAGGCAGATGCCGTGGCTGAGGCAGAACGTTTGCTGGCCGGGCGTTAGAAGCTGCAAg                                              <  2:111971/128‑1 (MQ=255)
                                                                                                          ccTTCAAGGCCAGTGAATGACCGACAGTCACGACTCACCAGCCTCGCGGGCCGAGTTTCATCGGCAGCACCAGGCAGATGCCGTGGCTGAGGCAGAACGTTTGCTGGCCGGGCGTGAGAAGCTGCAAg                                              >  1:111971/1‑128 (MQ=255)
                                                                                                           cTTCAAGGCCAGTGAATGACCGACAGTCACGACTCACCAGCCTCGCGGGCCGAGTTTCATCGGCAGCACCAGGCAGATGCCGTGGCTGAGGCAGAACGTTTGCTGGCCGGGCGTGAGAAGCTGCAAGGCGCCTGGTTGACTTGGGTAGc                        <  8:73124/149‑1 (MQ=255)
                                                                                                           cTTCAAGGCCAGTGAATGACCGACAGTCACGACTCACCAGCCTCGCGGGCCGAGTTTCATCGGCAGCACCAGGCAGATGCCGTGGCTGAGGCAGAACGTTTGCTGGCCGGGCGTGAGAAGCTGCAAGGCGCCTGGTTGACTTGGGTAGc                        >  8:244378/1‑149 (MQ=255)
                                                                                                            ttCAAGGCCAGTGAATGACCGACAGTCACGACTCACCAGCCTCGCGGGCCGAGTTTCATCGGCAGCACCAGGCAGATGCCGTGGCTGAGGCAGAACGTTTGCTGGCCGGGCGTGAGAAGCTGCAAGGCGCCTGGTTGACTTgg                             <  3:51415/143‑1 (MQ=255)
                                                                                                            ttCAAGGCCAGTGAATGACCGACAGTCACGACTCACCAGCCTCGCGGGCCGAGTTTCATCGGCAGCACCAGGCAGATGCCGTGGCTGAGGCAGAACGTTTGCTGGCCGGGCGTGAGAAGCTGCAAGGCGCCTGGTTGACTTgg                             >  4:51415/1‑143 (MQ=255)
                                                                                                                 ggCCAGTGAATGACCGACAGTCACGACTCACCAGCCTCGCGGGCCGAGTTTCATCGGCAGCACCAGGCAGATGCCGTGGCTGAGGCAGAACGTTTGCTGGCCGGGCGTGAGAAGCTGCAAGGCGCCTGGTTGACTTGGGTAGCAGCt                    >  4:93895/1‑147 (MQ=255)
                                                                                                                     aGTGAATGACCGACGGTCACGACTCACCAGCCTCGCGGGCCGAGTTTCATCGGCAGCACCAGGCAGATGCCGTGGCTGAGGCAGAAGGTTTGCTGGCCGGGCGTGAGAAGCTGCAAgg                                             >  5:148240/1‑118 (MQ=255)
                                                                                                                     aGTGAATGACCGACGGTCACGACTCACCAGCCTCGCGGGCCGAGTTTCATCGGCAGCACCAGGCAGATGCCGTGGCTGAGGCAGAAGGTTTGCTGGCCGGGCGTGAGAAGCTGCAAgg                                             >  5:143301/1‑118 (MQ=255)
                                                                                                                     aGTGAATGACCGACGGTCACGACTCACCAGCCTCGCGGGCCGAGTTTCATCGGCAGCACCAGGCAGATGCCGTGGCTGAGGCAGAACGTTTGCTGGCCGGGCGTGAGAAGCTGCAAgg                                             <  6:148240/118‑1 (MQ=255)
                                                                                                                     aGTGAATGACCGACGGTCACGACTCACCAGCCTCGCGGGCCGAGTTTCATCGGCAGCACCAGGCAGATGCCGTGGCTGAGGCAGAACGTTTGCTGGCCGGGCGTGAGAAGCTGCAAgg                                             <  6:143301/118‑1 (MQ=255)
                                                                                                                      gTGAATGACCGACAGTCACGACTCACCAGCCTCGCGGGCCGAGTTTCATCGGCAGCACCAGGCAGATGCCGTGGCTGAGGCAGAACGTTTGCTGGCCGGGCGTGAGAAGCTGCAAGGCGCCTGGTTGACTTGGGTAGCGGCTGAGctc              >  5:317871/1‑148 (MQ=255)
                                                                                                                         aaTGACCGACAGTCACGACTCACCAGCCTCGCGGGCCGAGTTTCATCGGCAGCACCAGGCAGATGCCGTGGCTGAGGCAGAACGTTTGCTGGCCgg                                                               >  8:57149/1‑96 (MQ=255)
                                                                                                                         aaTGACCGACAGTCACGACTCACCAGCCTCGCGGGCCGAGTTTCATCGGCAGCACCAGGCAGATGCCGTGGCTGAGGCAGAACGTTTGCTGGCCgg                                                               <  7:57149/96‑1 (MQ=255)
                                                                                                                         aaTGACCGACAGTCACGACTCACCAGCCTCGCGGGCCGAGTTTCATCGGCAGCACCAGGCAGATGCCGTGGCTGAGGCAGAACGTTTGCTGGCCGGGCGTGAGAAGCTGCAAGGCGCCTGGTTGAc                                 >  3:210577/1‑126 (MQ=255)
                                                                                                                         aaTGACCGACAGTCACGACTCACCAGCCTCGCGGGCCGAGTTTCATCGGCAGCACCAGGCAGATGCCGTGGCTGAGGCAGAACGTTTGCTGGCCGGGCGTGAGAAGCTGCAAGGCGCCTGGTTGAc                                 <  4:210577/126‑1 (MQ=255)
                                                                                                                         aaTGACCGACAGTCACGACTCACCAGCCTCGCGGGCCGAGTTTCATCGGCAGCACCAGGCAGATGCCGTGGCTGAGGCAGAACGTTTGCTGGCCGGGCGTGAGAAGCTGCAAGGCGCCTGGTTGACTTGGGTAGCGGCTGAGCTCTATc          <  5:9432/149‑1 (MQ=255)
                                                                                                                          aTGACCGACAGTCACGACTCACCAGCCTCGCGGGCCGAGTTTCATCGGCAGCACCAGGCAGATGCCGTGGCTGAGGCAGAACGTTTGCTGGCCGGGCGTGAGAAGCTGCAAGGCGCCTGGTTGACTTGGGTAGCGGCTGAGCTCTATc          <  8:270073/148‑1 (MQ=255)
                                                                                                                          aTGACCGACAGTCACGACTCACCAGCCTCGCGGGCCGAGTTTCATCGGCAGCACCAGGCAGATGCCGTGGCTGAGGCAGAACGTTTGCTGGCCGGGCGTGAGAAGCTGCAAGGCGCCTGGTTGACTTGGGTAGCGGCTGAGCTCTATc          <  2:78797/148‑1 (MQ=255)
                                                                                                                             aCCGACAGTCACGACTCACCAGCCTCGCGGGCCGAGTTTCATCGGCAGCACCAGGCAGATGCCGTGGCTGAGGCAGAACGTTTGCTGGCCGGGCGTGAGAAGCTGCAAGGCGCCTGGTTGACTTGGGTAGCGGCTGAGCTCTATCGCCt      >  5:189319/1‑149 (MQ=255)
                                                                                                                              ccGACAGTCACGACTCACCAGCCTCGCGGGCCGAGTTTCATCGGCAGCACCAGGCAGATGCCGTGGCTGAGGCAGAACGTTTGCTGGCCGGGCGTGAGAAGCTGCAAGGCGCCTGGTTGACTTGGGTAGCGGCTGAGCTCTATCGCCt      >  5:14760/1‑148 (MQ=255)
                                                                                                                              ccGACAGTCACGACTCACCAGCCTCGCGGGCCGAGTTTCATCGGCAGCACCAGGCAGATGCCGTGGCTGAGGCAGAACGTTTGCTGGCCGGGCGTGAGAAGCTGCAAGGCGCCTGGTTGACTTGGGTAGCGGCTGAGCTCTATCGCCTc     <  6:205474/149‑1 (MQ=255)
                                                                                                                                 aCAGTCACGACTCACCAGCCTCGCGGGCCGAGTTTCATCGGCAGCACCAGGCAGATGCCGTGGCTGAGGCAGAACGTTTGCTGGCCGGGCGTGAGAAGCTGCAAGGCGCCTGGTTGACTTGGGTAGCGGCTGAGCTCTATCGCCTCAGc  >  6:268548/1‑149 (MQ=255)
                                                                                                                                   |                                                                                                                                                  
GTGATGCTGGTATTCGCCGCGCAGCAGCAAGTGCAGGAAGGCCTGGGCAAGTTGCTTAAAGCAGTGACTGGCAAAGTCGCCAAATGGCGCAAGGGTGAACTGACCACCTTCAAGGCCAGTGAATGACCGACAGTCACGACTCACCAGCCTCGCGGGCCGAGTTTCATCGGCAGCACCAGGCAGATGCCGTGGCTGAGGCAGAACGTTTGCTGGCCGGGCGTGAGAAGCTGCAAGGCGCCTGGTTGACTTGGGTAGCGGCTGAGCTCTATCGCCTCAGC  >  NC_002947/2532634‑2532911

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: