Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,577,020 G→A 7.5% A102A (GCG→GCA)  PP_4057 → membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,577,0200GA 7.5% 153.8 / 2.5 53A102A (GCG→GCA) PP_4057membrane protein
Reads supporting (aligned to +/- strand):  ref base G (27/22);  new base A (2/2);  total (29/24)
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.52e-01

GCTCAGTGCCAGTGGGTTGAGTCTCGGCACCCTCAATGCGTTTTCGCAGCAGATGAGCGAACGCCACCAGCAGCCGGCACCGGACCTGATTGCCAGTCAGTGGGCGCAGTTCTTCCCCACCAGTGCGCGTGGCCGCGCGCTGCCCCCGGCACAGCTGGAGGCGCCCAGCCAGCAACTGATGATCGGCCCGGACCTGTTCGTGCGTGAAACTGCAGCGGGCGATGTGCACCGCGCAGGGATTTTCGTGGGGCACAACAACCTGCAAAGCAGCTTCAGAG  >  NC_002947/4576882‑4577159
                                                                                                                                          |                                                                                                                                           
gCTCAGTGCCAGTGGGTTGAGTCTCGGCACCCTCAATGCGTTTTCGCAGCAGATGAGCGAACGCCACCAGCAGCCGGCACCGGACCTGATTGCCAGTCAGTGGGCGCAGTTCTTCCCCACCAGTGCGCGTGGCCGCGCGCTGCCCCCgg                                                                                                                                   >  7:192586/1‑149 (MQ=255)
gCTCAGTGCCAGTGGGTTGAGTCTCGGCACCCTCAATGCGTTTTCGCAGCAGATGAGCGAACGCCACCAGCAGCCGGCACCGGACCTGATTGCCAGTCAGTGGGCGCAGTTCTTCCCCACCAGTGCGCGTGGCCGCGCGCTGCCCCCgg                                                                                                                                   <  6:42290/149‑1 (MQ=255)
   cAGTGCCAGTGGGTTGAGTCTCGGCACCCTCAATGCGTTTTCGCAGCAGATGAGCGAACGCCACCAGCAGCCGGCACCGGACCTGATTGCCAGTCAGTGGGCGCAGTTCTTCCCCACCAGTGCGCGTGGCCGCGCGCTGCCCCCGGca                                                                                                                                 >  2:152962/1‑148 (MQ=255)
           gTGGGTTGAGTCTCGGCACCCTCAATGCGTTTTCGCAGCAGATGAGCGAACGCCACCAGCAGCCGGCACCGGACCTGATTGCCAGTCAGTGGGCGCAGTTCTTCCCCACCAGTGCGCGTGGCCGCGCGCTGCCCCCGGCACAGCTGGAg                                                                                                                        <  2:265582/149‑1 (MQ=255)
               gTTGAGTCTCGGCACCCTCAATGCGTTTTCGCAGCAGATGAGCGAACGCCACCAGCAGCCGGCACCGGACCTGATTGCCAGTCAGTGGGCGCAGTTCTTCCCCACCAGTGCGCGTGGCCGCGCGCTGCCCCCGGCACAGCTGGAGgcgc                                                                                                                    >  4:62498/1‑149 (MQ=255)
               gTTGAGTCTCGGCACCCTCAATGCGTTTTCGCAGCAGATGAGCGAACGCCACCAGCAGCCGGCACCGGACCTGATTGCCAGTCAGTGGGCGCAGTTCTTCCCCACCAGTGCGCGTGGCCGCGCGCTGCCCCCGGCACAGCTGGAGgcgc                                                                                                                    <  8:319543/149‑1 (MQ=255)
                    gTCTCGGCACCCTCAATGCGTTTTCGCAGCAGATGAGCGAACGCCACCAGCAGCCGGCACCGGACCTGATTGCCAGTCAGTGGGCGCAGTTCTTCCCCACCAGTGCGCGTGGCCGCGCGCTGCCCCCGGCACAGCTGGAGGCGCccagc                                                                                                               <  1:152962/149‑1 (MQ=255)
                    gTCTCGGCACCCTCAATGCGTTTTCGCAGCAGATGAGCGAACGCCACCAGCAGCCGGCACCGGACCTGATTGCCAGTCAGTGGGCGCAGTTCTTCCCCACCAGTGCGCGTGGCCGCGCGCTGCCCCCGGCACAGCTGGAGGCGCccagc                                                                                                               >  7:102897/1‑149 (MQ=255)
                                  aaTGCGTTTTCGCAGCAGATGAGCGAACGCCACCAGCAGCCGGCACCGGACCTGATTGCCAGTCAGTGGGCGCAGTTCTTCCCCACCAGTGCGCGTGGCCGCGCGCTGCCCCCGGCACAGCTGGAGGCGCCCAGCCAGCAACTgatgat                                                                                                 >  4:111321/1‑149 (MQ=255)
                                  aaTGCGTTTTCGCAGCAGATGAGCGAACGCCACCAGCAGCCGGCACCGGACCTGATTGCCAGTCAGTGGGCGCAGTTCTTCCCCACCAGTGCGCGTGGCCGCGCGCTGCCCCCGGCACAGCTGGAGGCGCCCAGCCAGCAACTgatgat                                                                                                 >  5:102927/1‑149 (MQ=255)
                                   aTGCGTTTTCGCAGCAGATGAGCGAACGCCACCAGCAGCCGGCACGGGACCTGATTGCCAGTCAGTGGGCGCAGTTCTTCCCCACCAGTGCGCGTGGCCGCGCGCTGCCCCCGGCACAGCTGGAGGCGCCCAGCCAGCAACTGATGATc                                                                                                >  1:262937/1‑149 (MQ=255)
                                      cGTTTTCGCAGCAGATGAGCGAACGCCACCAGCAGCCGGCACCGGACCTGATTGCCAGTCAGTGGGCGCAGTTCTTCCCCACCAGTGCGCGTGGCCGCGCGCTGCCCCCGGCACAGCTGGAGGCGCCCAGCCAGCAACTGATGATCGGc                                                                                             <  2:242123/149‑1 (MQ=255)
                                         tttCGCAGCAGATGAGCGAACGCCACCAGCAGCCGGCACCGGACCTGATTGCCAGTCAGTGGGCGCAGTTCTTCCCCACCAGTGCGCGTGGCCGCGCGCTGCCCCCGGCACAGCTGGAGGCGCCCAGCCAGCAACTGATGATCGGCCCg                                                                                          >  5:63681/1‑149 (MQ=255)
                                              cagcagATGAGCGAACGCCACCAGCAGCCGGCACCGGACCTGATTGCCAGTCAGAGGGCGCAGTTCTTCCCCACCAGTGCGCGTGGCCGCGCACTGCCCCCGGCACAGCTGGAGGCGCCCAGCCAGCAACTGATGATCGGCCCGGACCt                                                                                     >  2:161689/1‑149 (MQ=255)
                                              cagcagATGAGCGAACGCCACCAGCAGCCGGCACCGGACCTGATTGCCAGTCAGAGGGCGCAGTTCTTCCCCACCAGTGCGCGTGGCCGCGCACTGCCCCCGGCACAGCTGGAGGCGCCCAGCCAGCAACTGATGATCGGCCCGGACCt                                                                                     >  2:167184/1‑149 (MQ=255)
                                                gcagATGAGCGAACGCCACCAGCAGCCGGCACCGGACCTGATTGCCAGTCAGTGGGCGCAGTTCTTCCCCACCAGTGCGCGTGGCCGCGCGCTGCCCCCGGCACAGCTGGAGGCGCCCAGCCAGCAACTGATGATCGGCCCGGACCTGt                                                                                   <  3:141110/149‑1 (MQ=255)
                                                   gATGAGCGAACTCCACCAGCAGCCGGCACCGGACCTGATTGCCAGTCAGTGGGCGCAGTTCTTCCCCACCAGTGCGCGTGGCCGCGCGCTGCCCCCGGCACAGCTGGAg                                                                                                                        <  1:209592/109‑1 (MQ=255)
                                                   gATGAGCGAACACCACCAGCAGCCGGCACCGGACCTGATTGCCAGTCAGTGGGCGCAGTTCTTCCCCACCAGTGCGCGTGGCCGCGCGCTGCCCCCGGCACAGCTGGAg                                                                                                                        >  2:209592/1‑109 (MQ=255)
                                                      gAGCGAACGCCACCAGCAGCCGGCACCGGACCTGATTGCCAGTCAGTGGGCGCAGTTCTTCCCCACCAGTGCGCGTGGCCGCGCGCTGCCCCCGGCACAGCTGGAGGCGCCCAGCCAGCAACTGATGATCGGCCCGGACCTGTTcgtgc                                                                             >  4:285369/1‑149 (MQ=255)
                                                       aGCGAACGCCACCAGCAGCCGGCACCGGACCTGATTGCCAGTCAGTGGGCGCAGTTCTTCCCCACCAGTGCGCGTGGCCGCGCGCTGCCCCCGGCACAGCTGGAGGCGCCCAGCCAGCAACTGATGATCGGCCCGGACCTGTTcgtgc                                                                             <  8:102897/148‑1 (MQ=255)
                                                        gCGAACGCCACCAGCAGCCGGCACCGGACCTGATTGCCAGTCAGTGGGCGCAGTTCTTCCCCACCAGTGCGCGTGGCCGCGCGCTGCCCCCGGCACAGCTGGAGGCGCCCAGCCAGCAACTGATGATCGGCCCGGACCTGTTcgtgcgt                                                                           <  2:198412/149‑1 (MQ=255)
                                                            aCGCCACCAGCAGCCGGCACCGGACCTGATTGCCAGTCAGTGGGCGCAGTTCTTCCCCACCAGTGCGCGTGGCCGCGCGCTGCCCCCGGCACAGCTGGAGGCGCCCAGCCAGCAACTGATGATCGGCCCGGACCTGTTCGTGCGTGaaa                                                                       <  3:62498/149‑1 (MQ=255)
                                                              gCCACCAGCAGCCGGCACCGGACCTGATTGCCAGTCAGTGGGCGCAGTTCTTCCCCACCAGTGCGCGTGGCCGCGCGCTGCCCCCGGCACAGCTGGAGGCGCCCAGCCAGCAACTGATGATCGGACCGGACCTGTTCGTGCGTGAAACt                                                                     >  3:238882/1‑149 (MQ=255)
                                                                caccaGCAGCCGGCACCGGACCTGATTGCCAGTCAGTGGGCGCAGTTCTTCCCCACCAGTGCGCGTGGCCGCGCGCTGCCCCCGGCACAGCTGGAGGCGCCCAGCCAGCAACTGATGATCGGCCCGGACCTGTTCGTGCGTGAAACTg                                                                    >  3:26765/1‑148 (MQ=255)
                                                                 accaGCAGCCGGCACCGGACCTGATTGCCAGTCAGTGGGCGCAGTTCTTCCCCACCAGTGCGCGTGGCCGCGCGCTGCCCCCGGCACAGCTGGAGGCGCCCAGCCAGCAACTGATGATCGGCCCGGACCTGTTCGTGCGTGAAACTg                                                                    >  7:99146/1‑147 (MQ=255)
                                                                 accaGCAGCCGGCACCGGACCTGATTGCCAGTCAGTGGGCGCAGTTCTTCCCCACCAGTGCGCGTGGCCGCGCGCTGCCCCCGGCACAGCTGGAGGCGCCCAGCCAGCAACTGATGATCGGCCCGGACCTGTTCGTGCGTGAAACTg                                                                    >  4:130400/1‑147 (MQ=255)
                                                                 accaGCAGCCGGCACCGGACCTGATTGCCAGTCAGAGGGCGCAGTTCTTCCCCACCAGTGCGCGTGGCCGCGCACTGCCCCCGGCACAGCTGGAGGCGCCCAGCCAGCAACTGATGATCGGCCCGGACCTGTTCGTGCGTGAAACTGCa                                                                  <  1:167184/149‑1 (MQ=255)
                                                                 accaGCAGCCGGCACCGGACCTGATTGCCAGTCAGAGGGCGCAGTTCTTCCCCACCAGTGCGCGTGGCCGCGCACTGCCCCCGGCACAGCTGGAGGCGCCCAGCCAGCAACTGATGATCGGCCCGGACCTGTTCGTGCGTGAAACTGCa                                                                  <  1:161689/149‑1 (MQ=255)
                                                                        gcCGGCACCGGACCTGATTGCCAGTCAGTGGGCGCAGTTCTTCCCCACCAGTGCGCGTGGCCGCGCGCTGCCCCCGGCACAGCTGGAGGCGCCCAGCCAGCAACTGATGATCGGCCCGGACCTGTTCGTGCGTGAAACTGCAGCGGGc                                                            >  5:189180/1‑148 (MQ=255)
                                                                          cGGCACCGGACCTGATTGCCAGTCAGTGGGCGCAGTTCTTCCCCACCAGTGCGCGTGGCCGCGCGCTGCCCCCGGCACAGCTGGAGGCGCCCAGCCAGCAACTGATGATCGGCCCGGACCTGTTCGTGCGTGAAACTGCAGCGGGCGAt                                                         <  1:122646/149‑1 (MQ=255)
                                                                           ggCACCGGACCTGATTGCCAGTCAGTGGGCGCAGTTCTTCCCCACCAGTGCGCGTGGCCGCGCGCTGCCCCCGGCACAGCTGGAGGCGCCCAGCCAGCAACTGATGATCGGCCCGGACCTGTTCGTGCGTGAAACTGCAGCGGGCGAtg                                                        >  1:255138/1‑149 (MQ=255)
                                                                              aCCGGACCTGATTGCCAGTCAGTTGGCGCAGTTCTTCCCCACCAGTGCGCGTGGCCGCGCGCTGCCCCCGGCACAGCTGGAGGCGCCCAGCCAGCAACTGATGATCGGCCCGGACCTGTTCGTGCGTGAAACTGCAGCGGGCGATGTGc                                                     <  8:192586/149‑1 (MQ=255)
                                                                              aCCGGACCTGATTGCCAGTCAGTGGGCGCAGTTCTTCCCCACCAGTGCGCGTGGCCGCGCGCTGCCCCCGGCACAGCTGGAGGCGCCCAGCCAGCAACTGATGATCGGCCCGGACCTGTTCGTGCGTGAAACTGCAGCGGGCGAtgtg                                                      >  2:7134/1‑148 (MQ=255)
                                                                                cGGACCTGATTGCCAGTCAGTGGGCGCAGTTCTTCCCCACCAGTGCGCGTGGCCGCGCGCTGCCCCCGGCACAGCTGGAGGCGCCCAGCCAGCAACTGATGATCGGCCCGGACCTGTTCGTGCGTGAAACTGCAGCGGGCGAtgtg                                                      <  2:1624/146‑1 (MQ=255)
                                                                                cGGACCTGATTGCCAGTCAGTGGGCGCAGTTCTTCCCCACCAGTGCGCGTGGCCGCGCGCTGCCCCCGGCACAGCTGGAGGCGCCCAGCCAGCAACTGATGATCGGCCCGGACCTGTTCGTGCGTGAAACTGCAGCGGGCGAtgtg                                                      >  1:1624/1‑146 (MQ=255)
                                                                                  gACCTGATTGCCAGTCAGTGGGCGCAGTTCTTCCCCACCAGTGCGCGTGGCCGCGCGCTGCCCCCGGCACAGCTGGAGGCGCCCAGCCAGCAACTGATGATCGGCCCGGACCTGTTCGTGCGTGAAACTGCAGCGGGCGATGTGCACcg                                                 >  3:332024/1‑149 (MQ=255)
                                                                                      tGATTGCCAGTCAGTGGGCGCAGTTCTTCCCCACCAGTGCGCGTGGCCGCGCGCTGCCCCCGGCACAGCTGGAGGCGCCCAGCCAGCAACTGATGATCGGCCCGGACCTGTTCGTGCGTGAAACTGCAGCGGGCGATGTGCACcgcg                                               >  1:267664/1‑147 (MQ=255)
                                                                                           gCCAGTCAGTGGGCGCAGTTCTTCCCCACCAGTGCGCGTGGCCGCGCGCTGCCCCCGGCACAGCTGGAGGCGCCCAGCCAGCAACTGATGATCGGCCCGGACCTGTTCGTGCGTGAAACTGCAGCGGGCGATGTGCACCGCGCAGGGa                                         <  3:285369/148‑1 (MQ=255)
                                                                                             cagtcagtGGGCGCAGTTCTTCCCCACCAGTGCGCGTGGCCGCGCGCTGCCCCCGGCACAGCTGGAGGCGCCCAGCCAGCAACTGATGATCGGCCCGGACCTGTTCGTGCGTGAAACTGCAGCGGGCGATGTGCACCGCGCAGGGAttt                                      >  5:136425/1‑149 (MQ=255)
                                                                                                    tGGGCGCAGTTCTTCCCCACCAGTGCGCGTGGCCGCGCGCTGCCCCCGGCACAGCTGGAGGCGCCCAGCCAGCAACTGATGATCGGCCCGGACCTGTTCGTGCGTGAAACTGCAGCGGGCGATGTGCACCGCGCAGGGATTTTCGTggg                               <  4:332024/149‑1 (MQ=255)
                                                                                                    tGGGCGCAGTTCTTCCCCACCAGTGCGCGTGGCCGCGCGCTGCCCCCGGCACAGCTGGAGGCGCCCAGCCAGCAACTGATGATCGGCCCGGACCTGTTCGTGCGTGAAACTGCAGCGGGCGATGTGCACCGCGCAGGGATTTTCGTggg                               <  3:111321/149‑1 (MQ=255)
                                                                                                      ggCGCAGTTCTTCCCCACCAGTGCGCGTGGCCGCGCGCTGCCCCCGGCACAGCTGGAGGCGCCCAGCCAGCAACTGATGATCGGCCCGGACCTGTTCGTGCGTGAAACTGCAGCGGGCGATGTGCACCGCGCAGGGATTTTCGTgggg                              >  3:120646/1‑148 (MQ=255)
                                                                                                      ggCGCAGTTCTTCCCCACCAGTGCGCGTGGCCGCGCGCTGCCCCCGGCACAGCTGGAGGCGCCCAGCCAGCAACTGATGATCGGCCCGGACCTGTTCGTGCGTGAAACTGCAGCGGGCGATGTGCACCGCGCAGGGATTTTCGTGGGGc                             <  6:102927/149‑1 (MQ=255)
                                                                                                       gcgcAGTTCTTCCCCACCAGTGCGCGTGGCCGCGCGCTGCCCCCGGCACAGCTGGAGGCGCCCAGCCAGCAACTGATGATCGGCCCGGACCTGTTCGTGCGTGAAACTGCAGCGGGCGATGTGCACCGCGCAGGGATTTTCGTGGGGca                            >  5:334052/1‑149 (MQ=255)
                                                                                                        cgcAGTTCTTCCCCACCAGTGCGCGTGGCCGCGCGCTGCCCCCGGCACAGCTGGAGGCGCCCAGCCAGCAACTGATGATCGGCCCGGACCTGTTCGTGCGTGAAACTGCAGCGGGCGATGTGCACCGCGCAGGGATTTTCGTGGGGcac                           <  8:99146/149‑1 (MQ=255)
                                                                                                         gcAGTTCTTCCCCACCAGTGCGCGTGGCCGCGCGCTGCCCCCGGCACAGCTGGAGGCGCCCAGCCAGCAACTGATGATCGGCCCGGACCTGTTCGTGCGTGAAACTGCAGCGGGCGATGTGCACCGCGCAGGGATTTTCGTGGGGcac                           >  3:1970/1‑148 (MQ=255)
                                                                                                           aGTTCTTCCCCACCAGTGCGCGTGGCCGCGCGCTGCCCCCGGCACAGCTGGAGGCGCCCAGCCAGCAACTGATGATCGGCCCGGACCTGTTCGTGCGTGAAACTGCAGCGGGCGATGTGCACCGCGCAGGGATTTTCGTGGGGCacaa                         >  5:23687/1‑148 (MQ=255)
                                                                                                                ttcCCCACCAGTGCGCGTGGCCGCGCGCTGCCCCCGGCACAGCTGGAGGCGCCCAGCCAGCAACTGATGATCGGCCCGGACCTGTTCGTGCGTGAAACTGCAGCGGGCGATGTGCACCGCGCAGGGATTTTCGTGGGGCACAACAACCt                   >  1:68906/1‑149 (MQ=255)
                                                                                                                      accaGTGCGCGTGGCCGCGCGCTGCCCCCGGCACAGCTGGAGGCGCCCAGCCAGCAACTGATGATCGGCCCGGACCTGTTCGTGCGTGAAACTGCAGCGGGCGATGTGCACCGCGCAGGGATTTTCGTGGGGCACAACAACCTGCAAag             >  4:2541/1‑149 (MQ=255)
                                                                                                                             cgcgTGGCCGCGCGCTGCCCCCGGCACAGCTGGAGGCGCCCAGCCAGCAACTGATGATCGGCCCGGACCTGTTCGTGCGTGAAACTGCAGCGGGCGATGTGCACCGCGCAGGGATTTTCGTGGGGCACAACAACCTGCAAAGCAGCtt       <  4:120646/148‑1 (MQ=255)
                                                                                                                             cgcgTGGCCGCGCGCTGCCCCCGGCACAGCTGGAGGCGCCCAGCCAGCAACTGATGATCGGCCCGGACCTGTTCGTGCGTGAAACTGCAGCGGGCGATGTGCACCGCGCAGGGATTTTCGTGGGGCACAACAACCTGCAAAGCAGCTTc      <  1:7134/149‑1 (MQ=255)
                                                                                                                                gTGGCCGCGCGCTGCCCCCGGCACAGCTGGAGGCGCCCAGCCAGCAACTGATGATCGGCCCGGACCTGTTCGTGCGTGAAACTGCAGCGGGCGATGTGCACCGCGCAGGGATTTTCGTGGGGCACAACAACCTGCAAAGCAGCTTCaga   <  3:2541/149‑1 (MQ=255)
                                                                                                                                 tGGCCGCGCGCTGCCCCCGGCACAGCTGGAGGCGCCCAGCCAGCAACTGATGATCGGCCCGGACCTGTTCGTGCGTGAAACTGCAGCGGGCGATGTGCACCGCGCAGGGATTTTCGTGGGGCACAACAACCTGCAAAGCAGCTTCagag  <  3:130400/149‑1 (MQ=255)
                                                                                                                                          |                                                                                                                                           
GCTCAGTGCCAGTGGGTTGAGTCTCGGCACCCTCAATGCGTTTTCGCAGCAGATGAGCGAACGCCACCAGCAGCCGGCACCGGACCTGATTGCCAGTCAGTGGGCGCAGTTCTTCCCCACCAGTGCGCGTGGCCGCGCGCTGCCCCCGGCACAGCTGGAGGCGCCCAGCCAGCAACTGATGATCGGCCCGGACCTGTTCGTGCGTGAAACTGCAGCGGGCGATGTGCACCGCGCAGGGATTTTCGTGGGGCACAACAACCTGCAAAGCAGCTTCAGAG  >  NC_002947/4576882‑4577159

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