Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,134,920 T→C 9.6% E7G (GAG→GGG)  PP_4520 ← hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,134,9200TC 9.6% 133.2 / 4.9 52E7G (GAG→GGG) PP_4520hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (24/23);  new base C (3/2);  total (27/25)
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

GGCATTGTTCAGGCCAAGGGCAGCCAGGCGCGCCGTGTTCAGGGCATGGGCAGTTTCACTGGGGTGCTGCAGCAGGATCACACGGGTGCGGCTGTCGAGCGCGGGGATCAGTGGGCATAGGCAATGGTCGAGCGGGCGCTGGCAGCGCTCGCAGCGGGGTCTGGGCATGGCGGTCTCCTTGGCGGGGACAGTTTGCCACA‑TAACCGGGGCCCGCTTTGCGGGCCCCGCGATTTCAGCGGTTGAACCTTTCCGCCAGGCTATGCAGGTACTCCGCCATGCGCTCGAGGTCCTGG  >  NC_002947/5134772‑5135064
                                                                                                                                                    |                                                                                                                                                 
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                                                                                                                                    <  1:149559/149‑1 (MQ=255)
               aaGGGCAGCCAGGCGCGCCGTGTTCAGGGCATGGGCAGTTTCACTGGGGTGCTGCAGCAGGATCACACGGGTGCGGCTGTCGAGCGCGGGGATCAGTGGGCATAGGCAATGGTCGAGCGGGCGCTGGCAGCGCTCGCAGCGGGGTCTgg                                                                                                                                    <  5:143582/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                                                                                                                   >  8:204686/1‑146 (MQ=255)
                                   tgtTCAGGGCATGGGCAGTTTCACTGGGGTGCTGCAGCAGGATCACACGGGTGCGGCTGTCGAGCGCGGGGATCAGTGGGCATAGGCAATGGTCGAGCGGGCGCTGGCAGCGCTCGCAGCGGGGTCTGGGCATGGCGGTCTCCTTg                                                                                                                   <  7:204686/146‑1 (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                                                                              >  2:99227/1‑128 (MQ=255)
                                                                                         ggCTGTCGAGCGCGGGGATCAGTGGGCATAGGCAATGGTCGAGCGGGCGCTGGCAGCGCTCGCAGCGGGGTCTGGGCATGGCGGTCTCCTTGGCGGGGACAGTTTGCCACA‑TAACCGGGGCCCGCttt                                                                              <  1:99227/128‑1 (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                                                  >  6:302123/1‑139 (MQ=255)
                                                                                                          aTCAGTGGGCATAGGCAATGGTCGAGCGGGCGCTGGCAGCGCTCGCAGCGGGGTCTGGGCATGGCGGTCTCCTTGGCGGGGACAGTTTGCCACA‑TAACCGGGGCCCGCTTTGCGGGCCCCGCGATTTCAGCGGTTGAAc                                                  <  5:302123/139‑1 (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                                                                 >  5:322257/1‑110 (MQ=255)
                                                                                                                        gCAATGGTCGAGCGGGCGCTGGCAGCGCTCGCAGCGGGGTCTGGGCATGGCGGTCTCCTTGGCGGGGACAGTTTGCCACA‑TAACCGGGGCCCGCTTTGCGGGCCCCGCGa                                                                 <  6:322257/110‑1 (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                                                                         <  8:118709/91‑1 (MQ=255)
                                                                                                                                 ctgCGAGCGCTGCCAGCGCCCGCAGCGGGGTCTGGGCATGGCGGTCTCCTTGGCGGGGACAGTTTGCCACA‑TAACCGGGGCCCGCTTTGCggg                                                                         >  7:118709/3‑93 (MQ=255)
                                                                                                                                   gCGAGCGCTGCCAGCGCCCGCAGCGGGGTCTGGGCATGGCGGTCTCCTTGGCGGGGACAGTTTGCCACA‑TAACCGGGGCACGCTTTGCGGGCCCCGCGATTTCAGCGGTTGAACCTTTCCGCCAGGCTATGCAGGTACTCCGCCATgcg               >  3:201080/1‑149 (MQ=255)
                                                                                                                                      cccgcTGCGAGCGCTACCAGCGGGGTCTGGGCATGGCGGTCTCCTTGGCGGGGACAGTTTGCCACA‑TAACCGGGGCCCGCTTTGCGGGCCCCGCGATTTCAGCGGTTGAACCTTTCCGCCAGGCTATGCAGGTACTCCGCCATGCGCTc            >  5:198990/3‑149 (MQ=255)
                                                                                                                                         gcTGGCAGCGCTCGCAGCGGGGTCTGGGCATGGCGGTCTCCTTGGCGGGGACAGTTTGCCACA‑TAACCGGGGCCCGCTTTGCGGGCCCCGCGATTTCAGCGGTTGAACCTTTCCGCCAGGCTATGCAGGTACTCCGCCATGCGCTCGAg         >  2:54548/1‑149 (MQ=255)
                                                                                                                                           tGGCAGCGCTCGCTGGGGGGTCTGGGCATGGCGGGCTCCTTGGCGGGGACTGTTTGCCACA‑TAACCGGGGCCCGCTTTGCGGGCCCCGCGATTTCAGCGGTTGAACCtt                                               <  4:37397/109‑1 (MQ=255)
                                                                                                                                           tGGCAGCGCTCGCAGCGGGGTCTGGGCATGGCGGTCTCCTTGGCGGGGACAGTTTGCCACA‑TAACCGGGGCCCGCTTTGCGGGCCCCGCGATTTAAGCGGTTGAACCtt                                               >  3:37397/1‑109 (MQ=255)
                                                                                                                                             gCAGCGCTCGCAGCGGGGTCTGGGCATGGCGGTCTCCTTGGCGGGGACAGTTTGCCACA‑TAACCGGGGCCCGCTTTGCGGGCCCCGCGATTTCAGCGGTTGAACCTTTCCGCCAGGCTATGCAGGTACTCCGCCATGCGCTCGAGGTcc     >  6:322646/1‑149 (MQ=255)
                                                                                                                                             gCAGCGCTCGCAGCGGGGTCTGGGCATGGCGGTCTCCTTGGCGGGGACAGTTTGCCACA‑TAACCGGGGCCCGCTTTGCGGGCCCCGCGATTTCAGCGGTTGAACCTTTCCGCCAGGCTATGCAGGTACTCCGCCATGCGCTCGAGGTcc     <  8:39095/149‑1 (MQ=255)
                                                                                                                                             gCAGCGCTCGCAGCGGGGTCTGGGCATGGCGGTCTCCTTGGCGGGGACAGTTTGCCACA‑TAACCGGGGCCCGCTTTGCGGGCCCCGCGATTTCAGCGGTTGAACCTTTCCGCCAGGCTATGCAGGTACTCCGCCATGCGCTCGAGGTcc     <  7:72538/149‑1 (MQ=255)
                                                                                                                                              cAGCGCTCGCAGCGGGGTCTGGGCATGGCGGTCTCCTTGGCGGGGACAGTTTGCCACA‑TAACCGGGGCCCGCTTTGCGGGCCCCGCGATTTCAGCGGTTGAACCTTTCCGCCAGGCTATGCAGGTACTCCGCCATGCGCTCGAGGTCCt    >  2:209293/1‑149 (MQ=255)
                                                                                                                                              cAGCGCTCGCAGCGGGGTCTGGGCATGGCGGTCTCCTTGGCGGGGACAGTTTGCCACA‑TAACCGGGGCCCGCTTTGCGGGCCCCGCGATTTCAGCGGTTGAACCTTTCCGCCAGGCTATGCAGGTACTCCGCCATGCGCTCGAGGTCCt    <  5:295098/149‑1 (MQ=255)
                                                                                                                                               aGCGCTCGCAGCGGGGTCTGGGCATGGCGGTCTCCTTGGCGGGGACAGTTTGCCACA‑TAACCGGGGCCCGCTTTGCGGGCCCCGCGATTTCAGCGGTTGAACCTTTCCGCCAGGCTATGCAGGTACTCCGCCATGCGCTCGAg         <  3:84797/143‑1 (MQ=255)
                                                                                                                                               aGCGCTCGCAGCGGGGTCTGGGCATGGCGGTCTCCTTGGCGGGGACAGTTTGCCACA‑TAACCGGGGCCCGCTTTGCGGGCCCCGCGATTTCAGCGGTTGAACCTTTCCGCCAGGCTATGCAGGTACTCCGCCATGCGCTCGAg         >  4:84797/1‑143 (MQ=255)
                                                                                                                                                gcgcTCGCAGCGGGGTCTGGGCATGGCGGTCTCCTTGGCGGGGACAGTTTGCCACA‑TAACCGGGGCCCGCTTTGCGGGCCCCGCGATTTCAGCGGTTGAACCTTTCCGCCAGGCTATGCAGGTACTCCGCCATGCGCTCGAGGTCCTgg  <  1:192781/149‑1 (MQ=255)
                                                                                                                                                    tCGCAGCGGGGTCTGGGCATGGCGGTCTCCTTGGCgg                                                                                                               <  6:307012/37‑1 (MQ=255)
                                                                                                                                                    tCGCAGCGGGGTCTGGGCATGGCGGTCTCCTTGGCgg                                                                                                               >  5:307012/1‑37 (MQ=255)
                                                                                                                                                    |                                                                                                                                                 
GGCATTGTTCAGGCCAAGGGCAGCCAGGCGCGCCGTGTTCAGGGCATGGGCAGTTTCACTGGGGTGCTGCAGCAGGATCACACGGGTGCGGCTGTCGAGCGCGGGGATCAGTGGGCATAGGCAATGGTCGAGCGGGCGCTGGCAGCGCTCGCAGCGGGGTCTGGGCATGGCGGTCTCCTTGGCGGGGACAGTTTGCCACA‑TAACCGGGGCCCGCTTTGCGGGCCCCGCGATTTCAGCGGTTGAACCTTTCCGCCAGGCTATGCAGGTACTCCGCCATGCGCTCGAGGTCCTGG  >  NC_002947/5134772‑5135064

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