Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,508,534	T→C	11.1%	F163L (TTC→CTC)	oprC →	copper receptor OprC

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	5,508,534	0	T	C	11.1%	88.4 / 4.6	36	F163L (TTC→CTC)	oprC	copper receptor OprC
Reads supporting (aligned to +/- strand):  ref base T (17/15);  new base C (2/2);  total (19/17)
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.99e-01

GGAGCCTGCCCCAATCGCATGGACGCACCCACGTCCTATATCTCGCCAGAAACCTACGACCGCCTGACCGTGATCAAAGGCCCGCAGAGCGTCATCTGGGGGCCCGGCGGCTCGGCCGGCACCATCCTCTTCGAGCGCGAACCTGAGAAATTCGGCACCCTCGGCAGCCGCGTCAACGCAAGCCTGCTGGCCGGCTCCAACGGCCGCTTCGACAAGGTGCTGGATGCCGCGGCCGGCAACAGTCA  >  NC_002947/5508405‑5508649                                                                                                                                  |                                                                                                                    ggAGCCTGCCCCAATCGCATGGACGCACCCACGTCCTATATCTCGCCAGAAACCTACGACCGCCTGACCGTGATCAAAGGCCCGCAGAGCGTCATCTGGGGGCCCGGCGGCTCGGCCGGCACCATCCTCTTCGAgc                                                                                                               <  2:282996/136‑1 (MQ=255)   aGCCTGCCCCAATCGCATGGACGCACCCACGTCCTATATCTCGCCAGAAACCTACGACCGCCTGACCGTGATCAAAGGCCCGCAGAGCGTCATCTGGGGGCCCGGCGGCTCGGCCGGCACCATCCTCTTCGAgcgc                                                                                                             >  2:379302/1‑136 (MQ=255)    gCCTGCGCCAATCGCTTGGACGCACCCACGTCCTATATCTCGCCAGAAACCTACGAGCGCCTGACCGTGATCAAAGGCCCGCAGAGCGTCATCTGGGGGCCCGGCGGCTCGGCCGGCACCATCCTCTTCGAGcgcg                                                                                                            <  1:126417/136‑1 (MQ=255)    gCCTGCCCCAATCGCATGGACGCACCCACGTCCTATATCTCGCCAGAAACCTACGACCGCCTGACCGTGATCAAAGGCCCGCAGAGCGTCATCTGGGGGCCCGGCGGCTCGGCCGGCACCATCCTCTTCGAGcgcg                                                                                                            <  1:214021/136‑1 (MQ=255)          cccAATCGCATGGACGCACCCACGTCCTATATCTCGCCAGAAACCTACGACCGCCTGACCGTGATCAAAGGCCCGCAGAGCGTCATCTGGGGGCCCGGCGGCTCGGCCGGCACCATCCTCCTCGAGCGCGAACCTg                                                                                                      >  2:193931/1‑136 (MQ=255)          cccAATCGCATGGACGCACCCACGTCCTATATCTCGCCAGAAACCTACGACCGCCTGACCGTGATCAAAGGCCCGCAGAGCGTCATCTGGGGGCCCGGCGGCTCGGCCGGCACCATCCTCCTCGAGCGCGAACCTg                                                                                                      >  2:196411/1‑136 (MQ=255)           ccAATCGCATGGACGCACCCACGTCCTATATCTCGCCAGAAACCTACGACCGCCTGACCGTGATCAAAGGCCCGCAGAGCGTCATCTGGGGGCCAGGCGGCTCGGCCGGCACCCTCCTCTTCGAGCGCGAACCTga                                                                                                     >  1:158458/1‑136 (MQ=255)            cAATCGCATGGACGCACCCACGTCCTATATCTCGCCAGAAACCTACGACCGCCTGACCGTGTTCAAAGGCCCGCAGAGCGTCATCTGGGGGCCCGGCGGCTCGGCCGGCACCATCCTCTTCGAGCGCGAACCTgag                                                                                                    >  1:134550/1‑136 (MQ=255)            cAATCGCATGGACGCACCCACGTCCTATATCTCGCCAGAAACCTACGACCGCCTGACCGTGATCAAAGGCCCGCAGAGCGTCATCTGGGGGCCCGGCGGCTCGGCCGGCACCATCCTCTTCGAGCGCGAACCTgag                                                                                                    >  1:456696/1‑136 (MQ=255)                   aTGGACGCACCCACGTCCTATATCTCGCCAGAAACCTACGACCGCCTGACCGTGATCAAAGGCCCGCAGAGCGTCATCTGGGGGCCCGGCGGCTCGGCCGGCACCATCCTCTTCGAGCGCGAACCTGAGAAATTCg                                                                                             <  1:379302/136‑1 (MQ=255)                              cACGTCCTATATCTCGCCAGAAACCTACGACCGCCTGACCGTGATCAAAGGCCCGCAGAGCGTCATCTGGGGGCCCGGCGGCTCGGCCGGCACCATCCTCTTCGAGCGCGAACCTGAGAAATTCGGCACCCTCGGc                                                                                  <  2:136297/136‑1 (MQ=255)                                    cTATATCTCGCCAGAAACCTACGACCGCCTGACCGTGATCAAAGGCCCGCAGAGCGTCATCTGGGGGCCCGGCGGCTCGGCCGGCACCATCCTCTTCGAGCGCGAACCTGAGAAATTCGGCACCCTCGGCAGCcgc                                                                            >  1:262656/1‑136 (MQ=255)                                        atCTCGCCAGAAACCTACGACCGCCTGACCGTGATCAAAGGCCCGCAGAGCGTCATCTGGGGGCCCGGCGGCTCGGCCGGCACCATCCTCTTCGGGCGCGAACCTGAGAAATTCGGAACCCTCGGCAGCCGCGTCa                                                                        >  2:290251/1‑136 (MQ=255)                                              ccAGAAACCTACGACCGCCTGACCGTGATCAAAGGCCCGCAGAGCGTCATCTGGGGGCCCGGCGGCTCGGCCGGCACCATCCTCTTCGAGCGCGAACCTGAGAAATTCGGCACCCTCGGCAGCCGCGTCAACGCaa                                                                  >  2:339307/1‑136 (MQ=255)                                              ccAGAAACCTACGACCGCCTGACCGTGATCAAAGGCCCGCAGAGCGTCATCTGGGGGCCCGGCGGCTCGGCCGGCACCATCCTCTTCGAGCGCGAACCTGAGAAATTCGGCACCCTCGGCAGCCGCGTCAACGCaa                                                                  >  2:505598/1‑136 (MQ=255)                                               cAGAAACCTACGACCGCCTGACCGTGATCAAAGGCCCGCAGAGCGTCATCTGGGGGCCCGGCGGCTCGGCCGGCACCATCCTCTTCGAGCGCGAACCTGAGAAATTCg                                                                                             <  1:183289/108‑1 (MQ=255)                                               cAGAAACCTACGACCGCCTGACCGTGATCAAAGGCCCGCAGAGCGTCATCTGGGGGCCCGGCGGCTCGGCCGGCACCATCCTCTTCGAGCGCGAACCTGAGAAATTCg                                                                                             >  2:183289/1‑108 (MQ=255)                                               cAGAAACCTACGACCGCCTGACCGTGATCAAAGGCCCGCAGAGCGTCATCTGGGGGCCCGGCGGCTCGGCCGGCACCATCCTCTTCGAGCGCGAACCTGAGAAATTCGGCTCCCTCGGCAGCCGCGTCAACGCAAg                                                                 >  2:140729/1‑136 (MQ=255)                                               cAGAAACCTACGACCGCCTGACCGTGATCAAAGGCCCGCAGAGCGTCATCTGGGGGCCCGGCGGCTCGGCCGGCACCATCCTCTTCGAGCGCGAACCTGAGAAATTCGGCACCCTCGGCAGCCGCGTCAACGCAAg                                                                 >  1:74232/1‑136 (MQ=255)                                                aGAAACCTACGACCGCCTGACCGTGATCAAAGGCCCGCAGAGCGTCATCTGGGGGCCCGGCGGCTCGGCCGGCACCATCCTCTTCGAGCGCGAACCTGAGAAATTCGGCACCCTCGGCAGCCGCGTCAACGCAAGc                                                                >  1:64765/1‑136 (MQ=255)                                                aGAAACCTACGACCGCCTGACCGTGATCAAAGGCCCGCAGAGCGTCATCTGGGGGCCCGGCGGCTCGGCCGGCACCATCCTCTTCGAGCGCGAACCTGAGAAATTCGGCACCCTCGGCAGCCGCGTCAACGCAAGc                                                                >  2:114199/1‑136 (MQ=255)                                                aGAAACCTACGACCGCCTGACCGTGATCAAAGGCCCGCAGAGCGTAATCTGGGGGCCCGGCGGCTCGGCCGGCACCATCCTCTTCGAGCGCGAACCTGAGAAATTCGGCACCCTCGGCAGCCGCGTCa                                                                        <  1:56542/128‑1 (MQ=255)                                                aGAAACCTACGACCGCCTGACCGTGATCAAAGGCCCGCAGAGCGTAATCTGGGGGCCCGGCGGCTCGGCCGGCACCATCCTCTTCGAGCGCGAACCTGAGAAATTCGGCACCCTCGGCAGCCGCGTCa                                                                        >  2:56542/1‑128 (MQ=255)                                                 gAAACCTACGACCGCCTGACCGTGATCAAAGGCCCGCAGAGCGTCATCTGGGGGCCCGGCGGCTCGGCCGGCACCATCCTCTTCGAGCGCGAACCTGAGAAATTCGGCACCCTCGGCAGCCGCGTCAACGCAAGcc                                                               >  2:39805/1‑136 (MQ=255)                                                    aCCTACGACCGCCTGACCGTGATCAAAGGCCCGCAGAGCGTCATCTGGGGGCCCGGCGGCTCGGCCGGCACCATCCTCTTCGAGCGCGAACCTGAGAAATTCGGCACCCTCGGCAGCCGCGTCAACGCAAGCctgc                                                            <  1:436109/136‑1 (MQ=255)                                                                    ccGTGATCAAAGGCCCGCAGAGCGTCATCTGGGGGCCCGGCGGCTCGGCCGGCACCATCCTCTTCGAGCGCGAACCTGAGAAATTCGGCACCCTCGGCAGCCGCGTCAACGCAAGCCTGCTGGCCGGCTCCAACgg                                            <  2:262656/136‑1 (MQ=255)                                                                            aaaGGCCCGCAGAGCGTCATCTGGGGGCCCGGCGGCTCGGCCGGCACCATCCTCTTCGAGCGCGAACCTGAGAAATTCGGCACCCTCGGCAGCCGCGTCAACGCAAGCCTGCTGGCCGGCTCCAACGGCCGCTTCg                                    <  1:39805/136‑1 (MQ=255)                                                                            aaaGGCCCGCAGAGCGTCATCTGGGGGCCCGGCGGCTCGGCCGGCACCATCCTCTTCGAGCGCGAACCTGAGAAATTCGGCACCCTCGGCAGCCGCGTCAACGCAAGCCTACTGGCCGGCTCCAACGGCCGCTTCg                                    <  2:456696/136‑1 (MQ=255)                                                                            aaaGGCCCGCAGAGCGTCATCTGGGGGCCCGGCGGCTCGGCCGGCACCATCCTCCTCGAGCGCGAACCTGAGAAATTCGGCACCCTCGGCAGCCGCGTCAACGCAAGCCTGCTGGCCGGCTCCAACGGCCGCTTCg                                    <  1:196411/136‑1 (MQ=255)                                                                            aaaGGCCCGCAGAGCGTCATCTGGGGGCCCGGCGGCTCGGCCGGCACCATCCTCCTCGAGCGCGAACCTGAGAAATTCGGCACCCTCGGCAGCCGCGTCAACGCAAGCCTGCTGGCCGGCTCCAACGGCCGCTTCg                                    <  1:193931/136‑1 (MQ=255)                                                                                             cATCTGGGGGCCCGGCGGCTCGGCCGGCACCATCCTCTTCGAGCGCGAACCTGAGAAATTCGGCACCCTCGGCAGCCGCGTCAACGCAAGCCTGCTGGCCGGCTCCAACGGCCGCTTCGCCAAGGTGCTGGATGcc                   >  2:127967/1‑136 (MQ=255)                                                                                                    gggCCCGGCGGCTCGGCCGGCACCATCCTCTTCGAGCGCGAACCTGAGAAATTCGGCACCCTCGGCAGCCGCGTCAACGCAAGCCTGCTGGCCGGCTCCAACGGCCGCTTCGACAAGGTGCTGGATGCCGcggccg            <  2:470417/136‑1 (MQ=255)                                                                                                      gCCCGGCGGCTCGGCCGGCACCATCCTCTTCGAGCGCGAACCTGAGAAATTCGGCACCCTCGGCAGCCGCGTCAACGCAAGCCTGCTGGCCGGCTCCAACGGCCGCTTCGACAAGGTGCTGGATGCCGcggccggc          <  1:505598/136‑1 (MQ=255)                                                                                                      gCCCGGCGGCTCGGCCGGCACCATCCTCTTCGAGCGCGAACCTGAGAAATTCGGCACCCTCGGCAGCCGCGTCAACGCAAGCCTGCTGGCCGGCTCCAACGGCCGCTTCGACAAGGTGCTGGATGCCGcggccggc          <  1:339307/136‑1 (MQ=255)                                                                                                          ggcggcTCGGCCGGCACCATCCTCTTCGAGCGCGAACCTGAGAAATTCGGCACCCTCGGCAGCCGCGTCAACGCAAGCCTGCTGGCCGGCTCCAACGGCCGCTTAGACAAGGTGCTGGATGCCGCGGCCGGCAACa      <  2:295277/136‑1 (MQ=255)                                                                                                              gcTCGGCCGGCACCATCCTCTTCGAGCGCGAACCTGAGAAATTCGGCACCCTCGGCAGCCGCGTCAAAGCAAGCCTGCTGGCCGGCTCCAACGGCCGCTTCGACAAGGTGCTGGATGCCGCGGCCGGCAACAGTCa  >  2:59932/1‑136 (MQ=255)                                                                                                                                  |                                                                                                                    GGAGCCTGCCCCAATCGCATGGACGCACCCACGTCCTATATCTCGCCAGAAACCTACGACCGCCTGACCGTGATCAAAGGCCCGCAGAGCGTCATCTGGGGGCCCGGCGGCTCGGCCGGCACCATCCTCTTCGAGCGCGAACCTGAGAAATTCGGCACCCTCGGCAGCCGCGTCAACGCAAGCCTGCTGGCCGGCTCCAACGGCCGCTTCGACAAGGTGCTGGATGCCGCGGCCGGCAACAGTCA  >  NC_002947/5508405‑5508649

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