Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,708,363	C→T	8.1%	A794V (GCG→GTG)	PP_5009 →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	5,708,363	0	C	T	8.1%	140.0 / 2.2	49	A794V (GCG→GTG)	PP_5009	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (25/20);  new base T (2/2);  total (27/22)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.73e-01

GGAGGCGCTGATCAAGCGCAGCGGTGACAACAAGATCGTCGCAGCCATAGACCTGATCACCAACCTGGCCTTGCTCCTTGCCCACCGCCCTGCAACGCCACCCACCAACAGCGAAGAGGCCGC‑TACGCTCCGTTTTTCAGGCCCTGCGCCGAAGGAAGATTCCCCGCCGCCAGCGCTCGAGCAGCCACAGGAGAAACCCTGGCAGGCCCCCGTCGAGCCGCCGCCGCCCAGGATCATGCGGGTAAATGCCTGGCGTGAGGAGCAGAGCCTGAGCAAGCTGTCGCCCCAGGTGCGT  >  NC_002947/5708217‑5708511                                                                                                                                                    |                                                                                                                                                     ggAGGCGCTGATCAAGCGCAGCGGTGACAACAAGATCGTCGCAGCCATAGACCTGATCACCAACCTGGCCTTGCTCCTTGCCCACCGCCCTGCAACGCCACCCACCAACAGCGAAGAGGCCGC‑TACGCTCCGTTTTTCAGGCCCTgcgc                                                                                                                                                    <  8:193756/149‑1 (MQ=255)     gcgcTGATCAAGCGCAGCGGTGACAACAAGATCGTCGCAGCCATAGACCTGATCACCAACCTGGCCTTGCTCCTTGCCCACCGCCCTGCAACGCCACCCACCAACAGCGAAGAGGCCGC‑TACGCTCCGTTTTTCAGGCCCTGCGCCgaa                                                                                                                                                <  5:77862/149‑1 (MQ=255)     gcgcTGATCAAGCGCAGCGGTGACAACAAGATCGTCGCAGCCATAGACCTGATCACCAACCTGGCCTTGCTCCTTGCCCACCGCCCTGCAACGCCACCCACCAACAGCGAAGAGGCCGC‑TACGCTCCGTTTTTCAGGCCCTGCGCCgaa                                                                                                                                                <  2:241352/149‑1 (MQ=255)     gcgcTGATCAAGCGCAGCGGTGACAACAAGATCGTCGCAGCCATAGACCTGATCACCAACCTGGCCTTGCTCCTTGCCCACCGCCCTGCAACGCCACCCACCAACAGCGAAGAGGCCGC‑TACGCTCCGTTTTTCAGGCCCTGCGCCgaa                                                                                                                                                <  4:336483/149‑1 (MQ=255)           aTCAAGCGCAGCGGTGACAACAAGATCGTCGCAGCCATAGACCTGATCACCAACCTGGCCTTGCTCCTTGCCCACCGCCCTGCAACGCCACCCACCAACAGCGAAGAGGCCGC‑TACGCTCCGTTTTTCAGGCCCTGCGCCgaagg                                                                                                                                              >  6:336993/1‑145 (MQ=255)                                   aTCGTCGCAGCCATAGACCTGATCACCAACCTGGCCTTGCTCCTTGCCCACCGCCCTGCAACGCCACCCACCAACAGCGAAGAGGCCGC‑TACGCTCCGTTTTTCAGGCCCTGCGCCGAAGGAAGATTCCCCGCCGCCAGCGCTCGagca                                                                                                                  <  8:101827/149‑1 (MQ=255)                                           aGCCATAGACCTGATCACCAACCTGGCCTTGCTCCTTGCCCACCGCCCTGCAACGCCACCCACCAACAGCGAAGAGGCCGC‑TACGCTCCGTTTTTCAGGCCCTGCGCCGAAGGAAGATTCCCCGCCGCCAGCGCTCGAGCAGCCACAgg                                                                                                          <  4:318021/149‑1 (MQ=255)                                                 aGACCTGATCACCAACCTGGCCTTGCTCCTTGCCCACCGCCCTGCAACGCCACCCACCAACAGCGAAGAGGCCGC‑TACGCTCCGTTTTTCAGGCCCTGTGCCGAAGGAAGATTCCCCGCCGCCAGCGCTCGAGCAGCCACAGGAGAAAc                                                                                                    >  6:187662/1‑149 (MQ=255)                                                  gACCTGATCACCAACCTGGCCTTGCTCCTTGCCCACCGCCCTGCAACGCCCCCCACCAACAGCGAAGAGGCCGC‑TACGCTCCGTTTTTCAGGCCCTGCGCCGAAGGAAGATTCCCCGCCGCTAGCGCTCGAGCAGCCACAGGAGAAAcc                                                                                                   >  8:2760/1‑149 (MQ=255)                                                     cTGATCACCAACCTGGCCTTGCTCCTTGCCCACCGCCCTGCAACGCCACCCACCAACAGCGAAGAGGCCGC‑TACGCTCCGTTTTTCAGGCCCTGCGCCGAAGGAAGATTCCCCGCCGCCAGCGCTCGAGCAGCCACAGGAGAAACCCTg                                                                                                >  6:64021/1‑149 (MQ=255)                                                        aTCACCAACCTGGCCTTGCTCCTTGCCCACCGCCCTGCAACGCCACCCACCAACAGCGAAGAGGCCGC‑TACGCTCCGTTTTTCAGGCCCTGCGCCGAAGGAAGATTCCCCGCCGCCAGCGCTCGAGCAGCCACAGGAGAAACCCTgg                                                                                               >  2:222018/1‑147 (MQ=255)                                                             cAACCTGGCCTTGCTCCTTGCCCACCGCCCTGCAACGCCACCCACCAACAGCGAAGAGGCCGC‑TACGCTCCGTTTTTCAGGCCCTGCGCCGAAGGAAGATTCCCCGCCGCCAGCGCTCGAGCAGCCACAGGAGAAACCCTGGCAGGccc                                                                                        >  8:208890/1‑149 (MQ=255)                                                                              ttGCCCACCGCCCTGCAACGCCACCCACCAACAGCGAAGAGGCCGC‑TACGCTCCGTTTTTCAGGCCCTGCGCCGAAGGAAGATTCCCCGCCGCCAGCGCTCGAGCAGGCACAGGAGAAACCCTgg                                                                                               >  5:310128/1‑125 (MQ=255)                                                                              ttGCCCACCGCCCTGCAACGCCACCCACCAACAGCGAAGAGGCCGC‑TACGCTCCGTTTTTCAGGCCCTGCGCCGAAGGAAGATTCCCCGCCGCCAGCGCTCGAGCAGCCACAGGAGAAACCCTgg                                                                                               <  6:310128/125‑1 (MQ=255)                                                                                gCCCACCGCCCTGCAACGCCACCCACCAACAGCGAAGAGGCCGC‑TACGCTCCGTTTTTCAGGCCCTGCGCCGAAGGAAGATTCCCCGCCGCCAGCGCTCGAGCAGCCACAGGAGAAACCCTGGCAGGCCCCCGTCGAgccgccgccgcc                                                                     >  8:108511/1‑149 (MQ=255)                                                                                  ccACCGCCCTGCAACGCCACCCACCAACAGCGAATAGGCAGC‑AACGCTCCGTTTTTCAGGCCCTGTGCCCAAGGAAGATTCCCCGCCGCCAGCGCTCGAGCAGCCACAGCAGAAACCCTGGCAGGCCCCCGTCGAGCCGCCGCCTCCCa                                                                   >  6:79958/1‑149 (MQ=255)                                                                                      cGCCCTGCAACGCCACCCACCAACAGCGAAGAGGCCGC‑TACGCTCCGTTTTTCAGGCCCTGCGCCGAAGGAAGATTCCCCGCCGCCAGCGCTCGAGCAGCCACAGGAGAAACCCTGGCAGGCCCCCGTCGAGCCGCCGCCGCCCAGGAt                                                               >  4:345747/1‑149 (MQ=255)                                                                                        cccTGCAACGCCACCCACCAACAGCGAAGAGGCCGCTTACGCTCCGTTTTTCAGGCCCTGCGCCGAAGGAAGATTCCCCGCCGCCAGCGCTCGAGCAGCCACAGGAGAAACCCTGGCAGGCCCCCGTCGAGCCGCCGCCGCCCAGGAt                                                               >  1:64198/1‑148 (MQ=255)                                                                                          cTGCAACGCCACCCACCAACAGCGAAGAGGCCGC‑TACGCTCCGTTTTTCAGGCCCTGCGCCGAAGGAAGATTCCCCGCCGCCAGCGCTCGAGCAGCCACAGGAGAAACCCTGGCAGGCCCCCGTCGAGCCGCCGCCGCCCAGGATCATg                                                           <  2:64198/149‑1 (MQ=255)                                                                                            gCAACGCCACCCACCAACAGCGAAGAGGCCGC‑TACGCTCCGTTTTTCAGGCCCTGCGCCGAAGGAAGAATCCCCGCCGCCAGCGCTCGAGCAGCCACAGGAGAAACCCTGGCAGGCCCCCGTCGAGCCGCCGCCGCCCAGGAGCATGCg                                                         >  7:261052/1‑149 (MQ=255)                                                                                            gCAAAGCCACCCACCAACAGCGAAGAGGCCGC‑TACGCTCCGTTTTTCAGGCCCTGCGCCGAAGGAAGATTCCCCGCCGCCAGCGCTCGAGCAGCCACAGGAGAAACCCTGGCAGGCCCCCGTCGAGCCGCCGCCGCCCAGGATCATGCg                                                         >  6:326729/1‑149 (MQ=255)                                                                                             cAACGCCACCCACCAACAGCGAAGAGGCCGC‑TACGCTCCGTTTTTCAGGCCCTGTGCCGAAGGAAGATTCCCCGCCGCCAGCGCTCGAGCAGCCACAGGAGAAACCCTGGCAGGCCCCCGTCGAGCCGCCGGCGCCCAGGATCATGCgg                                                        <  5:79958/149‑1 (MQ=255)                                                                                             cAACGCCACCCACCAACAGCGAAGAGGCCGC‑TACGCTCCGTTTTTCAGGCCCTGTGCCGAAGGAAGATTCCCCGCCGCCAGCGCTCGAGCAGCCACAGGAGAAACCCTGGCAGGCCCCCGTCGAGCCGCCGCCGCCCAGGATCATGCgg                                                        <  5:187662/149‑1 (MQ=255)                                                                                             cAACGCCACCCACCAACAGCGAAGAGGCCGC‑TACGCTCCGTTTTTCAGGCCCTGCGCCGAAGGAAGATTCCCCGCCGCCAGCGCTCGAGCAGCCACAGGAGAAACCCTGGCAGGCCCCCGTCGAGCCGCCGCCGCCCAGGATCATGCg                                                         <  5:147046/148‑1 (MQ=255)                                                                                                   cacccaccAACAGCGAAGAGGCCGC‑TACGCTCCGTTTTTCAGGCCCTGCGCCGAAGGAAGATTCCCCGCCGCCAGCGCTCGAGCAGCCACAGGAGAAACCCTGGCAGGCCCCCGTCGAGCCGCCGCCGCCCAGGATCATGCGGGTAAAt                                                  <  5:64021/149‑1 (MQ=255)                                                                                                   cacccaccAACAGCGAAGAGGCCGC‑TACGCTCCGTTTTTCAGGCCCTGCGCCGAAGGAAGATTCCCCGCCGCCAGCGCTCGAGCAGCCACAGGAGAAACCCTGGCAGGCCCCCGTCGAGCCGCCGCCGCCCAGGATCATGCGGGTAAAt                                                  >  8:106103/1‑149 (MQ=255)                                                                                                      ccaccaACAGCGAAGAGGCCGC‑TACGCTCCGTTTTTCAGGCCCTGCGCCGAAGGAAGATTCCCCGCCGCCAGCGCTCGAGCAGCCACAGGAGAAACCCTGGCAGGCCCCCGTCGAGCCGCCGCCGCCCAGGATCATGCGGGTAAATGcc                                               >  2:4329/1‑149 (MQ=255)                                                                                                        accaACAGCGAAGAGGCCGC‑TACGCTCCGTTTTTCAGGCCCTGCGCCGAAGGAAGATTCCCCGCCGCCAGCGCTCGAGCAGCCACAGGAGAAACCCTGGCAGGCCCCCGTCGAGCCGCCGCCGCCCAGGATCATGCGGGTAAATGCCTg                                             >  5:40378/1‑149 (MQ=255)                                                                                                        accaACAGCGAAGAGGCCGC‑TACGCTCCGTTTTTCAGGCCCTGCGCCGAAGGAAGATTCCCCGCCGCCAGCGCTCGAGCAGCCACAGGAGAAACCCTGGCAGGCCCCCGTCGAGCCGCCGCCGCCCAGGATCATGCGGGTAAATGCCTg                                             >  6:136339/1‑149 (MQ=255)                                                                                                          caACAGCGAAGAGGCCGC‑TACGCTCCGTTTTTCAGGCCCTGCGCCGAAGGAAGATTCCCCGCCGCCAGCGCTCGAGCAGCCACAGGAGAAACCCTGGCAGGCCCCCGTCGAGCCGCCGCCGCCCAGGATCATGCGGGTAAATGCCTGGc                                           >  5:230690/1‑149 (MQ=255)                                                                                                          caACAGCGAAGAGGCCGC‑TACGCTCCGTTTTTCAGGCCCTGCGCCGAAGGAAGATTCCCCGCCGCCAGCGCTCGAGCAGCCACAGGAGAAACCCTGGCAGGCCCCCGTCGAGCCGCCGCCGCCCAGGATCATGCGGGTAAATGCCTGGc                                           >  8:185349/1‑149 (MQ=255)                                                                                                          caACAGCGAAGAGGCCGC‑TACGCTCCGTTTTTCAGGCCCTGCGCCGAAGGAAGATTCCCCGCCGCCAGCGCTCGAGCAGCCACAGGAGAAACCCTGGCAGGCCCCCGTCGAGCCGCCGCCGCCCAGGATCATGCGGGTAAATGCCTGGc                                           <  7:208890/149‑1 (MQ=255)                                                                                                             cAGCGAAGAGGCCGC‑TACGCTCCGTTTTTCAGGCCCTGCGCCGAAGGAAGATTCCCCGCCGCCAGCGCTCGAGCAGCCACAGGAGAAACCCTGGCAGGCCCCCGTCGAGCCGCCGCCGCCCAGGATCATGCGGGTAAATGCCTGGCGTg                                        <  5:326729/149‑1 (MQ=255)                                                                                                               gCGAAGAGGCCGC‑TACGCTCCGTTTTTCAGGCCCTGCGCCGAAGGAAGATTCCCCGCCGCTAGCGCTCGAGCAGCCACAGGAGAAACCCTGGCAGGCCCCCGTCGAGCCGCCGCCGCCCAGGATCATGCGGGTAAATGCCTGGCGTgag                                      <  7:2760/149‑1 (MQ=255)                                                                                                               gCGAAGAGGCCGC‑TACGCTCCGTTTTTCAGGCCCTGCGCCGAAGGAAGATTCCCCGCCGCCAGCGCTCGAGCAGCCACAGGAGAAACCCTGGCAGGCCCCCGTCGAGCCGCCGCCGCCCAGGATCATGCGGGTAAATGCCTGGCGTgag                                      <  7:106103/149‑1 (MQ=255)                                                                                                                cGAAGAGGCCGC‑TACGCTCCGTTTTTCAGGCCCTGCGCCGAAGGAAGATTCCCCGCCGCCAGCGCTCGAGCAGCCACAGGAGAAACCCTGGCAGGCCCCCGTCGAGCCGCCGCCGCCCAGGATCATGCGGGTAAATGCCTGGCGTgag                                      >  6:228111/1‑148 (MQ=255)                                                                                                                         cGC‑TACGCTCCGTTTTTCAGGCCCTGCGCCGAAGGAAGATTCCCCGCCGCCAGCGCTCGAGCAGCCACAGGAGAAACCCTGGCAGGCCCCCGTCGAGCCGCCGCCGCCCAGGATCATGCGGGTAAATGCCTGGCGTGAGGAGCAGAGcc                            >  7:220237/1‑149 (MQ=255)                                                                                                                          gC‑TACGCTCCGTTTTTCAGGCCCTGCGCCGAAGGAAGATTCCCCGCCGCCAGCGCTCGAGCAGCCACAGGAGAAACCCTGGCAGGCCCCCGTCGAGCCGCCGCCGCCCAGGATCATGCGGGTAAATGCCTGGCGTGAGGAGCAGAGCCt                           <  7:185349/149‑1 (MQ=255)                                                                                                                               cGCTCCGTTTTTCAGGCCCTGCGCCGAAGGAAGATTCCCCGCCGCCAGCGCTCGAGCAGCCACAGGAGAAACCCTGGCAGGCCCCCGTCGAGCCGCCGCCGCCCAGGATCATGCGGGTAAATGCCTGGCGTGAGGAGCAGAGCCTGa                         <  1:4329/147‑1 (MQ=255)                                                                                                                                 cTCCGTTTTTCAGGCCCTGCGCCGAAGGAAGATTCCCCGCCGCCAGCGCTCGAGCAGCCACAGGAGAAACCCTGGCAGGCCCCCGTCGAGCCGCCGCCGCCCAGGATCATGCGGGTAAATGCCTGGCGTGAGGAGCAGAGCCTGAGCaa                     <  5:336993/149‑1 (MQ=255)                                                                                                                                  tCCGTTTTTCAGGCCCTGCGCCGAAGGAAGATTCCCCGCCGCCAGCGCTCGAGCAGCCACAGGAGAAACCCTGGCAGGCCCCCGTCGAGCCGCCGCCGCCCAGGATCATGCGGGTAAATGCCTGGCGTGAGGAGCAGAGCCTGAg                        >  5:322021/1‑145 (MQ=255)                                                                                                                                  tCCGTTTTTCAGGCCCTGCGCCGAAGGAAGATTCCCCGCCGCCAGCGCTCGAGCAGCCACAGGAGAAACCCTGGCAGGCCCCCGTCGAGCCGCCGCCGCCCAGGATCATGCGGGTAAATGCCTGGCGTGAGGAGCAGAGCCTGAg                        <  6:322021/145‑1 (MQ=255)                                                                                                                                     gTTTTTCAGGCCCTGCGCCGAAGGAAGATTCCCCGCCGCCAGCGCTCGAGCAGCCACAGGAGAAACCCTGGCAGGCCCCCGTCGAGCCGCCGCCGCCCAGGATCATGCGGGTAAATGCCTGGCGTGAGGAGCAGAGCCTGAGCAAGCTg                 >  2:359943/1‑149 (MQ=255)                                                                                                                                      tttttCAGGCCCTGCGCCGAAGGAAGATTCCCCGCCGCCAGCGCTCGAGCAGCCACAGGAGAAACCCTGGCAGGCCCCCGTCGAGCCGCCGCCGCCCAGGATCATGCGGGTAAATGCCTGGCGTGAGGAGCAGAGCCTGAGCAAGCTGt                >  3:112004/1‑149 (MQ=255)                                                                                                                                        tttCAGGCCCTGCGCCGAAGGAAGATTCCCCGCCGCCAGCGCTCGAGCAGCCACAGGAGAAACCCTGGCAGGCCCCCGTCGAGCCGCCGCCGCCCAGGATCATGCGGGTAAATGCCTGGCGTGAGGAGCAGAGCCTGAGCAAGCTGTc               <  4:112004/148‑1 (MQ=255)                                                                                                                                              gCCCTGCGCCGAAGGAAGATTCCCCGCCGCCAGCGCTCGAGCAGCCACAGGAGAAACCCTGGCAGGCCCCCGTCGAGCCGCCGCCGCCCAGGATCATGCGGGTAAATGCCTGGCGTGAGGAGCAGAGCCTGAGCa                      >  4:157951/1‑135 (MQ=255)                                                                                                                                              gCCCTGCGCCGAAGGAAGATTCCCCGCCGCCAGCGCTCGAGCAGCCACAGGAGAAACCCTGGCAGGCCCCCGTCGAGCCGCCGCCGCCCAGGATCATGCGGGTAAATGCCTGGCGTGAGGAGCAGAGCCTGAGCa                      <  3:157951/135‑1 (MQ=255)                                                                                                                                                ccTGCGCCGAAGGAAGATTCCCCGCCGCCAGCGCTCGAGCAGCCACAGGAGAAACCCTGGCAGGACTCCGTCGAGCCGCCGCCGCCCAGGATCATGCGGGTAAATGCCTGGCGTGAGGAGCAGAGCCTGAGCAAGCTGTCGCCCCAGGt      >  5:170422/1‑149 (MQ=255)                                                                                                                                                  tGCGCCGAAGGAAGATTCCCCGCCGCCAGCGCTCGAGCAGCCACAGGAGAAACCCTgg                                                                                               >  8:339714/1‑58 (MQ=255)                                                                                                                                                  tGCGCCGAAGGAAGATTCCCCGCCGCCAGCGCTCGAGCAGCCACAGGAGAAACCCTgg                                                                                               <  7:339714/58‑1 (MQ=255)                                                                                                                                                   gcgcCGAAGGAAGATTCCCCGCCGCCAGCGCTCGAGCAGCCACAGGAGAAACCCTGGCAGGCCCCCGTCGAGCCGCCGCCGCCCAGGATCATGCGGGTAAATGCCTGGCGTGAGGAGCAGAGCCTGAGCAAGCTGTCGCCCCAGGTGCg   >  4:142513/1‑149 (MQ=255)                                                                                                                                                    cgcCGAAGGAAGATTCCCCGCCGCCAGCGCTCGAGCAGCCACAGGAGAAACCCTGGCAGGCCCCCGTCGAGCCGCCGCCGCCCAGGATCATGCGGGTAAATGCCTGGCGTGAGTAGCAGAGCCTGAGCAAGCTGTCGCCCCAGGTGCgt  <  7:108511/149‑1 (MQ=255)                                                                                                                                                    |                                                                                                                                                     GGAGGCGCTGATCAAGCGCAGCGGTGACAACAAGATCGTCGCAGCCATAGACCTGATCACCAACCTGGCCTTGCTCCTTGCCCACCGCCCTGCAACGCCACCCACCAACAGCGAAGAGGCCGC‑TACGCTCCGTTTTTCAGGCCCTGCGCCGAAGGAAGATTCCCCGCCGCCAGCGCTCGAGCAGCCACAGGAGAAACCCTGGCAGGCCCCCGTCGAGCCGCCGCCGCCCAGGATCATGCGGGTAAATGCCTGGCGTGAGGAGCAGAGCCTGAGCAAGCTGTCGCCCCAGGTGCGT  >  NC_002947/5708217‑5708511

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