Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	2,612,160	C→A	22.2%	R560S (CGC→AGC)	PP_2284 →	tail tubular protein B

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	2,612,160	0	C	A	22.2%	35.0 / 3.5	18	R560S (CGC→AGC)	PP_2284	tail tubular protein B
Reads supporting (aligned to +/- strand):  ref base C (8/6);  new base A (2/2);  total (10/8)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.64e-02

TACTGGGACGGTGACACGAAGGTGCAAAGCTCCTGGAGCAAGTGGACCTTCCCGGACACCGATACCATCCTGAGCGCGGAGTTCCTGGACTCCGAAGTGTTCCTGGCGATCAACCGCGCCGATGGCTTGTACTTCGAAAAGCTCACTGTAGCCACGGACAGCTTGGGCACCAATGAGCCGTATCTGGTTCACCTGGACCGAAAGCAATACGTCACGAAAGACACCCTGAGTTACGCG  >  NC_002947/2612046‑2612282                                                                                                                   |                                                                                                                           tACTGGGACGGTGACACGAAGGTGCAAAGCTCCTGGAGCAAGTGGACCTTCCCGGACACCGATACCATCCTGAGCGCGGAGTTCCTGGACTCCGAAGTGTTCCTGGCGATCAACCGCGCCGATGGCTTGTACTTCg                                                                                                       >  1:217448/1‑136 (MQ=255)     gggACGGTGACACGAAGGTGCAAAGCTCCTGGAGCAAGTGGACCTTCCCGGACACCGATACCATCCTGAGCGCGGAGTTCCTGGACTCCGAAGTGTTCCTGGCGATCAACCGCGCCGATGGCTTGTACTTCGaaaa                                                                                                   <  1:26002/136‑1 (MQ=255)     gggACGGTGACACGAAGGTGCAAAGCTCCTGGAGCAAGTGGACCTTCCCGGACACCGATACCATCCTGAGCGCGGAGTTCCTGGACTCCGAAGTGTTCCTGGCGATCAACAGCGCCGATGGCTTGTACTTCGaaaa                                                                                                   >  1:230395/1‑136 (MQ=255)     gggACGGTGACACGAAGGTGCAAAGCTCCTGGAGCAAGTGGACCTTCCCGGACACCGATACCATCCTGAGCGCGGAGTTCCTGGACTCCGAAGTGTTCCTGGCGATCAACAGCGCCGATGGCTTGTACTTCGaaaa                                                                                                   >  1:286929/1‑136 (MQ=255)           gTGACACGAAGGTGCAAAGCTCCTGGAGCAAGTGGACCTTCCCGGACACCGATACCATCCTGAGCGCGGAGTTCCTGGACTCCGAAGTGTTCCTGGCGATCAACCGCGCCGATGGCTTGTACTTCGAAAAGCTCAc                                                                                             >  1:380135/1‑136 (MQ=255)           gTGACACGAAGGTGCAAAGCTCCTGGAGCAAGTGGACCTTCCCGGACACCGATACCATCCTGAGCGCGGAGTTCCTGGACTCCGAAGAGATCCTGTCGAACAACCGCGCCGATGGCATGTACTTCGAAAAGCTCAc                                                                                             >  1:335916/1‑136 (MQ=255)           gTGACACGAAGGTGCAAAGCTCCTGGAGCAAGTGGACCTTACCGGACACCGATACCATCCTGAGCGCGGAGTTCCTGGACTCCGAAGTGTTCCTGGCGATCAACCGCGCCGATGGCTTGTACTTCGAAAAGCTCAc                                                                                             >  1:420640/1‑136 (MQ=255)                acGAAGGTGCAAAGCTCCTGGAGCAAGTGGACCTTCCCGGACACCGATACCATCCTGAGCGCGGAGTTCCTGGACTCCGAAGTGTTCCTGGCGATCAACCGCGCCGATGGCTTGTACTTCGAAAAGCTCACTGTAg                                                                                        >  1:105692/1‑136 (MQ=255)                          aaaGCTCCTGGAGCAAGTGGACCTTCCCGGACACCGATACCATCCTGAGCGCGGAGTTCCTGGACTCCGAAGTGTTCCTGGCGATCAACCGCGCCGATGGCTTGTACTTCGAAAAGCTCACTGTAGCCACGGACAg                                                                              <  2:380135/136‑1 (MQ=255)                          aaaGCTCCTGGAGCAAGTGGACCTTACCGGACACCGATACCATCCTGAGCGCGGAGTTCCTGGACTCCGAAGTGTTCCTGGCGATCAACCGCGCCGATGGCTTGTACTTCGAAAAGCTCACTGTAGCCACGGACAg                                                                              <  2:420640/136‑1 (MQ=255)                             gCTCCTGGAGCAAGTGGACCTTCCCGGACACCGATACCATCCTGAGCGCGGAGTTCCTGGACTCCGAAGTGTTCCTGGCGATCAACCGCGCCGATGGCTTGTACTTCGAAAAGCTCACTGTAGCCACGGACAGCtt                                                                           <  2:294599/136‑1 (MQ=255)                                                        acacCGATACCATCCTGAGCGCGGAGTTCCTGGACTCCGAAGTGTTCCTGGCGATCAACCGCGCCGATGGCTTGTACTTCGAAAAGCTCACTGTAGCCACGGACAGCTTGGGCACCAATGAGCCGTATCTGGTTCa                                                <  1:223437/136‑1 (MQ=255)                                                                       tGAGCGCGGAGTTCCTGGACTCCGAAGTGTTCCTGGCGATCAACCGCGCCGATGGCTTGTACTTCGAAAAGCTCACTGTAGCCACGGACAGCTTGGGCACCAATGAGCCGTATCTGGTTCACCTGGACCGAAAGCa                                 >  1:79233/1‑136 (MQ=255)                                                                            gcgGAGTTCCTGGACTCCGAAGTGTTCCTGGCGATCAACCGCGCCGATGGCTTGTACTTCGAAAAGCTCACTGTAGCCACGGACAGCTTGGGCACCAATGAGCCGTATCTGGTTCACCTGGACCGAAAGCAATACg                            >  1:134895/1‑136 (MQ=255)                                                                             tggtgTTGCTTGACTCCGAAGTGTTCGTGGCGTTCTACCGCGGCGATGGCTTGTAGTTCGAAAAGCTCACTGTAGCCAGGGACAGCTTGGGCACCAATGAGCCGTATCTGGTTGACCTGGACCGAAAGCAATACGt                           <  2:335916/132‑1 (MQ=255)                                                                                 gTTCCTGGACTCCGAAGTGTTCCTGGCGATCAACCGCGCCGATGGCTTGTACTTCCAAAAGCTCACTGTAGCAACGGACAGCTTGGGCACCAATGAGCCGTATCTGCTTCACCTGGACCGAAAGCAATACGTCACg                       >  2:335736/1‑136 (MQ=255)                                                                                                      ccTGGCGATCAACAGCGCCGATGGCTTGTACTTCGAAAAGCTCACTGTAGCCACGGACAGCTTGGGCACCAATGAGCCGTATCTGGTTCACCTGGACCGAAAGCAATACGTCACGAAAGACACCCTGAGTTAcgcg  <  2:230395/136‑1 (MQ=255)                                                                                                      ccTGGCGATCAACAGCGCCGATGGCTTGTACTTCGAAAAGCTCACTGTAGCCACGGACAGCTTGGGCACCAATGAGCCGTATCTGGTTCACCTGGACCGAAAGCAATACGTCACGAAAGACACCCTGAGTTAcgcg  <  2:286929/136‑1 (MQ=255)                                                                                                                   |                                                                                                                           TACTGGGACGGTGACACGAAGGTGCAAAGCTCCTGGAGCAAGTGGACCTTCCCGGACACCGATACCATCCTGAGCGCGGAGTTCCTGGACTCCGAAGTGTTCCTGGCGATCAACCGCGCCGATGGCTTGTACTTCGAAAAGCTCACTGTAGCCACGGACAGCTTGGGCACCAATGAGCCGTATCTGGTTCACCTGGACCGAAAGCAATACGTCACGAAAGACACCCTGAGTTACGCG  >  NC_002947/2612046‑2612282

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