Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,799,626	C→T	7.4%	V79V (GTG→GTA)	PP_5077 ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	5,799,626	0	C	T	7.4%	159.0 / 2.3	54	V79V (GTG→GTA)	PP_5077	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (23/27);  new base T (2/2);  total (25/29)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.53e-01

CTTCCTGCCCGGTCAGTGCCAGCTGCACCACCTTCGACAGGATCGCCTGCACCTCGGGCTGGGCGACACCCAGCTCCTGCGCAACCTGCGCGAGCATACTGGAGGCATCTCTGGCGCTGCGGGCCGACACCACCACACTTTGCACAGCCTGCTTGTTGGTGCTGGCCACCAGGGCCTGGC‑GCAGCAGCGACTTGCCGCTGCCCAATGGGCCACTGACCACCAGCATCAACTGGCTATAGCGCGCCAGGTGGTGCAGTTGGCCAAGCACGGGCTTGCGCTGGGCCGGGAAG  >  NC_002947/5799484‑5799773                                                                                                                                               |                                                                                                                                                     cTTGCTGCCCTGTCAGTGCCTGCTGCACCACCTTCGACAGGATCTCCTGCACCTCGGTCTGTTAGACACCCAGCTCCTGCGCAACCTTCGCGAGCATACTGGAGGCATCTCTGGCTGTGCGGGCCGACAACACCACACTTTGCACAGcc                                                                                                                                                <  8:191130/149‑1 (MQ=255)        cccGGTCAGTGCCAGCTGCACCACCTTCGACAGGATCGCCTGCACCTCGGGCTGGGCGACACCCAGCTCCTGCGCAACCTGCGCGAGCATACTGGAGGCATCTCTGGCGCTGCGGGCCGACACCACCACACTTTGCACAGCCTGCttgt                                                                                                                                         >  5:181645/1‑149 (MQ=255)          cGGTCAGTGCCAGCTGCACCACCTTCGACAGGATCGCCTGCACCTCGGGCTGGGCGACACCCAGCTCCTGCGCAACCTGCGCGAGCATACTGGAGGCATCTCTGGCGCTGCGGGCCGACACCACCACACTTTGCACAGCCTGCttgttg                                                                                                                                       >  8:35639/1‑149 (MQ=255)          cGGTCAGTGCCAGCTGCACCACCTTCGACAGGATCGCCTGCACCTCGGGCTGGGCGACACCCAGCTCCTGCGCAACCTGCGCGAGCATACTGGAGGCATCTCTGGCGCTGCGGGCCGACACCACCACACTTTGCACAGCCTGCttgttg                                                                                                                                       <  6:181645/149‑1 (MQ=255)                   ccAGCTGCACCACCTTCGACAGGATCGCCTGCACCTCGGGCTGGGCGACACCCAGCTCCTGCGCAACCTGCGCGAGCATACTGGAGGCATCTCTGGCGCTGCGGGCCGACACCACCACACTTTGCACAGCCTGCTTGTTGGTGCTGGc                                                                                                                               >  2:164818/1‑148 (MQ=255)                     aGCTGCACCACCTTCGACAGGATCGCCTGCACCTCGGGCTGGGCGACACCCAGCTCCTGCGCAACCTGCGCGAGCATACTGGAGGCATCTCTGGCGCTGCGGGCCGACACCACCACACTTTGCACAGCCTGCTTGTTGGTGCTGGccac                                                                                                                            <  6:280866/149‑1 (MQ=255)                             caccTTCGACAGGATCGCCTGCACCTCGGGCTGGGCGACACCCAGCTCCTGCGCAACCTGCGCGAGCATACTGGAGGCATCTCTGGCGCTGCGGGCCGACACCACCACACTTTGCACAGCCTGCTTGTTGGTGCTGGccacc                                                                                                                           <  7:287171/142‑1 (MQ=255)                             caccTTCGACAGGATCGCCTGCACCTCGGGCTGGGCGACACCCAGCTCCTGCGCAACCTGCGCGAGCATACTGGAGGCATCTCTGGCGCTGCGGGCCGACACCACCACACTTTGCACAGCCTGCTTGTTGGTGCTGGccacc                                                                                                                           >  8:287171/1‑142 (MQ=255)                              accTTCGACAGGATCGCCTGCACCTCGGGCTGGGCGACACCCAGCTCCTGCGCAACCTGCGCGAGCATACTGGAGGCATCTCTGGCGCTGCGGGCCGACACCACCACACTTTGCACAGCCTGCTTGTTGGTGCTGGCCACCAGGGCCTg                                                                                                                   <  8:183871/149‑1 (MQ=255)                              accTTCGACAGGATCGCCTGCACCTCGGGCTGGGCGACACCCAGCTCCTGCGCAACCTGCGCGAGCATACTGGAGGCATCTCTGGCGCTGCGGGCCGACACCACCACACTTTGCACAGCCTGCTTGTTGGTGCTGGCCACCAGGGCCTg                                                                                                                   <  1:227106/149‑1 (MQ=255)                                   cGACAGGATCGCCTGCACCGCGGGCTGGGCGACACCCAGCTCCTGCGCAACCTCCGCGAGCATACTGGAGGCATCTCTGGCGCTGCGGGCCGACACCACCACACTTTGCACAGCCTGCTTGTTGGTGCTGGCCACCAGGGCCTGgc‑gc                                                                                                              <  5:221820/148‑1 (MQ=255)                                    gACAGGATCGCCTGCACCTCGGGCTGGGCGACACCCAGCTCCTGCGCAACCTGCGCGAGCATACTGGAGGCATCTCTGGCGCTGCGGGCCGACACCACCACACTTTGCACAGCCTGCTTGTTGGTGCTGGCCACCAGGGCCTGGC‑gca                                                                                                             <  2:44510/148‑1 (MQ=255)                                        ggATCGCCTGCACCTCGGGCTGGGCGACACCCAGCTCCTGCGCAACCTGCGCGAGCATACTGGAGGCATCTCTGGCGCTGCGGGCCGACACCACCACACTTTGCACAGCCTGCTTGTTGGTGCTGGCCACCAGGGCCTGGC‑GCagcagc                                                                                                        <  5:30225/149‑1 (MQ=255)                                              ccTGCACCTCGGGCTGGGCGACACCCAGCTCCTGCGCAACCTGCGCGAGCATACTGGAGGCATCTCTGGCGCTGCGGGCCGACACCACCACACTTTGCACAGCCTGCTTGTTGGTGCTGGCCACCAGGGCCTGGC‑GCAGCAGCGACTTg                                                                                                  >  5:273947/1‑149 (MQ=255)                                                  cACCTCGGGCTGGGCGACACCCAGCTCCTGCGCAACCTGCGCGAGCATACTGGAGGCATCTCTGGCGCTGCGGGCCGACACCACCACACTTTGCACAGCCTGCTTGTTGGTGCTGGCCACCAGGGCCTGGC‑GCAGCAGCGACTTGCCg                                                                                               <  7:35639/148‑1 (MQ=255)                                                   aCCTCGGGCTGGGCGACACCCAGCTCCTGCGCAACCTGCGCGAGCATACTGGAGGCATCTCTGGCGCTGCGGGCCGACACCACCACACTTTGCACAGCCTGCTTGTTGGTGCTGGCCACCAGGGCCTGGC‑GCAGCAGCGACTTGCCGCt                                                                                             >  7:297746/1‑149 (MQ=255)                                                     caggagCTGGGTGTCGCCCAGCCCCTGCGCAACCTGCGCGAGCATACTGGAGGCATCTCTGGCGCTGCGGGCCGACACCACCACACTTTGCACAGCCTGCttgt                                                                                                                                         <  1:316746/99‑1 (MQ=255)                                                     caggagCTGGGTGTCGCCCAGCCCCTGCGCAACCTGCGCGAGCATACTGGAGGCATCTCTGGCGCTGCGGGCCGACACCACCACACTTTGCACAGCCTGCttgt                                                                                                                                         >  2:316746/6‑104 (MQ=255)                                                     cTCGGGCTGGGCGACACCCAGCTCCTGCGCAACCTGCGCGAGCATACTGGAGGCATCTCTGGCGCTGCGGGCCGACACCACCACACTTTGCACAGCCTGCTTGTTGGTGCTGGCCACCAGGGCCTGGC‑GCAGCAGCGACTTGCCGCTGc                                                                                           >  4:210241/1‑149 (MQ=255)                                                       cGGGCTGGGCGACACCCAGCTCCTGCGCAACCTGCGCGAGCATACTGGAGGCATCTCTGGCGCTGCGGGCCGACACCACCACACTTTGCACAGCCTGCTTGTTGGTGCTGGCCACCAGGGCCTGGC‑GCAGCAGCGACTTGCCGCTGcc                                                                                          >  8:81735/1‑148 (MQ=255)                                                          gCTGGGCGACACCCAGCTCCTGCGCAACCTGCGCGAGCATACTGGAGGCATCTCTGGCGCTGCGGGCCGACACCACCACACTTTGCACAGCCTGCTTGTTGGTGCTGGCCACCAGGGCCTGGC‑GCAGCAGCGACTTGCCGCTGCCCAAt                                                                                      <  1:164818/149‑1 (MQ=255)                                                          gCTGGGCGACACCCAGCTCCTGCGCAACCTGCGCGAGCATACTGGAGGCATCTCTGGCGCTGCGGGCCGACACCACCACACTTTGCACAGCCTGCTTGTTGGTGCTGGCCACCAGGGCCTGGC‑GCAGCAGCGACTTGCCGCTGCCCAAt                                                                                      <  6:273947/149‑1 (MQ=255)                                                              ggCGACACCCAGCTCCTGCGCAACCTGCGCGAGCATACTGGAGGCATCTCTGGCGCTGCGGGCCGACACCACCACACTTTGCACAGCCTGCTTGTTGGTGCTGGCCACCAGGGCCTGGCAGCAGCAGCGACTTGCCGCTGCCCAATg                                                                                     <  1:172499/147‑1 (MQ=255)                                                               gCGACACCCAGCTCCTGCGCAACCTGCGCGAGCATACTGGAGGCATCTCTGGCGCTGCGGGCCGACACCACCACACTTTGCACAGCCTGCTTGTTGGTGCTGGCCACCAGGGc                                                                                                                      >  4:228498/1‑113 (MQ=255)                                                               gCGACACCCAGCTCCTGCGCAACCTGCGCGAGCATACTGGAGGCATCTCTGGCGCTGCGGGCCGACACCACCACACTTTGCACAGCCTGCTTGTTGGTGCTGGCCACCAGGGc                                                                                                                      <  3:228498/113‑1 (MQ=255)                                                                     cccAGCGCGTGCGCAAGGTGGGCGAGCATTCTGGAGGCATCTCTGACTCTGCGGGCCGACAGCACCACACTTTGCACAGCCTGCTTGTTGGTGCTGGCCACCAGGGCCTGGC‑GCAGCAGCGACTTGCCGCTGCCCATTGGGcgtgtgac                                                                           <  2:247393/149‑8 (MQ=255)                                                                      ccAGCTCCTGCGCAACCTGCGCGAGCATACTGGAGGCATCTCTGGCGCTGCGGGCCGACACCACCACACTTTGCACAGCCTGCTTGTTGGTGCTGGCCACCAGGGCCTGGC‑GCAGCAGCGACTTGCCGCTGCCCAATGGGCCACTGac                                                                           >  2:75155/1‑148 (MQ=255)                                                                        aGCTCCTGCGCAACCTGCGCGAGCATACTGGAGGCATCTCTGGCGCTGCGGGCCGACACCACCACACTTTGCACAGCCTGCTTGTTGGTGCTGGCCACCAGGGCCTGGC‑GCAGCAGCGACTTGCCGCTGCCCAATGGGc                                                                                  >  1:124558/1‑139 (MQ=255)                                                                        aGCTCCTGCGCAACCTGCGCGAGCATACTGGAGGCATCTCTGGCGCTGCGGGCCGACACCACCACACTTTGCACAGCCTGCTTGTTGGTGCTGGCCACCAGGGCCTGGC‑GCAGCAGCGACTTGCCGCTGCCCAATGGGc                                                                                  <  2:124558/139‑1 (MQ=255)                                                                             cTGCGCAACCTGCGCGAGCATACTGGAGGCATCTCTGGCGCTGCGGGCCGACACCACCACACTTTGCACAGCCTGCTTGTTGGTGCTGGCCACCAGGGCCTGGC‑GCAGCAGCGACTTGCCg                                                                                               <  5:286998/121‑1 (MQ=255)                                                                             cTGCGCAACCTGCGCGAGCATACTGGAGGCATCTCTGGCGCTGCGGGCCGACACCACCACACTTTGCACAGCCTGCTTGTTGGTGCTGGCCACCAGGGCCTGGC‑GCAGCAGCGACTTGCCg                                                                                               >  6:286998/1‑121 (MQ=255)                                                                               gcgcAACCTGCGCGAGCATACTGGAGGCATCTCTGGCGCTGCGGGCCGACACCACCACACTTTGCACAGCCTGCTTGTTGGTGCTGGCCACCAGGGCCTGGC‑GCAGCAGCGACTTGCCGCTGCCCAATGGGCCACTGACCACCAGCAt                                                                  >  8:279497/1‑148 (MQ=255)                                                                                   aaCCTGCGCGAGCATACTGGAGGCATCTCTGGCGCTGCGGGCCGACACCACCACACTTTGCACAGCCTGCTTGTTGGTGCTGGCCACCAGGGCCTGGC‑GCAGCAGCGACTTGCCGCTGCCCAATGGGCCACTGACCACCAGCATCAACt                                                             <  2:78231/149‑1 (MQ=255)                                                                                     ccTGCGCGAGCATACTGGAGGCATCTCTGGCGCTGCGGGCCGACACCACCACACTTTGTACAGCCTGCTTGTTGGTGCTGGCCACCAGGGCCTGGC‑GCAGCAGCGACTTGCCGCTGCCCAATGGGCCACTGAccacca                                                                      <  6:110345/138‑1 (MQ=255)                                                                                     ccTGCGCGAGCATACTGGAGGCATCTCTGGCGCTGCGGGCCGACACCACCACACTTTGTACAGCCTGCTTGTTGGTGCTGGCCACCAGGGCCTGGC‑GCAGCAGCGACTTGCCGCTGCCCAATGGGCCACTGAccacca                                                                      >  5:110345/1‑138 (MQ=255)                                                                                     ccTGCGCGAGCATACTGGAGGCATCTCTGGCGCTGCGGGCCGACACCACCACACTTTGCACAGCCTGCTTGTTGGTGCTGGCCACCAGGGCCTGGC‑GCAGCAGCGACTTGCCGCTGCCCAATGGGCCACTGACCACCAGCATCAACTgg                                                           <  8:297746/149‑1 (MQ=255)                                                                                         cgcgAGCATACTGGAGGCATCTCTGGCGCTGCGGGCCGACACCACCACACTTTGCACAGCCTGCTTGTTGGTGCTGGCCACCAGGGCCTGGC‑GCAGCAGCGACTTGCCGCTGCCCAATGGGCCACTGACCACCAGCATCAACTGGCta                                                        <  3:290986/148‑1 (MQ=255)                                                                                            gAGCATACTGGAGGCATCTCTGGCGCTGCGGGCCGACACCACCACACTTTGCACAGCCTGCTTGTTGGTGCTGGCCACCAGGGCCTGGC‑GCAGCAGCGACTTGCCGCTGCCCAATGGGCCACTGACCACCAGCATCAACTGGCTATAgc                                                    >  6:57485/1‑149 (MQ=255)                                                                                            gAGCATACTGGAGGCATCTCTGGCGCTGCGGGCCGACACCACCACACTTTGCACAGCCTGCTTGTTGGTGCTGGCCACCAGGGCCTGGC‑GCAGCAGCGACTTGCCGCTGCCCAATGGGCCACTGACCACCAGCATCAACTGGCTATAgc                                                    >  6:332128/1‑149 (MQ=255)                                                                                            gAGCATACTGGAGGCATCTCTGGCGCTGCGGGCCGACACCACCACACTTTGCACAGCCTGCTTGTTGGTGCTGGCCACCAGGGCCTGGC‑GCAGCAGCGACTTGCCGCTGCCCAATGGGCCACTGACCACCAGCATCAACTGGCTATAgc                                                    <  2:37525/149‑1 (MQ=255)                                                                                                       aGGCATCTCTGGCGCTGCGGGCCGACACCACCACACTTTGCACAGCCTGCTTGTTGGTGCTGGCCACCAGGGCCTGGC‑GCAGCAGCGACTTGCCGCTGCCCAATGGGCCACTGACCACCAGCATCAACTGGCTATAGCGCGCCAggtgg                                         >  5:141562/1‑149 (MQ=255)                                                                                                        ggCATCTCTGGCGCTGCGGGCCGACACCACCACACTTTGCACAGCCTGCTTGTTGGTGCTGGCCACCAGGGCCCGGC‑GCAGCAGCGACTTGCCGCTGCCCAATGGGCCACTGACCACCAGCATCAACTGGCTATAGCGCGCCAggtggt                                        >  8:140089/1‑149 (MQ=255)                                                                                                         gCATCTCTGGCGCTGCGGGCCGACACCACCACACTTTGTACAGCCTGCTTGTTGGTGCTGGCCACCAGGGCCTGGC‑GCAGCAGCGACTTGCCGCTGCCCAATGGGCCAc                                                                               >  7:127757/1‑109 (MQ=255)                                                                                                         gCATCTCTGGCGCTGCGGGCCGACACCACCACACTTTGTACAGCCTGCTTGTTGGTGCTGGCCACCAGGGCCTGGC‑GCAGCAGCGACTTGCCGCTGCCCAATGGGCCAc                                                                               <  8:127757/109‑1 (MQ=255)                                                                                                           aTCTCTGGCGCTGCGGGCCGACACCACCACACTTTGCACAGCCTGCTTGTTGGTGCTGGCCACCAGGGCCTGGC‑GCAGCAGCGACTTGCCGCTGCCCAATGGGCCACTGACCACCAGCATCAACTGGCTATAGCGCGCCa                                              <  3:249435/140‑1 (MQ=255)                                                                                                           aTCTCTGGCGCTGCGGGCCGACACCACCACACTTTGCACAGCCTGCTTGTTGGTGCTGGCCACCAGGGCCTGGC‑GCAGCAGCGACTTGCCGCTGCCAAATGGGCCACGGACCAGCCGCATCAACTGGCGATAGCGCGGCAgagcgcg                                       >  4:249435/1‑141 (MQ=255)                                                                                                                  gcgcTGCGGGCCGACACCACCACACTTTGCACAGCCTGCTTGTTGGTGCTGGCCACCAGGGCCTGGC‑GCAGCAGCGACTTGCCGCTGCCCAATGGGCCACTGACCACCAGCATCAACTGGCTATAGCGCGCCAGGTGGTGCAGTTgg                                <  3:210241/147‑1 (MQ=255)                                                                                                                        cGGGCCGACACCACCACACTTTGCACAGCCTGCTTGTTGGTGCTGGCCACCAGGGCCCGGC‑GCAGCAGCGACTTGCCGCTGCCCAATGGGCCACTGACCACCAGCATCAACTGGCTATAGCGCGCCAGGTGGTGCAGTTGGCCAAGc                          <  7:140089/147‑1 (MQ=255)                                                                                                                         gggCCAACACCACCACACTTTGCACAGCCTGCTTGTTGGTGCTGGCCACCAGGGCCTGGC‑GCAGCAGCGACTTGCCGCTGCCCAATGGGCCACTGACCACCAGCATCAACTGGCTATAGCGCGCCAGGTGGTGCAGTTGGCCAAGCACg                       >  7:286350/1‑149 (MQ=255)                                                                                                                          ggCCGACACCACCACACTTTGCACAGCCTGCTTGTTGGTGCTGGCCACCAGGGCCTGGC‑GCAGCAGCGACTTGCCGCTGCCCAATGGGCCACTGACCACCAGCATCAACTGGCTATAGCGCGCCAGGTGGTGCAGTTGGCCAAGCACgg                      >  1:248805/1‑149 (MQ=255)                                                                                                                            ccGACACCACCACACTTTGCACAGCCTGCTTGTTGGTGCTGGCCACCAGGGCCTGGC‑GCAGCAGCGACTTGCCGCTGCCCAATGGGCCACTGACCACCAGCATCAACTGGCTATAGCGCGCCAGGTGGTGCAGTTGGCCAAGCACGGGc                    <  1:139318/149‑1 (MQ=255)                                                                                                                               acacCACCACACTTTGCACAGCCTGCTTGTTGGTGCTGGCCACCAGGGCCTGGC‑GCAGCAGCGACTTGCCGCTGCCCAATGGGCCACTGACCACCAGCATCAACTGGCTATAGCGCGCCAGGTGGTGCAGTTGGCCAAGCACGGGCTTg                 >  5:286042/1‑149 (MQ=255)                                                                                                                                 accaccacACTTTGCACAGCCTGCTTGTTGGTGCTGGCCACCAGGGCCTGGC‑GCAGCAGCGACTTGCCGCTGCCCAATGGGCCACTGACCACCAGCATCAACTGGCTATAGCGCGCCAGGTGGTGCAGTTGGCCAAGCACGGGCTTgcg               >  2:202976/1‑149 (MQ=255)                                                                                                                                             tGCACAGCCTGCTTGTTGGTGCTGGCCACCAGGGCCTGGC‑GCAGCAGCGACTTGCCGCTGCCCAATGGGCCACTGACCACCAGCATCAACTGGCTATAGCGCGCCAGGTGGTGCAGTTGGCCAAGCACGGGCTTGCGCTGGGCCggg     <  7:81735/147‑1 (MQ=255)                                                                                                                                               cacaGCCTGCTTGTTGGTGCTGGCCACCAGGGCCTGGC‑GCAGCAGCGACTTGCCGCTGCCCAATGGGCCACTGACCACCAGCATCAACTGGCTATAGCGCGCCAGGTGGTGCAGTTGGCCAAGCACGGGCTTGCGCTGGGCCGGgaag  >  6:87857/1‑148 (MQ=255)                                                                                                                                               |                                                                                                                                                     CTTCCTGCCCGGTCAGTGCCAGCTGCACCACCTTCGACAGGATCGCCTGCACCTCGGGCTGGGCGACACCCAGCTCCTGCGCAACCTGCGCGAGCATACTGGAGGCATCTCTGGCGCTGCGGGCCGACACCACCACACTTTGCACAGCCTGCTTGTTGGTGCTGGCCACCAGGGCCTGGC‑GCAGCAGCGACTTGCCGCTGCCCAATGGGCCACTGACCACCAGCATCAACTGGCTATAGCGCGCCAGGTGGTGCAGTTGGCCAAGCACGGGCTTGCGCTGGGCCGGGAAG  >  NC_002947/5799484‑5799773

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