Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,791,818	C→A	16.0%	P76T (CCC→ACC)	trxC →	thioredoxin

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	5,791,818	0	C	A	16.0%	67.3 / 3.4	25	P76T (CCC→ACC)	trxC	thioredoxin
Reads supporting (aligned to +/- strand):  ref base C (11/10);  new base A (2/2);  total (13/12)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.07e-01

CGCTGCAAAAAGGATGTGCTGTTGAGCACCCCCTTCGAACTTACCGAAGCCAGCTACGCCAGCCAGATCAAAGGCGACCTGCCGCTGCTGGTCGACATCTGGGCCGACTGGTGCGGCCCGTGCAAGTCATTCGCCCCCACCTTCGAAGAGGCCGCCCGGAAGCTTGCTGGCCGCTGCCGCCTGGCCAAGCTCGACAGCGAAGCCAACCGTAACCTCGCAGGGCAACTGGGCATACG  >  NC_002947/5791683‑5791918                                                                                                                                        |                                                                                                     cGCTGCAAAAAGGATGGGCTGTTGAGCACCCCCTTCGAACTTACCGAAGCCAGCTACGCCAGCCAGATCAAAGGCGACCTGCCGCTGCTGGTCGACATCTGGGCCGACTGGTGCGGCCCGTGCAAGTCATTCGccc                                                                                                      <  2:647371/136‑1 (MQ=255)     gCAAAAAGGATGTGCTGTTGAGCACCCCCTTCGAACTTACCGAAGCCAGCTACGCCAGCCAGATCAAAGGCGACCTGCCGCTGCTGGTCGACATCTGGGCCGACTGGTGCGGCCCGTGCAAGTCATTCGCCCCCAc                                                                                                  >  1:103462/1‑136 (MQ=255)     gCAAAAAGGATGTGCTGTTGAGCACCCACTTCGAACTTACCGAAGCCAGCTACGCCAGCCAGATCAAAGGCGACCTGCCGCTGCTGGTCGACATCTGGGCCGACTGGTGCGGCCCGTGCAAGTCATTCGCCCCCAc                                                                                                  >  1:659248/1‑136 (MQ=255)              aTGTGCTGTTGAGCACCCCCTTCGAACTTACCGAAGCCAGCTACGCCAGCCAGATCAAAGGCGACCTGCCGCTGCTGGTCGACATCTGGGCCGACTGGTGCGGCCCGTGCAAGTCATTCGCCCCCACCTTCGAaga                                                                                         <  2:438779/136‑1 (MQ=255)                 tgCTGTTGAGCACCCCCTTCGAACTTACCGAAGCCAGCTACGCCAGCCAGATCAAAGGCGACCTGCCGCTGCTGGTCGACATCTGGGCCGACTGGTGCGGCCCGTGCAAGTCATTCGCCCCCACCTTCGAAGAGgc                                                                                      >  1:204619/1‑136 (MQ=255)                                ccTTCTAACTTACCGAAGCCAGCTACGCCAGCCAGATCAAAGGCGACCTGCCGCTGCTGGTCGACATCTGGGCCGACTGGTGCGGCCCGTGCAAGTCATTCGCCCCCACCTTCGAAGAGGCCGCCCGGAAGCTTGc                                                                       >  2:369208/1‑136 (MQ=255)                                 cTTCGAACTTACCGAAGCCAGCTACGCCAGCCAGATCAGAGGCGACCTGCCGCTGCTGGTCGACATCTGGGCCGACTGGTGCGGCCCGTGCAAGTCATTCGCCCCCACCTTCGAAGAGGCCGCCCGGAAGCTTGCt                                                                      >  2:125168/1‑136 (MQ=255)                                   tCGAACTTACCGAAGCCAGCTACGCCAGCCAGATCAAAGGCGACCTGCCGCTGCTGGTCGACATCTGGGCCGACTGGTGCGGCCCGTGCAAGTCATTCGCCCCCACCTTCGAAGAGGCCGCCCGGAAGCTTGCTgg                                                                    >  2:398267/1‑136 (MQ=255)                                           aCCGAAGCCAGCTACGCCAGCCAGATCAAAGGCGACCTGCCGCTGCTGGTCGACATCTGGGCCGACTGGTGCGGCCCGTGCAAGTCATTCGCCACCACCTTCGAAGAGGCCGCCCGGAAGCTTGCTGGCCGCTgcc                                                            >  1:669611/1‑136 (MQ=255)                                           aCCGAAGCCAGCTACGCCAGCCAGATCAAAGGCGACCTGCCGCTGCTGGTCGACATCTGGGCCGACTGGTGCGGCCCGTGCAAGTCATTCGCCACCACCTTCGAAGAGGCCGCCCGGAAGCTTGCTGGCCGCTgcc                                                            >  1:457191/1‑136 (MQ=255)                                              gAAGCCAGCTACGCCAGCCAGATCAAAGGCGACCTGCCGCTGCTGGTCGACATCTGGGCCGACTGGTGCGGCCCGTGCAAGTCATTCGCCCCCACCTTCGAAGAGGCCGCCCGGAAGCTTGCTGGCCGCTgccgcc                                                         <  1:319108/136‑1 (MQ=255)                                                            cagccagATCAAAGGCGACCTGCCGCTGCTGGTCGACATCTGGGCCGACTGGTGCGGCCCGTGCAAGTCATTCGCCCCCACCTTCGAAGAGGCCGCCCGGAAGCTTGCTGGCCGCTGCCGCCTGGCCAAGCTCGAc                                           <  1:609535/136‑1 (MQ=255)                                                                cagATCAAAGGCGACCTGCCGCTGCTGGTCGACATCTGGGCCGACTGGTGCGGCCCGTGCAAGTCATTCGCCCCCACCTTCGAAGAGGCCGCCCGGAAGCTTGCTGGCCGCTGCCGCCTGGCCAAGCTCGACAGCg                                       <  1:8122/136‑1 (MQ=255)                                                                cagATCAAAGGCGACCTGCCGCTGCTGGTCGACATCTGGGCCGACTGGTGCGGCCCGTGCAAGTCATTCGCCCCCACATTCGAAGAGGCCGCCCGGAAGCTTGCTGGCCGCTGCCGCCTGGCCAAGCTCGACAGCg                                       <  1:563409/136‑1 (MQ=255)                                                                 agATCAAAGGCGACCTGCCGCTGCTGGTCGACATCTGGGCCGACTGGTGCGGCCCGTGCAAGTCATTCGCCCCCACCTTCGAAGAGGCCGCCCGGAAGCTTGCTGGCCGCTGCCGCCTGGCCAAGCTCGACAGCGa                                      <  2:103462/136‑1 (MQ=255)                                                                  gATCAAAGGCGACCTGCCGCTGCTGGTCGACATCTGGGCCGACTGGTGCGGCCCGTGCAAGTCATTCGCCACCACCTTCGAAGAGGCCGCCCGGAAGCTTGCTGGCCGCTGCCGCCTGGCCAAGCTCGACAGCGaa                                     <  2:669611/136‑1 (MQ=255)                                                                  gATCAAAGGCGACCTGCCGCTGCTGGTCGACATCTGGGCCGACTGGTGCGGCCCGTGCAAGTCATTCGCCACCACCTTCGAAGAGGCCGCCCGGAAGCTTGCTGGCCGCTGCCGCCTGGCCAAGCTCGACAGCGaa                                     <  2:457191/136‑1 (MQ=255)                                                                      aaaGGCGACCTGCCGCTGCTGGTCGACATCTGGGCCGACTGGTGCGGCCCGTGCAAGTCATTCGCCCCCACCTTCGAAGAGGCCGCCCGGAAGCTTGCTGGCCGCTGCCGCCTGGCCAAGCTCGACAGCGAAGCCa                                 >  1:334686/1‑136 (MQ=255)                                                                      aaaGGCGACCTGCCGCTGCTGGTCGACATCTGGGCCGACTGGTGCGGCCCGTGCAAGTCATTCGCCCCCACCTTCGAAGAGGCCGCCCGGAAGCTTGCTGGCCGCTGCCGCCTGGCCAAGCTCGACAGCGAAGCCa                                 >  1:331364/1‑136 (MQ=255)                                                                       aaGGCGACCTGCCGCTGCTGGTCGACATCTGGGCCGACTGGTGCGGCCCGTGCAAGTCATTCGCCCCCACCTTCGAAGAGGCCGCCCGGAAGCTTGCTGGCCGCTGCCGCCTGGCCAAGCTCGACAGCGAAGCCaa                                <  1:272334/136‑1 (MQ=255)                                                                                tGCCGCTGCTGGTCGACATCTGGGCCGACTGGTGCGGCCCGTGCAAGTCATTCGCCCCCACCTTCGAAGAGGCCGCCCGGAAGCTTGCTGGCCGCTGCCGCCTGGCCAAGCTCGACAGCGAAGCCAACCGTAACCt                       >  2:374030/1‑136 (MQ=255)                                                                                tGCCGCTGCTGGTCGACATCTGGGCCGACTGGTGCGGCCCGTGCAAGTCATTCGCCCCCACCTTCGAAGAGGCCGCCCGGAAGCATGCTGGCCGCTGCCGCCTGGCCAAGCTCGACAGCGAAGCCAACCgccgtcg                       >  1:64301/1‑130 (MQ=255)                                                                                       gctgGTCGACATCTGGGCCGACTGGTGCGGCCCGTGCAAGTCATTCGCCCCCACCTTCGAAGAGGCCGCCCGGAAGCTTGCTGGCCGCTGCCGCCTGGCCAAGCTCGACAGCGAAGCCAACCGTAACCTCGCAggg                <  1:398267/136‑1 (MQ=255)                                                                                                     gggCCGACTGGTGCGGCCCGTGCAAGTCATTCGCCCCCACCTTCGAAGAGGCCGCCCGGAAGCTTGCTGGCCGCTGCCGCCTGGCCAAGCTCGACAGCGAAGCCAACCGTAACCTCGCAGGGCAACTGGGCATACg  <  2:512763/136‑1 (MQ=255)                                                                                                                        gTGCAAGTCATTCGCCCCCACCTTCGAAGAGGCCGCCCGGAAGCTTGCTGGCCGCTgccgcc                                                         >  2:642997/1‑62 (MQ=255)                                                                                                                        gTGCAAGTCATTCGCCCCCACCTTCGAAGAGGCCGCCCGGAAGCTTGCTGGCCGCTgccgcc                                                         <  1:642997/62‑1 (MQ=255)                                                                                                                                        |                                                                                                     CGCTGCAAAAAGGATGTGCTGTTGAGCACCCCCTTCGAACTTACCGAAGCCAGCTACGCCAGCCAGATCAAAGGCGACCTGCCGCTGCTGGTCGACATCTGGGCCGACTGGTGCGGCCCGTGCAAGTCATTCGCCCCCACCTTCGAAGAGGCCGCCCGGAAGCTTGCTGGCCGCTGCCGCCTGGCCAAGCTCGACAGCGAAGCCAACCGTAACCTCGCAGGGCAACTGGGCATACG  >  NC_002947/5791683‑5791918

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