Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 874,538 C→T 9.5% V39V (GTC→GTT PP_0759 → hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947874,5380CT 9.5% 143.6 / 3.9 52V39V (GTC→GTTPP_0759hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (23/24);  new base T (3/2);  total (26/26)
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

AGCCACTGGCGACTCGGCCTAGGATACGGCGCATGTGGATAGGGAGTGTTTTCATGCGTAAGTCTTTTGTCATCAGCCTGCTGAGTGCTGGCATCCTGCTGGCCGGCTGCCAGGCGGTGAATACCACCAGCGGCGGTGCTGTTGGCGTCGAACGCCAACAGTACATGTTCAGCATGCTCTCGACCGATGAGGTCAACCAGATGTACGCGCAGTCGTATCAGCAGACCCTGGGTGAAGCATCGAGCAAGGGTGTGCTCGACAAGTCCAGCGCCGATGCCAAGCGTGTGCAGGCCATT  >  NC_002947/874390‑874685
                                                                                                                                                    |                                                                                                                                                   
aGCCACTGGCGACTCGGCCTAGGATACGGCGCATGTGGATAGGGAGTGTTTTCATGCGTAAGTCTTTTGTCATCAGCCTGCTGAGTGCTGGCATCCTGCTGGCCGGCTGCCAGGCGGTGAATACCACCAGCGGCGGTGCTGTTGGCGTc                                                                                                                                                     >  1:289845/1‑149 (MQ=255)
         cGACTCGGCCTAGGATACGGCGCATGTGGATAGGGAGTGTTTTCATGCGTAAGTCTTTTGTCATCAGCCTGCTGAGTGCTGGCATCCTGCTGGCCGGCTGCCAGGCGGTGAATACCACCAGCGGCGGTGCTGTTGGCGTCGAACGCCaa                                                                                                                                            <  5:88689/149‑1 (MQ=255)
           aCTCGGCCTAGGATACGGCGCATGTGGATAGGGAGTGTTTTCATGCGTAAGTCTTTTGTCATCAGCCTGCTGAGTGCTGGCATCCTGCTGGCCGGCTGCCAGGCGGTGAATACCACCAGCGGCGGTGCTGTTGGCGTTCGACGCCAACa                                                                                                                                          <  7:220489/149‑1 (MQ=255)
             tCGGCCTAGGATACGGCGCATGTGGATAGGGAGTGTTTTCATGCGTAAGTCTTTTGTCATCAGCCTGCTGAGTGCTGGCATCCTGCTGGCCGGCTGCCAGGCGGTGAATACCACCAGCGGCGGTGCTGTTGGCGTTCGACGCCAAc                                                                                                                                           >  3:642/1‑146 (MQ=255)
             tCGGCCTAGGATACGGCGCATGTGGATAGGGAGTGTTTTCATGCGTAAGTCTTTTGTCATCAGCCTGCTGAGTGCTGGCATCCTGCTGGCCGGCTGCCAGGCGGTGAATACCACCAGCGGCGGTGCTGTTGGCGTTCGACGCCAAc                                                                                                                                           <  4:642/146‑1 (MQ=255)
                  cTAGGATACGGCGCATGTGGATAGGGAGTGTTTTCATGCGTAAGTCTTTTGTCATCAGCCTGCTGAGTGCTGGCATCCTGCTGGCCGGCTGCCAGGCGGTGAATACCACCAGCGGCGGTGCTGTTGGCGTCGAACGCCAACAGt                                                                                                                                        <  8:274477/144‑1 (MQ=255)
                  cTAGGATACGGCGCATGTGGATAGGGAGTGTTTTCATGCGTAAGTCTTTTGTCATCAGCCTGCTGAGTGCTGGCATCCTGCTGGCCGGCTGCCAGGCGGTGAATACCACCAGCGGCGGTGCTGTTGGCGTCGAACGCCAACAGt                                                                                                                                        >  7:274477/1‑144 (MQ=255)
                     ggATACGGCGCATGTGGATAGGGAGTGTTTTCATGCGTAAGTCTTTTGTCATCAGCCTGCTGAGTGCTGGCATCCTGCTGGCCGGCTGCCAGGCGGTGAATACCACCAGCGGCGGTGCTGTTGGCGTCGAACGCCAACAGTACATGt                                                                                                                                  <  7:129187/147‑1 (MQ=255)
                     ggATACGGCGCATGTGGATAGGGAGTGTTTTCATGCGTAAGTCTTTTGTCATCAGCCTGCTGAGTGCTGGCATCCTGCTGGCCGGCTGCCAGGCGGTGAATACCACCAGCGGCGGTGCTGTTGGCGTCGAACGCCAACAGTACATGt                                                                                                                                  <  8:204708/147‑1 (MQ=255)
                                  gtgGATAGGGAGTGTTTTCATGCGTAAGTCTTTTGTCATCAGCCTGCTGAGTGCTGGCATCCTGCTGGCCGGCTGCCAGGCGGTGAATACCACCAGCGGCGGTGCTGTTGGCGTCGAACGCCAACAGTACATGTTCAGCATGCTCTCGa                                                                                                                   >  8:23738/1‑149 (MQ=255)
                                        aGGGAGTGTTTTCATGCGTAAGTCTTTTGTCATCAGCCTGCTGAGTGCTGGCATCCTGCTGGCCGGCTGCCAGGCGGTGAATACCACCAGCGGCGGTGCTGTTGGCGTCGAACGCCAACAGTACATGTTCAGCATGATCTCGACCGAt                                                                                                              >  6:291346/1‑148 (MQ=255)
                                           gAGTGTTTTCATGCGTAAGTCTTTTGTCATCAGCCTGCTGAGTGCTGGCATCCTGCTGGCCGGCTGCCAGGCGGTGAATACCACCAGCGGCGGTGCTGTTGGCGTCGAACGCCAACAGTACATGTTCAGCATGCTCTCGACCGATGAgg                                                                                                          >  1:312164/1‑149 (MQ=255)
                                                     aTGCGTCAGTCTTTTGTCATCAGCCTGCTGAGTGCTGGCATCCTGCTGGCCGGCTGCCAGGCGGTGAATACCACCAGCGGCGGTGCTGTTGGCGTCGAACGCCAACAGTACATGTTCAGCATGCTCTCGACCGATGAGGTCAACCAGAt                                                                                                <  4:215850/149‑1 (MQ=255)
                                                     aTGCGTAAGTCTTTTGTCATCAGCCTGCTGAGTGCTGGCATCCTGCTGGCCGGCTGCCAGGCGGTGAATACCACCAGCGGCGGTGCTGTTGGCGTCGAACGCCAACAGTACATGTTCAGCATGCTCTCGACCGATGAGGTCAACCAGa                                                                                                 <  3:17887/148‑1 (MQ=255)
                                                         gTAAGTCTTTTGTCATCAGCCTGCTGAGTGCTGGCATCCTGCTGGCCGGCTGCCAGGCGGTGAATACCACCAGCGGCGGTGCTGTTGGCGTCGAACGCCAACAGTACATGTTCAGCATGCTCTCGACCGATGAGGTCAACCAGATGTAc                                                                                            >  4:79180/1‑149 (MQ=255)
                                                                ttttGTCATCAGCCTGCTGAGTGCTGGCATCCTGCTGGCCGGCTGCCAGGCGGTGAATACCACCAGCGGCGGTGCTGTTGGCGTCGAACGCCAACAGTACATGTTCAGCATGCTCTCGACCGATGa                                                                                                            <  3:125515/126‑1 (MQ=255)
                                                                ttttGTCATCAGCCTGCTGAGTGCTGGCATCCTGCTGGCCGGCTGCCAGGCGGTGAATACCACCAGCGGCGGTGCTGTTGGCGTCGAACGCCAACAGTACATGTTCAGCATGCTCTCGACCGATGa                                                                                                            >  4:125515/1‑126 (MQ=255)
                                                                ttttGTCATCAGCCTGCTGAGTGCTGGCATCCTGCTGGCCGGCTGCCAGGCGGTGAATACCACCAGCGGCGGTGCTGTTGGCGTCGAACGCCAACAGTACATGTTCAGCATGCTCTCGACCGATGAg                                                                                                           >  6:20049/1‑127 (MQ=255)
                                                                ttttGTCATCAGCCTGCTGAGTGCTGGCATCCTGCTGGCCGGCTGCCAGGCGGTGAATACCACCAGCGGCGGTGCTGTTGGCGTCGAACGCCAACAGTACATGTTCAGCATGCTCTCGACCGATGAg                                                                                                           <  5:20049/127‑1 (MQ=255)
                                                                ttttGTCATCAGCCTGCTGAGTGCTGGCATCCTGCTGGCCGGCTGCCAGGCGGTGAATACCACCAGCGGCGGTGCTGTTGGCGTCGAACGCCAACAGTACATGTTCAGCATGCTCTCGACCGATGAGGTCAACCAGATGTACGCGCAGt                                                                                     >  5:84228/1‑149 (MQ=255)
                                                                         caGCCTGCTGAGTGCTGGCATCCTGCTGGCCGGCTGCCAGGCGGTGAATACCACCAGCGGCGGTGCTGTTGGCGTCGAACGCCAACAGTACATGTTCAGCATGCTCTCGACCGATGAGGTCAACCAGATGTACGCGCAGTCGTAt                                                                                >  4:194657/1‑145 (MQ=255)
                                                                         caGCCTGCTGAGTGCTGGCATCCTGCTGGCCGGCTGCCAGGCGGTGAATACCACCAGCGGCGGTGCTGTTGGCGTCGAACGCCAACAGTACATGTTCAGCATGCTCTCGACCGATGAGGTCAACCAGATGTACGCGCAGTCGTAt                                                                                <  3:194657/145‑1 (MQ=255)
                                                                                ctgAGTGCTGGCATCCTGCTGGCCGGCTGCCAGGCGGTGAATACCACCAGCGGCGGTGCTGTTGGCGTCGAACGCCAACAGTACATGTTCAGCATGCTCTCGACCGATGAGGTCAACCAGATGTACGCGCAGTCGTATCAGCAGAcc                                                                       <  2:289845/147‑1 (MQ=255)
                                                                                ctgAGTGCTGGCATCCTGCTGGCCGGCTGCCAGGCGGTGAATACCACCAGCGGCGGTGCTGTTGGCGTCGAACGCCAACAGTACATGTTCAGCATGCTCTCGACCGATGAGGTCAACCAGATGTACGCGCAGTCGTATCAGCAGACCCt                                                                     >  2:24862/1‑149 (MQ=255)
                                                                                 tgAGTGCTGGCATCCTGCTGGCCGGCTGCCAGGCGGTGAATACCACCAGCGGCGGTGCTGTTGGCGTCGAACGCCAACAGTACATGTTCAGCATGCTCTCGACCGATGAGGTCAACCAGATGTACGCGCAGTCGTATCAGCAGACCCTg                                                                    >  2:149013/1‑149 (MQ=255)
                                                                                 tgAGTGCTGGCATCCTGCTGGCCGGCTGCCAGGCGGTGAATACCACCAGCGGCGGTGCTGTTGGCGTCGAACGCCAACAGTACATGTTCAGCATGCTCTCGACCGATGAGGTCAACCAGATGTACGCGCAGTCGTATCAGCAGACCCTg                                                                    >  1:1397/1‑149 (MQ=255)
                                                                                   aGTGCTGGCATCCTGCTGGCCGGCTGCCAGGCGGTTAATACCACCAGCGGCGGTGCTGTTGGCGTCGAACGCCAACAGTACATGTTCAGCATGCTCTCGACCGATGAGGTCAACCAGATGTACGCTCAGTCGTATCAGCAGACCCTgg                                                                   <  2:252376/148‑1 (MQ=255)
                                                                                    gTGCTGGCATCCTGCTGGCCGGCTGCCAGGCGGTGAATACCACCAGCGGCGGTGCTGTTGGCGTCGAACGCCAACAGTACATGTTCAGCATGCTCTCGACCGATGAGGTCAACCAGATGTACGCGCAGTCGTATCAGCAGACCCTggg                                                                  <  7:23738/148‑1 (MQ=255)
                                                                                       cTGGCAGCCGGCCAGCAGGATGCCAGGCGGTGAATACCACCAGCGGCGGTGCTGTTGGCGTCGAACGCCAACAGTACATGTTCAGCATGCTCTCGACCGATGAGGTCAACCAGATGTACGCGCAGTCGTATCAGCAGACCCTGGGTg                                                                >  8:267800/1‑147 (MQ=255)
                                                                                        tGGCATCCTGCTGGCCGGCTGCCAGGCGGTGAATACCACCAGCGGCGGTGCTGTTGGCGTCGAACGCCAACAGTACATGTTCAGCATGCTCTCGACCGATGAGGTCAACCAGATGTACGCGCAGTCGTATCAGCAGACCCTGGGt                                                                 >  6:302545/1‑145 (MQ=255)
                                                                                        tGGCATCCTGCTGGCCGGCTGCCAGGCGGTGAATACCACCAGCGGCGGTGCTGTTGGCGTCGAACGCCAACAGTACATGTTCAGCATGCTCTCGACCGATGAGGTCAACCAGATGTACGCGCAGTCGTATCAGCAGACCCTGGGt                                                                 <  5:302545/145‑1 (MQ=255)
                                                                                        tGGCATCCTGCTGGCCGGCTGCCAGGCGGTGAATACCACCAGCGGCGGTGCTGTTGGCGTCGAACGCCAACAGTACATGTTCAGCATGCTCTCGACCGATGAGGTCAACCAGATGTACGCGCAGTCGTATCAGCAGACCCTGGGTGAAg                                                             >  8:86445/1‑149 (MQ=255)
                                                                                        tGGCATCCTGCTGGCCGGCTGCCAGGCGGTGAATACCACCAGCGGCGGTGCTGTTGGCGTCGAACGCCAACAGTACATGTTCAGCATGCTCTCGACCGATGAGGTCAACCAGATGTACGCGCAGTCGTATCAGCAGACCCTGGGTGAAg                                                             >  2:69042/1‑149 (MQ=255)
                                                                                            aTCCTGCTGGCCGGCTGCCAGGCGGTGAATACCACCAGCGGCGGTGCTGTTGGCGTCGAACGCCAACAGTACATGTTCAGCATGCTCTCGACCGATGAGGTCAACCAGATGTACGCGCAGTCGTATCAGCAGACCCTGGGTGAAGCATc                                                         <  2:312164/149‑1 (MQ=255)
                                                                                                         gCTGCCAGGCGGTGAATACCACCAGCGGCGGTGCTGTTGGCGTCGAACGCCAACAGTACATGTTCAGCATGCTCTCGACCGATGAGGTCAACCAGATGTACGCGCAGTCGTATCAGCAGACCCTGGGTGAAGCATCGAGCAAGGGtgtg                                            >  7:172843/1‑149 (MQ=255)
                                                                                                          cTGCCAGGCGGTGAATACCACCAGCGGCGGTGCTGTTGGCGTCGAACGCCAACAGTACATGTTCAGCATGCTCTCGACCGATGAGGACAACCAGATGTACGCGCAGTCGTATCAGCAGACCCTGGGTGAAGCATCGAGCAAGGGTGTGc                                           <  8:15069/149‑1 (MQ=255)
                                                                                                              cAGGCGGTGAATACCACCAGCGGCGGTGCTGTTGGCGTCGAACGCCAACAGTACATGTTCAGCATGCTCTCGACCGATGAGGTCAACCAGATGTACGCGCAGTCGTATCAGCAGACCCTGGGTGAAGCATCGAGCAAGGGTGTGCt                                          >  8:210380/1‑146 (MQ=255)
                                                                                                              cAGGCGGTGAATACCACCAGCGGCGGTGCTGTTGGCGTCGAACGCCAACAGTACATGTTCAGCATGCTCTCGACCGATGAGGTCAACCAGATGTACGCGCAGTCGTATCAGCAGACCCTGGGTGAAGCATCGAGCAAGGGTGTGCt                                          <  7:210380/146‑1 (MQ=255)
                                                                                                                     tGAATACCACCAGTGGCGGTGCTGTTGGCGTCGAACGCCAACAGTACATGTTCAGCATGCTCTCGACCGATGAGGTCAACCAGATGTACGCGCAGTCGTATCAGCAGACCCTGGGTGAAGCATCGAGCAAGGGTGTGCTCGACAAGTc                                 >  7:19640/1‑148 (MQ=255)
                                                                                                                         tACCACCAGCGGCGGTGCTGTTGGCGTCGAACGCCAACAGTACATGTTCAGCATGCTCTCGACCGATGAGGTCAACCAGATGTACGCGCAGTCGTATCAGCAGACCCTGGGTGAAGCATCGAGCAAGGGTGTGCTCGACAAGTCCAgcg                            <  2:1397/149‑1 (MQ=255)
                                                                                                                            caccaGCGGCGGTGCTGTTGGCGTCGAACGCCAACAGTACATGTTCAGCATGCTCTCGACCGATGAGGTCAACCAGATGTACGCGCAGTCGTATCAGCAGACCCTGGGTGAAGCATCGAGCAAGGGTGTGCTCGACAAGTCCAGCGcc                          >  3:250601/1‑148 (MQ=255)
                                                                                                                             accaGCGGCGGTGCTGTTGGCGTCGAACGCCAACAGTACATGTTCAGCATGCTCTCGACCGATGAGGTCAACCAGATGTACGCGCAGTCGTATCAGCAGACCCTGGGTGAAGCATCGAGCAAGGGTGTGCTCGACAAGTCCAGCGCCGa                        <  3:79180/149‑1 (MQ=255)
                                                                                                                                     cggTGCTGTTGGCGTCGAACGCCAACAGTACATGTTCAGCATGCTCTCGACCGATGAGGTCAACCAGATGTACGCGCAGTCGTATCAGCAGACCCTGGGTGAAGCATCGAGCAAGGGTGTGCTCGACAAGTCCAGCGCCGATGCCAAGc                <  8:19640/149‑1 (MQ=255)
                                                                                                                                      tgtacTGTTGGCGTTCGACGCCAACGGTACATGTTCAGCATGCTCTCGACCGATGAGGTCAACCAGATGTACGCGCAGTCGTATCAGCAGACCCTGGGTGAAGCATCGAGCAAGGGTGTGCTCGAc                                      >  7:232526/5‑126 (MQ=255)
                                                                                                                                      ggTGCTGTTGGCGTCGAACGCCAACAGTACATGTTCAGCATGCTCTCGACCGATGAGGTCAACCAGATGTACGCGCAGTCGTATCAGCAGACCCTGGGTGAAGCATCg                                                        >  4:109877/1‑108 (MQ=255)
                                                                                                                                      ggTGCTGTTGGCGTCGAACGCCAACAGTACATGTTCAGCATGCTCTCGACCGATGAGGTCAACCAGATGTACGCGCAGTCGTATCAGCAGACCCTGGGTGAAGCATCg                                                        <  3:109877/108‑1 (MQ=255)
                                                                                                                                       gTGCTGTTGGCGTCGAACGCCAACAGTACATGTTCAGCATGCTCTCGACCGATGAGGTCAACCAGATGTACGCGCAGTCGTATCAGCAGACCCTGGGTGAAGCATCGAGCAAGGGTGTGCTCGACAAGTCCAGCGCCGATGCCAAGCgt              <  4:91344/149‑1 (MQ=255)
                                                                                                                                          cTGTTGGCGTTCGACGCCAACAGTACATGTTCAGCATGCTCTCGACCGATGAGGTCAACCAGATGTACGCGCAGTCGTATCAGCAGACCCTGGGTGAAGCATCGAGCAAGGGTGTGCTCGACAAGTCCAGCGCCGATGCCAAGCGtgtg           >  1:168370/1‑149 (MQ=255)
                                                                                                                                                gCGTCGAACGCCAACAGTACATGTTCAGCATGCTCTCGACCGATGAGGTCAACCAGATGTACGCGCAGTCGTATCAGCAGACCCTGGGTGAAGCATCGAGCaa                                                   >  6:94695/1‑103 (MQ=255)
                                                                                                                                                gCGTCGAACGCCAACAGTACATGTTCAGCATGCTCTCGACCGATGAGGTCAACCAGATGTACGCGCAGTCGTATCAGCAGACCCTGGGTGAAGCATCGAGCaa                                                   <  5:94695/103‑1 (MQ=255)
                                                                                                                                                 cGTCGAACGCCAACAGTACATGTTCAGCATGCTCTCGACCGATGAGGTCAACCAGATGTACGCGCAGTCGTATCAGCAGACCCTGGGTGAAGCATCGAGCAAGGGTGTGCTCGACAAGTCCAGCGCCGATGCCAAGCGTGTGCAGGCCa    >  7:8065/1‑149 (MQ=255)
                                                                                                                                                  gTCGAACGCCAACAGTACATGTTCAGCATGCTCTCGACCGATGAGGTCAACCAGATGTACGCGCAGTCGTATCAGCAGACCCTGGGTGAAGCATCGAGCAAGGGTGTGCTCGACAAGTCCAGCGCCGATGCCAAGCGTGTGCAGGCCAt   <  4:275676/149‑1 (MQ=255)
                                                                                                                                                   tCGAACGCCAACAGTACATGTTCAGCATGCTCTCGACCGATGAGGTCAACCAGATGTACGCGCAGTCGTATCAGCAGACCCTGGGTGAAGCATCGAGCAAGGGTGTGCTCGACAAGTCCAGCGCCGATGCCAAGCGTGTGCAGGCCAtt  <  6:84228/149‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                   
AGCCACTGGCGACTCGGCCTAGGATACGGCGCATGTGGATAGGGAGTGTTTTCATGCGTAAGTCTTTTGTCATCAGCCTGCTGAGTGCTGGCATCCTGCTGGCCGGCTGCCAGGCGGTGAATACCACCAGCGGCGGTGCTGTTGGCGTCGAACGCCAACAGTACATGTTCAGCATGCTCTCGACCGATGAGGTCAACCAGATGTACGCGCAGTCGTATCAGCAGACCCTGGGTGAAGCATCGAGCAAGGGTGTGCTCGACAAGTCCAGCGCCGATGCCAAGCGTGTGCAGGCCATT  >  NC_002947/874390‑874685

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: