Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,134,906 G→A 10.6% P12S (CCG→TCG)  PP_4520 ← hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,134,9060GA10.6% 124.0 / 5.7 47P12S (CCG→TCG) PP_4520hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (21/21);  new base A (3/2);  total (24/23)
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.55e-01

AATCCGCAGCTCGGCATTGTTCAGGCCAAGGGCAGCCAGGCGCGCCGTGTTCAGGGCATGGGCAGTTTCACTGGGGTGCTGCAGCAGGATCACACGGGTGCGGCTGTCGAGCGCGGGGATCAGTGGGCATAGGCAATGGTCGAGCGGGCGCTGGCAGCGCTCGCAGCGGGGTCTGGGCATGGCGGTCTCCTTGGCGGGGACAGTTTGCCACA‑TAACCGGGGCCCGCTTTGCGGGCCCCGCGATTTCAGCGGTTGAACCTTTCCGCCAGGCTATGCAGGTACTCCGCCATGCG  >  NC_002947/5134760‑5135051
                                                                                                                                                  |                                                                                                                                                  
aaTCCGCAGCTCGGCATTGTTCAGGCCAAGGGCAGCCAGGCGCGCCGTGTTCAGGGCATGGGCAGTTTCACTGGGGTGCTGCAGCAGGATCACACGGGTGCGGCTGTCGAGCGCGGGGATCAGTGGGCATAGGCAATGGTCGAGCGGgc                                                                                                                                                  <  3:180976/149‑1 (MQ=255)
  tCCGCAGCTCGGCATTGTTCAGGCCAAGGGCAGCCAGGCGCGCCGTGTTCAGGGCATGGGCAGTTTCACTGGGGTGCTGCAGCAGGATCACACGGGTGCGGCTGTCGAGCGCGGGGATCAGTGGGCATAGGCAATGGTCGAGCGGgcgc                                                                                                                                                >  2:149559/1‑149 (MQ=255)
     gCAGCTCGGCATTGTTCAGGCCAAGGGCAGCCAGGCGCGCCGTGTTCAGGGCATGGGCAGTTTCACTGGGGTGCTGCAGCAGGATCACACGGGTGCGGCTGTCGAGCGCGGGGATCAGTGGGCATAGGCAATGGTCGAGCGGGCGCTgg                                                                                                                                             >  4:169756/1‑149 (MQ=255)
        gCTCGGCATTGTTCAGGCCAAGGGCAGCCAGGCGCGCCGTGTTCAGGGCATGGGCAGTTTCACTGGGGTGCTGCAGCAGGATCACACGGGTGCGGCTGTCGAGCGCGGGGATCAGTGGGCATAGGCAATGGTCGAGCGGGCGCTGGCAg                                                                                                                                          <  4:282081/149‑1 (MQ=255)
          tCGGCATTGTTCAGGCCAAGGGCAGCCAGGCGCGCCGTGTTCAGGGCATGGGCAGTTTCACTGGGGTGCTGCAGCAGGATCACACGGGTGCGGCTGTCGAGCGCGGGGATCAGTGGGCATAGGCAATGGTCGAGCGGgc                                                                                                                                                  >  1:138936/1‑139 (MQ=255)
            ggCATTGTTCAGGCCAAGGGCAGCCAGGCGCGCCGTGTTCAGGGCATGGGCAGTTTCACTGGGGTGCTGCAGCAGGATCACACGGGTGCGGCTGTCGAGCGCGGGGATCAGTGGGCATAGGCAATGGTCGAGCGGGCGCTGGCAgcg                                                                                                                                        <  6:162063/147‑1 (MQ=255)
            ggCATTGTTCAGGCCAAGGGCAGCCAGGCGCGCCGTGTTCAGGGCATGGGCAGTTTCACTGGGGTGCTGCAGCAGGATCACACGGGTGCGGCTGTCGAGCGCGGGGATCAGTGGGCATAGGCAATGGTCGAGCGGGCGCTGGCAGCGCt                                                                                                                                      <  8:291393/149‑1 (MQ=255)
               aTTGTTCAGGCCAAGGGCAGCCAGGCGCGCCGTGTTCAGGGCATGGGCAGTTTCACTGGGGTGCTGCAGCAGGATCACACGGGTGCGGCTGTCGAGCGCGGGGATCAGTGGGCATAGGCAATGGTCGAGCGGGCGCTGGCAGCGCTCg                                                                                                                                    >  3:101166/1‑148 (MQ=255)
                      aGGCCAAGGGCAGCCAGGCGCGCCGTGTTCAGGGCATGGGCAGTTTCACTGGGGTGCTGCAGCAGGATCACACGGGTGCGGCTGTCGAGCGCGGGGATCAGTGGGCATAGGCAATGGTCGAGCGGGCGCTGGCAGCGCTCGCAGCggg                                                                                                                             <  3:72655/148‑1 (MQ=255)
                       ggCCAAGGGCAGCCAGGCGCGCCCTGTTCAGGGCATGGGCAGTTTCACTGCGGAGCTGCAGCAGGATCACACGGGTGCGGCTGTCGAGCGCGGGGATCAGTGGGCATAGGCAATGGTCGAGCGGGCGCTGGCAGCGCTCGCAGCgggc                                                                                                                            >  5:175454/1‑147 (MQ=255)
                         ccAAGGGCAGCCAGGCGCGCCGTGTTCAGGGCATGGGCAGTTTCACTGGGGTGCTGCAGCAGGATCACACGGGTGCGGCTGTCGAGCGCGGAGATCAGTGGGCATAGGCAATGGTCGAGCGGGCGCTGGCAGCGCTCGCAGCGGGGTCt                                                                                                                         >  1:179116/1‑149 (MQ=255)
                           aaGGGCAGCCAGGCGCGCCGTGTTCAGGGCATGGGCAGTTTCACTGGGGTGCTGCAGCAGGATCACACGGGTGCGGCTGTCGAGCGCGGGGATCAGTGGGCATAGGCAATGGTCGAGCGGGCGCTGGCAGCGCTCGCAGCGGGGTCTgg                                                                                                                       <  5:143582/149‑1 (MQ=255)
                           aaGGGCAGCCAGGCGCGCCGTGTTCAGGGCATGGGCAGTTTCACTGGGGTGCTGCAGCAGGATCACACGGGTGCGGCTGTCGAGCGCGGGGATCAGTGGGCATAGGCAATGGTCGAGCGGGCGCTGGCAGCGCTCGCAGCGGGGTCTgg                                                                                                                       <  1:149559/149‑1 (MQ=255)
                              ggCAGCCAGGCGCGCCGTGTTCAGGGCATGGGCAGTTTCACTGGGGTGCTGCAGCAGGATCACACGGGTGCGGCTGTCGAGCGCGGGGATCAGTGGGCATAGGCAATGGTCGAGCGAGCGCTGCCAGCGCCCGctcg                                                                                                                                <  3:293528/137‑4 (MQ=255)
                              ggCAGCCAGGCGCGCCGTGTTCAGGGCATGGGCAGTTTCACTGGGGTGCTGCAGCAGGATCACACGGGTGCGGCTGTCGAGCGCGGGGATCAGTGGGCATAGGCAATGGTCGAGCGAGCGCTGCCAGCGCCCGctcg                                                                                                                                >  4:293528/1‑134 (MQ=255)
                                 aGCCAGGCGCGCCGTGTTCAGGGCATGGGCAGTTTCACTGGGGTGCTGCAGCAGGATCACACGGGTGCGGCTGTCGAGCGCGGGGATCAGTGGGCATAGGCAATGGTCGAGCGGGCGCTGGCAGCGCTCGCAGCGGGGTCTGGGCATgg                                                                                                                 <  4:101166/149‑1 (MQ=255)
                                              gtgtTCAGGGCATGGGCAGTTTCACTGGGGTGCTGCAGCAGGATCACACGGGTGCGGCTGTCGAGCGCGGGGATCAGTGGGCATAGGCAATGGTCGAGCGGGCGCTGGCAGCGCTCGCAGCGGGGTCTGGGCATGGCGGTCTCCTTGGc                                                                                                    <  6:162821/149‑1 (MQ=255)
                                               tgtTCAGGGCATGGGCAGTTTCACTGGGGTGCTGCAGCAGGATCACACGGGTGCGGCTGTCGAGCGCGGGGATCAGTGGGCATAGGCAATGGTCGAGCGGGCGCTGGCAGCGCTCGCAGCGGGGTCTGGGCATGGCGGTCTCCTTg                                                                                                      <  7:204686/146‑1 (MQ=255)
                                               tgtTCAGGGCATGGGCAGTTTCACTGGGGTGCTGCAGCAGGATCACACGGGTGCGGCTGTCGAGCGCGGGGATCAGTGGGCATAGGCAATGGTCGAGCGGGCGCTGGCAGCGCTCGCAGCGGGGTCTGGGCATGGCGGTCTCCTTg                                                                                                      >  8:204686/1‑146 (MQ=255)
                                                 ttCAGGGCATGGGCAGTTTCACTGGGGTGCTGCAGCAGGATCACACGGGTGCGGCTGTCGAGCGCGGGGATCAGTGGGCATAGGCAATGGTCGAGCGGGCGCTGGCAGCGCTCGCAGCGGGGTCTGGGCATGGCGGTCTCCTTGGCggg                                                                                                 <  3:345789/149‑1 (MQ=255)
                                                                     aCTGGGGTGCTGCAGCAGGATCACACGGGTGCGGCTGTCGAGCGCGGGGATCAGTGGGCATAGGCAATGGTCGAGCGGGCGCTGGCAGCGCTCGCAGCGGGGTCTGGGCATGGCGGTCTCCTTGGCGGGGACAGTTTGCCACA‑TAACCg                                                                            >  2:192781/1‑149 (MQ=255)
                                                                                gcagcaGGATCACACGGGTGCGGCTGTCGAGCGCGGGGATCAGTGGGCATAGGCAATGGTCGAGCGGGCGCTGGCAGCGCTCGCAGCGGGGTCTGGGCATGGCGGTCTCCTTGGCGGGGACAGTTTGCCACAGGGTCCGGGGCCCGCtt                                                                  >  1:119472/1‑149 (MQ=255)
                                                                                  agcagGATCACACGGGTGCGGCTGTCGAGCGCGGGGATCAGTGGGCATAGGCAATGGTCGAGCGGGCGCTGGCAGCGCTCGCAGCGGGGTCTGGGCATGGCGGTCTCCTTGGCGGGGACAGTTTGCCGCA‑TAACCGGGGCCCGCTTTGc                                                               <  2:82061/149‑1 (MQ=255)
                                                                                    cagGATCACACGGGTGCGGCTGTCGAGCGCGGGGATCAGTGGGCATAGGCAATGGTCGAGCGGGCGCTGGCAGCGCTCGCAGCGGGGTCTGGGCATGGCGGTCTCCTTGGCGGGGACAGTTTGCCACA‑TAACCGGGGCCCGCTTTGCga                                                             >  5:29555/1‑148 (MQ=255)
                                                                                    cagGATCACACGGGTGCGGCTGTCGAGCGAGGGGATCAGTGGGCATAGGCAATGGTCGAGCGGGCGCTGGCAGCGCTCGCAGCGGGGTCTGGGCATGGCGGTCTCCTTGGCGGGGACAGTTTGCCACA‑TAACCGGGGCCCGCTTTGc                                                               <  6:29555/147‑1 (MQ=255)
                                                                                         tCACACGGGTGCGGCTGTCGAGCGCGGAGATCAGTGGGCATAGGCAATGGTCGAGCGGGCGCTGGCAGCGCTCGCAGCGGGGTCTGGGCATGGCGGTCTCCTTGGCGGGGACAGTTTGCCACA‑TAACCGGGGCCCGCAAAGCGGGcccc                                                        <  2:179116/149‑1 (MQ=255)
                                                                                                     ggCTGTCGAGCGCGGGGATCAGTGGGCATAGGCAATGGTCGAGCGGGCGCTGGCAGCGCTCGCAGCGGGGTCTGGGCATGGCGGTCTCCTTGGCGGGGACAGTTTGCCACA‑TAACCGGGGCCCGCttt                                                                 <  1:99227/128‑1 (MQ=255)
                                                                                                     ggCTGTCGAGCGCGGGGATCAGTGGGCATAGGCAATGGTCGAGCGGGCGCTGGCAGCGCTCGCAGCGGGGTCTGGGCATGGCGGTCTCCTTGGCGGGGACAGTTTGCCACA‑TAACCGGGGCCCGCttt                                                                 >  2:99227/1‑128 (MQ=255)
                                                                                                      gCTGTCGAGCGCGGGGATCAGTGGGCATAGGCAATGGTCGAGCGGGCGCTGGCAGCGCTCGCAGCGGGGTCTGGGCATGGCGGTCTCCTTGGCGGGGACAGTTTGCCACA‑TAACCGGGGCCCGCTTTGCGGGCCCCGCGATTTCAGCgg                                           >  8:72538/1‑149 (MQ=255)
                                                                                                       cTGTCGAGCGCGGGGATCAGTGGGCATAGGCAATGGTCGAGCGGGCGCTGGCAGCGCTCGCAGCGGGGTCTGGGCATGGCGGTCTCCTTGGCGGGGACAGTTTGCCACA‑TAACCGGGGCCCGCTTTGCGGGCCCCGCGATTTCAGCgg                                           >  4:301098/1‑148 (MQ=255)
                                                                                                       cTGTCGAGCGCGGGGATCAGTGGGCATAGGCAATGGTCGAGCGGGCGCTGGCAGCGCTCGCAGCGGGGTCTGGGCATGGCGGTCTCCTTGGCGGGGACAGTTTGCCACA‑TAACCGGGGCCCGCTTTGCGGGCCCCGCGATTTCAGCgg                                           <  3:301098/148‑1 (MQ=255)
                                                                                                                gcgGGGATCAGTGGGCATAGGCAATGGTCGAGCGGGCGCTGGCAGCGCTCGCAGCGGGGTCTGGGCATGGCGGTCTCCTTGGCGGGGACAGTTTGCCACA‑TAACCGGGGCCCGCTTTGCGGGCCCCGCGATTTCAGCGGTTGAACCttt                                 >  7:39095/1‑149 (MQ=255)
                                                                                                                      aTCAGTGGGCATAGGCTTTGGTAGTGCGGGCGCTGGCAGCGGTGGCAGCGGGGTCTGGGCATGGCGGTCTCGTTGGCGGTGAGAGTTTGGCACA‑TAACGGGGGGGCGGTTTGGGGGCGCGGCGATTTCAGCGGTTGAACCTTTCCGCCa                           <  6:175454/149‑1 (MQ=255)
                                                                                                                      aTCAGTGGGCATAGGCAATGGTCGAGCGGGCGCTGGCAGCGCTCGCAGCGGGGTCTGGGCATGGCGGTCTCCTTGGCGGGGACAGTTTGCCACA‑TAACCGGGGCCCGCTTTGCGGGCCCCGCGATTTCAGCGGTTGAAc                                     <  5:302123/139‑1 (MQ=255)
                                                                                                                      aTCAGTGGGCATAGGCAATGGTCGAGCGGGCGCTGGCAGCGCTCGCAGCGGGGTCTGGGCATGGCGGTCTCCTTGGCGGGGACAGTTTGCCACA‑TAACCGGGGCCCGCTTTGCGGGCCCCGCGATTTCAGCGGTTGAAc                                     >  6:302123/1‑139 (MQ=255)
                                                                                                                              gCATAGGCAATGGTCGAGCGGGCGCTGGCAGCGCTCGCAGCGGGGTCTGGGCATGGCGGTCTCCTTGGCGGGGACAGTTTGCCACA‑TAACCGGGGCCCGCTTTGCGGGCCCCGCGATTTCAGCGGt                                          >  4:231305/1‑126 (MQ=255)
                                                                                                                              gCATAGGCAATGGTCGAGCGGGCGCTGGCAGCGCTCGCAGCGGGGTCTGGGCATGGCGGTCTCCTTGGCGGGGACAGTTTGCCACA‑TAACCGGGGCCCGCTTTGCGGGCCCCGCGATTTCAGCGGt                                          <  3:231305/126‑1 (MQ=255)
                                                                                                                              gCATAGGCAATGGGCGAGCGGGCGCTGGCAGCGCTCGCAGCGGGGGCTGGGCATGGCGGGCTCCTTGGCGGGGACAGTTTGCCACA‑TAACCGGGGCCCGCTTTGCGGGCCCCGCGATTTCAGCGGTTGAACCTTTCCGCCAGGCTATGc                   <  3:169756/149‑1 (MQ=255)
                                                                                                                                   ggCAATGGTCGAGCGGGCGCTGGCAGCGCTCGCAGCGGGGTCTGGGCATGGCGGTCTCCTTGGCGGGGACAGTTTGCCACA‑TAACCGGGGCCCGCTTTGCGGGCCCCGCGATTTCAGCGGTTGAACCTTTCCGCCAGGCTATGCAGGTa              >  3:216468/1‑149 (MQ=255)
                                                                                                                                    gCAATGGTCGAGCGGGCGCTGGCAGCGCTCGCAGCGGGGTCTGGGCATGGCGGTCTCCTTGGCGGGGACAGTTTGCCACA‑TAACCGGGGCCCGCTTTGCGGGCCCCGCGa                                                    <  6:322257/110‑1 (MQ=255)
                                                                                                                                    gCAATGGTCGAGCGGGCGCTGGCAGCGCTCGCAGCGGGGTCTGGGCATGGCGGTCTCCTTGGCGGGGACAGTTTGCCACA‑TAACCGGGGCCCGCTTTGCGGGCCCCGCGa                                                    >  5:322257/1‑110 (MQ=255)
                                                                                                                                    gCAATGGTCGAGCGGGCGCTGGCAGCGCTCGCAGCGGGGTCTGGGCATGGCGGTCTCCTTGGCGGGGACAGTTTGCCACA‑TAACCGGGGCCCGCTTTGCGGGCCCCGCGATTTCAGCGGTTGAACCTTTCCGCCAGGCTATGCAGGTAc             >  6:273468/1‑149 (MQ=255)
                                                                                                                                           tCGAGCGGGCGCTGGCAGCGCTCGCAGCGGGGTCTGGGCATGGCGGTCTCCTTGGCGGGGACAGTTTGCCACA‑TAACCGGGGCCCGCTTTGCGGGCCCCGCGATTTCAGCGGTTGAACCTTTCCGCCAGGCTATGCAGGTACTCCGCCa      >  6:295098/1‑149 (MQ=255)
                                                                                                                                           tCGAGCGGGCGCTGGCAGCGCTCGCAGCGGGGTCTGGGCATGGCGGTCTCCTTGGCGGGGACAGTTTGCCACA‑TAACCGGGGCCCGCTTTGCGGGCCCCGCGAGTTCAGCGGTTGAACCTTTCCGCCAGGCTATGCAGGTACTCCGCCa      >  5:255669/1‑149 (MQ=255)
                                                                                                                                             ctgCGAGCGCTGCCAGCGCCCGCAGCGGGGTCTGGGCATGGCGGTCTCCTTGGCGGGGACAGTTTGCCACA‑TAACCGGGGCCCGCTTTGCggg                                                            >  7:118709/3‑93 (MQ=255)
                                                                                                                                             ctgCGAGCGCTGCCAGCGCCCGCAGCGGGGTCTGGGCATGGCGGTCTCCTTGGCGGGGACAGTTTGCCACA‑TAACCGGGGCCCGCTTTGCggg                                                            <  8:118709/91‑1 (MQ=255)
                                                                                                                                               gCGAGCGCTGCCAGCGCCCGCAGCGGGGTCTGGGCATGGCGGTCTCCTTGGCGGGGACAGTTTGCCACA‑TAACCGGGGCACGCTTTGCGGGCCCCGCGATTTCAGCGGTTGAACCTTTCCGCCAGGCTATGCAGGTACTCCGCCATgcg  >  3:201080/1‑149 (MQ=255)
                                                                                                                                                  |                                                                                                                                                  
AATCCGCAGCTCGGCATTGTTCAGGCCAAGGGCAGCCAGGCGCGCCGTGTTCAGGGCATGGGCAGTTTCACTGGGGTGCTGCAGCAGGATCACACGGGTGCGGCTGTCGAGCGCGGGGATCAGTGGGCATAGGCAATGGTCGAGCGGGCGCTGGCAGCGCTCGCAGCGGGGTCTGGGCATGGCGGTCTCCTTGGCGGGGACAGTTTGCCACA‑TAACCGGGGCCCGCTTTGCGGGCCCCGCGATTTCAGCGGTTGAACCTTTCCGCCAGGCTATGCAGGTACTCCGCCATGCG  >  NC_002947/5134760‑5135051

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