Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	1,751,350	G→T	13.4%	P36Q (CCG→CAG)	PP_1544 ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	1,751,350	0	G	T	13.4%	107.6 / 2.7	35	P36Q (CCG→CAG)	PP_1544	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (17/13);  new base T (2/2);  total (20/15)
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.56e-01

CAAGGCGACGCGCCCGCTCTTCAGGATCCTGCGCGATCCATGACAGCGCCAGGTCCTGCCAGAGTTCCTGCACCTGGTCGTAGCCGTGCCGCTCGATCTCCGCAGCCATGGCAGCCTTGATCCCTTCCGGCATGTTGAAGAAAACCTTCTCGATCCCCAGTTTCGCGGCTGCCTTCTTCTGGCGGTCTCGGTAGTCCGCTGAGTGCTTCGCGGCTCCGGTCTTCTTCTCGGCCATGACCGATACCTCCCAAGCCGCTTGGCGGC  >  NC_002947/1751222‑1751485                                                                                                                                 |                                                                                                                                        cAAGGCGACGCGCCCGCTCTTCAGGATCCTGCGCGATCCATGACAGCGCCAGGTCCTGCCAGAGTTCCTGCACCTGGTCGTAGCCGTGCCGCTCGATCTCCGCAGCCATGGCAGCCTTGATCCCTTCCGGCATGtt                                                                                                                                  >  2:195510/1‑136 (MQ=255)              ccGCTCTTCAGGATCCTGCGCGATCCATGACAGCGCCAGGTCCTGCCAGAGTTCCTGCACCTGGTCGTAGCCGTGCCGCTCGATCTCCGCAGCCATGGCAGCCTTGATCCCTTCCGGCATGTTGAAGAAAACCTtc                                                                                                                     <  1:472922/136‑1 (MQ=255)                  tctTCAGGATCCTGCGCGATCCATGACAGCGCCAGGTCCTGCCAGAGTTCCTGCACCTGGTCGTAGCCGTGCCGCTCGATCTCCGCAGCCATGGCAGCCTTGATCCCTTCCGGCATGTTGAAGAAAACCTTCTCGa                                                                                                                 >  1:163771/1‑136 (MQ=255)                   ctTCAGGATCCTGCGCGATCCATGACAGCGCCAGGTCCTGCCAGAGTTCCTGCACCTGGTCGTAGCCGTGCCGCTCGATCTCCGCAGCCATGGCAGCCTTGATCCCTTCCGGCATGTTGAAGAAAACCTTCTCGAt                                                                                                                <  1:230020/136‑1 (MQ=255)                              tGCGCGATCCATGACAGCGCCAGGTCCTGCCAGAGTTCCTGCACCTGGTCGTAGCCGTGCCGCTCGATCTCCGCAGCCATGGCAGCCTTGATCCCTTCCGGCATGTTGAAGAAAACCTTCTCGATCCCCAGTTTcg                                                                                                     >  2:170304/1‑136 (MQ=255)                                 gcgATCCATGACAGCGCCAGGTCCTGCCAGAGTTCCTGCACCTGGTCGTAGCCGTGCCGCTCGATCTCCGCAGCCATGGCAGCCTTGATCCCTTCCGGCATGTTGAAGAAAACCTTCTCGATCCCCAGTTTCGCgg                                                                                                  >  1:288593/1‑136 (MQ=255)                                   gATCCATGACAGCGCCAGGTCCTGCCAGAGTTCCTGCACCTGGTCGTAGCCGTGCCGCTCGATCTCCGCAGCCATGGCAGCCTTGATCCCTTCCGGCATGTTGAAGAAAACCTTCTCGATCCCCAGTTTCGCGGCt                                                                                                >  1:268121/1‑136 (MQ=255)                                   gATCCATGACAGCGCCAGGTCCTGCCAGAGTTCCTGCACCTGGTCGTAGCCGTGCCGCTCGATCTCCGCAGCCATGGCAGCCTTGATCCCTTCCGGCATGTTGAAGAAAACCTTCTCGATCCCCAGTTTCGCGGCt                                                                                                >  1:101715/1‑136 (MQ=255)                                   gATCCATGACAGCGCCAGGTCCTGCCAGAGTTCCTGCACCTGGTCGTAGCCGTGCCGCTCGATCTCCGCAGCCATGGCAGCCTTGATCCCTTCCGGCATGTTGAAGAAAACCTTCTCGATCCCCAGTTTCGCGGCt                                                                                                >  1:77323/1‑136 (MQ=255)                                   gATCCATGACAGCGCCAGGTCCTGCCAGAGTTCCTGCACCTGGTCGTAGCCGTGCCGCTCGATCTCCGCAGCCATGGCAGCCTTGATCCCTTCCGGCATGTTGAAGAAAACCTTCTCGATCCCCAGTTTCGCGGCt                                                                                                >  1:377494/1‑136 (MQ=255)                                   gATCCATGACAGCGCCAGGTCCTGCCAGAGTTCCTGCACCTGGTCGTAGCCGTGCCGCTCGATCTCCGCAGCCATGGCAGCCTTGATACCTTCCGGCATGTTGAAGAAAACCTTCTCGATCCCCAGTTTCGCGGCt                                                                                                >  1:76641/1‑136 (MQ=255)                                   gATCCATGACAGCGCCAGGTCCTGCCAGAGTTCCTGCACCTAGTCGTAGCCGTGCCGCTCGATCTCCGCAGCCATGGCAGCCTTGATCCCTTCCGGCATGTTGAAGAAAACCTTCTCGATCCCCAGTTTCGCGGCt                                                                                                <  1:245159/136‑1 (MQ=255)                                       cATGACAGCGCCAGGTCCTGCCAGAGTTCCTGCACCTGGTCGTAGCCGTGCCGCTCGATCTCCGCAGCCATGGCAGCCTTGATCCCTTCCGGCATGTTGAAGAAAACCTTCTCGATCCCCAGTTTCGCGGCTGCct                                                                                            >  2:102458/1‑136 (MQ=255)                                             cgcgcCAGGTCCTGCCAGAGTTCCTGCACCTGGTCGTAGCCGTGCCGCTCGATCTCCGCAGCCATGGCAGCCTTGATCCCTTCCGGCATGTTGAAGAAAACCTTCTCGATCCCCAGTTTCGCGGCTGCCTtcttct                                                                                      <  2:379913/135‑1 (MQ=255)                                             aGCGCCAGGTCCTGCCAGAGTTCCTGCACCTGGTCGTAGCCGTGCCGCTCGATCTCCGCAGCCATGGCAGCCTTGATCCCTTCCGGCATGTTGAAGAAAACCTTCTCGATCCCCAGTTTCGCGGCTGCCTtcttct                                                                                      <  1:411676/136‑1 (MQ=255)                                             aGCGCCAGGTCCTGCCAGAGTTCCTGCACCTGGTCGTAGCCGTGCCGCTCGATCTCCGCAGCCATGGCAGCCTTGATCCCTTCCGGCATGTTGAAGAAAACCTTCTCGATCCCCAGTTTCGCGGCTGCCTtcttct                                                                                      <  2:467006/136‑1 (MQ=255)                                                  cAGGTCCTGCCAGAGTTCCTGCACCTGGTCGTAGCCGTGCCGCTCGATCTCCGCAGCCATGGCAGCCTTGATCCCTTCCTGCATGTTGAAGAAAACCTTCTCGATCCCCAGTTTCGCGGCTGCCTTCTTCTGGCgg                                                                                 >  2:326135/1‑136 (MQ=255)                                                  cAGGTCCTGCCAGAGTTCCTGCACCTGGTCGTAGCCGTGCCGCTCGATCTCCGCAGCCATGGCAGCCTTGATCCCTTCCTGCATGTTGAAGAAAACCTTCTCGATCCCCAGTTTCGCGGCTGCCTTCTTCTGGCgg                                                                                 >  2:328710/1‑136 (MQ=255)                                                  cAGGTCCTGCCAGAGTTCCTGCACCTGGTCGTAGCCGTGCCGCTCGATCTCCGCAGCCATGGCAGCCTTGATCCCTTCCGGCATGTTGAAGAAAACCTTCTCGATCCCCAGTTTCGCGGCTGCCTTCTTCTGGCgg                                                                                 >  1:455493/1‑136 (MQ=255)                                                            cAGAGTTCCTGCACCTGGTCGTAGCCGTGCCGCTCGATCTCCGCAGCCATGGCAGCCTTGATCCCTTCCGGCATGTTGAAGAAAACCTTCTCGATCCCCAGTTTCGCGGCTGCCTTCTTCTGGCGGTCTCGGTAGt                                                                       >  2:312362/1‑136 (MQ=255)                                                            cAGAGTTCCTGCACCTGGTCGTAGCCGTGCCGCTCGATCTCCGCAGCCATGGCAGCCTTGATCCCTTCCAGCATGTTGAAGAAAACCTTCTCGATCCCCAGTTTCGCGGCTGCCTTCTTCTGGCGGTCTCGGTAGt                                                                       >  2:289290/1‑136 (MQ=255)                                                                              tCGTAGCCGTGCCGCTCGATCTCCGCAGCCATGGCAGCCTTGATCCCTTCCGGCATGTTGAAGAAAACCTTCTCGATCCCCAGTTTCGCGGCTGCCTTCTTCTGGCGGTCTCGGTAGTCCGCTGAGTGCTTCGCgg                                                     >  2:290623/1‑136 (MQ=255)                                                                              gcGTAGCCGTGCCGCTCGATCTCCGCAGCCATGGCAGCCTTGATCCCTTCCGGCATGTTGAAGAAAACCTTCTCGATCCCCAGTTTCGCGGCTGCCTTCTTCTGGCGGTCTCGGTAGTCCGCTGAGTGCTTCGCgg                                                     <  1:78161/135‑1 (MQ=255)                                                                                    ccGTGCCGCTCGATCTCCGCAGCCATGGCAGCCTTGATCCCTTCCGGCATGTTGAAGAAAACCTTCTCGATCCCCAGTTTCGCGGCTGCCTTCTTCTGGCGGTCTCGGTAGTCCGCTGAGTGCTTCGCGGCTCCgg                                               <  2:467771/136‑1 (MQ=255)                                                                                      gTGCCGCTCGATCTCCGCAGCCATGGCAGCCTTGATCCCTTCCGGCATGTTGAAGAAAACCTTCTCGATCCCCAGTTTCGCGGCTGCCTTCTTCTGGCGGTCTCGGTAGTCCGCTGAGTGCTTCGCGGCTCCGGtc                                             >  2:432570/1‑136 (MQ=255)                                                                                          cGCTCGATCTCCGCAGCCATGGCAGCCTTGATCCCTTCCGGCATGTTGAAGAAAACCTTCTCGATCCCCAGTTTCGCGGCTGCCTTCTTCTGGCGGTCTCGGTAGTCCGCTGAGTGCTTCGCGGCTCCGGtcttct                                         >  1:344774/1‑136 (MQ=255)                                                                                               gATCTCCGCAGCCATGGCAGCCTTGATCCCTTCCGGCATGTTGAAGAAAACCTTCTCGATCCCCAGTTTCGCGGCTGCCTTCTTCTGGCGGTCTCGGTAGTCCGCTGAGTGCTTCGCGGCTCCGGTCTTCTTCTCg                                    <  2:328398/136‑1 (MQ=255)                                                                                               caTCTCCGCAGCCATGGCAGCCTTGATCCCGTCCGGCATGTTGAAGAAAACCTTCTCGATCCCCAGTTTCGCGGCTGCCTTCTTCTGGCGGTCTCGGTAGTCCGCTGAGTGCTTCGCGGCTCCGGTCTTCTTCTCg                                    <  2:79075/135‑1 (MQ=255)                                                                                                        aGCCATGGCAGCCTTGATCCCTTCCTGCATGTTGAAGAAAACCTTCTCGATCCCCAGTTTCGCGGCTGCCTTCTTCTGGCGGTCTCGGTAGTCCGCTGAGTGCTTCGCGGCTCCGGTCTTCTTCTCGGCCATGAcc                           <  1:326135/136‑1 (MQ=255)                                                                                                        aGCCATGGCAGCCTTGATCCCTTCCTGCATGTTGAAGAAAACCTTCTCGATCCCCAGTTTCGCGGCTGCCTTCTTCTGGCGGTCTCGGTAGTCCGCTGAGTGCTTCGCGGCTCCGGTCTTCTTCTCGGCCATGAcc                           <  1:328710/136‑1 (MQ=255)                                                                                                         gCCATGGCAGCCTTGATCCCTTCCGGCATGTTGAAGAAAACCTTCTCGATCCCCAGTTTCGCGGCTGCCTTCTTCTGGCGGTCTCGGTAGTCCGCTGAGTGCTTCGCGGCTCCGGTCTTCTTCTCGGCCATGACCg                          <  2:377494/136‑1 (MQ=255)                                                                                                                cAGCCTTGATCCCTTCCGGCATGTTGAAGAAAACCTTCTCGATCCCCAGTTTCGCGGCTGCCTTCTTCTGGCGGTCTCGGTAGTCCGCTGAGTGCTTCGCGGCTCCGGTCTTCTTCTCGGCCATGACCGATACCTc                   <  2:363275/136‑1 (MQ=255)                                                                                                                cAGCCTTGATCCCTTCCGGCATGTTGAAGAAAACCTTCTCGATCCCCAGTTTCGCGGCTGCCTTCTTCTGGCGGTCTCGGTAGTCCGCTGAGTGCTTCGCGGCTCCGGTCTTCTTCTCGGCCATGACCGATACCTc                   <  2:91090/136‑1 (MQ=255)                                                                                                                       gATCCCTTCCGGCATGTTGAAGAAAACCTTCTCGATCCCCAGTTTCGCGGCTGCCTTCTTCTGGCGGTCTCGGTAGTCCGCTGAGTGCTTCGCGGCTCCGGTCTTCTTCTCGGCCATGACCGATACCTCCCAAGcc            >  2:470591/1‑136 (MQ=255)                                                                                                                                cGGCATGTTGAAGAAAACCTTCTCGATCCCCAGTTTCGCGGCTGCCTTCTTCTGGCGGTCTCGGTAGTCCGCTGAGTGCTTCGCGGCTCCGGTCTTCTTCTCGGCCATGACCGATACCTCCCAAGCCGCTTggcgg   >  1:183132/1‑136 (MQ=255)                                                                                                                                 ggCATGTTGAAGAAAACCTTCTCGATCCCCAGTTTCGCGGCTGCCTTCTTCTGGCGGTCTCGGTAGTCCGCTGAGAGCTTCGCGGCTCCGGTCTTCTTCTCGGCCATGACCGATACCTCCCAAGCCGCTTggcggc  >  1:247063/1‑136 (MQ=255)                                                                                                                                 |                                                                                                                                        CAAGGCGACGCGCCCGCTCTTCAGGATCCTGCGCGATCCATGACAGCGCCAGGTCCTGCCAGAGTTCCTGCACCTGGTCGTAGCCGTGCCGCTCGATCTCCGCAGCCATGGCAGCCTTGATCCCTTCCGGCATGTTGAAGAAAACCTTCTCGATCCCCAGTTTCGCGGCTGCCTTCTTCTGGCGGTCTCGGTAGTCCGCTGAGTGCTTCGCGGCTCCGGTCTTCTTCTCGGCCATGACCGATACCTCCCAAGCCGCTTGGCGGC  >  NC_002947/1751222‑1751485

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