Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,134,913 G→C 12.8% C9W (TGC→TGG)  PP_4520 ← hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,134,9130GC12.8% 128.1 / 11.0 47C9W (TGC→TGG) PP_4520hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (20/21);  new base C (4/2);  total (24/23)
Fisher's exact test for biased strand distribution p-value = 6.66e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.71e-01

GCAGCTCGGCATTGTTCAGGCCAAGGGCAGCCAGGCGCGCCGTGTTCAGGGCATGGGCAGTTTCACTGGGGTGCTGCAGCAGGATCACACGGGTGCGGCTGTCGAGCGCGGGGATCAGTGGGCATAGGCAATGGTCGAGCGGGCGCTGGCAGCGCTCGCAGCGGGGTCTGGGCATGGCGGTCTCCTTGGCGGGGACAGTTTGCCACA‑TAACCGGGGCCCGCTTTGCGGGCCCCGCGATTTCAGCGGTTGAACCTTTCCGCCAGGCTATGCAGGTACTCCGCCATGCGCTCGAGGTCC  >  NC_002947/5134765‑5135061
                                                                                                                                                    |                                                                                                                                                     
gCAGCTCGGCATTGTTCAGGCCAAGGGCAGCCAGGCGCGCCGTGTTCAGGGCATGGGCAGTTTCACTGGGGTGCTGCAGCAGGATCACACGGGTGCGGCTGTCGAGCGCGGGGATCAGTGGGCATAGGCAATGGTCGAGCGGGCGCTgg                                                                                                                                                       >  4:169756/1‑149 (MQ=255)
   gCTCGGCATTGTTCAGGCCAAGGGCAGCCAGGCGCGCCGTGTTCAGGGCATGGGCAGTTTCACTGGGGTGCTGCAGCAGGATCACACGGGTGCGGCTGTCGAGCGCGGGGATCAGTGGGCATAGGCAATGGTCGAGCGGGCGCTGGCAg                                                                                                                                                    <  4:282081/149‑1 (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                                                                                                                                 <  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                                                    <  3:231305/126‑1 (MQ=255)
                                                                                                                         gCATAGGCAATGGTCGAGCGGGCGCTGGCAGCGCTCGCAGCGGGGTCTGGGCATGGCGGTCTCCTTGGCGGGGACAGTTTGCCACA‑TAACCGGGGCCCGCTTTGCGGGCCCCGCGATTTCAGCGGt                                                    >  4:231305/1‑126 (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                                                                      >  7:118709/3‑93 (MQ=255)
                                                                                                                                        ctgCGAGCGCTGCCAGCGCCCGCAGCGGGGTCTGGGCATGGCGGTCTCCTTGGCGGGGACAGTTTGCCACA‑TAACCGGGGCCCGCTTTGCggg                                                                      <  8:118709/91‑1 (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  <  7:72538/149‑1 (MQ=255)
                                                                                                                                                    gCAGCGCTCGCAGCGGGGTCTGGGCATGGCGGTCTCCTTGGCGGGGACAGTTTGCCACA‑TAACCGGGGCCCGCTTTGCGGGCCCCGCGATTTCAGCGGTTGAACCTTTCCGCCAGGCTATGCAGGTACTCCGCCATGCGCTCGAGGTcc  >  6:322646/1‑149 (MQ=255)
                                                                                                                                                    gCAGCGCTCGCAGCGGGGTCTGGGCATGGCGGTCTCCTTGGCGGGGACAGTTTGCCACA‑TAACCGGGGCCCGCTTTGCGGGCCCCGCGATTTCAGCGGTTGAACCTTTCCGCCAGGCTATGCAGGTACTCCGCCATGCGCTCGAGGTcc  <  8:39095/149‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
GCAGCTCGGCATTGTTCAGGCCAAGGGCAGCCAGGCGCGCCGTGTTCAGGGCATGGGCAGTTTCACTGGGGTGCTGCAGCAGGATCACACGGGTGCGGCTGTCGAGCGCGGGGATCAGTGGGCATAGGCAATGGTCGAGCGGGCGCTGGCAGCGCTCGCAGCGGGGTCTGGGCATGGCGGTCTCCTTGGCGGGGACAGTTTGCCACA‑TAACCGGGGCCCGCTTTGCGGGCCCCGCGATTTCAGCGGTTGAACCTTTCCGCCAGGCTATGCAGGTACTCCGCCATGCGCTCGAGGTCC  >  NC_002947/5134765‑5135061

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