Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 197,089 T→C 14.3% 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,0890TC14.3% 78.7 / 9.3 45N865N (AAT→AACPP_0168putative surface adhesion protein
Reads supporting (aligned to +/- strand):  ref base T (17/19);  new base C (2/4);  total (21/24)
Fisher's exact test for biased strand distribution p-value = 6.73e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.83e-01

CGTTGAAGCCGGCAAGACCCAGGGCAGCGTCGATTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACCGTCAGCGTCACCATCGAGAACGCCACAGGCGGTAATTTCGAGCAACTGACCCCGAATCCGACCCCGGCTCAGACCACGATCAACGACTCGGTCGATGCCACCACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTAATCACCTACACCGTGACCCTGAGCAATCCTGCCCAGACG  >  NC_002947/196957‑197218
                                                                                                                                    |                                                                                                                                 
cGTTGAAGCCGGCAAGACCCAGGGCAGCGTCGATTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACCGTCAGCGTCACCATCGAGAACGCCACAGGCGGTAATTTCGAGCAACTGACCCCGAATCCg                                                                                                                                >  2:377785/1‑136 (MQ=34)
   tGAAGCCGGCAAGACCCAGGGCAGCGTCGATTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACCGTCAGCGTCACCATCGAGAACGCCACAGGCGGTAATTTCGAGCAACTGACCCCGAATCCGAcc                                                                                                                             <  2:170683/136‑1 (MQ=35)
     aaGCCGGCAAGACCCAGGGCAGCGTCGATTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACCGTCAGCGTCACCATCGAGAACGCCACAGGCGGTAATTTCGAGCAACTGACCCCGAATCCGAcccc                                                                                                                           <  1:464455/136‑1 (MQ=37)
     aaGCCGGCAAGACCCAGGGCAGCGTCGATTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACAGTCAGCGTCACCATCGAGAACGCCACAGGCGGTAATTTCGAGCAACTGACCCCGAAGCCGAcccc                                                                                                                           >  2:14746/1‑136 (MQ=17)
       gCCGGCAAGACCCAGGGCAGCGTCGATTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACCGTCAGCGTCACCATCGAGAACGCCACAGGCGGTAATTTCGAGCAACTGACCCCGAATCCGACCCCgg                                                                                                                         >  1:469732/1‑136 (MQ=37)
       gCCGGCAAGACCCAGGGCAGCGTCGATTTCCAGACCCCGGCGAATGACGTCTACAACAACGGATCGACCGTCAGCGTCACCATCGAGAACGCCACAGGCGGAAATTTCGAGCAACTGACCCCGAATCCCACCCCgc                                                                                                                         >  1:310256/1‑135 (MQ=255)
             aaGACCCAGGGCAGCGTCGATTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACCGTCAGCGTCACCATCGAGAACGCCACAGGCGGTAATTTCGAGCAACTGACCCCGAATCCGACCCCGGCTCAgc                                                                                                                   >  1:38005/1‑135 (MQ=37)
             aaGACCCAGGGCAGCGTCGATTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACCGTCAGCGTCACCATCGAGAACGCCACAGGCGGTAATTTCGAGCAACTGACCCCGAATCCGACCCCGGCTCAGa                                                                                                                   <  1:52865/136‑1 (MQ=37)
               gACCCAGGGCAGCGTCGTTTTCCATACCCCGTCGAATGAGGTCTACAACTTCGGTTCGTCCGTCTGCGTCACCATCTAGAAGGCCACTGGCGGTAATTTCGAGCAACTGACCCCGAATCCGACCCCGGCTCAGAcc                                                                                                                 <  2:310256/136‑1 (MQ=11)
               gACCCAGGGCAGCGTCGATTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACCGTCAGCGTCACCATCGAGAACGCCACAGGCGGTAATTTCGAGCAACTGACCCCGAATCCGACCCCGGCTCAGAcc                                                                                                                 >  1:417639/1‑136 (MQ=35)
               gACCCAGGGCAGCGTCGATTTCCAGACCCAGGCGAATGACGTCTACAACAACGGTTCGACCGTCAGCGTCACCATCGAGAACGCCACAGGCGGTAATTTCGAGCAACTGACCCCGAATCCGACCCCGGCTCAGAcc                                                                                                                 <  2:115825/136‑1 (MQ=25)
                           ggTCGATTTCCACACCCCGGCTAATGACGTCTACAACAACGGTTCGACCGTCAGCGTAACCATCGTTAACGCCACAGGCGTTATTTTCGAGCAACTGACCCCTAATCCGACCCCGGCTCACACCACGATCAACGAc                                                                                                     <  2:89150/135‑1 (MQ=11)
                              cGATTTCCATACCCCGGCGAATGACGTCTACAACAACGGTTCGACCGTCAGCGTCACCATCGAGAACGCCACAGGCGGTAATTTCGAGCAACTGACCCCGAATCCGACCCCGGCTCAGACCACGATCAACGACTcc                                                                                                  >  2:308954/1‑135 (MQ=21)
                              cGATTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACCGTCAGCGTCACCATCGAGAACGCCACAGGCGGTAATTTCGAGCAACTGACCCCGAATCCGACCCCGGCTCAGACCACGATCAACGACTCg                                                                                                  <  2:417639/136‑1 (MQ=34)
                              cGATTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACCGTCAGCGTCACCATCGAGAACGCCACAGGCGGTAATTTCGAGCAACTGACCCCGAATCCGACCCCGGCTCAGACCACGATCAACGACTCg                                                                                                  >  1:123219/1‑136 (MQ=34)
                                aTTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACCGTCAGCGTCACCATCGAGAACGCCACAGGCGGTAATTTCGAGCAACTGACCCCGAATCCGACCCCGGCTCAGACCACGATCAACGACTCGGt                                                                                                <  1:199946/136‑1 (MQ=34)
                                      aGACCCCGGCGAATGACGTCTACAACAACGGTTCGACCGTCAGCGTCACCATCGAGAACGCCACAGGCGGTAATTTCGAGCAACTGACCCCGAAACCGACCCCGGCTCAGACCACGATCAACGACTCGGTCGATGc                                                                                          >  2:174510/1‑136 (MQ=37)
                                        aCCCCGGCGAATGACGTCTACAACAACGGTTCGACCGTCAGCGTCACCATCGAGAACGCCACAGGCGGTAATTTCGAGCAACTGACCCCGAATCCGACGCCGGCTCAGACCACGATCAACGACTCGGTCGATGcca                                                                                        >  2:127149/1‑136 (MQ=21)
                                             ggCGAATGACGTCTACAACAACGGTTCGACTGTCAGCGTCACCATCGAGAACGCCACAGGCGGTAATTTCGAGCAACTGACCCCGAACCCGACCCCGGCTCAGACCACGATCAACGACTCGGTCGATGCCaccacc                                                                                   >  2:503096/1‑136 (MQ=21)
                                             ggCGAATGACGTCTACAACAACGGTTCGACCGTCAGCGTCACCATCGAGAACGCCACAGGCGGTAATTTCGAGCAACTGACCCCGAATCCGACCCCGGCTCAGACCACGATCAACGACTCGGt                                                                                                <  1:294196/123‑1 (MQ=35)
                                             ggCGAATGACGTCTACAACAACGGTTCGACCGTCAGCGTCACCATCGAGAACGCCACAGGCGGTAATTTCGAGCAACTGACCCCGAATCCGACCCCGGCTCAGACCACGATCAACGACTCGGt                                                                                                >  2:294196/1‑123 (MQ=35)
                                                  ttGAGGTCTTCATCAACTGTTCGACCGTCAGGGTAACCATCGAGAAGGCCACAGGCGTAATTTTCGAGCACGTGACCGCGAATCCGACGCCGGCCCATACCAGGTTCAACTACTCGGTCGATGCCACCACCGCGAc                                                                              <  2:484917/135‑1 (MQ=2)
                                                  aTGACGTCTACAACAACGGTTCGACCGTCAGCGTCACCATCGAGAACGCCACAGGCGGTAATTTCGAGCAACTGACCCCGAATCCGACCCCGGCTCAGACCACGATCAACGACTCGGTCGATGCCACCACCGCGAc                                                                              >  1:432028/1‑136 (MQ=37)
                                                           acaacaACGGTTCGACCGTCAGCGTCACCATCGAGAACGCCACAGGCGGTAATTTCGAGCAACTGACCCCGAATCCGACCCCGGCTCAGACCACGATCAACGACTCGGTCGATGCCACCACCGCGACCCTGACGGc                                                                     >  2:348631/1‑136 (MQ=37)
                                                            caacaaCGGTTCGACGGTCAGCGTCACCATCGAGAACGCCACAGGCGGTAATTTCGAGCAACTGACCCCGAATCCGa                                                                                                                               <  2:97544/77‑1 (MQ=21)
                                                                aacGGTTCGACCGTCAGCGTCACCATCGAGAACGCCACAGGCGGTAATTTCGAGCAACTGACCCCGAATCCGACCCCGGCTCAGACCACGATCAACGACTCGGTCGATGCCACCACCGCGACCCTGACGGCGAGcc                                                                >  1:267167/1‑136 (MQ=37)
                                                                       cGACCGTCAGCGTCACCATCGAGAACGCCACAGGCGGTAATTTCGAGCAACTGACCCCGAATCCGACCCCGGCTCAGACCACGATCAACGACTCGGTCGATGCCACCACCGCGACCCTGACGGCGAGCCCGTCGGt                                                         >  2:431908/1‑136 (MQ=37)
                                                                               aGCGTCACCATCGAGAACGCCACAGGCGTTAATTTCTAGCCACTGACCCCGAATCAGACCCGGGCTCATACCAGGTTCAACGACTCGGTCGATGCCACCACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAg                                                 <  2:396077/136‑1 (MQ=11)
                                                                                  gTCACCATCGAGAACGCCACAGGCGGTAATTTCGAGCAACTGACCCCGAATCCGACCCCGGCTCAGACCACGATCAACGACTCGGTCGATGCCACCACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAggcg                                              >  2:193195/1‑136 (MQ=34)
                                                                                    cACCATCGAGAACGCCACAGGCGTTAATTTCGAGCAACTGACCCCGAATCCGACCCCGGCTCAGACCACGATCAACGACTCGGTCGATGCCACCACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAggcggc                                            <  1:46483/136‑1 (MQ=21)
                                                                                      ccATCGAGAACGCCACTGGCGGTAATTTCGAGCAACTGACCCCCAAACCGACCCCGGCTCAGACCACGATCAACGACTCGGTCGATGCCACCACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGt                                          <  1:156810/136‑1 (MQ=12)
                                                                                      ccATCGACACCGCCACAGGCGTTAATTTCGAGCAACTGACCCCGACCCCGCCCCCGGCTCAGACCACGAGCAACGACTCGGTCGATGCCACCACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGt                                          <  1:11693/136‑1 (MQ=11)
                                                                                            ctaaCGCCACTGGGGTTAATTTCGAGCAAGTCACCCCGAATCCGACCCCGGCACACACCACCACCAACGACTCGGTCGATGCCACCACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTAATCAc                                    <  2:298488/134‑1 (MQ=11)
                                                                                            agaACGCCACAGGCGTTAATTTCGAGCAACTGACCCCGAATCCGACCCCGGCTCAGACCACGATCAACGACTCGGTCGATGCCACCACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTAATCAc                                    <  1:141303/136‑1 (MQ=37)
                                                                                            agaACGCCACAGGCGGTAATTTCGAGCAACTGACCCCGAATCCGACCCCGGCTCAGACCACGATCAACGACTCGGTCGATGCCACCACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTAATCAc                                    <  1:203799/136‑1 (MQ=37)
                                                                                                      cggCGTTAATTTCGAGCAACTGACCCCGACCCCGACCCCGGCTCAGACCACGATCAACGACTCGGTCGATGCCACCACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTAATCACCTACACCGTg                          <  1:14746/135‑1 (MQ=255)
                                                                                                           gTAATTTCGAGCAACTGACCCCGAACCCGACCCCGGCTCAGACCACGATCAACGACTCGGTCGATGCCACCACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTAATCACCTACACCGTGACCCt                     >  1:309202/1‑136 (MQ=21)
                                                                                                              aTTTCGAGCAACTGACCCCGAATCCGACCCCGGCTCAGACCACGATCAACGACTCGGTCGATGCCACCACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTAATCAc                                    >  2:45517/1‑118 (MQ=35)
                                                                                                              aTTTCGAGCAACTGACCCCGAATCCGACCCCGGCTCAGACCACGATCAACGACTCGGTCGATGCCACCACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTAATCAc                                    <  1:45517/118‑1 (MQ=35)
                                                                                                                ttCGAGCAACTGACCCCGAATCCGACCCCGGCTCAGACCACGATCAACGACTCGGTCGATGCCACCACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTAATCACCTACACCGTGACCCTGAGCa                <  1:431908/136‑1 (MQ=34)
                                                                                                                           gACCCCGAATCCGACCCCGGCTCAGACCACGATCAACGACTCGGTCGATGCCACCACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTAATCACCTACACCGTGACCCTGAGCAATCCTGCCCAg     >  2:387033/1‑136 (MQ=34)
                                                                                                                            cccccGAACCCGACCCCGGCTCAGACCACGATCAACGACTCGGTCGATGCCACCACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTAATCACCTACACCGTGACCCTGAGCAATCCTGCCCAGa    <  2:190284/135‑1 (MQ=17)
                                                                                                                              ccccAACCCGACCCCGGATCAGACCACTATCAACGACTCGGTCGATGCCACCACCGCGACCCTGACGGCGAGCCCGTCGTTCACCGAAGGCGGCGTAATCACATACACCGTGACCCTGAGCAATCCTGCCCAGACg  <  2:356690/136‑1 (MQ=11)
                                                                                                                                  aaTCCGACCCCGGCTCAGACCACGATCAACGACTCGGTCGATGCCACCACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAggcggc                                            >  2:393935/1‑90 (MQ=33)
                                                                                                                                  aaTCCGACCCCGGCTCAGACCACGATCAACGACTCGGTCGATGCCACCACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAggcggc                                            <  1:393935/90‑1 (MQ=33)
                                                                                                                                    |                                                                                                                                 
CGTTGAAGCCGGCAAGACCCAGGGCAGCGTCGATTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACCGTCAGCGTCACCATCGAGAACGCCACAGGCGGTAATTTCGAGCAACTGACCCCGAATCCGACCCCGGCTCAGACCACGATCAACGACTCGGTCGATGCCACCACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTAATCACCTACACCGTGACCCTGAGCAATCCTGCCCAGACG  >  NC_002947/196957‑197218

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: