Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	196,879	G→A	23.5%	V795V (GTG→GTA)	PP_0168 →	putative surface adhesion protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	196,879	0	G	A	23.5%	56.8 / 14.1	34	V795V (GTG→GTA)	PP_0168	putative surface adhesion protein
Reads supporting (aligned to +/- strand):  ref base G (15/11);  new base A (4/4);  total (19/15)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.12e-01

CGTCAGCGTCACCATCGAGAACGCCACTGGCGGCAACTTCGAACAGCTGACCCCGAACCCGACGCCGGCTCAGACCACGATCAACGACTCGGTCGATGCCACCACCGCG‑ACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACCTACACCGTGACCCTGAGCAATCCTGCCCAGACGCCGGTGACAGTGACCCT  >  NC_002947/196732‑196935                                                                                                                                                     |                                                         cGTCAGCGTCACCATCGAGAACGCCACTGGCGGCAACTTCGAACAGCTGACCCCGAACCCGACGCCGGCTCAGACCACGATCAACGACTCGGTCGATGCCACCACCGCG‑ACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGa                                                         >  5:92097/1‑149 (MQ=32)     aGCGTCAACATCGAGAACGCCACTGGCGGCAACTTCGAACAGCTGACCCCGAACCCGACGCCGGCTCAGACCACGATCAACGACTCGGTCGATGCCACCACCGCG‑ACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCa                                                      <  6:333263/148‑1 (MQ=17)      gCGTCACCATCGAGAACGCCACTGGCGGCAACTTCGAACAGCTGACCCCGAACCCGACGCCGGCTCAGACCACGATCAACGACTCGGTCGATGCCACCACCGCG‑ACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCAcc                                                    <  3:27621/149‑1 (MQ=32)             cATCGAGAACGCCACTGGCGGCAACTTCGAACAGCTGACCCCGAACCCGACGCCGGCTCAGACCACGATCAACGACTCGGTCGATGCCACCACCGCG‑ACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACCTACACCg                                             >  7:41094/1‑149 (MQ=32)             cATCGAGAACGCCACTGGCGGCAACTTCGAACAGCTGACCCCGAACCCGACGCCGGCTCAGACCACGATCAACGACTCGGTCGATGCCACCACCGCG‑ACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACCTACACCg                                             >  6:269746/1‑149 (MQ=32)             cATCGAGAACGCCACTGGCGGCAACTTCGAACAGCTGACCCCGAACCCGACGCCGGCTCAGACCACGATCAACGACTCGGTCGATGCCACCACCGCG‑ACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACCTACACCg                                             >  2:384979/1‑149 (MQ=32)              aTCGAGAACGCCACTGGCGGCAACTTCGAACAGCTGACCCCGAACCCGACGCCGGCTCAGACCACGATCAACGACTCGGTCGATGCCACCACCGCG‑ACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACCTACACCGt                                            <  8:41094/149‑1 (MQ=32)                        ccACTGGCGGCAACTTCGAACAGCTGACCCCGAACCCGACGCCGGCTCAGACCACGATCAACGCCTCGGTCGATGCCACCACCGCG‑ACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACCTACACCGTGACCCTGAg                                   >  6:2497/1‑148 (MQ=21)                                ggcAACTTCGAACAGCTGACCCCGAACCCGACGCCGGCTCAGACCACGATCAACGACTCGGTCGATGCCACCACCGCG‑ACCCTGACGGCGAGCGCGTCGGTCACCGAAGGCGGCGTGATCACCTACACCGTGACCCTGAGCAATCCTGc                          <  4:112301/149‑1 (MQ=17)                                ggcAACTTCGAACAGCTGACCCCGAACCCGACGCCGGCTCAGACCACGATCAACGACTCGGTCGATGCCACCACCGCG‑ACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACCTACACCGTGACCCTGAGCAATCCTGc                          >  6:276728/1‑149 (MQ=32)                                 gcAACTTCGAACAGCTGACCCCGAACCCGACGCCGGCTCAGACCACGATCAACGACTCGGTCGATGCCACCACCGCG‑ACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACCTACACCGTGACCCTGAGCAATCCTGcc                         <  3:180160/149‑1 (MQ=32)                                   aaCTTCGAACAGCTGACCCCGAACCCGACGCCGGCTCAGACCACGATCAACGACTCGGTCGATGCCACCACCGCG‑ACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACCTACACCGTGACCCTGAGCAATCCTGcc                         >  2:152757/1‑147 (MQ=32)                                     cTTCGAACAGCTGACCCCGAACCCGACGCCGGCTCAGACCACGATCAACGACTCGGTCGATGCCACCACCGCG‑ACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACCTACACCGTGACCCTGAGCAATCGTGCCCAGa                     <  5:115658/149‑1 (MQ=17)                                      ttCGAACAGCTGACCCCGAACCCGACGCCGGCTCAGACCACGATCAACGACTCGGTCGATGCCACCACCGCG‑ACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACCTACACCGTTACCCTGa                                    >  8:150549/1‑133 (MQ=17)                                      ttCGAACAGCTGACCCCGAACCCGACGCCGGCTCAGACCACGATCAACGACTCGGTCGATGCCACCACCGCG‑ACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACCTACACCGTGACCCTGa                                    <  7:150549/133‑1 (MQ=32)                                      ttCGAACAGCTGACCCCGAACCCGACGCCGGCTCAGACCACGATCAACGACTCGGTCGATGCCACCACCGCG‑ACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACCTACACCGTGACCCTGAGCAATCCTGCCCAGAc                    >  4:251878/1‑149 (MQ=32)                                         gAACAGCTGACCCCGAACCCGACGCCGGCTCAGACCACGATCAACGACTCGGTCGATGCCACCACCGCG‑ACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACCTACACCGTGACCCTGAGCAATCCTGCCCAGACGc                  <  1:384979/148‑1 (MQ=32)                                          aaCAGCTGACCCCGAACCCGACGCCGGCTCAGACCACGATCAACGACTCGGTCGATGCCACCACCGCG‑ACCCTGACGGCGAGTCCGTCGGTCACCGAAGGCGGCGTGATCACCTACACCGTGACCCTGAGCAATCCTGCCCAGACGCCg                <  1:114096/149‑1 (MQ=17)                                          aaCAGCTGACCCCGAACCCGACGCCGGCTCAGACCACGATCAACGACTCGGTCGATGCCACCACCGCG‑ACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACCTACACCGTGACCCTGAGCAATCCTGCCCAGACGCCg                >  3:232773/1‑149 (MQ=32)                                          aaCAGCTGACCCCGAACCCGACGCCGGCTCAGACCACGATCAACGACTCGGTCGATGCCACCACCGCG‑ACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACCTACACCGTGACCCTGAGCAATCCTGCCCAGACGCCg                >  7:276339/1‑149 (MQ=32)                                          aaCAGCTGACCCCGAACCCGACGCCGGCTCAGACCACGATCAACGACTCGGTCGATGCCACCACCGCG‑ACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACCTACACCGTGACCCTGAGCAATCCTGCCCAGACGCCg                >  5:63092/1‑149 (MQ=32)                                              gCTGACCCCGAACCCGACGCCGGCTCAGACCACGATCAACGACTCGGTCGATGCCACCACCGCG‑ACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACCTACACCGTGACCCTGAGCAATCCTGCCCAGACGCCGGTGa            <  4:232773/149‑1 (MQ=32)                                                  aCCCCGAACCCGACGCCGGCTCAGACCACGATCAACGACTCGGTCGATGCCACCACCGCG‑ACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACCTACACCGTGACCCTGAGCAATCCTGCCCAGACGCCGGTGACAGt        >  1:275946/1‑149 (MQ=32)                                                  aCCCCGAACCCGACCCCGGCTCAGACCACGATCAACGACTCGGTCGATGCCACCACCGCG‑ACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTAATCACCTACACCGTGACCCTGAGCAATCCTGCCCAGACGCCGGTGACa          <  5:298595/147‑1 (MQ=11)                                                  aCCCCGAACCCGACCCCGGCTCAGACCACGATCAACGACTCGGTCGATGCCACCACCGCG‑ACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTAATCACCTACACCGTGACCCTGAGCAATCCTGCCCAGACGCCGGTGACa          <  7:108991/147‑1 (MQ=11)                                                    cccGAACCCGACGCCGGCTCAGACCACGATCAACGACTCGGTCGATGCCACCACCGCG‑ACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACCTACACCGTGACCCTGAGCAATCCTGCCCAGACGCCGGTGACAGTGa      <  3:362202/149‑1 (MQ=32)                                                    cccGAACCCGACGCCGGCTCAGACCACGATCAACGACTCGGTCGATGCCACCACCGCG‑ACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACCTACACCGTGACCCTGAGCAATCCTGCCCAGACGCCGGTGACAGTGa      >  1:209846/1‑149 (MQ=32)                                                    cccGAACCCGACGCCGGCTCAGACCACGATCAACGACTCGGTCGATGCCACCACCGCG‑ACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACCTACACCGTGACCCTGAGCAATCCTGCCCAGACGCCGGTGACAGTGa      >  4:362202/1‑149 (MQ=32)                                                    cccGAACCCGACCCCGGCTCAGACCACGATCAACGACTCGGTCGATGCCACCACCGCGTACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTAATCACCTACACCGTGACCCTGAGCAATCCTGCCCAGACGCCGGTGACAGt        >  1:323727/1‑148 (MQ=255)                                                    cccGAACCCGACCCCGGCTCAGACCACGATCAACGACTCGGTCGATGCCACCACCGCG‑ACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTAATCACCTACACCGTGACCCTGAGCAATCCTGCCCAGACGCCGGTGACAGTGa      >  1:141527/1‑149 (MQ=17)                                                     ccGAACCCGACCCCGGCTCAGACCACGATCAACGACTCGGTCGATGCCACCACCGCG‑ACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTAATCACCTACACCGTGACCCTGAGCAATCCTGCCCAGACGCCGGTGACAGTGAc     <  5:96547/149‑1 (MQ=17)                                                      cGAACCCGACCCCGGCTCAGACCACGATCAACGACTCGGTCGATGCCACCACCGCG‑ACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTAATCACCTACACCGTGACCCTGAGCAATCCTGCCCAGACGCCGGTGACAGTGAcc    >  6:181133/1‑149 (MQ=17)                                                      cGAACCCGACCCCGGCTCAGACCACGATCAACGACTCGGTCGATGCCACCACCGCG‑ACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTAATCACCTACACCGTGACCCTGAGCAATCCTGCCCAGACGCCGGTGACAGTGAcc    >  5:373094/1‑149 (MQ=17)                                                        aaCCCGACCCCGGCTCAGACCACGATCAACGACTCGGTCGATGCCACCACCGCG‑ACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGTCGTAATCACCTACACCGTGACCCTGAGCAATCCTGCCCAGACGCCGGTGACAGTGACCCt  <  3:218008/149‑1 (MQ=12)                                                                                                                                                     |                                                         CGTCAGCGTCACCATCGAGAACGCCACTGGCGGCAACTTCGAACAGCTGACCCCGAACCCGACGCCGGCTCAGACCACGATCAACGACTCGGTCGATGCCACCACCGCG‑ACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACCTACACCGTGACCCTGAGCAATCCTGCCCAGACGCCGGTGACAGTGACCCT  >  NC_002947/196732‑196935

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 21 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑