Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,332,994 T→A 7.8% P266P (CCT→CCA)  PP_3804 → hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,332,9940TA 7.8% 144.8 / 3.2 51P266P (CCT→CCA) PP_3804hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (24/23);  new base A (2/2);  total (26/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 = 8.08e-01

GCAGCCCTCAAGGCGCTGGGGGAGAACCTGAAGGCGCAGGACGTGGCGCTGGTGCTGGACCATCGTCAGCCTGAATCGGCAGTGGCCGAGGTGATCAAGGCGGCCGGGGCCAAGCTGGTGGTGGTGGAGAGTGACCCTGAAGATGCGATTGCGGGGCTCATAGCCAGCGTTGACCAAGTGGTTGGGGCATTGGGCGGGCGCTGATCGGACAGGCCCGCGCGGTGGATGGCACGGGCTTCGCCCGTGTTCGCGGGTAAACCCGCTCAAGGATCGCGGTTGCTGC  >  NC_002947/4332857‑4333139
                                                                                                                                         |                                                                                                                                                 
gCAGCCCTCAAGGCGCTGGGGGAGACCCTGAAGGCGCAGGACGTGGCGCTGGTGCTGGACCATCGTCAGCCTGAATCGGCAGTGGCCGAGGTGATCAAGGCGGCCGGGGCCAAGCTGGTGGTGGTGGAGAGTGACCCTGAAGATGCGAt                                                                                                                                        >  5:279913/1‑149 (MQ=255)
    agcTCAATGCGCGGGGGGAGAACCTGAAGGCGCAGGACGTGGCGCTGGTGCTGGACCTTCGTCAGCCTGAATCGGCAGTGGCCGAGGTGATCAAGGCGGCCGGGTCCAAGCTGGTGGTGGTGGAGAGTGACCCTGAAGATGCGATTGCg                                                                                                                                    <  4:130196/147‑1 (MQ=255)
     ccTCAAGGCGCTGGGGGAGAACCTGAAGGCGCAGGACGTGGCGCTGGTGCTGGACCATCGTCAGCCTGAATCGGCAGTGGCCGAGGTGATCAAGGCGGCCGGGGCCAAGCTGGTGGTGGTGGAGAGTGACCCAGAAGATGCGATTGCgg                                                                                                                                   >  7:246066/1‑149 (MQ=255)
     ccTCAAGGCGCTGGGGGAGAACCTGAAGGCGCAGGACGTGGCGCTGGTGCTGGACCATCGTCAGCCTGAATCGGCAGTGGCCGAGGTGATCAAGGCGGCCGGGGCCAAGCTGGTGGTGGTGGAGAGTGACCCAGAAGATGCGATTGCgg                                                                                                                                   >  3:325164/1‑149 (MQ=255)
      cTCAAGGCGCTGGGGGAGAACCTGAAGGCGCAGGACGTGGCGCTGGTGCTGGACCATCGTCAGCCTGAATCGGCAGTGGCCGAGGTGATCAAGGCGGCCGGGGCCAAGCTGGTGGTGGTGGAGAGTGACCCTGAAGATGCGATTGCggg                                                                                                                                  >  6:126409/1‑149 (MQ=255)
      cTCAAGGCGCTGGGGGAGAACCTGAAGGCGCAGGACGTGGCGCTGGTGCTGGACCATCGTCAGCCTGAATCGGCAGTGGCCGAGGTGATCAAGGCGGCCGGGGCCAAGCTGGTGGTGGTGGAGAGTGACCCTGAAGATGCGATTGCggg                                                                                                                                  >  8:103341/1‑149 (MQ=255)
                  ggggAGAACCTGAAGGCGCAGGACGTGGCGCTGGTGCTGGACCATCGTCAGCCTGAATCGGCAGTGGCCGAGGTGATCAAGGCGGCCGGGGCCAAGCTGGTGGTGGTGGAGAGTGACCCTGAAGATg                                                                                                                                            >  6:135080/1‑127 (MQ=255)
                  ggggAGAACCTGAAGGCGCAGGACGTGGCGCTGGTGCTGGACCATCGTCAGCCTGAATCGGCAGTGGCCGAGGTGATCAAGGCGGCCGGGGCCAAGCTGGTGGTGGTGGAGAGTGACCCTGAAGATg                                                                                                                                            <  5:135080/127‑1 (MQ=255)
                     gagaACCTGAAGGCGCAGGACGTGGCGCTGGTGCTGGACCATCGTCAGCCTGAATCGGCAGTGGCCGAGGTGATCAAGGCGGCCGGGGCCAAGCTGGTGGTGGTGGAGAGTGACCCTGAAGATGCGATTGCGGGGCTCATAGCCAGCGt                                                                                                                   >  1:148187/1‑149 (MQ=255)
                          ccTGAAGGCGCAGGACGTGGCGCTGGTGCTGGACCATCGTCAGCCTGAATCGGCAGTGGCCGAGGTGATCAAGGCGGCCGGGGCCAAGCTGGTGGTGGTGGAGAGTGACCCTGAAGATGCGATTGCGGGGCTCATAGCCAGCGTTGAc                                                                                                               >  5:242988/1‑148 (MQ=255)
                                   gcATGACGTGGCGCTGGTGCTGGACAATCGTCAGCCTGAATCGGCAGTGGACGAGGTGATCAAGGCGGCCGGGGCCAAGCTGGTGGTTGTGTAGAGTGACCCTGAAGATGCGATTGCGGGGCTCATAGCCAGCGTTg                                                                                                                 >  4:205409/1‑137 (MQ=255)
                                   gcAGGACGTGGCGCTGGTGCTGGACCATCGTCAGCCTGAATCGGCAGTGGCCGAGGTGATCAAGGCGGCCGGGGCCAAGCTGGTGGTGGTGGAGAGTGACCCTGAAGATGCGATTGCGGGGCTCATAGCCAGCGTTGACCAAGTGGTTg                                                                                                     >  7:51570/1‑149 (MQ=255)
                                   gcAGGACGTGGCGCTGGTGCTGGACCATCGTCAGCCTGAATCGGCAGTGGCCGAGGTGATCAAGGCGGCCGGGGCCAAGCTGGTGGTGGTGGAGAGTGACCCTGAAGATGCGATTGCGGGGCTCATAGCCAGCGTTGACCAAGTGGTTg                                                                                                     >  2:213702/1‑149 (MQ=255)
                                   gcAGGACGTGGCGATGGTGCTGGACAATCGTCAGCCTGACTCGGCAGTGGACGAGGTGATCAAGGCGGCCGGGGCCAAGATGTTGTTGGTGGAGAGTGACACTGAAGATGCGATTGCGGGGCTCATATCCAGCGTTg                                                                                                                 <  3:205409/137‑1 (MQ=255)
                                     aGGACGTGGCGCTGGTGCTGGACCATCGTCAGCCTGAATCGGCAGTGGCCGAGGTGATCAAGGCGGCCGGGGCCAAGCTGGTGGTGGTGGAGAGTGACCCTGAAGATGCGATTGCGGGGCTCATAGCCAGCGTTGACCAAGTGGTTggg                                                                                                   <  8:51570/149‑1 (MQ=255)
                                          gTGGCGCTGGTGCTGGACCATCGTCAGCCTGAATCGGCAGTGGCCGAGGTGATCAAGGCGGCCGGGGCCAAGCTGGTGGTGGTGGAGAGTGACCCTGAAGATGCGATTGCGGGGCTc                                                                                                                              <  3:26822/117‑1 (MQ=255)
                                          gTGGCGCTGGTGCTGGACCATCGTCAGCCTGAATCGGCAGTGGCCGAGGTGATCAAGGCGGCCGGGGCCAAGCTGGTGGTGGTGGAGAGTGACCCTGAAGATGCGATTGCGGGGCTc                                                                                                                              >  4:26822/1‑117 (MQ=255)
                                          gTGGCGCTGGTGCTGGACCATCGTCAGCCTGAATCGGCAGTGGCCGAGGTGATCAAGGCGGCCGGGGCCAAGCTGGTGGTGGTGGAGAGTGACCCTGAAGATGCGATTGCGGGGCTCATAGCCAGCGTTg                                                                                                                 >  2:227866/1‑130 (MQ=255)
                                          gTGGCGCTGGTGCTGGACCATCGTCAGCCTGAATCGGCAGTGGCCGAGGTGATCAAGGCGGCCGGGGCCAAGCTGGTGGTGGTGGAGAGTGACCCTGAAGATGCGATTGCGGGGCTCATAGCCAGCGTTGACCAAGTGGTTGGGGCAtt                                                                                              >  6:327276/1‑149 (MQ=255)
                                               gcTGGTGCTGGACCATCGTCAGCCTGAATCGGCAGTGGCCGAGGTGATCAAGGCGGCCGGGGCCAAGCTGGTGGTGGTGGAGAGTGACCCTGAAGATGCGATTGCGGGGCTCATAGCCAGCGTTGACCAAGTGGTTGGGGCATTgggcg                                                                                         >  5:107171/1‑149 (MQ=255)
                                                cTGGTGCTGGACCATCGTAAGCCTGAATCGGCAGTGGCCGAGGTGATCAAGGCGGCCGGGGCCAAGCTGGTGGTGGTGGAGAGTGACCCTGAAGATGCGATTGCGGGGCTCATAGCCAGCGTTGACCAAGTGGTTGGGGCATTgggcgg                                                                                        <  6:264358/149‑1 (MQ=255)
                                                   gTGCTGGACCATCGTCAGCCTGAATCGGCAGTGGCCGAGGTGATCAAGGCGGCCGGGGCCAAGCTGGTGGTGGTGGAGAGTGACCCTGAAGATGCGATTGCGGGGCTCATAGCCAGCGTTGACCAAGTGGTTGGGGCATTgggcgggc                                                                                      <  5:126409/148‑1 (MQ=255)
                                                   gTGCTGGACCATCGTCAGCCTGAATCGGCAGTGGCCGAGGTGATCAAGGCGGCCGGGGCCAAGCTGGTGGTGGTGGAGAGTGACCCTGAAGATGCGATTGCGGGGCTCATAGCCAGCGTTGACCAAGTGGTTGGGGCATTgggcggg                                                                                       <  7:103341/147‑1 (MQ=255)
                                                   gTGCTGGACCATCGTCAGCCTGAATCGGCAGTGGCCGAGGTGATCAAGGCGGCCGGGGCCAAGCTGGTGGTGGTGGAGAGTGACCCAGAAGATGCGATTGCGGGGCTCATAGCCAGCGTTGACCAAGTGGTTGGGGCATTgggcggg                                                                                       <  4:325164/147‑1 (MQ=255)
                                                   gTGCTGGACCATCGTCAGCCTGAATCGGCAGTGGCCGAGGTGATCAAGGCGGCCGGGGCCAAGCTGGTGGTGGTGGAGAGTGACCCAGAAGATGCGATTGCGGGGCTCATAGCCAGCGTTGACCAAGTGGTTGGGGCATTgggcggg                                                                                       <  8:246066/147‑1 (MQ=255)
                                                         gACCATCGTCAGCCTGAATCGGCAGTGGCCGAGGTGATCAAGGCGGCCGGGGCCAAGCTGGTGGTGGTGGAGAGTGACCCTGAAGATGCGATTGCGGGGCTCATAGCCAGCGTTGACCAAGTGGTTGGGGCATTGGGCGGGCGCTGATc                                                                               >  7:284546/1‑149 (MQ=255)
                                                            cATCGTCAGCCTGAATCGGCAGTGGCCGAGGTGATCAAGGCGGCCGGGGCCAAGCTGGTGGTGGTGGAGAGTGACCCTGAAGATGCGATTGCGGGGCTCATAGCCAGCGTTGACCAAGTGGTTGGGGCATTGGGCGGGCGCTGATCGGa                                                                            >  6:177631/1‑149 (MQ=255)
                                                             aTCGTCAGCCTGAATCGGCAGTGGCCGAGGTGATCAAGGCGGCCGGGGCCAAGCTGGTGGTGGTGGAGAGTGACCCTGAAGATGCGATTGCGGGGCTCATAGCCAGCGTTGACCAAGTGGTTGGGGCATTGGGCGGGCGCTGATCGGAc                                                                           <  7:27779/149‑1 (MQ=255)
                                                                   aGCCTGAATCGGCAGTGGCCGAGGTGATCAAGGCGGCCGGGGCCAAGCTGGTGGTGGTGGAGAGTGACCCTGAAGATGCGATTGCGGGGCTCATAGCCAGCGTTGACCAAGTGGTTGGGGCATTGGGCGGGCGCTGATCGGACAGGccc                                                                     >  3:293726/1‑149 (MQ=255)
                                                                     ccTGAATCGGCAGTGGCCGAGGTGATCAAGGCGGCCGGGGCCAAGCTGGTGGTGGTGGAGAGTGACCCTGAAGATGCGATTGCGGGGCTCATAGCCAGCGTTGACCAAGTGGTTggg                                                                                                   <  2:100926/117‑1 (MQ=255)
                                                                     ccTGAATCGGCAGTGGCCGAGGTGATCAAGGCGGCCGGGGCCAAGCTGGTGGTGGTGGAGAGTGACCCTGAAGATGCGATTGCGGGGCTCATAGCCAGCGTTGACCAAGTGGTTGGGGCATTGGCCGGGCGCTGATCGGACAGGCCcgc                                                                   >  7:261069/1‑149 (MQ=255)
                                                                         aaTCGGCAGTGGCCGAGGTGATCAAGGCGGCCGGGGCCAAGCTGGTGGTGGTGGAGAGTGACCCTGAAGATGCGATTGCGGGGCTCATAGCCAGCGTTGACCAAGTGGTTGGGGCATTGGGCGGGCGCTGATCGGACAGGCCcgcgc                                                                 <  6:107171/147‑1 (MQ=255)
                                                                          aTCGGCAGTGGCCGAGGTGATCAAGGCGGCCGGGGCCAAGCTGGTGGTGGTGGAGAGTGACCCTGAAGATGCGATTGCGGGGCTCATAGCCAGCGTTGACCAAGTGGTTGGGGCATTGGGCGGGCGCTGATCGGACAGGCCCGCGCGGt                                                              <  5:327276/149‑1 (MQ=255)
                                                                              gCAGTGGCCGAGGTGATCAAGGCGGCCGGGGCCAAGCTGGTGGTGGTGGAGAGTGACCCTGAAGATGCGATTGCGGGGCTCATAGCCAGCGTTGACCAAGTGGTTGGGGCATTGGGCGGGCGCTGATCGGACAGGCCCGCGCGGTGGAt                                                          >  8:190830/1‑149 (MQ=255)
                                                                              gCAGTGGCCGAGGTGATCAAGGCGGCCGGGGCCAAGCTGGTGGTGGTGGAGAGTGACCCTGAAGATGCGATTGCGGGGCTCATAGCCAGCGTTGACCAAGTGGTTGGGGCATTGGGCGGGCGCTGATCGGACAGGCCCGCGCGGTGGAt                                                          >  3:169854/1‑149 (MQ=255)
                                                                                      cGAGGTGATCAAGGCGGCCGGGGCCAAGCTGGTGGTGGTGGAGAGTGACCCTGAAGATGCGATTGCGGGGCTCATAGCCAGCGTTGACCAAGTGGTTGGGGCATTGGGCGGGCGCTGATCGGACAGGCCCGCGCGGt                                                              >  1:306445/1‑137 (MQ=255)
                                                                                      cGAGGTGATCAAGGCGGCCGGGGCCAAGCTGGTGGTGGTGGAGAGTGACCCTGAAGATGCGATTGCGGGGCTCATAGCCAGCGTTGACCAAGTGGTTGGGGCATTGGGCGGGCGCTGATCGGACAGGCCCGCGCGGt                                                              <  2:306445/137‑1 (MQ=255)
                                                                                           tGATCAAGGCGGCCGGGGCCAAGCTGGTGGTGGTGGAGAGTGACCCTGAAGATGCGATTGCGGGGCTCATAGCCAGCGTTGACCAAGTGGTTGGGGCATTGGGCGGGCGCTGATCgg                                                                             >  4:333850/1‑117 (MQ=255)
                                                                                           tGATCAAGGCGGCCGGGGCCAAGCTGGTGGTGGTGGAGAGTGACCCTGAAGATGCGATTGCGGGGCTCATAGCCAGCGTTGACCAAGTGGTTGGGGCATTGGGCGGGCGCTGATCgg                                                                             <  3:333850/117‑1 (MQ=255)
                                                                                             aTCAAGGCGGCCGGGGCCAAGCTGGTGGTGGTGGAGAGTGACCCTGAAGATGCGATTGCGGGGCTCATAGCCAGCGTTGACCAAGTGGTTGGGGCATTGGGCGGGCGCTGATCGGACAGGCCCGCGCGGTGGATGGCACGGGCTTCGc                                            <  6:242988/148‑1 (MQ=255)
                                                                                                 aGGCGGCCGGGGCCAATCTGGTGGTGGTGGAGAGTGACCCTGAAGATGCGATTGCGGGGCTCATAGCCAGCGTTGACCAAGTGGTTGGGGCATTGGGCGGGCGCTGATCGGACAGGCCCGCGCGGTGGATGGCACGGGCTTCGCCCgtg                                       <  6:279913/149‑1 (MQ=255)
                                                                                                 aGGCGGCCGGGGCCAAGCTGGTGGTGGTGGAGAGTGACCCTGAAGATGCGATTGCGGGGCTCATAGCCAGCGTTGACCAAGTGGTTGGGGCATTGGGCGGGCGCTGATCGGACAGGCCCGCGCGGTGGATGGCACGGGCTTCGCCCgtg                                       >  7:249759/1‑149 (MQ=255)
                                                                                                    cggcCGGGGCCAAGCTGGTGGTGGTGGAGAGTGACCCTGAAGATGCGATTGCGGGGCTCATAGCCAGCGTTGACCAAGTGGTTGGGGCATTGGGCGGGCGCTGATCGGACAGGCCCGCGCGGTGGATGGCACGGGCGAAGCCCgtgcca                                    <  1:213702/149‑4 (MQ=255)
                                                                                                        cGGGGCCAAGCTGGTGGTGGTGGAGAGTGACCCTGAAGATGCGATTGCGGGGCTCATAGCCAGCGTTGACCAAGTGGTTGGGGCATTGGGCGGGCGCTGATCGGACAGGCCCGCGCGGTGGATGGCACGGGCTTCGCCCGTGTTCGCgg                                <  5:177631/149‑1 (MQ=255)
                                                                                                                 gCTGGTGGTGGTGGAGAGTGACCCTGAAGATGCGATTGCGGGGCTCATAGCCAGCGTTGACCAAGTGGTTGGGGCATTGGGCGGGCGCTGATCGGACAGGCCCGCGCGGTGGATGGCACGGGCTTCGCCCgt                                        >  4:150286/1‑132 (MQ=255)
                                                                                                                 gCTGGTGGTGGTGGAGAGTGACCCTGAAGATGCGATTGCGGGGCTCATAGCCAGCGTTGACCAAGTGGTTGGGGCATTGGGCGGGCGCTGATCGGACAGGCCCGCGCGGTGGATGGCACGGGCTTCGCCCgt                                        <  3:150286/132‑1 (MQ=255)
                                                                                                                        gtggtggAGAGTGACCCTGAAGATGCGATTGCGGGGCTCATAGCCAGCGTTGACCAAGTGGTTGGGGCATTGGGCGGGCGCTGATCGGACAGGCCCGCGCGGTGGATGGCACGGGCTTCGCCCGTGTTCGCGGGTAAACCCGCTCAAg                 <  7:190830/148‑1 (MQ=255)
                                                                                                                                      cccTGAAGATGCGATTGCGGGGCTCATAGCCAGCGTTGACCATGTGGTTGGGGCATTGGGCGGGCGCTGATCGGACAGGCCCGCGCGGTGGATGGCACGGGCTTCGCCCGTGTTCGCGGGTAAACCCGCTCAAGGATCGCGGTtgctgc  <  4:293726/149‑1 (MQ=255)
                                                                                                                                      cccTGAAGATGCGATTGCGGGGCTCATAGCCAGCGTTGACCAAGTGGTTGGGGCATTGGGCGGGCGCTGATCGGACAGGCCCGCGCGGTGGATGGCACGGGCTTCGCCCGTGTTCGCGGGTAAACCCGCTCAAGGATCGCGGTtgctgc  <  2:148187/149‑1 (MQ=255)
                                                                                                                                       ccTGAAGATGCGATTGCGGGGCTCATAGCCAGCGTTGACCAAGTGGTTGGGGCATTGGGCGGGCGCTGATCGGACAGGCCCGCGCGGTGGATGGCACGGGCTTCGCCCGTGTTCGCGGGTAAACCCGCTc                    <  6:333601/130‑1 (MQ=255)
                                                                                                                                       ccTGAAGATGCGATTGCGGGGCTCATAGCCAGCGTTGACCAAGTGGTTGGGGCATTGGGCGGGCGCTGATCGGACAGGCCCGCGCGGTGGATGGCACGGGCTTCGCCCGTGTTCGCGGGTAAACCCGCTc                    >  5:333601/1‑130 (MQ=255)
                                                                                                                                         |                                                                                                                                                 
GCAGCCCTCAAGGCGCTGGGGGAGAACCTGAAGGCGCAGGACGTGGCGCTGGTGCTGGACCATCGTCAGCCTGAATCGGCAGTGGCCGAGGTGATCAAGGCGGCCGGGGCCAAGCTGGTGGTGGTGGAGAGTGACCCTGAAGATGCGATTGCGGGGCTCATAGCCAGCGTTGACCAAGTGGTTGGGGCATTGGGCGGGCGCTGATCGGACAGGCCCGCGCGGTGGATGGCACGGGCTTCGCCCGTGTTCGCGGGTAAACCCGCTCAAGGATCGCGGTTGCTGC  >  NC_002947/4332857‑4333139

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 21 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑