Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 197,572 G→A 27.3% K1026K (AAG→AAA PP_0168 → putative surface adhesion protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947197,5720GA27.3% 45.1 / 19.2 33K1026K (AAG→AAAPP_0168putative surface adhesion protein
Reads supporting (aligned to +/- strand):  ref base G (12/12);  new base A (5/4);  total (17/16)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.95e-01

CCACCACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACCTACACCGTGATCCTGAGCAATCCTGCCCAGACGCCGGTAACCGTGACCCTGTCCAACGGCCAAACCATTACCGTTGAAGCCGGCAAGAC‑‑CCAGGGCAGCGTCGATTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACTGTCAGCGTCACCATCGAGAACGCCACAGGCGGTAATTTCGAGCAACTGACCCCGAACCCGACGCCGGCTCAGACCACGATCACCG  >  NC_002947/197430‑197717
                                                                                                                                               |                                                                                                                                                   
ccaccaCCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACCTACACCGTGATCCTGAGCAATCCTGCCCAGACGCCGGTAACCGTGACCCTGTCCAACGGCCAAACCATTACCGTTGAAGCCGGCAAGAC‑‑CCAg                                                                                                                                             >  2:340996/1‑149 (MQ=39)
ccaccaCCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACCTACACCGTGATCCTGAGCAATCCTGCCCAGACGCCGGTAACCGTGACCCTGTCCAACGGCCAAACCATTACCGTTGAAGCCGGCAAGAC‑‑CCAg                                                                                                                                             <  1:340996/149‑1 (MQ=39)
  accaccGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACCTACACCGTGATCCTGAGCAATCCTGCCCAGACGCCGGTAACCGTGACCCTGTCCAACGGCCAAACCATTACCGTTGAAGCCGGCAAGAC‑‑CCAGGGc                                                                                                                                          >  2:51716/1‑150 (MQ=38)
          gACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACCTACACCGTGATCTTGGGCAATCTTGCCCAGACGCCGGTAACCGTGACCCTGTCCAACGGCCAAATCATTACCGTTGAAGCCGGCAAGAC‑‑CCAGGGCGGCGTCGa                                                                                                                                  >  1:88701/1‑150 (MQ=11)
             ccTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACCTACACCGTGATCCTGAGCAATCCTGCCCAGACGCCGGTAACCGTGACCCTGTCCAACGGCCAAACCATTACCGTTGAAGCCGGCAAGAC‑‑CCAGGGCAGCGTCGAttt                                                                                                                               >  1:418499/1‑150 (MQ=39)
              cTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACCTACACCGTGATCCTGAGCAATCCTGCCCCGACGGCGGTAACCGTGACCCTGTCCAACGGCCAAACCATTACCGTTGAAGCCGGCAAGAC‑‑CCAGGGCAGCGTCGATTTc                                                                                                                              <  2:463099/150‑1 (MQ=255)
                      gAGCCCGTCGGTCACCGAAGGCGGCGTGATCACCTACACCGTGATCCTGAGCAATCCTGCCCAGACGCCGGTAACCGTGACCCTGTCCAACGGCCAAACCATTACCGTTGAAGCCGGCAAGAC‑‑CCAGGGCAGCGTCGATTTCCAGAcccc                                                                                                                      <  1:477300/150‑1 (MQ=38)
                            gTCGGTCACCGAAGGCGGCGTGATCACCTACACCGTGATCCTGAGCAATCCTGCCCAGACGCCGGTAACCGTGACCCTGTCCAACGGCCAAACCATTACCGTTGAAGCCGGCAAGAC‑‑CCAGGGCAGCGTCGATTTCCAGACCCCGGCGaa                                                                                                                <  1:351638/150‑1 (MQ=39)
                                           cggcgTGATCACCTACACCGTGATCCTGAGCAATCCTGCCCAGACGCCGGTAACCGTGACCCTGTCCAACGGCCAAACCATTACCGTTGAAGCCGGCAAGAC‑‑CCAGGGCAGCGTCGATTTCCAGACCCCGGCGAATGACGTCTAcaacaa                                                                                                 <  2:476454/150‑1 (MQ=39)
                                              cgTGATCACCTACACCGTGATCCTGAGCAATCCTGCCCAGACGCCGGTAACCGTGACCCTGTCCAACGGCCAAACCATTACCGTTGAAGCCGGCAAGAC‑‑CCAGGGCAGCGTCGATTTCCAGACCCCGGCGAATGACGTCTACAACAACgg                                                                                              >  2:62754/1‑150 (MQ=38)
                                                               tGATCCTGAGCAATCCTGCCCAGACGCCGGTAACCGTGACCCTGTCCAACGGCCAAACCATTACCGTTGAAGCCGGCAAGAC‑‑CCAGGGCAGCGTCGATTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACCGTCAGCGTCa                                                                             <  2:286858/150‑1 (MQ=18)
                                                                   ccTGAGCAATCCTGCCCAGACGCCGGTAACCGTGACCCTGTCCAACGGCCAAACCATTACCGTTGAAGCCGGCAAGAC‑‑CCAGGGCAGCGTCGATTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACCGTCAGCGTCACCAt                                                                         <  1:51716/150‑1 (MQ=17)
                                                                   ccTGAGCAATCCTGCCCAGACGCCGGTAACCGTGACCCTGTCCAACGGCCAAACCATTACCGTT‑‑‑GCCGTAACAGTCGTGCACTTCAGCGTCGATTTCCAGACCCCGGCGAATGACGTCTACAACAACGGCTCGACTGTCAGCGTCACCAt                                                                         >  2:142114/1‑150 (MQ=11)
                                                                     tGAGCAATCCTGCCCAGACGCCGGTAACCGTGACCCTGTCCAACGGCCAAACCATTACCGTTGAAGCCGGCAAGAC‑‑CCAGGGCAGCGTCGATTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACTGTCAGCGTCACCATCg                                                                       >  1:128307/1‑150 (MQ=38)
                                                                     tGAGCAATCCTGCCCAGACGCCGGTAACCGTGACCCTGTCCAACGGCCAAACCATTACCGTTGAAGCCGGCAAGAC‑‑CCAGGGCAGCGTCGATTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACCGTCAGCGTCACCATCg                                                                       >  1:423440/1‑150 (MQ=17)
                                                                      gAGCAATCCTGCCCAGACGCCGGTAACCGTGACCCTGTCCAACGGCCAAACCATTACCGTTGAAGCCGGCAAGAC‑‑CCAGGGCAGCGTCGATTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACTGTCAGCGTCACCATCga                                                                      <  1:62754/150‑1 (MQ=39)
                                                                                   cAGACGCCGGTAACCGTGACCCTGTCCAACGGCCAAACCATTACCGTTGAAGCCGGCAAGAC‑‑CCAGGGCAGCGTCGATTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACTGTCAGCGTCACCATCGAGAACGCCACAGGc                                                         >  2:459625/1‑150 (MQ=39)
                                                                                    aGACGCCGGTGACCGTGACCCTGTCCAACGGCCAAACGATTACCGTTGAAGCCGGCAAAAC‑‑CCAGGGCAGCGTCGACTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACTGTCAGCGTCACCATCGAGAACGCCACAGGCg                                                        <  2:373933/150‑1 (MQ=11)
                                                                                    aGACGCCGGTAACCGTGACCCTGTCCAACGGCCAAACCATTACCGTTGAAGCCGGCAAGAC‑‑CCAGGGCAGCGTCGATTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACCGTCAGCGTCACCATCGAGAACGCCACAGGCg                                                        <  1:416698/150‑1 (MQ=17)
                                                                                     gACGCCGGTAACCGTGACCCTGTCCAACGGCCAAACCATTACCGTTGAAGCCGGCAAGAC‑‑CCAGGGCAGCGTCGATTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACTGTCAGCGTCACCATCGAGAACGCCACAGGCgg                                                       <  2:418499/150‑1 (MQ=39)
                                                                                              gaCCGTGACCCTGTCCAACGGCCAAACGATTACCGTTGCAGCCGGCGAAAC‑‑GCAGGTCAGCGTCGACTTCCAGACCCCGGCGAATGACGTCTACAACACCGGTTCGACTGTCAGCGTAACAATCGAGAACGCCACAGGCGGTATTTTCGa                                              <  2:399172/149‑1 (MQ=11)
                                                                                                       ccTGTCCAACGGCCAAACGATTACCGTTGAAGCCGGCAAAAC‑‑CCAGGGCAGCGTCGACTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACTGTCAGCGTCACCATCGAGAACGCCACAGGCGGTAATTTCGAGCAACTGAc                                     >  2:282579/1‑150 (MQ=12)
                                                                                                        cTGTCCAACGGCCAAACGATTACCGTTGAAGCCGGCAAAAC‑‑CCAGGGCAGCGTCGACTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACTGTCAGCGTCACCATCGAGAACGCCACAGGCGGTAATTTCGAGCAACTGAcc                                    >  1:48356/1‑150 (MQ=12)
                                                                                                         tGTCCAACGGCCAAACCATTACCGTTGAAGCCGGCAAGAC‑‑CCAGGGCAGCGTCGATTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACTGTCAGCGTCACCATCGAGAACGCCACAGGCGGTAATTTCGAGCAACTGAc                                     >  1:405224/1‑148 (MQ=32)
                                                                                                         tGTCCAACGGCCAAACCATTACCGTTGAAGCCGGCAAGAC‑‑CCAGGGCAGCGTCGATTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACTGTCAGCGTCACCATCGAGAACGCCACAGGCGGTAATTTCGAGCAACTGAc                                     <  2:405224/148‑1 (MQ=32)
                                                                                                                  gCCAAACGATTACCGTTGAAGCCGGCAAAAC‑‑CCAGGGCAGCGTCGACTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACTGTCAGCGTCACCATCGAGAACGCCACAGGCGGTAATTTCGa                                              <  2:265389/130‑1 (MQ=255)
                                                                                                                  gCCAAACGATTACCGTTGAAGCCGGCAAAAC‑‑CCAGGGCAGCGTCGACTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACTGTCAGCGTCACCATCGAGAACGCCACAGGCGGTAATTTCGa                                              >  1:265389/1‑130 (MQ=255)
                                                                                                                                gTTGAAGCCGGCAAGAC‑‑CCAGGGCAGCGTCGATTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACTGTCAGCGTCACCATCGAGAACGCCACAg                                                           <  2:120249/103‑1 (MQ=33)
                                                                                                                                gTTGAAGCCGGCAAGAC‑‑CCAGGGCAGCGTCGATTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACTGTCAGCGTCACCATCGAGAACGCCACAg                                                           >  1:120249/1‑103 (MQ=33)
                                                                                                                                    aaGCCGGCAAGAC‑‑CCAGGGCAGCGTCGATTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACTGTCAGCGTCACCATCGAGAACGCCACAGGCGGTAATTTCGAGCAACTGACCCCGAACCCGACGCCGGCTCAGACCACGa        >  2:440031/1‑150 (MQ=38)
                                                                                                                                      gCCGGCAAAAC‑‑CCAGGGCAGCGTCGACTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACTGTCAGCGTCACCATCGAGAACGCCACAGGCGGTAATTTCGAGCAACTGACCCCGAAc                              <  2:266883/126‑1 (MQ=255)
                                                                                                                                      gCCGGCAAAAC‑‑CCAGGGCAGCGTCGACTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACTGTCAGCGTCACCATCGAGAACGCCACAGGCGGTAATTTCGAGCAACTGACCCCGAAc                              >  1:266883/1‑126 (MQ=255)
                                                                                                                                         ggCAAAAC‑‑CCAGGGCAGCGTCGACTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACTGTCAGCGTCACCTCGAGAACGCCACAGGCGGTAATTTCGAGCAACTGACCCCGAACCCGACCCCGGCTCAGACCACGATCATCg  >  1:112218/1‑150 (MQ=11)
                                                                                                                                               |                                                                                                                                                   
CCACCACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACCTACACCGTGATCCTGAGCAATCCTGCCCAGACGCCGGTAACCGTGACCCTGTCCAACGGCCAAACCATTACCGTTGAAGCCGGCAAGAC‑‑CCAGGGCAGCGTCGATTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACTGTCAGCGTCACCATCGAGAACGCCACAGGCGGTAATTTCGAGCAACTGACCCCGAACCCGACGCCGGCTCAGACCACGATCACCG  >  NC_002947/197430‑197717

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