Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 195,559 C→A 24.7% T355T (ACC→ACA PP_0168 → putative surface adhesion protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947195,5590CA24.7% 86.4 / 16.4 41T355T (ACC→ACAPP_0168putative surface adhesion protein
Reads supporting (aligned to +/- strand):  ref base C (14/16);  new base A (6/4);  total (21/20)
Fisher's exact test for biased strand distribution p-value = 7.16e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CGGTGACCGTGACCCTGTCCAACGGCCAGGTCATCACTGTTGAAGCCGGCAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCCAATGACGTCTACAACAACGGTTCGACCGTCAGCGTCACCATCGAGAACGCCACCGGCGGTAATTTCGAGCAACTGACCCCGAATCCGACCCCGGCTCAGACCACGATCAACGACTCGGTCGACACCACCACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACTTACACCGTGACCCTGAGCAA  >  NC_002947/195420‑195705
                                                                                                                                           |                                                                                                                                                   
cGGTGACCGTGACCCTGTCCAACGGCCAGGTCATCACTGTTGAAGCCGGCAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCCAATGACGTCTACAACAACGGTTCGACCGTCAGCGTCACCATCGAGAACGCCACCGGCGGTAAtt                                                                                                                                           >  1:601558/1‑150 (MQ=37)
            cccTGTCCAACGGCCAGGTCATCACTGTTGAAGCCGGCAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCCAATGACGTCTACAACAACGGTTCGACCGTCAGCGTCACCATCGAGAACGCCACCGGCGGTAATTTCGAGCAACTGa                                                                                                                               <  2:701938/150‑1 (MQ=255)
                    aaCGGCCAGGTCATCACTGTTGAAGCCGGCAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCCAATGACGTCTACAACAACGGTTCGACCGTCAGCGTCACCATCGAGAACGCCACCGGCGGTAATTTCGAGCAACTGACCCCGAAt                                                                                                                       <  2:634495/150‑1 (MQ=255)
                                atcaCTGTTGAAGCCGGCAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCCAATGACGTCTACAACAACGGTTCGACCGTCAGCGTCACCATCGAGAACGCCACCGGCGGTAATTTCGAGCAACTGACCCCGAATCCGACCCCGGCt                                                                                                           <  1:550115/150‑1 (MQ=255)
                                atcaCTGTTGAAGCCGGCAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCCAATGACGTCTACAACAACGGTTCGACCGTCAGCGTCACCATCGAGAACGCCACCGGCGATAATTTCGAGCAACTGACCCCGAATCCGACCCCGGCt                                                                                                           >  2:518213/1‑150 (MQ=255)
                                               ggCAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCCAATGACGTCTACAACAACGGTTCGACCGTCAGCGTCACCATCGAGAACGCCACCGGCGGTAATTTCGAGCAACTGACCCCGAATCCGACCCCGGCTCAGACCACGATCAAc                                                                                            <  2:583406/150‑1 (MQ=38)
                                               ggCAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCCAATGACGTCTACAACAACGGTTCGACCGTCAGCGTCACCATCGAGAACGCCACCGGCGATAATTTCGAGCAACTGACCCCGAATCCGACCCCGGCTCAGACCACGATCAAc                                                                                            <  1:518213/150‑1 (MQ=255)
                                                  acaaCCCAGGCGAGCGTCGAATGCCAGACCCCGGCCAATGAGGTCTATAACAACGGTTCGACCGTCAGCGTCACCATCGAGAACGCCACCGGCGGTAATTTCGAGCAACTGACCCCGAATCCGACCCCGGCTCAGACCACGATCAACGAc                                                                                        <  1:170912/148‑1 (MQ=11)
                                                               gCGTCGATTTCCAGACCCCGGCCAATGACGTCTACAACAACGGTTCGACCGTCAGCGTCACCATCGAGAACGCCACCGGCGGTAATTTCGAGCaa                                                                                                                                   <  1:194810/95‑1 (MQ=38)
                                                               gCGTCGATTTCCAGACCCCGGCCAATGACGTCTACAACAACGGTTCGACCGTCAGCGTCACCATCGAGAACGCCACCGGCGGTAATTTCGAGCaa                                                                                                                                   >  2:194810/1‑95 (MQ=38)
                                                                      tttCCAGACCCCGGCCAATGACGTCTACAACAACGGTTCGACCGTCAGCGTCACCATCGAGAACGCCACCGGCGGTAATTTCGAGCAAACCTCGCCCACTCCGACCCCGGCTCAGACCACGATCAACGACTCGGTCGACACCACCACcgc                                                                    <  2:592689/150‑1 (MQ=14)
                                                                          cAGACCCCGGCCAATGACGTCTACAACAACGGTTCGACCGTCAGCGTCACCATCGAGAACGCCACCGGCGGTAATTTCGAGCAACTGACCCCGAATCCGACCCCGGCTCAGACCACGATCAACGACTCGGTCGACACCACCACCGCGAcc                                                                >  1:81570/1‑150 (MQ=255)
                                                                          cAGACCACGGCAAATAACGTCTACAACAAAGGTTCGACCGTCAGCGTCACCATCGAGAACGCCACAGGCGGTAATTTCGAGCAACTGACCCCGAATCCGACCCCGGCTCAGACCACGATCAACGACTCGGTCGACACCACTACCGCGAcc                                                                >  2:628684/1‑150 (MQ=11)
                                                                           aGACCCCGGCCAATGACGTCTACAACAACGGTTCGACCGTCAGCGTCACCATCGAGAACGCTACGGGCGGCAATTTCGAGCAACTGACTCCGAATCCGACCCCGGGTCAGACACCGACAACAGGATCGGTCGACACCACCACCGCGAccc                                                               >  1:123878/1‑150 (MQ=11)
                                                                           aGACCCCGGCCAATGACGTCTACAACAACGGTTCGACCGTCAGCGTCACCATCGAGAACGCCACCGGCGGTAATTTCGAg                                                                                                                                      >  2:270378/1‑80 (MQ=255)
                                                                           aGACCCCGGCCAATGACGTCTACAACAACGGTTCGACCGTCAGCGTCACCATCGAGAACGCCACCGGCGGTAATTTCGAg                                                                                                                                      <  1:270378/80‑1 (MQ=255)
                                                                               cccGGCCAATGACGTCTACAACAACGGTTCGACCGTCAGCGTCACCATCGAGAACGCCACCGGCGGTAATTTCGAGCAACTGACCCCCAATCCGACCCCGGCTCAGACCACGATCAACGACTCGGTCGACACCACCACCGAGACCCTGAc                                                           >  2:6310/1‑150 (MQ=255)
                                                                                           cGTCTACAACAACGGTTCGACCGTCAGCGTCACCATCGAGAACGCCACCGGCGGTAATTTCGAGCAACTGACCCCGAATCCGACCCCGGCTCAGAACACGAACAACGACTCGGTCGACACCACCACCGCGACCCTGACGGCTAGCCCGTc                                               >  1:320609/1‑150 (MQ=255)
                                                                                            gTCTACAACAACGGTTCGACCGTCAGCGTCACCATCGAGAACGCCACAGGCGGTAATTTCGAGCAACTGACCCCGAATCCGACCCCGGCTCAGACCACGATCAACGACTCGGTCGACACCACTACCGCGACCCTGACGGCGAGCCCGTCg                                              <  1:628684/150‑1 (MQ=18)
                                                                                              cTACAACAACGGTTCGACCGTCAGCATCACCATCGAGAACGCCACCGGCGGTAATTTCGAGCAACTGACCCCGAATCCGACCCCGGCTCAGACCACGATCAACGACTCGGTCGACACCACCACCGCGACCCTg                                                             <  2:693103/133‑1 (MQ=255)
                                                                                              cTACAACAACGGTTCGACCGTCAGCATCACCATCGAGAACGCCACCGGCGGTAATTTCGAGCAACTGACCCCGAATCCGACCCCGGCTCAGACCACGATCAACGACTCGGTCGACACCACCACCGCGACCCTg                                                             >  1:693103/1‑133 (MQ=255)
                                                                                                  aacaacGGTTCGACCGTCAGCGTCACCATCGAGAACGCCACCGGCGGTAATTTCGAGCAACTGACCCCGAATCCGACCCCGGCTCAGACCACGATCAACGACTCGGTCGACACCACCACCGCGACCCTGACGGCGAGCCCGTCGGTCAcc                                        >  1:23351/1‑150 (MQ=39)
                                                                                                   acaacGGTTCGACCGTCAGCGTCACCATCGAGAACGCCACCGGCGGTAATTTCGAGCAAc                                                                                                                                  >  1:648702/1‑60 (MQ=38)
                                                                                                   acaacGGTTCGACCGTCAGCGTCACCATCGAGAACGCCACCGGCGGTAATTTCGAGCAAc                                                                                                                                  <  2:648702/60‑1 (MQ=38)
                                                                                                   acaacGGTTCGACCGTCAGCGTCACCAACGAGAACGCCACCGCCGGTAATTCCGAGCAACTGACCCCGAATCCGACCCCGGCTCAGACCACGATCAACAAATCGGTCGACAGCCGCACCGCGACCCTGACAGCGAGCCAGTCgagcacaa                                       >  2:379757/1‑143 (MQ=11)
                                                                                                         gTTCGACCGTCAGCGTCACCATCGAGAACGCCACCGGCGGTAATTTCGAGCAACTGACCCCGAATCCGACCCCGGCTCAGACCACGATCAACGACTCGGTCGACACCACCACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAggcg                                 <  2:601558/150‑1 (MQ=38)
                                                                                                         gTTCGACCGTCAGCGTCACCATCGAGAACGCCACCGGCGGTAATATCGCGCAACAGACCCCGAATCCGACCCCGGCTCAGACCACGATCAACGACGCGGTCGACACCACCACCACGGCCCTGACGGCGAGCCCGTCGCTCACCGAAggcg                                 >  2:29321/1‑150 (MQ=11)
                                                                                                          ttCGACCGTCAGCGTCACCATCGAGAACGCCACAGGCGGTAATTTCGAGCAACTGACCCCGAATCCGACCCCGGCTCAGACCACGATCAACGACTCGGTCGACACCACTACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAggcgg                                >  1:348121/1‑150 (MQ=12)
                                                                                                                cGTCAGCGTCACCATCGAGAACGCCACCGGCGGTAATTTCGAGCAACTGACCCCGAATCCGACCCCGGCTCAGACCACGATCAACGACTCGGTCGACACCACCACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGAt                          <  2:81570/150‑1 (MQ=40)
                                                                                                                     gCGTCACCATCGAGAACGCCACAGGCGGTAATTTCGAGCAACTGACCCCGAATCCGACCCCGGCTCAGACCACGATCAACGACTCGGTCGACACCACTACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACtt                     <  2:480960/150‑1 (MQ=12)
                                                                                                                          aCCATCGAGAACGCCACAGGCGGTAATTTCGAGCAACTGACCCCGAATCCGACCCCGGCTCAGACCACGATCAACGACTCGGTCGACACCACTACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACTTACAcc                >  2:321561/1‑150 (MQ=17)
                                                                                                                          aCCATCGAGAACGCCACAGGCGGTAATTTCGAGCAACTGACCCCGAATCCGACCCCGGCTCAGACCACGATCAACGACTCGGCCGACACCACTACCGGCACCCTGACGCCGAGCCCGTCGGCCACCACAGCCGGCGAGATCACTTAcccc                >  2:110299/1‑150 (MQ=11)
                                                                                                                             aTCGAGAACGCCACCGGCGGTAATTTCGAGCAACTGACCCCGAATCCGACCCCGGCTCAGACCACGATCAACGACTCGGTCGAcaccaccac                                                                       >  2:429237/1‑92 (MQ=255)
                                                                                                                             aTCGAGAACGCCACCGGCGGTAATTTCGAGCAACTGACCCCGAATCCGACCCCGGCTCAGACCACGATCAACGACTCGGTCGAcaccaccac                                                                       <  1:429237/92‑1 (MQ=255)
                                                                                                                              tCGAGAACGCCACCGGCGGTAATTTCGAGCAACTGACCCCGAATCCGACCCCGGCTCAGACCACGATCAACGACTCGGTCGACACCACCACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACTTACACCGTGa            <  2:23351/150‑1 (MQ=255)
                                                                                                                              tCGAGAACGCCACAGGCGGTAATTTCGAGCAACTGACCCCGAATCCGACCCCGGCTCAGACCACGATCAACGACTCGGTCGACACCACTACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACTTACACCGTGa            >  1:470444/1‑150 (MQ=18)
                                                                                                                                gagaACGCCACAGGCGGTAATTTCGAGCAACTGACCCCGAATCCGACCCCGGCTCAGACCACGATCAACGACTCGGTCGACACCACTACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACTTACACCGTGAcc          <  2:420816/150‑1 (MQ=12)
                                                                                                                                  gaACGCCACAGGCGGTAATTTCGAGCAACTGACCCCGAATCCGACCCCGGCTCAGACCACGATCAACGACTCGGTCGACACCACTACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACTTACACCGTGACCCt        >  1:132429/1‑150 (MQ=12)
                                                                                                                                    aCGCCACAGGCGGTAATTTCGAGCAACTGACCCCGAATCCGACCCCGGCTCAGACCACGATCAACGACTCGGTCGACACCACTACCGCGACCCTGACGGCGGGCCCGTCGGTCACCGAAGGCGGCGTGATCACTTACACCGTGACCCTGa      <  1:321561/150‑1 (MQ=12)
                                                                                                                                        cacaGGCGGTAATTTCGAGCAACTGACCCCGAATCCGACCCCGGCTCAGACCACGATCAACGACTCGGTCGACACCACTACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACTTACACCGTGACCCTGAGCaa  >  1:228979/1‑150 (MQ=14)
                                                                                                                                        cACCGGCGGTAATTTCGAGCAACTGACCCCGAATCCGACCCCGGCTCAGACCACGATCAACGACTCGGTCGACACCACCACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACTTACACCGTGACCCTGa      <  2:168989/146‑1 (MQ=255)
                                                                                                                                        cACCGGCGGTAATTTCGAGCAACTGACCCCGAATCCGACCCCGGCTCAGACCACGATCAACGACTCGGTCGACACCACCACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACTTACACCGTGACCCTGa      >  1:168989/1‑146 (MQ=255)
                                                                                                                                           |                                                                                                                                                   
CGGTGACCGTGACCCTGTCCAACGGCCAGGTCATCACTGTTGAAGCCGGCAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCCAATGACGTCTACAACAACGGTTCGACCGTCAGCGTCACCATCGAGAACGCCACCGGCGGTAATTTCGAGCAACTGACCCCGAATCCGACCCCGGCTCAGACCACGATCAACGACTCGGTCGACACCACCACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACTTACACCGTGACCCTGAGCAA  >  NC_002947/195420‑195705

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