Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 202,309 C→G 9.6% V2605V (GTC→GTG PP_0168 → putative surface adhesion protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947202,3090CG 9.6% 134.4 / 8.2 52V2605V (GTC→GTGPP_0168putative surface adhesion protein
Reads supporting (aligned to +/- strand):  ref base C (28/19);  new base G (3/2);  total (31/21)
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.63e-01

GCTGGCCAGATCAGCCTGGGCATCACGTCCGCTGTGGATGTAGACGGTCACACCTTCGAGAACCTGCAACTGGGCGGCAACGCTTCGGTTCAAGTGACCGACACCCTCGACGAGGTGGTGGCGAAGCTGACCGCGACCCCTTCGGTCACCGAAGGCGGTGAGATCACCTACACCATCACCCTGACCAACAAAGACGGTCTGCCGATCAACAACCACTCGGAGCTGTACTTCAAGCTGACCGACGGCACCACCGTCGTCGTGGCAGCCAACAGCACCACCGGTTCGGCCACTGCAA  >  NC_002947/202163‑202457
                                                                                                                                                  |                                                                                                                                                    
gCTGGCCAGATCAGCCTGGGCATCACGTCCGCTGTGGATGTAGACGGTCACACCTTCGAGAACCTGCAACTGGGCGGCAACGCTTCGGTTCAAGTGACCGACACCCTCGACGAGGTGGTGGCGAAGCTGACCGCGACCCCTTCGGTCAcc                                                                                                                                                   <  2:1084264/150‑1 (MQ=38)
                        aCGTCCGCTGTGGATGTAGACGGTCACACCTTCGAGAACCTGCAACTGGGCGGCAACGCTTCGGTTCAAGTGACCGACACCCTCGACGAGGTGGTGGCGAAGCTGACCGCGACCCCTTCGGTCAc                                                                                                                                                    >  1:1080792/1‑125 (MQ=33)
                        aCGTCCGCTGTGGATGTAGACGGTCACACCTTCGAGAACCTGCAACTGGGCGGCAACGCTTCGGTTCAAGTGACCGACACCCTCGACGAGGTGGTGGCGAAGCTGACCGCGACCCCTTCGGTCAc                                                                                                                                                    <  2:1080792/125‑1 (MQ=33)
                              gCTGTGGATGTAGACTGTCACACCTTCGAGAACCTGCAACTGGGCGGCAACGCTTCGGTTCAAGTGACCGACACCCTCGACGAGGTGGTGGCGAAGCTGACCGCGACCCCTTCGGTCACCGAAGGCGGTGAGATAACCTACACCATCAcc                                                                                                                     >  2:1257720/1‑150 (MQ=17)
                              gCTGTGGATGTAGACGGTCACACCTTCGACAACCTGCAACTGGGCGGCAACGCTTCGGTTCAAGTGACCGACACCCTCGACGAGGCGGTGGCGAAGCTGACCGCGACCCCTTCGGTCACCGAAGGCGGTGAGATCACCTACACCATCAcc                                                                                                                     >  1:1460203/1‑150 (MQ=17)
                                tgtgGATGTAGACGGTCACACCTTCGAGAACCTGCAACTGGGCGGCAACGCTTCGGTTCAAGTGACCGACACCCTCGACGAGGTGGTGGCGAAGCTGACCGCGACCCCTTCGGTCACCGAAGGCGGTGAGATCACCTACACCATCACCCt                                                                                                                   >  1:622845/1‑150 (MQ=32)
                                   ggATGTAGACGGTCACACCTTCGAGAACCTGCAACTGGGCGGCAACGCTTCGGTTCAAGTGACCGACACCCTCGACGAGGTGGTGGCGAAGCTGACCGCGACCCCTTCGGTCACCGAAGGCGGTGAGATCACCTACACCATCACCCTGAc                                                                                                                <  1:1625934/150‑1 (MQ=32)
                                    gATGTAGACGGTCACACCTTCGAGAACCCGCAACTGGGCGGCAACGCTTCGGTTCAAGTGACCGACACCCTCGACGAGGTGGTGGCGAAGCTGACCGCGACCCCTTCGGTCACCGAAGGCGGTGAGATCACCTACACCATCa                                                                                                                       >  2:1161396/1‑142 (MQ=17)
                                    gATGTAGACGGTCACACCTTCGAGAACCCGCAACTGGGCGGCAACGCTTCGGTTCAAGTGACCGACACCCTCGACGAGGTGGTGGCGAAGCTGACCGCGACCCCTTCGGTCACCGAAGGCGGTGAGATCACCTACACCATCa                                                                                                                       <  1:1161396/142‑1 (MQ=17)
                                               tCACACCTTCGAGAACCTGCAACTGGGCGGCAACGCTTCGGTTCAAGTGACCGACACCCTCGACGAGGTGGTGGCGAAGCTGACCGCGACCCCTTCGGTCACCGAAGGCGGTGAGATCACCTACACCATCACCCTGACCAACAAAGACgg                                                                                                    >  1:830378/1‑150 (MQ=35)
                                                  cacCTTCGAGAACCTGCAACTGGGCGGCAACGCTTCGGTTCAAGTGACCGACACCCTCGACGAGGTGGTGGCGAAGCTGACCGCGACCCCTTCGGTCACCGAAGGCGGTgaga                                                                                                                                      >  1:247586/1‑113 (MQ=33)
                                                  cacCTTCGAGAACCTGCAACTGGGCGGCAACGCTTCGGTTCAAGTGACCGACACCCTCGACGAGGTGGTGGCGAAGCTGACCGCGACCCCTTCGGTCACCGAAGGCGGTgaga                                                                                                                                      <  2:247586/113‑1 (MQ=33)
                                                          agaACCTGCAACTGGGCGGCAACGCTTCGGTTCAAGTGACCGACACCCTCGACGAGGTGGTGGCGAAGCTGACCGCGACCCCTTCGGTCACCGAAGGCGGTGAGATCACCTACACCATCACCCTGACCAACAAAGACGGTCTGCCGATca                                                                                         >  2:1477171/1‑150 (MQ=35)
                                                               cTGCAACTGGGCGGCAACGCTTCGGTTCAAGTGACCGACACCCTCGACGAGGTGGTGGCGAAGCTGACCGCGACCCCTTCGGTCACCGAAGGCGGTGAGATCACCTACACCATCACCCTGACCAACAAAGACGGTCTGCCGATCaacaac                                                                                    <  1:879512/150‑1 (MQ=35)
                                                               cTGCAACTGGGCGGCAACGCTTCGGTTCAAGTGACCGACACCCTCGACGAGGTGGTGGCGAAGCTGACCGCGACCCCTTCGGTCACCGAAGGCGGTGAGATCACCTACACCATCACCCTGACCAACAAAGACGGTCTGCCGATCaacaac                                                                                    <  2:1406371/150‑1 (MQ=35)
                                                                      tGGGCGGCAACGCTTCGGTTCAAGTGACCGACACCCTCGACGAGCTGGTGGCGAAGCTGACCGCGACCCCTTCGGTCACCGAAGGCGGTGAGATCACATACACCCTGACCCTGACGAACAAAGACGGTCTGCCGATCAACAGCCAGTCgg                                                                             >  1:385560/1‑150 (MQ=11)
                                                                                        ttCAAGTGACCGACACCCTCGACGAGGTGGTGGCGAAGCTGACCGCGACCCCTTCGGTCACCGAAGGCg                                                                                                                                            >  1:1623006/1‑69 (MQ=37)
                                                                                        ttCAAGTGACCGACACCCTCGACGAGGTGGTGGCGAAGCTGACCGCGACCCCTTCGGTCACCGAAGGCg                                                                                                                                            <  2:1623006/69‑1 (MQ=37)
                                                                                           aaGTGACCGACACCCCCGACGAGGTGGTGGCGAAGCTGACCGCGACCCCTTCGGTCACCGAAGGCGGTGAGATCACCTACACCATCACCCTGACCAACAAAGACGGTCTGCCGATCAACAACCACTCGGAGCTGTACTTCAAGCTGACCg                                                        >  1:265162/1‑150 (MQ=21)
                                                                                                  cGACACCCTCGACGAGGTGGTGGCGAAGCTGACCGCGACCCCTTCGGTCACCGAAGGCGGt                                                                                                                                          >  2:1322527/1‑61 (MQ=37)
                                                                                                  cGACACCCTCGACGAGGTGGTGGCGAAGCTGACCGCGACCCCTTCGGTCACCGAAGGCGGt                                                                                                                                          <  1:1322527/61‑1 (MQ=37)
                                                                                                  cGACACCCTCGACGAGGTGGTGGCGAAGCTGACCGCGACCCCTTCGGTCACCGAAGGCGGTGAGATCACCTa                                                                                                                               <  1:737821/72‑1 (MQ=35)
                                                                                                  cGACACCCTCGACGAGGTGGTGGCGAAGCTGACCGCGACCCCTTCGGTCACCGAAGGCGGTGAGATCACCTa                                                                                                                               >  2:737821/1‑72 (MQ=35)
                                                                                                     cacCCTCGACGAGGTGGTGGCGAAGCTGACCGCGACCCCTTCGGTCACCGAAGGCGGTGAGATCACCTACACCATCACCCTGACCAACAAAGACGGTCTGCCGATCAACAACCACTCGGAGCTGTACTTCAAGCTGACCGACggtggtgc                                              >  1:413568/1‑144 (MQ=25)
                                                                                                        ggtCGACGAGGTGGTGGCGAAGCTGACCGCGACCCCTTCGGTCACCGAAGGCGGTGAGATCACCTACACCATCACCCTGACCAACAAAGACGGTCTGCCGATCAACAACCACTCGGAGCTGTACTTCAAGCTGACCGACGGCACCACcgt                                           >  1:1056496/3‑150 (MQ=38)
                                                                                                        ccTCGACGAGGTGGTGGCGAAGCTGACCGCGACCCCTTCGGTCACCGAAGGCGGTGAGATCACCTACACCATCACCCTGACCAACAAAGACGGTCTGCCGATCaacaac                                                                                    <  1:647406/109‑1 (MQ=38)
                                                                                                        ccTCGACGAGGTGGTGGCGAAGCTGACCGCGACCCCTTCGGTCACCGAAGGCGGTGAGATCACCTACACCATCACCCTGACCAACAAAGACGGTCTGCCGATCaacaac                                                                                    >  2:647406/1‑109 (MQ=38)
                                                                                                        ccTCGACGAGGTGGTGGCGAAGCTGACCGCGACCCCTTCGGTCACCGAAGGCGGTGAGATCACCTACACCATCACCCTGACCAACAAAGACGGTCTGCCGATCAACAACCACTCGGAGCTGTACTTCAAGCTGACCGACGGCACCACcgt                                           >  2:885985/1‑150 (MQ=38)
                                                                                                          tCGACGAGGTGGTGGCGAAGCTGACCGCGACCCCTTCGGTCACCGAAGGCGGTGAGATCACCTACACCATCACCCTGACCAACAAAGACGGTCTGCCGATCAACACCCACCCGGAGCTGTACTTCAAGCTGACCGACGGCACCACcgtcg                                         >  1:438234/1‑150 (MQ=255)
                                                                                                             acgaGGTGGTGGCGAAGCTGACCGCGACCCCTTCGGTCACCGAAGGCGGTGAGATCACCTACACCATCACCCTGACCAACAAAGACGGTCTGCCGATCAACAACCACTTGGAGCTGTACTTCAAGCTGACCGACGGCACCACcgtcgtcg                                      >  1:249341/1‑150 (MQ=255)
                                                                                                             acgaGGTGGTGGCGAAGCTGACCGCGACCCCTTCGGTCACCGAAGGCGGTGAGATCACCTACACCATCACCCTGACCAACAAAGACGGTCTGCCGATCAACAACCACTCGGAGCTGTACTTCAAGCTGACCGACGGCACCACcgtcgtcg                                      <  2:1460203/150‑1 (MQ=39)
                                                                                                              cgaaggggtggcagtcagcttggccaccacttcgGTGACCGAAGGCGGTGAGATCACCTACACCATCACCCTGACCAACAAAGACGGTCTGCCGATCAACAACCACTCGGAGCTGTa                                                                      <  2:18095‑M2/83‑1 (MQ=255)
                                                                                                              cgaaggggtggcagtcagcttggccaccacttcgGTGACCGAAGGCGGTGAGATCACCTACACCATCACCCTGACCAACAAAGACGGTCTGCCGATCAACAACCACTCGGAGCTGTa                                                                      >  1:18095‑M2/35‑117 (MQ=255)
                                                                                                               gaaggggtggcagtcagcttggccaccacttcgGTGACCGAAGGCGGTGAGATCACCTACACCATCACCCTGACCAACAAAGACGGTCTGCCGATCAATAACCACTCGGa                                                                            >  2:90759‑M2/34‑110 (MQ=255)
                                                                                                               gaaggggtggcagtcagcttggccaccacttcgGTGACCGAAGGCGGTGAGATCACCTACACCATCACCCTGACCAACAAAGACGGTCTGCCGATCAATAACCACTCGGa                                                                            <  1:90759‑M2/77‑1 (MQ=255)
                                                                                                               gaaggggtggcagtcagcttggccaccacttcgGTGACCGAAGGCGGTGAGATCACCTACACCATCACCCTGACCAACAAAGACGGTCTGCCGATCAACAACCACTCGGAGCTGTACTTCAAGCTGACCGACGGCACCACCGTCGTCGTg                                    >  2:588097‑M2/34‑150 (MQ=255)
                                                                                                                         cGAAGCTGACCGCGACCCCTTCGGTCACCGAAGGCGGTGAGATCACCTACACCATCACCCTGACCAACAAAGACGGTCTGCCGATCAACAACCACTCGGAGCTGTACTTCAAGCTGACCGACGGCACCACCGTCTTGTTGGCAGCCAACa                          <  2:1056496/150‑1 (MQ=255)
                                                                                                                          gAAGCTGACCGCGACCCCTTCGGTCACCGAAGGCGGTGAGATCACCTACACCATCACCCTGACCAACAAAGACGGTCTGCCGATCAACAACCACTCGGAGCTGTACTTCAAGCTGACCGACGGCACCACCGTCGTCGTGGCAGCCAACAg                         >  2:1242213/1‑150 (MQ=37)
                                                                                                                            aGCTGACCGCGACCCCTTCGGTCACCGAAGGCGGTGAGATCACCTACACCATCACCCTGACCAACAAAGACGGTCTGCCGATCAACAACCACTCGGAGCTGTACTTCAAGCTGACCGACGGCACCACCGTCGTCGTGGCAGCCAACAGca                       <  1:1477171/150‑1 (MQ=39)
                                                                                                                               tGACCGCGACCCCTTCGGTCACCGAAGGCGGTGAGATCACCTACACCATCACCCTGACCAACAAAGACGGTCTGCCGATCAACAACCACTCGGAGCTGTACTTCAAGCTGACCGACGGCACCACCGTCGTCGTGGCAGCCAACAGcacca                    <  2:265162/150‑1 (MQ=37)
                                                                                                                                gACCGCGACCCCTTCGGTCACCGAAGGCGGTGAGATCACCTACACCATCACCCTGACCAACAAAGACGGTCTGCCGATCAACAACCACTCGGAGCTGTACTTCAAGCTGACCGACGGCACCACCGTCGTCGTGGCAGCCAACAGcaccac                   >  1:236196/1‑150 (MQ=38)
                                                                                                                                    tcaccgccTTCGGTCACCGAAGGGGGTGAGATCACCTACACCATCACCCTGACCAACAAAGACGGTCTGCCGATCAACAACCACTCGGAGCTGTACTTCAAGCTGACCGACGGCACCACCGTCGTCGTGGCAGCCAACAGCACCACCGGt               >  1:607606/7‑150 (MQ=12)
                                                                                                                                    tcaccgccTTCGGTCACCGAAGGGGTGAGATCACCTAAAGCAAAACAATCAACAGTAAAGACAGTCTGCCGAGCAACAACCACTCGGAGCTGTACTTCAAGCTGACCGACGGCACCACCGTCGTCGTGGCAGCCAATAGCACCACCGGtt              >  1:30506/7‑150 (MQ=11)
                                                                                                                                        ccccTTCGGTCACCGAAGGCGGTGAGATCACCTACACCATCACCCTGACCAACAAAGACGGTCTGCCGATCAACAACCCCGCGGAGCTGTACTTCAAGCTGACCGACGGCACCACCGTCGTCGTGGCAGCCAACAGCACCACCGGTTCgg           <  2:310077/150‑1 (MQ=17)
                                                                                                                                         cccTTGGGTCACCCAAGGCGGTGAGATCTCCTACACCCTCACCCTGACCAACAAAGACGGTCTGCCGTTCACCAACCACTCGCAGCTGTACTACCAGCTTACCGACGGCACAACCGTAGTCGTGGCTACCAACAGCCCCAGCgccacggc          >  2:300385/1‑143 (MQ=11)
                                                                                                                                         cccTTCGTTCACCGAAGGCGGTGCGATCACCTACACCACCACCCTGACCAACAAAGACGGTCTGCCGATCAAGATCCACTCGGAGCTGTACTTCAAGATGACCGACTCCACCCCCGGCTTCGGGGCAAGCAACAGCACCACCGGTTCGGc          <  2:438234/150‑1 (MQ=11)
                                                                                                                                         cccTTCGGTCACCGAAGGCGGTGAGATCACCTACACCATCACCCTGACCAACAAAGACGGTCTGCCGATCAACAACCACTCGGAGCTGTACTTCAAGCTGACCGACGGCACCACCGTCGTCGTGGCAGCCAACAGCACCACCGGTTCGGc          >  2:469104/1‑150 (MQ=32)
                                                                                                                                         cccTTCGGTCACCGAAGGCGGTGAGATCACCTACACCATCACCCTGACCAACAAAGACGGTCTGCCGATCAACAACCACTCGGAGCTGTACTTCAAGCTGACCGACGGCACCACCGTCGTCGTGGCAGCCAACAGCACCACCGGTTCGGc          <  2:1037125/150‑1 (MQ=32)
                                                                                                                                            ttcgGTCACCGAAGGCGGTGAGATCACCTACACCATCACCCTGACCAACAAAGACGGTCTGc                                                                                               <  2:227773‑M2/58‑1 (MQ=255)
                                                                                                                                            ttcgGTCACCGAAGGCGGTGAGATCACCTACACCATCACCCTGACCAACAAAGACGGTCTGc                                                                                               >  1:227773‑M2/5‑62 (MQ=255)
                                                                                                                                                gTCACCGAAGGCGGTGAGATCACCTACACCATCACCCTGACCAACAAAGACGGTCTGCCGATCAACAACCACTCGGAGCTGTACTTCAAGCTGACCGACGGCACCACCGTCGTCGTGGCAGCCAACAGCACCACCGGTTCGGCCACTGCa   >  1:678864/1‑150 (MQ=32)
                                                                                                                                                 tCACCGAAGGCGGTGAGATCACCTACACCATCACCCTGACCAACAAAGACGGTCTGCCGATCAACAACCACTCGGAGCTGTACTTCAAGCTGACCGACGGCACCACCGTCGTCGTGGCAGCCAACAGCACCACCGGTTCGGCCACTGCaa  >  2:1472019/1‑150 (MQ=31)
                                                                                                                                                  |                                                                                                                                                    
GCTGGCCAGATCAGCCTGGGCATCACGTCCGCTGTGGATGTAGACGGTCACACCTTCGAGAACCTGCAACTGGGCGGCAACGCTTCGGTTCAAGTGACCGACACCCTCGACGAGGTGGTGGCGAAGCTGACCGCGACCCCTTCGGTCACCGAAGGCGGTGAGATCACCTACACCATCACCCTGACCAACAAAGACGGTCTGCCGATCAACAACCACTCGGAGCTGTACTTCAAGCTGACCGACGGCACCACCGTCGTCGTGGCAGCCAACAGCACCACCGGTTCGGCCACTGCAA  >  NC_002947/202163‑202457

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: