Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 197,089 T→C 27.4% 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,0890TC27.4% 76.7 / 32.6 59N865N (AAT→AACPP_0168putative surface adhesion protein
Reads supporting (aligned to +/- strand):  ref base T (24/15);  new base C (3/12);  total (31/28)
Fisher's exact test for biased strand distribution p-value = 1.35e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.40e-01

CGTTGAAGCCGGCAAGACCCAGGGCAGCGTCGATTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACCGTCAGCGTCACCATCGAGAACGCCACAGGCGGTAATTTCGAGCAACTGACCCCGAATCCGACCCCGGCTCAGACCACGATCAACGACTCGGTCGATGCCACCACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTAATCACCTACACCGTGACCCTGAGCA  >  NC_002947/196957‑197204
                                                                                                                                    |                                                                                                                   
cGTTGAAGCCGGCAAGACCCAGGGCAGCGTCGATTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACCGTCAGCGTCACCATCGAGAACGCCACAGGCGGTAATTTCGAGCAACTGACCCCGAATCCg                                                                                                                  <  1:253820/136‑1 (MQ=34)
cGTTGAAGCCGGCAAGACCCAGGGCAGCGTCGATTTCCAGACCCCGGCGAATGACGTCTACAACAACGGATCGACCGTCAGCGTCACCATCGAGAACGCCACAGGCGGTAAGTTCGAGCAACTGCCCCCGAAACCg                                                                                                                  >  1:248386/1‑136 (MQ=255)
     aaGCCGGCAAGACCCAGGGCAGCGTCGATTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACCGTCAGCGTCACCATCGAGAACGCCACAGGCGGTAATTTCGAGCAACTGACCCCGAATCCGAcccc                                                                                                             >  1:241954/1‑136 (MQ=37)
          ggCAAGACCCAGGGCAGCGTCGATTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACCGTCAGCGTCACCATCGAGAACGCCACAGGCGGAAAGTTCGAGCAACTGACCCCGAAACCGACCCCGGCTc                                                                                                        >  2:278774/1‑136 (MQ=255)
              aGACCCTGGGCAGCGTGGATTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACCGTCTGCGTCACCATCGATAACGCCACAGGCGGTATTTTCGAGCAACTGACCCCTAATCCGACCCCGGCTCAGAc                                                                                                    <  2:105864/136‑1 (MQ=11)
                       gCAGCGTCGATTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACCGTCAGCGTCACCATCGAGAACGCCACAGGCGGTAATTTCGAGCAACTGACCCCGAATCCGACCCCGGCACAGACCACGATCaa                                                                                           >  2:342488/1‑136 (MQ=21)
                         aGCGTCGATTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACCGTCAGCGTCACCATCGAGAACGCCACAGGCGGTAATTTCGAGCAACTGACCCCGAATCCGACCCCGGCTCAGACCACGATCAACg                                                                                         >  2:399055/1‑136 (MQ=34)
                             tCGATTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACCGTCAGCGTCACCATCGAGAACGCCACAGGCGGTAATTTCGAGCAACTGACCCCGAATCCGACCCCGGCTCAGACCACGATCAACg                                                                                         >  1:444891/1‑132 (MQ=34)
                             tCGATTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACCGTCAGCGTCACCATCGAGAACGCCACAGGCGGTAATTTCGAGCAACTGACCCCGAATCCGACCCCGGCTCAGACCACGATCAACg                                                                                         <  2:444891/132‑1 (MQ=34)
                              cGATTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACCGTCAGCGTCACCATCGAGAACGCCACAGGCGGTAATTTCGAGCAACTGACCCCGAATCCGACCCCGGCTCAGACCACGATCAACGACTCg                                                                                    >  1:60225/1‑136 (MQ=34)
                              cGATTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACCGTCAGCGTCACCATCGAGAACGCCACAGGCGGTAATTTCGAGCAACTGACCCCGAATCCGACCCCGGCTCAGACCACGATCAACGACTCg                                                                                    >  1:139615/1‑136 (MQ=34)
                              cGATTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACCGTCAGCGTCACCATCGAGAACGCCACAGGCGGTAATTTCGAGCAACTGACCCCGAATCCGACCCCGGCTCAGACCACGATCAACGACTCg                                                                                    >  1:436581/1‑136 (MQ=34)
                              cGATTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACCGTCAGCGTCACCATCGAGAACGCCACAGGCGGTAATTTCGAGCAACTGACCCCGAATCCGACCCCGGCTCAGACCACGATCAACGACTCg                                                                                    >  1:231112/1‑136 (MQ=34)
                              cGATTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACCGTCAGCGTCACCATCGAGAACGCCACAGGCGGTAATTTCGAGCAACTGACCCCGAATCCGACCCCGGCTCAGACCACGATCAACGACTCg                                                                                    >  1:66413/1‑136 (MQ=34)
                              cGATTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACAGTCAGCGTCACCATCGAGAACGCCACAGGCGGTAATTTCGAGCAACTGACCCCGAAACCGACCCCGGCTCAGACCACGATCAACGACTCg                                                                                    >  1:432350/1‑136 (MQ=21)
                                aTTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACCGTCAGCGTCACCATCGAGAACGCCACAGGCGGTAATTTCGAGCAACTGACCCCGAATCCGACCCCGGCTCAGACCACGATCAACGACTCGGt                                                                                  <  1:342656/136‑1 (MQ=34)
                                      aGACCCCGGCGAATGACGTCTACAACAACGGTTCGACTGTCAGCGCCACCATCGAGAACGCCACAGGCGGTAATTTCGAGCAACTGACCCCGAACCCGCCCCCGGCTCAGACCACGATCAACGACTCGGTCGATGc                                                                            <  2:177278/136‑1 (MQ=255)
                                      aGACCCCGGCGAATGACGTCTACAACAACGGTTCGACCGTCAGCGTCACCATCGAGAACGCCACAGGCGGTAATTTCGAGCAACTGACCCCGAATCCGACCCCGGCTCAGACCACGATCAACGACTCGGTCGATGc                                                                            <  1:399055/136‑1 (MQ=37)
                                        aCCCCGGCGAATGACGTCTACAACAACGGTTCGACCGTCAGCGTCACCATCGAGAACGCCACAGGCGGTAATTTCGAGCAACTGACCCCGAATCCGACCCCGGCTCAGACCACGATCAACGACTCGGTCGATGcca                                                                          >  1:59425/1‑136 (MQ=39)
                                              gCGAATGACGTCTACAACAACGGTTCGACCGTCAGCGTCACCATCGAGAACGCCACAGGCGGTAATTTCGAGCAACTGACCCCGAATCCGACCCCGGCTCAGACCACGATCAACGACTCGGTCGATGCCaccaccc                                                                    >  2:3627/1‑135 (MQ=40)
                                               cGAATGACGTCTACAACAACGGTTCGACCGTCAGCGTCACCATCGAGAACGCCACAGGCGGTAATTTCGAGCAACTGACCCCGAATCCGACCCCGGCTCAGACCACGATCAACGACTCGGTCGATGCCACCACcgc                                                                   >  1:25701/1‑136 (MQ=39)
                                                gAATGACGTCTACAACAACGGTTCGACCGTCAGCGTCACCATCGAGAACGCCACAGGCGGTAATTTCGAGCAACTGACCCCGAATCCGACCCCGGCTCAGCCCACGATCAACGACTCGGTCGATGCCACCACcgcg                                                                  >  1:187296/1‑136 (MQ=38)
                                                                aacGGTTCGACCGTCAGCGTCACCATCGAGAACGCCACAGGCGGTAATTTCGAGCAACTGACCCCGAATCCGACCCCGGCTCAGACCACGATCAACGACTCGGTCGATGCCACCACCGCGACCCTGACGGCGAGcc                                                  >  2:462618/1‑136 (MQ=37)
                                                                  cGGTTCGACTGTCAGCGTCACCATCGAGAACGCCACAGGCGGTAATTTCGAGCAACTGACCCCGAACCCGACCCCGGCTCAGACCACGATCAACGACTCGGTCGATGCCACCACCGCGACCCTGACGGCGAGCCCg                                                >  2:422463/1‑136 (MQ=21)
                                                                  cGGTTCGACTGTCAGCGTCACCATCGAGAACGCCACAGGCGGTAATTTCGAGCAACTGACCCCGAACACGACCCCGGCTCAGCCCACGATCAACGCCTCGGTCGATGCCACCACCGCGACCCTGACGGCGAGcccc                                                >  2:427166/1‑135 (MQ=11)
                                                                   ggTTCGACGGTCAGCGTCACCATCGACACCCCCACCGGCGTTAATTTCGACCAACTCACCCCGAACCCCACCCCGGCGCAGACCACGATCAACGACTCGGTCGATGCCACCACCGCGACCCTGACGGCGAGCCCGt                                               <  1:3627/136‑1 (MQ=11)
                                                                   ggTTCGACCGTCAGCGTCACCATCGAGAACGCCACAGGCGGTAATTTCGAGCAACTGACCCCGAACCCGACCCCGGCTCAGACCACGATCAACGACTCGGTCGATGCCACCACCGCGACCCTGACGGCGAGCCCGt                                               <  2:128207/136‑1 (MQ=37)
                                                                    gTTCGACCGTCAGCGTCACCATCGAGAACGCCACAGGCGGTAATTTCGAGCAACTGACCCTGAATCCg                                                                                                                  >  1:24693/1‑68 (MQ=25)
                                                                    gTTCGACCGTCAGCGTCACCATCGAGAACGCCACAGGCGGTAATTTCGAGCAACTGACCCTGAATCCg                                                                                                                  <  2:24693/68‑1 (MQ=25)
                                                                     ttCGACTGTCAGCGTCACCATCGAGAACGCCACAGGCGGTAATTTCGAGCAACTGACCCCGAACCCGACCCCGGCTCAGACCACGATCAACGACTCGGTCGATGCCACCACCGCGACCCTGACGGCGAGCCCGTCg                                             >  1:83461/1‑136 (MQ=21)
                                                                       cGACCGTCAGCGTCACCATCGAGAACGCCACAGGCGGTAATTTCGAGCAACTGACCCCGAATCCGACCCCGGCTCAGACCACGATCAACGACTCGGTCGATGCCACCACCGCGACCCTGACGGCGAGCCCGTCGGt                                           >  1:414069/1‑136 (MQ=37)
                                                                       cGACCGTCAGCGTCACCATCGAGAACGCCACAGGCGGGAAGTTCGAGCAACTGACCCCGAAACCGACCCCGGCTCAGACCACGATCAACGACTCGGTCGATGCCACCACCGCGACCCTGACGGCGAGCCCGTCGGt                                           >  2:4178/1‑136 (MQ=255)
                                                                         cccGTCTCCGTCACCACCGACAACGCCACCGGCTGTAATTTCCACCATCTCACCCCTAACCCCACCCCGGCTCACACCACGACCAACGACTCGGTCGATGCCACCACCGCGACCCTGACGGCGAGCCCGTCGGTCa                                         <  2:414069/135‑1 (MQ=11)
                                                                          ccGTCAGCGTCACCATCGAGAACGCCACAGGCGGTAATTTCGAGCAACTGACCCCGAATCCGACCCCGGCTCAGACCACGATCAACGACTCGGTCGATGCCACCACCGCGACCCTGACGGCGAGCCCGTCGGTCAc                                        >  1:233116/1‑136 (MQ=37)
                                                                                gCGTCACCATCGAGAACGCCACAGGCGGTAATTTCGAGCAACTGACCCCGAATCCGACCCCGGCTCAGACCACGATCAACGACTCGGTCGATGCCACCACCGCGAc                                                                >  1:447483/1‑106 (MQ=38)
                                                                                gCGTCACCATCGAGAACGCCACAGGCGGTAATTTCGAGCAACTGACCCCGAACCCGACCCCGGCTCAGACCACGATCAACGACTCGGTCGATGCCACCACCGCGAc                                                                <  2:447483/106‑1 (MQ=255)
                                                                                  gTCACCATCGATAACGCCACAGGCGTAAATTTCGAGCACCTGCCCCCGAACCCGACCCCGGGTCTGACCACTATCAACGACTCGGTCGATGCCACCACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAggcg                                <  1:465716/136‑1 (MQ=11)
                                                                                  cccACCATCCACACCCCCACTGGCGGTAATTTCGAGCACCTCACCCCGAACCCCACCCCGGCCCAGACCACCAGCACCGACTCGGTCGATGCCACCACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAggcg                                <  1:265507/134‑1 (MQ=11)
                                                                                    cACCATCGATAACGCCACAGGCGGAAATTTCGAGCAACTGACCCCGAACCCGACCCCGGCTCAGACCACGATCAACGACTCGGTCGATGCCACCACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAggcggc                              <  2:432350/136‑1 (MQ=255)
                                                                                    aaCCATCGTCAACGCCACAGGCGTTAATTTCTAGCAACTGACCCCCAAACCGACCCCGGCACAGACCACTATCAACGACTCGGTCGATGCCACCACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAggcggc                              <  1:342488/135‑1 (MQ=11)
                                                                                      ccATCGTCAACGCCACAGGCGGTAATTTCGACCAACTCACCCCGAATCCCACCCCGGCTCAGACCACGATCAACGACTCGGTCGATGCCACCACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGt                            <  1:723/136‑1 (MQ=11)
                                                                                      ccATCGAGAACGCCACAGGCGGTAATTTCGAGCAACTGACCCCGAATCCGACCCCGGCTCAGACCACGATCAACGACTCGGTCGATGCCACCACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGt                            <  2:294512/136‑1 (MQ=34)
                                                                                      ccATCGAGAACCCCACAGGCGGTAATTTCGCGCAACTCCCCCCGAACCCGACCCCGGCTCACACCACGATCAACGACTCGGTCGATGCCACCACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGt                            <  2:169676/136‑1 (MQ=11)
                                                                                      ccATCGACAACGCCACAGGCGGAAATTTCGAGCAACTGACCCCCAACCCGACCCCGGCTCACACCACGATCAACGACTCGGTCGATGCCACCACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGt                            <  2:167417/136‑1 (MQ=11)
                                                                                          cGAGAACGCCACAGGCGGTAATTTCGAGCAACTGACCCCGAATCCGACCCCGGCTCAGACCACGATCAACGACTCGGTCGATGCCACCACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTAATc                        <  2:25701/136‑1 (MQ=37)
                                                                                           gagaACGCCACAGGCGGTAATTTCGAGCAACTGACCCCGAACCCGACCCCGGCTCAGACCACGATCAACGACTCGGTCGATGCCACCACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTAATCa                       <  2:333160/136‑1 (MQ=37)
                                                                                            ccaccGCCACTGGCGTTAATTTCTAGCAACTCACCCCGAACCCGCCCCCGGCTCAGACCACGTTCAACGACTCGGTCGATGCCACCACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTAATCAc                      <  2:233116/132‑1 (MQ=11)
                                                                                            agaACGCCACAGGCGTTAATTTCGAGCAACTGACCCCGAATCCGACCCCGGCTCAGACCACGATCAACGACTCGGTCGATGCCACCACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTAATCAc                      <  2:231112/136‑1 (MQ=37)
                                                                                            agaACGCCACAGGCGGTAATTTCGAGCAACTGACCCCGAATCCGACCCCGGCTCAGACCACGATCAACGACTCGGTCGATGCCACCACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCAc                      <  2:458417/136‑1 (MQ=17)
                                                                                             gaACGCCACAGGCGGTAATTTCGAGCAACTGACCCCGAATCCGACCCCGGCTCAGACCACGATCAACGACTCGGTCGATGCCACCACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCAcc                     >  1:169676/1‑136 (MQ=17)
                                                                                             gaACGCCACAGGCGGTAATTTCGAGCAACTGACCCCGAATCCGACCCCGGCTCAGACCACGATCAACGACTCGGTCGATGCCACCACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCAcc                     >  1:167417/1‑136 (MQ=17)
                                                                                                 gCCACAGGGGGTAATTTCGAGAAACTGACCCCGAATCCGACCCCGGCTCAGACCACGATCAACGACTCGGTCGATGCCACCACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACTTaca                 <  2:60225/136‑1 (MQ=255)
                                                                                                 gCCACAGGCGGTAATTTCGAGCAACTGACCCCGAATCCGACCCCGGCTCAGACCACGATCAACGACTCGGTCGATGCCACCACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTAATCACCTaca                 <  2:436581/136‑1 (MQ=37)
                                                                                                      aGGCGGTAATTTCGAGCAACTGACCCCGAATCCGACCCCGGCTCAGACCACGATCAACGACTCGGTCGATGCCACCACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTAATCACCTACACCGTg            >  2:319210/1‑136 (MQ=35)
                                                                                                      aGGCGGTAATTTCGAGCAACTGACCCCGAATCCGACCCCGGCTCAGACCACGATCAACGACTCGGTCGATGCCACCACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTAATCACCTACACCGTg            <  2:241954/136‑1 (MQ=35)
                                                                                                      aGGCGGTAATTTCGAGCAACTGACCCCGAATCCGACCCCGGCTCAGACCACGATCAACGACTCGGTCGATGCCACCACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTAATCACCTACACCGTg            >  1:352207/1‑136 (MQ=35)
                                                                                                              aTTTCGAGCAACTGACCCCGAATCCGACCCCGGCTCAGACCACGATCAACGACTCGGTCGATGCCACCACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTAATCACCTACACCGTGACCCTGAg    >  2:102127/1‑136 (MQ=34)
                                                                                                                ttCGAGCAACTGACCCCGAATCCGACCCCGGCTCAGACCACGATCAACGACTCGGTCGATGCCACCACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTAATCACCTACACCGTGACCCTGAGCa  >  2:426916/1‑136 (MQ=34)
                                                                                                                ttCGAGAAATTTACCCCTAATCGCACCCCGGGTCAGACAACCTTAAACGTCTCGGTCGTTGCCACCACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTAATCACCTACACCGTGACCCTGAGCa  <  1:462618/136‑1 (MQ=11)
                                                                                                                                    |                                                                                                                   
CGTTGAAGCCGGCAAGACCCAGGGCAGCGTCGATTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACCGTCAGCGTCACCATCGAGAACGCCACAGGCGGTAATTTCGAGCAACTGACCCCGAATCCGACCCCGGCTCAGACCACGATCAACGACTCGGTCGATGCCACCACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTAATCACCTACACCGTGACCCTGAGCA  >  NC_002947/196957‑197204

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: