Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	4,093,303	C→T	12.1%	H178H (CAC→CAT)	PP_3598 →	peptidase C26

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	4,093,303	0	C	T	12.1%	89.7 / 4.0	33	H178H (CAC→CAT)	PP_3598	peptidase C26
Reads supporting (aligned to +/- strand):  ref base C (12/17);  new base T (2/2);  total (14/19)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.33e-01

TGGGTGGCGACATTTACCAGAAGGTGTATGCCGAACCTGGGTTCAATGACCACCGGGAAAACCCGGAAGACCCGGTCGAAGTGCAGTACGCCCAGGTCCATGGGGTGAAGATCCAGCCGGGCAGCTGGCTGCACGACACCCTGGGTACCGATGAGATTCGCGTCAACTCGCTGCATGGCCAAGGCCTGCACAAGCTTGGCGCCGGTATTGAAGCGATCGCCCGTGCCGAAGATGGCCTGGTCGAAGCGATTCACGCGCCAAGCATTTCG  >  NC_002947/4093170‑4093438                                                                                                                                      |                                                                                                                                        tGGGTGGCGACATTTACCAGAAGGTGTATGCCGAACCTGGGTTCAATGACCACCGGGAAAACCCGGAAGACCCGGTCGAAGTGCAGTACGCCCAGGTCCATGGGGTGAAGATCCAGCCGGGCAGCTGGCTGCACGa                                                                                                                                       >  1:507256/1‑136 (MQ=255)  gggTGGCGACATTTACCAGAAGGTGTATGCCGAACCTGGGTTCAATGACCACCGGGAAAACCCGGAAGACCCGGTCGAAGTGCAGTACGCCCAGGTCCATGGGGTGAAGATCCAGCCGGGCAGCTGGCTGCACGac                                                                                                                                      <  2:40117/136‑1 (MQ=255)         taCATTTACCAGAAGGTGTATGCCGAACCTGGGTTCAATGACCACCGGGAAAACCCGGAAGACCCGGTCGAAGTGCAGTACGCCCAGGTCCATGGGGTGAAGATCCAGCCGGGCAGCTGGCTGCACGACACCCTgg                                                                                                                               <  1:368329/135‑1 (MQ=255)           cATTTACCAGAAGGTGTATGCCGAACCTGGGTTCAATGACCACCGGGAAAACCCGGAAGACCCGGTCGAAGTGCAGTACGCCCAGGTCCATGGGGTGAAGATCCAGCCGGGCAGCTGGCTGCATGACACCCTGGGt                                                                                                                             >  1:478579/1‑136 (MQ=255)           cATTTACCAGAAGGTGTATGCCGAACCTGGGTTCAATGACCACCGGGAAAACCCGGAAGACCCGGTCGAAGTGCAGTACGCCCAGGTCCATGGGGTGAAGATCCAGCCGGGCAGCTGGCTGCATGACACCCTGGGt                                                                                                                             >  1:481975/1‑136 (MQ=255)           cATTTACCAGAAGGTGTATGCCGAACCTGGGTTCAATGACCACCGGGAAAACCCGGAAGACCCGGTCGAAGTGCAGTACGCCCAGGTCCATGGGGTGAAGATCCAGCCGGGCAGCTGGCTGCACGACACCCTGGGt                                                                                                                             >  2:617520/1‑136 (MQ=255)                aCCAGAAGGTGTATGCCGAACCTGGGTTCAATGACCACCGGGAAAACCCGGAAGACCCGGTCGAAGTGCAGTACGCCCAGGTCCATGGGGTGAAGATCCAGCCGGGCAGCTGGCTGCACGACACCCTGGGTACCGa                                                                                                                        >  2:649086/1‑136 (MQ=255)                        gtgtATGCCGAACCTGGGTTCAATGACCACCGGGAAAACCCGGAAGACCCGGTCGAAGTGCAGTACGCCCAGGTCCATGGGGTGAAGATCCAGCCGGGCAGCTGGCTGCACGACACCCTGGGTACCGATGAGATTc                                                                                                                >  1:649137/1‑136 (MQ=255)                        gtgtATGCCGAACCTGCGTTCAATGACCACCGGGAAAACACGGAAGACCCGGTCGAAGTGCAGTACGCCCAGGTCCATGCGGTGAAGATCAAGCCGGGCAGCTGGCTGCACGACACCCTGGGTACCGATGAGATTc                                                                                                                >  1:138478/1‑136 (MQ=255)                              cccGAACCTGGGTTCAACGACCACCGGGAAAACCCGGAAGACCCGGTCGAAGTGCAGTACGCCCAGGTCCATGGGGTGAAGATCCAGCCGGGCAGCTGGCTGCACGACACCCTGGGTACCGATGAGATTCGCGTCa                                                                                                          <  2:308927/135‑1 (MQ=255)                                cGAACCTGGGTTCAATGACCACCGGGAAAACCCGGAAGACCCGGTCGAAGTGCAGTACGCCCAGGTCCATGGGGTGAAGATCCAGCCGGGCAGCTGGCTGCACGACACCCTGGGTACCGATGAGATTCGCGTCAAc                                                                                                        <  1:246789/136‑1 (MQ=255)                                  aaCCTGGGTTCAATGACCACCGGGAAAACCCGGAAGACCCGGTCGAAGTGCAGTACGCCCAGGTCCATGGGGTGAAGATCCAGCCGGGCAGCTGGCTGCACGACACCCTGGGTACCGATGAGATTCGCGTCAACTc                                                                                                      >  1:87549/1‑136 (MQ=255)                                     cTGGGTTCAATGACCACCGGGAAAACCCGGAAGACCCGGTCGAAGTGCAGTACGCCCAGGTCCATGGGGTGAAGATCCAGCCGGGCAGCTGGCTGCACGACACCCTGGGTACCGATGAGATTCGCGTCAACTCGCt                                                                                                   <  1:289210/136‑1 (MQ=255)                                     cTGGGTTCAATGACCACCGGGAAAACCCGGAAGACCCGGTCGAAGTGCAGTACGCCCAGGTCCATGGGGTGAAGATCCAGCCGGGCAGCTGGCTGCACGACACCCTGGGTACCGATGAGATTCGCGTCAACTCGCt                                                                                                   <  1:310910/136‑1 (MQ=255)                                      tGGGTTCAATGACCACCGGGAAAACCCGGAAGACCCGGTCGAAGTGCAGTACGCCCAGGTCCATGGGGTGAAGATCCAGCCGGGCAGCTGGCTGCACGACACCCTGGGTACCGATGAGATTCGCGTCAACTCGCTg                                                                                                  >  1:589063/1‑136 (MQ=255)                                               tGACCACCGGGAAAACCCGGAAGACCCGGTCGAAGTGCAGTACGCCCAGGTCCATGGGGTGAAGATCCAGCCGGGCAGCTGGCTGCACGACACCCTGGGTACCGATGAGATTCGCGTCAACTCGCTGCATGGCCaa                                                                                         <  1:617520/136‑1 (MQ=255)                                                     ccGGGAAAACCCGGAAGACCCGGTCGAAGTGCAGTACGCCCAGGTCCATGGGGTGAAGATCCAGCCGGGCAGCTGGCTGCACGACACCCTGGGTACCGATGAGATTCGCGGCAACTCGCTGCATGGCCAAGGCCTg                                                                                   <  1:19917/136‑1 (MQ=255)                                                       gggAAAACCCGGAAGACCCGGTCGAAGTGCAGTACGCCCAGGTCCATGGGGTGAAGATCCAGCCGGGCAGCTGGCTGCACGACACCCTGGGTACCGATGAGATTCGCGTCAACTCGCTGCATGGCCAAGGCCTGca                                                                                 >  1:395953/1‑136 (MQ=255)                                                                cGGAAGACCCGGTCGAAGTGCAGTACGCCCAGGTCCATGGGGTGAAGATCCAGCCGGGCAGCTGGCTGCATGACACCCTGGGTACCGATGAGATTCGCGTCAACTCGCTGCATGGCCAAGGCCTGCACAAGCTTgg                                                                        <  2:481975/136‑1 (MQ=255)                                                                cGGAAGACCCGGTCGAAGTGCAGTACGCCCAGGTCCATGGGGTGAAGATCCAGCCGGGCAGCTGGCTGCATGACACCCTGGGTACCGATGAGATTCGCGTCAACTCGCTGCATGGCCAAGGCCTGCACAAGCTTgg                                                                        <  2:478579/136‑1 (MQ=255)                                                                cGGAAGACCCGGTCGAAGTGCAGTACGCCCAGGTCCATGGGGTGAAGATCCAGCCGGGCAGCTGGCTGCACGACACCCTGGGTACCGATGAGATTCGCGTCAACTCGCTGCATGGCCAAGGCCTGCACAAGCTTgg                                                                        <  1:599623/136‑1 (MQ=255)                                                                cGGAAGACCCGGTCGAAGTGCAGTACGCCCAGGTCCATGGGGTGAAGATCCAGCCGGGCAGCTGGCTGCACGACACCCTGGGTACCGATGAGATTCGCGTCAACTCGCTGCATGGCCAAGGCCTGCACAAGCTTgg                                                                        <  1:193015/136‑1 (MQ=255)                                                                           gTCGAAGTGCAGTACGCCCAGGTCCATGGGGTGAAGATCCAGCCGGGCAGCTGGCTGCACGACACCCTGGGTACCGATGAGATTCGCGTCAACTCGCTGCATGGCCAAGGCCTGCACAAGCTTGGCGCCGGTATTg                                                             >  2:455198/1‑136 (MQ=255)                                                                                     aGTACGCCCAGGTCCATGGGGTGAAGTTCCAGCCGGGCAGCTGGCTGCACGACACCCTGGGTACCGATGAGATTCGCGTCAACTCGCTGCATGGCCAAGGCCTGCACAAGCTTgg                                                                        <  2:261373/115‑1 (MQ=255)                                                                                     aGTACGCCCAGGTCCATGGGGTGAAGATCCAGCCGGGCAGCTGGCTGCACGACACCCTGGGTACCGATGAGATTCGCGTCAACTCGCTGCATGGCCAAGGCCTGCACAAGCTTgg                                                                        >  1:261373/1‑115 (MQ=255)                                                                                       tACGCCCAGTTCCTTGGGGTGAAGATCCAGCGGGGCAGCTGGCTGCACGACACCCTGGGTACCGATGAGATTCGCGTCAACTCGCTGCATGGCCAAGGCCTGCAGAAGCTTGGCGCCGGTATTGAAGCGATCGccc                                                 <  2:225704/136‑1 (MQ=255)                                                                                       tACGCCCAGGTCCATGGGGTGAAGATCCAGCCGGGCAGCTGGCTGCACGACACCCTGGGTACCGATGAGATTCGCGTCAACTCGCTGCATGGCCAAGGCCTGCACAAGCTTGGCGCCGGTATTGAAGCGATCGccc                                                 >  1:627417/1‑136 (MQ=255)                                                                                               ggTCCATGGGGTGAAGATCCAGCCGGGCAGCTGGCTGCACGACACCCTGGGTACCGATGAGATTCGCGTCAACTCGCTGCATGGCCAAGGCCTGCACAAGCTTGGCGCCGGTATTGAAGCGATCGCCCGTGCCGaa                                         >  1:654910/1‑136 (MQ=255)                                                                                                  ccATGGGGTGAAGATCCAGCGGGGCAGCTGGCTGCACGACACCCTGGGTACCGATGAGATTCGCGTCAACTCGCTGCATGGCGAAGGCCTGCACAAGCTTGGCGCCGGTATTGAAGCGATCGCCCGTGCCGAAGAt                                      <  1:480990/136‑1 (MQ=255)                                                                                                         gTGAAGAGCCAGTCGGGCAGCTGGCTGCACGACACCCTGGGTACCGATGAGATTCGCGTCAACTCGCTGCATGGCCAAGGCCTGCACAAGCTTGGCGCCGGTATTGAAGCGATCGCCCGTGCCGAAGATGGCCTgg                               <  1:50648/136‑1 (MQ=255)                                                                                                                   aGCCGGGCAGCTGGCTGCACGACACCCTGGGTACCGATGAGATTCGCGTCAACTCGCTGCATGGCCAAGGCCTGCACAAGCTTGGCGCCGGTATTGAAGCGATCGCCCGTGCCGAAGATGGCCTGGTCGAAGCGAt                     <  1:649086/136‑1 (MQ=255)                                                                                                                                 ctgCACGACACCCTGGGTACCGATGAGATTCGCGTCAACTCGCTGCATGGCCAAGGCCTGCACAAGCTTGGCGCCGGTATTGAAGCGATCGCCCGTGCCGAAGATGGCCTGGTCGAAGCGATTCACGCGCCAAGCa       <  2:477851/136‑1 (MQ=255)                                                                                                                                    cACGACACCCTGGGTACCGATGAGATTCGCGTCAACTCGCTGCATGGCCAAGGCCTGCACAAGCTTGGCGCCGGTATTGAAGCGATCGCCCGTGCCGAAGATGGCCTGGTCGAAGCGATTCACGCGCCAAGCAttt    <  2:589063/136‑1 (MQ=255)                                                                                                                                      cGACACCCTGGGTACCGATGAGATTCGCGTCAACTCGCTGCATGGCCAAGGCCTGCACAAGCTTGGCGCCGGTATTGAAGCGATCGCCCGTGCCGAAGATGGCCTGGTCGAAGCGATTCACGCGCCAAGCATTTCg  <  1:165580/136‑1 (MQ=255)                                                                                                                                      |                                                                                                                                        TGGGTGGCGACATTTACCAGAAGGTGTATGCCGAACCTGGGTTCAATGACCACCGGGAAAACCCGGAAGACCCGGTCGAAGTGCAGTACGCCCAGGTCCATGGGGTGAAGATCCAGCCGGGCAGCTGGCTGCACGACACCCTGGGTACCGATGAGATTCGCGTCAACTCGCTGCATGGCCAAGGCCTGCACAAGCTTGGCGCCGGTATTGAAGCGATCGCCCGTGCCGAAGATGGCCTGGTCGAAGCGATTCACGCGCCAAGCATTTCG  >  NC_002947/4093170‑4093438

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: ‑