Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 197,632 T→C 17.8% T1046T (ACT→ACC PP_0168 → putative surface adhesion protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947197,6320TC17.8% 69.5 / 4.5 33T1046T (ACT→ACCPP_0168putative surface adhesion protein
Reads supporting (aligned to +/- strand):  ref base T (17/10);  new base C (3/3);  total (20/13)
Fisher's exact test for biased strand distribution p-value = 6.59e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ACCTACACCGTGATCCTGAGCAATCCTGCCCAGACGCCGGTAACCGTGACCCTGTCCAACGGCCAAACCATTACCGTTGAAGCCGGCAAGACCCAGGGCAGCGTCGATTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACTGTCAGCGTCACCATCGAGAACGCCACAGGCGGTAATTTCGAGCAACTGACCCCGAACCCGACGCCGGCTCAGACCACGATCACCGACTCGGTCGACACCACTACTGCGACCCTGACCGCCTCGCCAAGCGTGACCGAAGGCGGCGTGAT  >  NC_002947/197483‑197781
                                                                                                                                                     |                                                                                                                                                     
aCCTACACCGTGATCCTGAGCAATCCTGCCCAGACGCCGGTAACCGTGACCCTGTCCAACGGCCAAACCATTACCGTTGAAGCCGGCAAGACCCAGGGCAGCGTCGATTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACt                                                                                                                                                       >  1:287871/1‑150 (MQ=38)
      acCGTGATCCTGAGCAATCCTGCCCAGACGCCGGTAACCGTGACCCTGTCCAACGGCCAAACCATTACCGTTGAAGCCGGCAAGACCCAGGGCAGCGTCGATTTCCAGACCCCNGCGAATGACGTCTACAACAACGGTTCGACCGTCAGc                                                                                                                                                 >  2:647826/1‑150 (MQ=25)
                tGAGCAATCCTGCCCAGACGCCGGTAACCGTGACCCTGTCCAACGGCCAAACCATTACCGTTGAAGCCGGCAAGACCCAGGGCAGCGTCGATTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACCGTCAGCGTCACCATCg                                                                                                                                       >  1:242684/1‑150 (MQ=18)
                 gAGCAATCCTGCCCAGACGCCGGTAACCGTGACCCTGTCCAACGGCCAAACCATTACCGTTGAAGCCGGCAAGACCCAGGGCAGCGTCGATTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACTGTCAGCGTCACCATCga                                                                                                                                      >  2:277188/1‑150 (MQ=38)
                 gAGCAATCCTGCCCAGACGCCGGTAACCGTGACCCTGTCCAACGGCCAAACCATTACCGTTGAAGCCGGCAAGACCCAGGGCAGCGTCGATTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACCGTCAGCGTCACCATCga                                                                                                                                      <  2:260554/150‑1 (MQ=17)
                 gAGCAATCCTGCCCAGACGCCGGTAACCGTGACCCTGTCCAACGGCCAAACCATTACCGTTGAAGCCGGCAAGACCCAGGGCAGCGTCGATTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTGCGACTGTCAGCGTCACCATCga                                                                                                                                      <  2:669628/150‑1 (MQ=38)
                  aGCAATCCTGCCCAGACGCCGGTAACCGTGACCCTGTCCAACGGCCAAACCATTACCGTTGAAGCCGGCAAGACCCAGGGCAGCGTCGATTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACTGTCAGCGTCACCATCgag                                                                                                                                     <  2:310421/150‑1 (MQ=39)
                     aaTCCTGCCCAGACGCCGGTAACCGTGACCCTGTCCAACGGCCAAACCATTACCGTTGAAGCCGGCAAGACCCAGGGCAGCGTCGATTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACCGTCAGCGTCACCATCGAGAAc                                                                                                                                  <  1:794521/150‑1 (MQ=17)
                              cAGACGCCGGTGACCGTGACCCTGTCCAACGGCCAAACGATTACCGTTGAAGCCGGCAAAACCCAGGGCAGCGTCGACTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACTGTCAGCGTCACCATCGAGAACGCCACAGGc                                                                                                                         <  1:229821/150‑1 (MQ=11)
                                    ccGGTAACCGTGACCCTGTCCAACGGCCAAACCATTACCGTTGAAGCCGGCAAGACCCAGGGCAGCGTCGATTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACTGTCAGCGTCACCATc                                                                                                                                        >  1:199142/1‑129 (MQ=40)
                                    ccGGTAACCGTGACCCTGTCCAACGGCCAAACCATTACCGTTGAAGCCGGCAAGACCCAGGGCAGCGTCGATTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACTGTCAGCGTCACCATc                                                                                                                                        <  2:199142/129‑1 (MQ=40)
                                      ggTAACCGTGACCCTGTCCAACGGCCAAACCATTACCGTTGAAGCCGGCAAGACCCAGGGCAGCGTCGATTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACCGTCAGCGTCACCATCga                                                                                                                                      <  1:153796/129‑1 (MQ=17)
                                      ggTAACCGTGACCCTGTCCAACGGCCAAACCATTACCGTTGAAGCCGGCAAGACCCAGGGCAGCGTCGATTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACCGTCAGCGTCACCATCga                                                                                                                                      >  2:153796/1‑129 (MQ=18)
                                       gtgaCCGTGACCCTGTCCAACGGCCAAACGATTACCGTTGAAGCCGGCAAAACCCAGGGCAGCGTCGACTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACTGTCAGCGTCACCATCGAGAACGCCACAGGCGGTAATTTc                                                                                                                <  2:365480/147‑1 (MQ=255)
                                                 cccTGTCCAACGGCCAAACCATTACCGTTGAAGCCGGCAAGACCCAGGGCAGCGTCGATTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGCCTGTCAGCGTCACCATCGAGAACGCCACAGGCGGTAATTTCGAGCAACTGa                                                                                                      >  1:617331/1‑150 (MQ=17)
                                                             gCCAAACGATTACCGTTGAAGCCGGCAAAACCCAGGGCAGCGTCGACTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACTGTCAGCGTCACCATCGAGAACGCCACAGGCGGTAATTTCGAGCAACTGACCCCGAACCCGa                                                                                          >  1:778161/1‑150 (MQ=255)
                                                               cAAACGATTACCGTTGAAGCCGGCAAAACCCAGGGCAGCGTCGACTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACTGTCAGCGTCACCATCGAGAACGCCACAGGCGGTAATTTCGa                                                                                                              >  2:550813/1‑128 (MQ=255)
                                                               cAAACGATTACCGTTGAAGCCGGCAAAACCCAGGGCAGCGTCGACTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACTGTCAGCGTCACCATCGAGAACGCCACAGGCGGTAATTTCGa                                                                                                              <  1:550813/128‑1 (MQ=255)
                                                                       tACCGTTGAAGCCGGCAAGACCCAGGGCAGCGTCGATTTCCAGACCCCGGTGAATGACGTCTACAACAACGGTTCGACTGTCAGCGTCACCATCGAGAACGCCACAGGCGGTAATTTCGAGCAACTGACCCCGAACCCGACGCCGGCTcg                                                                                >  2:313024/1‑149 (MQ=37)
                                                                       tACCGTTGAAGCCGGCAAAACCCAGGGCAGCGTCGACTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACTGTCAGCGTCACCATCGAGAACGCCACAGGCGGTAATTTCGAGCAACTGACCCCGAACCCGACCCCGGCTCa                                                                                <  1:591679/150‑1 (MQ=255)
                                                                                 gCCGGCAAGACCCAGGGCAGCGTCGATTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACTGTCAGCGTCACCATCGAGAACGCCACAGGCGGTAATCTCGAGCAACTGACCCCGAACCCGACGCCGGCTCAGACCACGATc                                                                      >  2:575864/1‑150 (MQ=38)
                                                                                     gCAAGACCCAGGGCAGCGTCGATTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACTGTCAGCGTCACCATCGAGAACGCCACAGGCGGTAATTTCGAGCAACTGACCCCGAACCCGACGCCGGCTCAGACCACGATCACCg                                                                  >  1:387270/1‑150 (MQ=40)
                                                                                      cAAGACCCAGGGCAGCGTCGATTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACTGTCAGCGTCACCATCGAGAACGCCACAGGCGGTAATTTCGAGCAACTGACCCCGAACCCGACGCCGGCTCAGACCACGATCACCGa                                                                 >  1:298412/1‑150 (MQ=255)
                                                                                                    gCGTCGATTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACTGTCAGCGTCACCATCGAGAACGCCACAGGCGGTAATTTCGAGCAACTGACCCCGAACCCGACGCCGGCTCAGACCACGATCACCGACTCGGTCGAcacca                                                   <  2:298412/150‑1 (MQ=255)
                                                                                                         gacttCCACACCCCGGCGACTGCCGTCTACAACCACGGTTCTATTGTCAGCGTCACCATCGAGAACGCCACAGGCGGTAATTTCGAGCAACTGACCCCGAACCCGCCGCCGGCTCAGACCACGATCACCGACTCGGTCGACACCactact                                              <  2:173806/147‑1 (MQ=14)
                                                                                                                     ccGGCGAATGACGTCTACAACAACGGTTCGACTGTCAGCGTCACCATCGAGAACGCCACAGGCGGTAATTTCGAGCAACTGACCCCGAACCCGACGCCGGCTCAGACCACGATCACCGACTCGGTCGACACCACTACTGCGACCCTGAcc                                  >  2:184117/1‑150 (MQ=255)
                                                                                                                       ggCGAATGACGTCTACAACAACGGTTCGACTGTCAGCGTCACCATCGAGAACGCCACAGGCGGTAATTTCGAGCAACTGACCCCGAACCCGACGCCGGCTCAGACCACGATCACCGACTCGGTCGACACCACTACTGCGACCCTGACCGc                                >  1:341919/1‑150 (MQ=255)
                                                                                                                                cGTCTACAACAACGGTTCGACTGTCAGCGTCACCATCGAGAACGCCACAGGCGGTAATTTCGAGCAACTGACCCCGAACCCGACGCCGGCTCAGACCACGATCACCGACTCGGTCGACACCACTACTGCGACCCTGACCGCCTCGCCAAg                       >  1:552437/1‑150 (MQ=255)
                                                                                                                                 gTCTACAACAACGGTTCGACTGTCAGCGTCACCATCGAGAACGCCACAGGCGGTAATTTCGAGCAACTGACCCCGAACCCGACGCCGGCTCAGACCACCATCACCCACTCGGTCGACACCACTACTGCGACCCTTACCGCCTCGCCactc                      >  2:182500/1‑147 (MQ=255)
                                                                                                                                    tACAACAACGGTTCGACTGTCAGCGTCACCATCGAGAACGCCACAGGCGGTAATTTCGAGCAACTGACCCCGAACCCGACGCCGGCTCAGACCACGATCACCGACTCGGTCGACACCACTACTGCGACCCTGACCGCCTCGCCAAGCGTg                   <  1:277188/150‑1 (MQ=255)
                                                                                                                                     acaacaACGGTTCGACTGTCAGCGTCACCATCGAGAACGCCACAGGCGGTAATTTCGAGCAACTGACCCCGAACCCGACGCCGGCTCAGACCACGATCACCGACTCGGTCGACACCACTACTGCGACCCTGACCGCCTCGCCAAGCGTGa                  >  2:415330/1‑150 (MQ=255)
                                                                                                                                       aacaacGGTTCGACTGTCAGCGTCACCATCGAGAACGCCACAGGCGGTAATTTCGAGCAACTGACCCCGAACCCGACGCCGGCTCAGACCACGATCACCGACTCGGTCGACACCACTACTGCGACCCTGACCGCCTCGCCAAGCGTGAcc                >  1:435714/1‑150 (MQ=255)
                                                                                                                                           acGGTTCGACTGTCAGCGTCACCATCGAGAACGCCACAGGCGGTAATTTCGAGCAACTGACCCCGAACACGACGCCGGCTCAGACCACGATCACCGACTCGGTCGACACCACTACTGCGACCATGACCGCCTCGCCAAGCGTGACCGAAg            >  1:521365/1‑150 (MQ=255)
                                                                                                                                                  gACTGTCAGCGTCACCATCGAGAACGCCACAGGCGGTAATTTCGAGCAACTGACCCCGAACCCGACGCCGGCTCAGACCACGATCACCGACTCGGTCGACACCACTACTGCGACCCTGACCGCCTCGCCAAGCGTGACCGAAGGCGGCGt     >  1:411173/1‑150 (MQ=255)
                                                                                                                                                     tGTCAGCGTCACCATCGAGAACGCCACAGGCGGTAATTTCGAGCAACTGACCCCGAACCCGACGCCGGCTCAGACCACGATCACCGACTCGGTCGACACCACTACTGCGACCCTGACCGCCTCGCCAAGCGTGACCGAAGGCGGCGTGAt  <  2:496609/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
ACCTACACCGTGATCCTGAGCAATCCTGCCCAGACGCCGGTAACCGTGACCCTGTCCAACGGCCAAACCATTACCGTTGAAGCCGGCAAGACCCAGGGCAGCGTCGATTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACTGTCAGCGTCACCATCGAGAACGCCACAGGCGGTAATTTCGAGCAACTGACCCCGAACCCGACGCCGGCTCAGACCACGATCACCGACTCGGTCGACACCACTACTGCGACCCTGACCGCCTCGCCAAGCGTGACCGAAGGCGGCGTGAT  >  NC_002947/197483‑197781

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