Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	197,590	T→C	31.9%	D1032D (GAT→GAC)	PP_0168 →	putative surface adhesion protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	197,590	0	T	C	31.9%	29.4 / 17.6	28	D1032D (GAT→GAC)	PP_0168	putative surface adhesion protein
Reads supporting (aligned to +/- strand):  ref base T (8/11);  new base C (5/4);  total (13/15)
Fisher's exact test for biased strand distribution p-value = 6.89e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.85e-01

CCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACCTACACCGTGATCCTGAGCAATCCTGCCCAGACGCCGGTAACCGTGACCCTGTCCAACGGCCAAACCATTACCGTTGAAGCCGGCA‑AGAC‑‑CCAGGGCAGCGTCGATTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACTGTCAGCGTCACCATCGAGAACGCCACAGGCGGTAATTTCGAGCAACTGACCCCGAACCCGACGCCGGCTCAGACCACGATCACCG  >  NC_002947/197443‑197717                                                                                                                                                       |                                                                                                                                ccTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACCTACACCGTGATCCTGAGCAATCCTGCCCAGACGCCGGTAACCGTGACCCTGTCCAACGGCCAAACCATTACCGTTGAAGCCGGCA‑AGAC‑‑CCAGGGCAGCGTCGAttt                                                                                                                               >  1:418499/1‑150 (MQ=39)  cTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACCTACACCGTGATCCTGAGCAATCCTGCCCCGACGGCGGTAACCGTGACCCTGTCCAACGGCCAAACCATTACCGTTGAAGCCGGCA‑AGAC‑‑CCAGGGCAGCGTCGATTTc                                                                                                                              <  2:463099/150‑1 (MQ=255)          gAGCCCGTCGGTCACCGAAGGCGGCGTGATCACCTACACCGTGATCCTGAGCAATCCTGCCCAGACGCCGGTAACCGTGACCCTGTCCAACGGCCAAACCATTACCGTTGAAGCCGGCA‑AGAC‑‑CCAGGGCAGCGTCGATTTCCAGAcccc                                                                                                                      <  1:477300/150‑1 (MQ=38)                gTCGGTCACCGAAGGCGGCGTGATCACCTACACCGTGATCCTGAGCAATCCTGCCCAGACGCCGGTAACCGTGACCCTGTCCAACGGCCAAACCATTACCGTTGAAGCCGGCA‑AGAC‑‑CCAGGGCAGCGTCGATTTCCAGACCCCGGCGaa                                                                                                                <  1:351638/150‑1 (MQ=39)                               cggcgTGATCACCTACACCGTGATCCTGAGCAATCCTGCCCAGACGCCGGTAACCGTGACCCTGTCCAACGGCCAAACCATTACCGTTGAAGCCGGCA‑AGAC‑‑CCAGGGCAGCGTCGATTTCCAGACCCCGGCGAATGACGTCTAcaacaa                                                                                                 <  2:476454/150‑1 (MQ=39)                                  cgTGATCACCTACACCGTGATCCTGAGCAATCCTGCCCAGACGCCGGTAACCGTGACCCTGTCCAACGGCCAAACCATTACCGTTGAAGCCGGCA‑AGAC‑‑CCAGGGCAGCGTCGATTTCCAGACCCCGGCGAATGACGTCTACAACAACgg                                                                                              >  2:62754/1‑150 (MQ=38)                                                   tGATCCTGAGCAATCCTGCCCAGACGCCGGTAACCGTGACCCTGTCCAACGGCCAAACCATTACCGTTGAAGCCGGCA‑AGAC‑‑CCAGGGCAGCGTCGATTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACCGTCAGCGTCa                                                                             <  2:286858/150‑1 (MQ=18)                                                       ccTGAGCAATCCTGCCCAGACGCCGGTAACCGTGACCCTGTCCAACGGCCAAACCATTACCGTTGAAGCCGGCA‑AGAC‑‑CCAGGGCAGCGTCGATTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACCGTCAGCGTCACCAt                                                                         <  1:51716/150‑1 (MQ=17)                                                       ccTGAGCAATCCTGCCCAGACGCCGGTAACCGTGACCCTGTCCAACGGCCAAACCATTACCGTT‑‑‑GCCGTAACAGTCGTGCACTTCAGCGTCGATTTCCAGACCCCGGCGAATGACGTCTACAACAACGGCTCGACTGTCAGCGTCACCAt                                                                         >  2:142114/1‑150 (MQ=11)                                                         tGAGCAATCCTGCCCAGACGCCGGTAACCGTGACCCTGTCCAACGGCCAAACCATTACCGTTGAAGCCGGCA‑AGAC‑‑CCAGGGCAGCGTCGATTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACTGTCAGCGTCACCATCg                                                                       >  1:128307/1‑150 (MQ=38)                                                         tGAGCAATCCTGCCCAGACGCCGGTAACCGTGACCCTGTCCAACGGCCAAACCATTACCGTTGAAGCCGGCA‑AGAC‑‑CCAGGGCAGCGTCGATTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACCGTCAGCGTCACCATCg                                                                       >  1:423440/1‑150 (MQ=17)                                                          gAGCAATCCTGCCCAGACGCCGGTAACCGTGACCCTGTCCAACGGCCAAACCATTACCGTTGAAGCCGGCA‑AGAC‑‑CCAGGGCAGCGTCGATTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACTGTCAGCGTCACCATCga                                                                      <  1:62754/150‑1 (MQ=39)                                                                       cAGACGCCGGTAACCGTGACCCTGTCCAACGGCCAAACCATTACCGTTGAAGCCGGCA‑AGAC‑‑CCAGGGCAGCGTCGATTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACTGTCAGCGTCACCATCGAGAACGCCACAGGc                                                         >  2:459625/1‑150 (MQ=39)                                                                        aGACGCCGGTGACCGTGACCCTGTCCAACGGCCAAACGATTACCGTTGAAGCCGGCA‑AAAC‑‑CCAGGGCAGCGTCGACTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACTGTCAGCGTCACCATCGAGAACGCCACAGGCg                                                        <  2:373933/150‑1 (MQ=11)                                                                        aGACGCCGGTAACCGTGACCCTGTCCAACGGCCAAACCATTACCGTTGAAGCCGGCA‑AGAC‑‑CCAGGGCAGCGTCGATTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACCGTCAGCGTCACCATCGAGAACGCCACAGGCg                                                        <  1:416698/150‑1 (MQ=17)                                                                         gACGCCGGTAACCGTGACCCTGTCCAACGGCCAAACCATTACCGTTGAAGCCGGCA‑AGAC‑‑CCAGGGCAGCGTCGATTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACTGTCAGCGTCACCATCGAGAACGCCACAGGCgg                                                       <  2:418499/150‑1 (MQ=39)                                                                                  gaCCGTGACCCTGTCCAACGGCCAAACGATTACCGTTGCAGCCGGCG‑AAAC‑‑GCAGGTCAGCGTCGACTTCCAGACCCCGGCGAATGACGTCTACAACACCGGTTCGACTGTCAGCGTAACAATCGAGAACGCCACAGGCGGTATTTTCGa                                              <  2:399172/149‑1 (MQ=11)                                                                                           ccTGTCCAACGGCCAAACGATTACCGTTGAAGCCGGCA‑AAAC‑‑CCAGGGCAGCGTCGACTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACTGTCAGCGTCACCATCGAGAACGCCACAGGCGGTAATTTCGAGCAACTGAc                                     >  2:282579/1‑150 (MQ=12)                                                                                            cTGTCCAACGGCCAAACGATTACCGTTGAAGCCGGCA‑AAAC‑‑CCAGGGCAGCGTCGACTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACTGTCAGCGTCACCATCGAGAACGCCACAGGCGGTAATTTCGAGCAACTGAcc                                    >  1:48356/1‑150 (MQ=12)                                                                                             tGTCCAACGGCCAAACCATTACCGTTGAAGCCGGCA‑AGAC‑‑CCAGGGCAGCGTCGATTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACTGTCAGCGTCACCATCGAGAACGCCACAGGCGGTAATTTCGAGCAACTGAc                                     >  1:405224/1‑148 (MQ=32)                                                                                             tGTCCAACGGCCAAACCATTACCGTTGAAGCCGGCA‑AGAC‑‑CCAGGGCAGCGTCGATTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACTGTCAGCGTCACCATCGAGAACGCCACAGGCGGTAATTTCGAGCAACTGAc                                     <  2:405224/148‑1 (MQ=32)                                                                                                      gCCAAACGATTACCGTTGAAGCCGGCA‑AAAC‑‑CCAGGGCAGCGTCGACTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACTGTCAGCGTCACCATCGAGAACGCCACAGGCGGTAATTTCGa                                              >  1:265389/1‑130 (MQ=255)                                                                                                      gCCAAACGATTACCGTTGAAGCCGGCA‑AAAC‑‑CCAGGGCAGCGTCGACTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACTGTCAGCGTCACCATCGAGAACGCCACAGGCGGTAATTTCGa                                              <  2:265389/130‑1 (MQ=255)                                                                                                                    gTTGAAGCCGGCA‑AGAC‑‑CCAGGGCAGCGTCGATTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACTGTCAGCGTCACCATCGAGAACGCCACAg                                                           <  2:120249/103‑1 (MQ=33)                                                                                                                    gTTGAAGCCGGCA‑AGAC‑‑CCAGGGCAGCGTCGATTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACTGTCAGCGTCACCATCGAGAACGCCACAg                                                           >  1:120249/1‑103 (MQ=33)                                                                                                                        aaGCCGGCA‑AGAC‑‑CCAGGGCAGCGTCGATTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACTGTCAGCGTCACCATCGAGAACGCCACAGGCGGTAATTTCGAGCAACTGACCCCGAACCCGACGCCGGCTCAGACCACGa        >  2:440031/1‑150 (MQ=38)                                                                                                                          gCCGGCA‑AAAC‑‑CCAGGGCAGCGTCGACTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACTGTCAGCGTCACCATCGAGAACGCCACAGGCGGTAATTTCGAGCAACTGACCCCGAAc                              <  2:266883/126‑1 (MQ=255)                                                                                                                          gCCGGCA‑AAAC‑‑CCAGGGCAGCGTCGACTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACTGTCAGCGTCACCATCGAGAACGCCACAGGCGGTAATTTCGAGCAACTGACCCCGAAc                              >  1:266883/1‑126 (MQ=255)                                                                                                                             ggCA‑AAAC‑‑CCAGGGCAGCGTCGACTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACTGTCAGCGTCACC‑TCGAGAACGCCACAGGCGGTAATTTCGAGCAACTGACCCCGAACCCGACCCCGGCTCAGACCACGATCATCg  >  1:112218/1‑150 (MQ=11)                                                                                                                                                       |                                                                                                                                CCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACCTACACCGTGATCCTGAGCAATCCTGCCCAGACGCCGGTAACCGTGACCCTGTCCAACGGCCAAACCATTACCGTTGAAGCCGGCA‑AGAC‑‑CCAGGGCAGCGTCGATTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACTGTCAGCGTCACCATCGAGAACGCCACAGGCGGTAATTTCGAGCAACTGACCCCGAACCCGACGCCGGCTCAGACCACGATCACCG  >  NC_002947/197443‑197717

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