Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	3,793,692	0	T	G	34.0%	68.7 / 10.9	52	S531R (AGC→CGC)	PP_3354	acyl‑CoA dehydrogenase
Reads supporting (aligned to +/- strand):  ref base T (22/11);  new base G (0/17);  total (24/28)
Fisher's exact test for biased strand distribution p-value = 2.40e-06
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

CCGCCACCTGCCGGTCGAGGTCGGGCAGCAGGTAGCCGAGGAAGAACCCAACGCTTTCCTGCTTTTCCCGGCGCTGCGCCGCCCCCACTGGCAACGCCTGGATGCAGCGCGCGCTGCGGGCCCAGGCGAAGGCCTGCAGCACCGTTGCGCACAACTGCAGGAAGTCGCCTGCGGCGCGGTACGGAAATTCGGGATCGCCCGCGGTCATCGACAGCACCTTGGCCAGCGCTGCCTCCAGCGTGTCGCAC  >  NC_002947/3793578‑3793825                                                                                                                   |                                                                                                                                      ccgccCCCGGCCGGTCGTGGGCGGGCAGCAGGTCGCCGAGGAACAACCCAACGCTTTCCTGCTTTTCCCGGCGCGGCGCCGCCCCCACTGCCACCGCCTGGATGCAGCGCGCGCGGCGGGCCCAGGCGAAGGCCTg                                                                                                                  <  2:635044/136‑1 (MQ=255) ccgccACCTGCCGGTCGAGGGCGGGCAGCAGGTCGCCGAGGAAGAACCCAACGCTTTCCTGCTTTCCCCGGCGCTGCGCCGCCCCCACTGCCAACGCCTGTATGCAGCGCGCGCTGCGGGCCCAGGCGAAGGCCTg                                                                                                                  <  2:258727/136‑1 (MQ=255)       ccTGCCGGTCGAGGTCGGGCAGCAGGTAGCCGAGGAAGAACCCAACGCTTTCCTGCTTTTCCCGGCGCTGCGCCGCCCCCACTGGCAACGCCTGGATGCAGCGCGCGCTGCGGGCCCAGGCGAAGGCCTGCAGCAc                                                                                                            >  1:777261/1‑136 (MQ=255)       ccTGCCGGTCGAGGTCGGGCAGCAGGTAGCCGAGGAAGAACCCAACGCTTTCCTGCTTTTCCCGGCGCTGCGCCGCCCCCACTGGCAACGCCTGGATGCAGCGCGCGCTGCGGGCCCAGGCGAAGGCCTGCAGCAc                                                                                                            >  2:581208/1‑136 (MQ=255)              ggcGAGGGCGGGCAGCCGGTCGCCGAGTAAGAACCCACCCCTTTCGTGCTTTCCCCGGCGCGGCGCCGCCCCCACTGGCACCCCCTGGATGCCGCGCGCGCGGCGGGCCCAGGCGAAGGCCTGCAGCACCGTTgcg                                                                                                     <  2:863863/134‑1 (MQ=255)                         gcagcaGGTAGCCGCGGAAGAACCCAACGCTTACCTGATTTTCCCGGCGCTGCGCCGCCCCCACTGGCAACGCCTGGATGCAGCGCGCGCTGCGGGCCCAGGCGAAGGCCTGCAGCACCGTTGCGCACAACTGCAg                                                                                          <  2:673339/136‑1 (MQ=255)                         ccagcagGTACCCGAGGAACAACCCAACGCTTTCCTGCTTTACCCGGCGCTGCGCCGCCCCCACTGGCACCGCCTGGATGCAGCGCGCGCGGCGGGCCCAGGCGAAGGCCTGCAGCACCGTTGCGCACAACTGCAg                                                                                          <  1:501959/135‑1 (MQ=255)                           agcagGTAGCCGAGGAAGAACCCAACGCTTTCCTGCTTTTCCCGGCGCTGCGCCGCCCCCACTGGCAACGCCTGGATGCAGCGCGCGCTGCGGGCCCAGGCGAAGGCCTGCAGCACCGTTGCGCACAACTGCAGGa                                                                                        >  1:471980/1‑136 (MQ=255)                                      gAGGAAGAACCCAACGCTTTCCTGCTTTTCCCGGCGCTGCGCCGCCCCCACTGGCACCGCCTGGATGCAGCGCGCGCTGCGGGCCCAGGCGAAGGCCTGCAGCACCGTTGCGCACAACTGCAGGAAGTCGCCTgcg                                                                             <  2:374547/136‑1 (MQ=255)                                          aagaaCCCAACGCTTTCCTGCTTTTCCCGGCGCTGCGCCGCCCCCACTGGCAACGCCTGGATGCAGCGCGCGCTCCGGGCCCAGGCGAAGGCCTGCAGCACCGTTGCGCACAACAGCAGGAAGTCCCCAGCGgcgc                                                                         >  1:946084/1‑136 (MQ=255)                                                 cAACGCTTTCCTGCTTTTCCCGGCGCTGCGCCGCCCCCACTGGCAACCCCAGGATGCAGCGCGCGCTGCCGGCCCAGGCGAACGCCTGCAGCACCGATGCGCACAACTGCAGGAAGACGCCTGCGGCGCGGGACgg                                                                  >  2:648664/1‑136 (MQ=255)                                                       cttCCTGCTTTGCCGGGCGCGCCCCCGCCCCCCCTGGCTCCGCTGGCTTGCCGCGCGCGCTGCGGGCCCAGGCGAAGGCCTGCAGCACCGTTGCGCACAACTGCAGGAAGTCGCCTGCGGCGCGGTACGGAAATTc                                                            <  1:583178/135‑1 (MQ=255)                                                            tGCTATTCCCGGCGCTGCGCCGCCCCCACTGGCAACGCCTGGATGCAGCGCGCGCTGCGGGCCCAGGCGAAGGCCAGCAGCACCGTTGCGCACAACTCCAGGAACTCCCCGGCGCCGCGGTACGGGACATCGGGag                                                       >  1:877694/1‑135 (MQ=255)                                                             gCTTTCCCCGGCGCGGCGCCGCCCCCACTCGCACCGCCTGGATGCAGCGCGCGCGGCGGGCCCAGGCGAAGGCCTGCAGCACCGTTGCGCACAACTGCAGGAAGTCGCCTGCGGCGCGGTACGGAAATTCGGGATc                                                      <  1:485082/136‑1 (MQ=255)                                                               ttttCCCGGCGCTGCGCCGCCCCCCCTGGCAACGCCTGGACGCAGCGCGCGCTGCGGGCCCAGGCGAAGGCCTGCAGCACCGTTGCGCACAACTGCAGGAAGTCGCCTGCGGCGCGGTACGGAAATTCGGGATCGc                                                    <  2:50225/136‑1 (MQ=255)                                                               tttacgcGGGGCGGCGCCGCCGCCACTGGCACCCCCTGGTGGCCGCGCGCGCGGCGGGCCCAGGCGAAGGCCTGCAGCACCGTTGCGCACAACTGCAGGAAGTCGCCTGCGGCGCGGTACGGAAATTCGGTATCGc                                                    <  1:529156/130‑1 (MQ=255)                                                               gtttCCGGGCGCGGCGCCGCCCCCTCTGCCCCGGGCTGGCTGCCGCGCGCGCGGCGGGCCCAGGCGAAGGCCTGCAGCACCGTTGCGCACAACTGCAGGAAGTCGCCTGCGGCGCGGTACGGAAATTCGGGATCGc                                                    <  1:903057/135‑1 (MQ=255)                                                                 ttCCCGGCGCTGCGCCGCCCCCACTGGCAACGCCTGGATGCAGCGCGCGCAGCGGGCCCAGGCGAAGGCCTGCAGCACCGTTGCGCACAACTGCAGGAAGTCGCCTGCGGCGCGGGACGGAAAATAGTGATAGccc                                                  >  1:796070/1‑136 (MQ=255)                                                                    cccgcgcGGCGCCGCCCCCACTGCCACCGCGGGGAGGCAGCGCGCGCGGCGGGCCCAGGCGAAGGCCTGCAGCACCGTTGCGCACAACTGCAGGAAGTCGCCTGCGGCGCGGTACGGAAATTCGGGATCGCCcgcg                                               <  1:825016/133‑1 (MQ=255)                                                                    ccGGCGCTGCGCCGCCCCCACTGGCAAAGCCTGGATGCAGCGCGCGCAGCGGGCCCAGGCGAAGGCCTGCAGCACCGTTGCGCACAAATGCAGGAAGTCGCCTGCGGCGCGGTACGGAAAGTAGGGATCGCCcgcg                                               >  2:608074/1‑136 (MQ=255)                                                                          cTGCGCCGCCCCCACTGGCAACCCCTGGATCCCGCGCGCGCGGCGGGCCCAGGCGAAGGCCTGCAGCACCGTTGCGCACAACTGCAGGAAGTCGCCTGCGGCGCGGTACGGAAATTCGGGATCGCCCGCGGTCATc                                         <  1:709379/136‑1 (MQ=255)                                                                           ggcgcCGCCCCCACTGGCAACGCGTGGATGCAGCGCGCGCGGCGGGCCCAGGCGAAGGCCTGCAGCACCGTTGCGCACAACTGCAGGAAGTCGCCTGCGGCGCGGTACGGAAATTCGGGATCGCCCGCGGTCATCg                                        <  1:753744/135‑1 (MQ=255)                                                                           ggcgcCGCCCCCACTGGCAACGCCTGGATGCCGCGCGCGCGGCGGGCCCAGGCGAAGGCCTGCAGCACCGTTGCGCACAACTGCAGGAAGTCGCCTGCGGCGCGGTACGGAAATTCGGGATCGCCCGCGGTCATCg                                        <  1:375356/135‑1 (MQ=255)                                                                                     ccACTGGCAACGCCTGGATGCAGCGCGCGCTGCGGGCCCAGGCGAAGGCCTGCAGCACCGTTGCGCACAACTGCAGGAAGTCGCCTGCGGCGCGGTACGGAAATTCGGGATCGCCCGCGGTCATCGACAGCACCtt                              >  2:907743/1‑136 (MQ=255)                                                                                     ccACTGGCAACGCCTGGATGCAGCGCGCGCTGCGGGCCCAGGCGAAGGCCTGCAGCACCGTTGCGCACAACTGCAGGAAGTCGCCGGCGGCGCGGGTCGGCACTTGGGGATCGCCCGCGGGCATGGGCAGCAccat                              >  2:897641/1‑134 (MQ=255)                                                                                      cccTGCCACCGCCGGGATGCAGCGCGCGCGGCGGGCCCAGGCGAAGGCCTGCAGCACCGTTGCGCACAACTGCAGGAAGTCGCCTGCGGCGCGGTACGGAAATTCGGGATCGCCCGCGGTCATCGACAGCACCTTg                             <  1:85579/134‑1 (MQ=255)                                                                                          ggCACCGCCTGGATGCAGCGCGCGCTGCGGGCCCAGGCGAAGGCCTGCAGCACCGTTGCGCACAACTGCAGGAAGTCGCCTGCGGCGCGGTACGGAAATTCGGGATCGCCCGCGGTCATCGACAGCACCTTGGCCa                         <  1:122984/136‑1 (MQ=255)                                                                                          ggCAACGCCTGGATGCAGCGCGCGCTGCGGGCCCAGGCGAAGGCCTGCAGCACCGTTGCGCACAACTGCAGGAAGTCGCCTGCGGCGCGGTACGGAAATTCGGGATCGCCCGCGGTCATCGACAGCACCTTGGCCa                         >  2:647691/1‑136 (MQ=255)                                                                                           gCACCGCCTGGATGCAGCGCGCGCGGCGGGCCCAGGCGAAGGCCTGCAGCACCGTTGCGCACAACTGCAGGAAGTCGCCTGCGGCGCGGTACGGAAATTCGGGATCGCCCGCGGTCATCGACAGCACCTTGGCCAg                        <  1:177413/136‑1 (MQ=255)                                                                                           ccaccGCCTGGATGCCGCGCGCGCGGCGGGCCCAGGCGAAGGCCTGCAGCACCGTTGCGCACAACTGCAGGAAGTCGCCTGCGGCGCGGTACGGAAATTCGGGATCGCCCGCGGTCATCGACAGCACCTTGGCCAg                        <  1:640285/132‑1 (MQ=255)                                                                                           ccAACGCCGGGATGCAGCGCGCGCGGCGGGCCCAGGCGAAGGCCTGCAGCACCGTTGCGCACAACTGCAGGAAGTCGCCTGCGGCGCGGTACGGAAATTCGGGATCGCCCGCGGTCATCGACAGCACCTTGGCCAg                        <  2:185562/135‑1 (MQ=255)                                                                                            cAACGCCTGGATGCAGCGCGCGCGGCGGGCCCAGGCGAAGGCCTGCAGCACCGTTGCGCACAACTGCGGGAAGTCGCCTGCGGCGCGGTACGGAAATTCGGGATCGCCCGCGGTCATCGACAGCACCTTGGCCAgc                       <  2:495790/136‑1 (MQ=255)                                                                                             aaCGCCTGGATGCAGCGCGCGCTGCGGGCCCAGGCGAAGGCCTGCAGCACCGTTGCGCACAACTGCAGGAAGTCGCCAGCGGCGCGGTACGGAAATTCGGGATCGCCCGCGGTCATCGACAGCACCATGGCCAgcc                      >  1:398223/1‑135 (MQ=255)                                                                                              ccGCCTGGCGGCCGGGCGCGCGGCGGGCCCAGGCGAAGGCCTGCAGCACCGTTGCGCACAACTGCAGGAAGTCGCCTGCGGCGCGGTACGGAAATTCGGGATCGCCCGCGGTCATCGACAGCACCTTGGCCAgcgc                     <  1:581208/135‑1 (MQ=255)                                                                                                    ggATGCAGCGCGCGCTGCGGGCCCAGGCGAAGGCCTGCAGCCCCGTTGCGCACAACTGCAGGAAGTCGCCTGCGGTGCGGTACGGAAATTCGGGATCGCCCGCGGTCATCGACAGCACCTTGGCCAGCGCTGCCTc               >  2:876698/1‑136 (MQ=255)                                                                                                    ggATGCAGCGCGCGCTGCGGGCCCAGGCGAAGGCCTGCAGCACCGTTGCGCACAACTGCAGGAAGTCGCCTGCGGTGCGGTACGGAAATTCGGGATCGCCCGCGGTCATCGACAGCACCTTGGCCAGCGCTGCCTc               >  2:881738/1‑136 (MQ=255)                                                                                                       tGCAGCGCGCGCTGCGGGCCCAGGCGAAGGCCTGCAGCACCGTTGCGCACAACTGCAGGAAGTCGCCTGCGGCGCGGTACGGAAATTCGGGATCGCCCGCGGTCATCGACAGCACCTTGGCCAGCGCTGCCTCCAg            <  2:662314/136‑1 (MQ=255)                                                                                                       cgCAGCGCGCGCTGCGGGCCCAGGCGAAGGCCTGCAGCACCGTTGCGCACAACTGCAGGAAGTCGCCTGCGGCGCGGTACGGAAATTCGGGATCGCCCGCGGTCATCGACAGCACCTTGGCCAGCGCTGCCTCCAg            <  2:753853/135‑1 (MQ=255)                                                                                                        gCAGCGCGCGCTGCGGGCCCAGGCGAAGGCCTGCAGCACCGTTGCGCACAACTGCAGGAAGTCGCCAGCGGCGCGGGACGGAAATTCGGGCTCGCCCGCGGGCATCGACCGCCCCGTGGCCCGCGCGgccgccccc           >  1:581552/1‑133 (MQ=255)                                                                                                            ggcgcggGGCGGGCCCAGGCGAAGGCCTGCAGCACCGTTGCGCACAACTGCAGGAAGTCGCCTGCGGCGCGGTACGGAAATTCGGGATCGCCCGCGGTCATCGACAGCACCTTGGCCAGCGCTGCCTCCAGCgtgt       <  2:581552/135‑1 (MQ=255)                                                                                                            cgcgcgcTGCGGGCCCAGGCGAAGGCCTGCAGCACCGTTGCGCACAACTGCAGGAAGTCGCCTGCGGCGCGGTACGGAAATTCGGGATCGCCCGCGGTCATCGACAGCACCTTGGCCAGCGCTGCCTCCAGCgtgt       <  2:540431/136‑1 (MQ=255)                                                                                                            cgcgcgcTGCGGGCCCAGGCGAAGGCCTGCAGCACCGTTGCGCACAACTGCAGGAAGTCGCCTGCGGCGCGGTACGGAAATTCGGGATCGCCCGCGGTCATCGACAGCACCTTGGCCAGCGCTGCCTCCAGCgtgt       >  2:242293/1‑136 (MQ=255)                                                                                                            cgcgcgcTGCGGGCCCAGGCGAAGGCCTGCAGCACCGTTGCGCACAACTGCAGGAAGTCGCCTGCGGCGCGGTACGGAAATTCGGGATCGCCCGCGGTCATCGACAGCACCTTGGCCAGCGCTGCCTCCAGCgtgt       <  2:777261/136‑1 (MQ=255)                                                                                                             gcgcgcTGCGGGCCCAGGCGAAGGCCTGCAGCACCGTTGCGCACAACTGCAGGAAGTCGCCTGCGGCGCGGTACGGAAATTCGGGATCGCCCGCGGTGATCGACAGCACCTTGGCCAGCGCTGCCTCCAGCGTGTc      >  2:697045/1‑136 (MQ=255)                                                                                                             gcgcgcTGCGGGCCCAGGCGAAGGCCTGCAGCACCGTTGCGCACAACTGCAGGAAGTCGCCTGCGGCGCGGTACGGAAATTCGGGATCGCCCGCGGTCATCGGCAGCACCTTGGCCCGCGCTGCCTCCAGCgtgtg      >  1:154256/1‑135 (MQ=255)                                                                                                             gcgcgcTGCGGGCCCAGGCGAAGGCCTGCAGCACCGTTGCGCACAACTGCAGGAAGTCGCCTGCGGCGCGGTACGGAAATTCGGGATCGCCCGCGGTCATCGACCGCCCCTTGGCCCGCGCTGCCTCCCGCGTGTc      >  1:107959/1‑136 (MQ=255)                                                                                                             gcgcgcTGCGGGCCCAGGCGAAGGCCTGCAGCACCGTTGCGCACAACTGCAGGAAGTCGCCTGCGGCGCGGTACGGAAATTCGGGATCGCCCGCGGTCATCGACAGCACCTTGGCCAGCGCTGCCTCCAGCGTGTc      >  1:247212/1‑136 (MQ=255)                                                                                                             gcgcgcTGCGGGCCCAGGCGAAGGCCTGCAGCACCGTTGCGCACAACTGCAGGAAGTCGCCTGCGGCGCGGTACGGAAATTCGGGATCGCCCGCGGTCATCGACAGCACCTTGGCCAGCGCTGCCTCCAGCGTGTc      >  1:282625/1‑136 (MQ=255)                                                                                                             gcgcgcTGCGGGCCCAGGCGAAGGCCTGCAGCACCGTTGCGCACAACTGCAGGAAGTCGCCTGCGGCGCGGGACGGAAATTCGGGCTCGCCCGCGGTCATCGACAGCACCTTGGCCCGCGCTGCCTCCAGCGTGTc      >  2:571692/1‑136 (MQ=255)                                                                                                               gcgcTGCGGGCCCAGGCGAAGGCCTGCAGCACCGTTGCGCACAACTGCAGGAAGTCGCCTGCGGCGCGGGACGGAAATTCGGGATCGCCCGCGGTCATCGACCGCACCTTGGCCCGCGCTGCCTCCCGCGTGTcgc    >  1:196081/1‑136 (MQ=255)                                                                                                                 gcTGCGGGCCCAGGCGAAGGCCTGCAGCACCGTTGCGCACAACTGCAGGAAGTCGCCTGCGGCGCGGTACGGAAATTCGGGATCGCCCGCGGTCATCGACAGCACCTTGGCCAGCGCTGCCTCCAGCGTGTcgcac  >  1:516473/1‑136 (MQ=255)                                                                                                                 gcTGCGGGCCCAGGCGAAGGCCTGCAGCACCGTTGCGCACAACTGCAGGAAGTCGCCTGCGGCGCGGTACGGAAATTCGGGATCGCCCGCGGTCATCGACAGCACCTTGGCCAGCGCTGCCTCCAGCGTGTcgcac  <  1:907743/136‑1 (MQ=255)                                                                                                                   |                                                                                                                                      CCGCCACCTGCCGGTCGAGGTCGGGCAGCAGGTAGCCGAGGAAGAACCCAACGCTTTCCTGCTTTTCCCGGCGCTGCGCCGCCCCCACTGGCAACGCCTGGATGCAGCGCGCGCTGCGGGCCCAGGCGAAGGCCTGCAGCACCGTTGCGCACAACTGCAGGAAGTCGCCTGCGGCGCGGTACGGAAATTCGGGATCGCCCGCGGTCATCGACAGCACCTTGGCCAGCGCTGCCTCCAGCGTGTCGCAC  >  NC_002947/3793578‑3793825

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