Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 197,089 T→C 11.0% N865N (AAT→AAC PP_0168 → putative surface adhesion protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947197,0890TC11.0% 94.6 / 3.3 49N865N (AAT→AACPP_0168putative surface adhesion protein
Reads supporting (aligned to +/- strand):  ref base T (22/18);  new base C (2/3);  total (25/24)
Fisher's exact test for biased strand distribution p-value = 6.52e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CGTTGAAGCCGGCAAGACCCAGGGCAGCGTCGATTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACCGTCAGCGTCACCATCGAGAACGCCACAGGCGGTAATTTCGAGCAACTGACCCCGAATCCGACCCCGGCTCAGACCACGATCAACGACTCGGTCGATGCCACCACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTAATCACCTACACCGTGACCCTGAGCAATCCTGCCCAGACGCCGGTA  >  NC_002947/196957‑197224
                                                                                                                                    |                                                                                                                                       
cGTTGAAGCCGGCAAGACCCAGGGCAGCGTCGATTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACCGTCAGCGTCACCATCGAGAACGCCACAGGCGGTAATTTCGAGCAACTGACCCCGAATCCg                                                                                                                                      >  1:116801/1‑136 (MQ=34)
cGTTGAAGCCGGCAAGACCCAGGGCAGCGTCGATTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACCGTCAGCGTCACCATCGAGAACGCCACAGGCGGTAATTTCGAGCAACTGACCCCGAATCCg                                                                                                                                      >  1:138796/1‑136 (MQ=34)
cGTTGAAGCCGGCAAGACCCAGGGCAGCGTCGATTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACCGTCAGCGTCACCATCGAGAACGCCACAGGCGGTAATTTCGAGCAACTGACCCCGAATCCg                                                                                                                                      <  1:163405/136‑1 (MQ=34)
cGTTGAAGCCGGCAAGACCCAGGGCAGCGTCGATTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACCGTCAGCGTCACCATCGAGAACGCCACAGGCGGTAATTTCGAGCAACTGACCCCGAATCCg                                                                                                                                      <  1:205274/136‑1 (MQ=34)
     aaGCCGGCAAGACCCAGGGCAGCGTCGATTTCCAGACCCCGGCGAATGACGTCTACAACAACCGTTCGACCGTCAGCGTCACCATCGAGAACGCCACACGAGGTAATTTCGAGCAACTGACCCCGAAACCGAcacc                                                                                                                                 >  1:337262/1‑136 (MQ=11)
     aaGCCGGCAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACCGTCAGCGTCACCATCGAGAACGCCACAGGCGGTAATTTCGAGCAACTGACCCCGAATCCGAcccc                                                                                                                                 >  2:323440/1‑136 (MQ=17)
      aGCCGGCAAGACCCAGGGCAGCGTCGATTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACCGTCAGCGTCACCATCGAGAAGGCCACAGGCGGTAATTTCGAGCAACTGACCCCGAATCCGACCCCg                                                                                                                                <  1:357811/136‑1 (MQ=37)
      aGCCGGCAAGACCCAGGGCAGCGTCGATTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACCGTCAGCGTCACCATCGAGAACGCCACAGGCGGTAATTTCGAGCAACTGACCCCGAATCCGACCCCg                                                                                                                                <  1:65955/136‑1 (MQ=37)
      aGCCGGCAAGACCCAGGGCAGCGTCGATTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACCGTCAGCGTCACCATCGAGAACGCCACAGGCGGTAATTTCGAGCAACTGACCCCGAATCCGACCCCg                                                                                                                                <  1:7538/136‑1 (MQ=37)
       ccccccAAGACCCAGGGCAGCCACGAATTCCACACCCCGGCGAATCACGTCAACAACAACGGAACGACCCTCACCGTCACCATCGACAACCCCACAGGCGCTAAATACGAGCAACAGACCCCGAATCCGACCCCgc                                                                                                                               >  2:70612/6‑135 (MQ=2)
          ggCAAGACCCAGGGCAGCGTCGATTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACCGTCAGCGTCACCATCGAGAACGCCACAGGCGGTAATTTCGAGCAACTGACCCCGAATCCGACCCCGGCTc                                                                                                                            <  2:406682/136‑1 (MQ=37)
             aaGACCCAGGGCAGCGTCGATTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACCGTCAGCGTCACCATCGAGAACGCCACAGGCGGTAATTTCGAGCAACTGACCCCGAATCCGACCCCGGCTCAgc                                                                                                                         >  1:33754/1‑135 (MQ=37)
             aaGACCCAGGGCAGCGTCGATTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACCGTCAGCGTCACCATCGAGAACGCCACAGGCGGTAATTTCGAGCAACTGACCCCGAATCCGACCCCGGCTCAGa                                                                                                                         >  2:101462/1‑136 (MQ=37)
              aGACCCAGGGCAGCGTCGATTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACCGTCAGCGTCACCATCGATAACGCCACAGGCGGTAATTTCGAGCAACTGACCCCTAATCCGACCCCGGCTCAGAc                                                                                                                        <  1:364088/136‑1 (MQ=255)
                    aGGGCTGCGTCGATTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACCGTCAGCGTCACCATCGAGAACGCCACAGGCGGTAATTTCGAGCAACTGACCCCGAATCCGACCCCGGCTCTGACCACGAt                                                                                                                  <  1:132608/136‑1 (MQ=21)
                         aGCTTCGATTTCCTGACCCCGGCTAATTACGTCTACAACAACGTTTCGACCGTCAGCGTCACAATCGATAACGCCACAGTCGGTAATTTCTAGAAACTGACCCCTAATCGGACCCCGTCTCAGACCACGATCAACg                                                                                                             <  2:43931/136‑1 (MQ=11)
                         aGCGTCGATTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACCGTCAGCGTCACCATCGAGAACGCCACAGGCGGTAATTTCGAGCAACTGACCCCGAATCCGACCCCGGCTCAGACCACGATCAACg                                                                                                             <  1:323440/136‑1 (MQ=34)
                            gTCGATTTCCAGCCCCCGGCTAATGACGTCTACAACAACGGTTCGACCGTCAGCGTCACCATCGAGAACGCCACAGGCGGTAATTTCGAGCAACTGACCCCGAATCCGACCCCGGCTCAGACCACGATCAACGACt                                                                                                          <  1:67791/136‑1 (MQ=17)
                              ctaTTTCCAGACCCCGGCTAATGACGTCTACAACAACGGTTCGACCGTCAGCGTCACCATCGAGAACGCCACAGGCGGTAATTTCGAGCAACTCACCCCGAATCCGACCCCGGCTCAGACCACGATCAACGACTCg                                                                                                        <  1:200781/134‑1 (MQ=12)
                              cGATTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACCGTCAGCGTCACCATCGAGAACGCCACAGGCGGTAATTTCGAGCAACTGACCCCGAATCCGACCCCGGCTCAGACCACGATCAACGCCTCg                                                                                                        >  2:23959/1‑136 (MQ=21)
                              cGATTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACCGTCAGCGTCACCATCGAGAACGCCACAGGCGGTAATTTCGAGCAACTGACCCCGAATCCGACCCCGGCTCAGACCACGATCAACGACTcc                                                                                                        >  1:194000/1‑135 (MQ=21)
                              cGATTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACCGTCAGCGTCACCATCGAGAACGCCACAGGCGGTAATTTCGAGCAACTGACCCCGAATCCGACCCCGGCTCAGACCACGATCAACGACTCg                                                                                                        >  1:88586/1‑136 (MQ=34)
                              cGATTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACCGTCAGCGTCACCATCGAGAACGCCACAGGCGGTAATTTCGAGCAACTGACCCCGAATCCGACCCCGGCTCAGACCACGATCAACGACTCg                                                                                                        >  2:294839/1‑136 (MQ=34)
                              cGATTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACCGTCAGCGTCACCATCGAGAACGCCACAGGCGGTAATTTCGAGCAACTGACCCCGAATCCGACCCCGGCTCAGACCACGATCAACGACTCg                                                                                                        <  1:176921/136‑1 (MQ=34)
                                   tCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACCGTCAGCGTCACCATCGAGAACGCCACAGGCGGTAATTTCGAGCAACTGACCCCGAATCCGACCCCg                                                                                                                                <  1:50180/107‑1 (MQ=35)
                                   tCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACCGTCAGCGTCACCATCGAGAACGCCACAGGCGGTAATTTCGAGCAACTGACCCCGAATCCGACCCCg                                                                                                                                >  2:50180/1‑107 (MQ=35)
                                   tCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACCGTCAGCGTCACCATCGAGAACGCCACAGGCGGTAATTTCGAGCAACTGACCCCGAATCCGACCCCGGCTCAGACCACGATCAACGACTCGGTCGa                                                                                                   >  2:227715/1‑136 (MQ=34)
                                   tCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACCGTCAGCGTCACCATCGAGAACGCCACAGGCGGTAATTTCGAGCAACTGACCCCGAATCCGACCCCGGCTCAGACCACGATCAACGACTCGGTCGa                                                                                                   >  2:88882/1‑136 (MQ=34)
                                        aCCCCGGCCAATGACGTCTACAACAACGGTTCGACCGTCAGCGTCACCATCGAGAACGCCACAGGCGGTAATTTCGAGCAACTGACCCCGAACCCGACCCCGGCTCAGACCACGATCAACGACTCCGTCGATGcca                                                                                              >  2:257130/1‑136 (MQ=255)
                                          cccGGCGAATGACGTCTACAACAACGGTTCGACCGTCAGCGTCACCATCGAGAACGCCACAGGCGGTAATTTCGAGCAACTGACCCCGAATCCGACCCCGGCTCAGACCACGATCAACGACTCGGTCGATGccacc                                                                                            >  2:221442/1‑136 (MQ=40)
                                           ccGGCGAATGACGTCAACAACAACCCTTCGACCGTCAGCGTCACCATCGAGAACGACACAGGCGGTAAATTCGAGCAACTGACCCCGAATCCGACCCCGGCTCAGACAACGATCAACGACTCGGTCGATGccacca                                                                                           >  1:403144/1‑136 (MQ=14)
                                                       gTCTACAACAACGGTTCGACGGTCAGCGTCACCATCGAGAACGCCACAGGCGGTAATTTCGAGCAACTGACCCCGAACCCGACCCCGGCTCAGACCACGATCAACGACTCGGTCGAt                                                                                                  <  1:84812/117‑1 (MQ=2)
                                                       gTCTACAACAACGGTTCGACCGTCAGCGTCACCATCGAGAACGCCACAGGCGGTAATTTCGAGCAACTGACCCCGAACCCGACCCCGGCTCAGACCACGATCAACGACTCGGTCGAt                                                                                                  >  2:84812/1‑117 (MQ=18)
                                                        tCTACAACAAGGTTGCGACCGTCTGCGTCACCATCGTCAAGCCCTCTGGGGTTATTTTCCAGCAACTGACCCGCTATCCGCCCCCGGCTCACACCACTATCAACGACTCGGTCGATGCCACCACCGCGACCCTGAc                                                                              <  2:21275/136‑1 (MQ=2)
                                                                     ttCGACCGTCAGCGTCACCATCGAGAACGCCACAGGCGGTAATTTCGAGCAACTGACCCCGAATCCGACCCCGGCTCAGACCACGATCAACGACTCGGTCGATGCCACCACCGCGACCCTGACGGCGAGCCCGTcc                                                                 >  1:239117/1‑135 (MQ=37)
                                                                                gCGTCACCATCGAGAACGCCACAGGCGGTAATTTCGAGCAACTGACCCCGAATCCGACCCCGGCTCAGACCACGATCAACGACTCGGTCGATGCCACCACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAgg                                                      >  1:116169/1‑136 (MQ=37)
                                                                                    cACCATCGACACCGCCACAGGCGTTAATTTCGACCAACTCACCCCGAACCCCCCCCCGGCCCACCCCACCATCAACGACTCGGTCGATGCCACCACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAggcggc                                                  <  2:252019/136‑1 (MQ=11)
                                                                                    cACCATCGACAACGCCACAGGCGGTAATTTCGAGCAACTGACCCCGAATCCGACCCCGGCTCAGACCACGATCAACGACTCGGTCGATGCCACCACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAggcggc                                                  <  1:117373/136‑1 (MQ=21)
                                                                                            agaACGCCACAGGCGGTAATTTCGAGCAACTGACCCCGAATCCGACCCCGCCTCAGACCACGATCAACGACTCGGTCGATGCCACCACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTAATCAc                                          <  1:221442/136‑1 (MQ=37)
                                                                                            acaccGCCACCGGCGGTAATTTCGAGCAACTCCCCCCGAAGCCGACCCCGGCTCAGACCACCATCAACGACTCGGTCGATGCCACCACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTAATCAc                                          <  2:194000/132‑1 (MQ=11)
                                                                                            acaaCGCCACAGGCGGTAATTTCGAGCAACTGACCCCGAATCCGACCCCGGCTCAGACCACGATCAACGACTCGGTCGATGCCACCACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCAc                                          <  1:138285/134‑1 (MQ=17)
                                                                                              aaCGCCACAGGCGGTAATTTCGAGCAACTGACCCCGAATCCGACCCCGGCTCAGACCACGATCAACGACTCGGTCGATGCCACCACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCActt                                        >  2:365816/1‑134 (MQ=17)
                                                                                                    acaGGCGGTAATTTCGAGCAACTGACCCCGAATCCGACCCCGGCTCAGACCACGATCAACGACTCGGTCGATGCCACCACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTAATCACCTACACCg                                  >  1:93983/1‑136 (MQ=37)
                                                                                                    acaGGCGGTAATTTCGAGCAACTGACCCCGAATCCGACCCCGGCTCAGACCACGATCAACGACTCGGTCGATGCCACCACCGCGACCCTGACGGCGAGCCAGTCGGTCACCGAAGGCGGCGTAAACACCTACACCg                                  >  2:200185/1‑136 (MQ=255)
                                                                                                      tggctttAATTGCTATCAACTCCCCCCTAAGCCCACCCGGGCTCAGACAACCATCAACGACTCGGTCGATGCCACCACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTAATCACCTACACCGTg                                <  1:355404/130‑1 (MQ=11)
                                                                                                      aGGCGGTAATTTCGAGCAACTGACCCCGAATCCGACCCCGGCTCAGACCACGATCAACGACTCGGTCGATGCCACCACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTAATCACCTACACCGTg                                >  2:73037/1‑136 (MQ=35)
                                                                                                              aTTTCGAGCAACTGACCCCGAATCCGACCCCGGCTCAGACCACGATCAACGACTCGGTCGATGCCACCACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTAATCACCTACACCGTGACCCTGAg                        >  2:106365/1‑136 (MQ=34)
                                                                                                                        attGACCCCCATCCCGACCCGGGTTCACACCAGTATCTACTACTCGGGCGATGCCACCAGCGCGACCGTGACGGCGAGGCCGTGGTTCACCTAAGGCGGCGTAATCACCTACACCGTTACCATGAGCAATCCTGcc              <  2:239117/134‑1 (MQ=255)
                                                                                                                          tcaCCCCGAAGCCGACCCCGGCTCAGCCCAGTATCAACGACTCGGGCGATGCCACCACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTAATCACCTACACCGTGACCCTGAGCAATCCTGCCCa            <  1:101462/134‑1 (MQ=11)
                                                                                                                                    tCCGACCCCGGCTCAGACCACGATCAACGACTCGGTCGATGCCACCACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTAATCACTTACACCGTGACCCTGAGCAACCCTGCCCAGACGCCGGtg  >  1:163809/1‑135 (MQ=12)
                                                                                                                                    |                                                                                                                                       
CGTTGAAGCCGGCAAGACCCAGGGCAGCGTCGATTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACCGTCAGCGTCACCATCGAGAACGCCACAGGCGGTAATTTCGAGCAACTGACCCCGAATCCGACCCCGGCTCAGACCACGATCAACGACTCGGTCGATGCCACCACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTAATCACCTACACCGTGACCCTGAGCAATCCTGCCCAGACGCCGGTA  >  NC_002947/196957‑197224

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