Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	4,917,689	A→G	10.6%	I30V (ATC→GTC) ‡	PP_4328 →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	4,917,689	0	A	G	10.6%	123.6 / 5.4	47	I30V (ATC→GTC) ‡	PP_4328	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (21/21);  new base G (3/2);  total (24/23)
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.95e-01

ATTATTACACGCGCCGCCGCGCTGCCAAGAGGCCGGGCTCGACAGGTTGTATACAATCATGACTGAAATCAATAGTCTCGGCGCACAAACCGCGATCAACAGCCAAGCCATGAAGGCGGGCATGACCGGCGATCTGCTGCGCCTGATCCAGGCGCAGCCTGGGTTGCTAGCTCCTGGCGAAACCGCGCAAGCCGAGGTTTTATCGTTGCGTCAGACGGGCCAGGACTTTCAGCTGGTGCTACGGCTGATGCAGGCCAACGGCGTGCAGACCCAGCTTCAGGCCAACGC  >  NC_002947/4917544‑4917831                                                                                                                                                  |                                                                                                                                               attattACACGCGCCGCCGCGCTGCCAAGAGGCCGGGCTCGACAGGTTGTATACAATCATGACTGAAATCAATAGTCTCGGCGCACAAACCGCGATCAACAGCCAAGCCATGAAGGCGGGCATGACCGGCGATCTGCTGCGCCTGATc                                                                                                                                              >  6:111305/1‑148 (MQ=255)              ccgccgcGCTGCCAAGAGGCCGGGCTCGACAGGTTGTATACAATCATGACTGAAATCAATAGTCTCGGCGCACAAACCGCGATCAACAGCCAAGCCATGAAGGCGGGCATGACCGGCGATCTGCTGCGCCTGATCCAGGCGCAGCCTgg                                                                                                                                <  4:263671/149‑1 (MQ=255)                 ccgcGCTGCCAAGAGGCCGGGCTCGACAGGTTGTATACAATCATGACTGAAATCAATAGTCTCGGCGCACAAACCGCGATCAACAGCCAAGCCATGAAGGCGGGCATGACCGGCGATCTGCTGCGCCTGATCCAGGCGCAGCCTGGGtt                                                                                                                             <  2:149161/149‑1 (MQ=255)                   gcgcTGCCAAGAGGCCGGGCTCGACAGGTTGTATACAATCATGACTGAAATCAATAGTCTCGGCGCACAAACCGCGATCAACAGCCAAGCCATGAAGGCGGGCATGACCGGCGATCTGCTGCGCCTGGATCAGGCGCAg                                                                                                                                     <  2:252735/139‑1 (MQ=255)                   gcgcTGCCAAGAGGCCGGGCTCGACAGGTTGTATACAATCATGACTGAAATCAATAGTCTCGGCGCACAAACCGCGATCAACAGCCAAGCCATGAAGGCGGGCATGACCGGCGATCTGCTGCGCCTGGATCAGGCGCAg                                                                                                                                     >  1:252735/1‑139 (MQ=255)                    cgcTGCCAAGAGGCCGGGCTCGACAGGTTGTATACAATCATGACTGAAATCAATAGTCTCGGCGCACAAACCGCGATCAACAGCCAAGCCATGAAGGCGGGCATGACCGGCGATCTGCTGCGCCTGATCCAGGCGCAGCCTGGGTTGCt                                                                                                                          <  2:239264/149‑1 (MQ=255)                    cgcTGCCAAGAGGCCGGGCTCGACAGGTTGTATACAATCATGACTGAAATCAATAGTCTCGGCGCACAAACCGCGATCAACAGCCAAGCCATGAAGGCGGGCATGACCGGCGATCTGCTGCGCCTGATCCAGGCGCAGCCTGGGTTGCt                                                                                                                          <  5:111305/149‑1 (MQ=255)                           aaGAGGCCGGGCTCGACAGGTTGTATACAATCATGACTGAAATCAATAGTCTCGGCGCACAAACCGCGATCAACAGCCAAGCCATGAAGGCGGGCATGACCGGCGATCTGCTGCGCCTGATCCAGGCGCAGCCTGGGTTGCTAGCTCCt                                                                                                                   <  5:76840/149‑1 (MQ=255)                           aaGAGGCCGGGCTCGACAGGTTGTATACAATCATGACTGAAATCAATAGTCTCGGCGCACAAACCGCGATCAACAGCCAAGCCATGAAGGCGGGCATGACCGGCGATCTGCTGCGCCTGATCCAGGCGCAGCCTGGGTTGCTAGCTCCt                                                                                                                   >  5:249777/1‑149 (MQ=255)                            agagGCCGGGCTCGACAGGTTGTATACAATCATGACTGAAATCAATAGTCTCGGCGCACAAACCGCGATCAACAGCCAAGCCATGAAGGCGGGCATGACCGGCGATCTGCTGCGCCTGATCCAGGCGCAGCCTGGGTTGCTAGCTCCTg                                                                                                                  <  8:112687/149‑1 (MQ=255)                                gCCGGGCTCGACAGGTTGTATACAATCATGACTGAAATCAATAGTCTCGGCGCACAAACCGCGATCAACAGCCAAGCCATGAAGGCGGGCATGACCGGCGATCTGCTGCGCCTGATCCAGGCGCAGCCTGGGTTGCt                                                                                                                          >  1:209830/1‑137 (MQ=255)                                gCCGGGCTCGACAGGTTGTATACAATCATGACTGAAATCAATAGTCTCGGCGCACAAACCGCGATCAACAGCCAAGCCATGAAGGCGGGCATGACCGGCGATCTGCTGCGCCTGATCCAGGCGCAGCCTGGGTTGCt                                                                                                                          <  2:209830/137‑1 (MQ=255)                                gCCGGGCTCGACAGGTTGTATACAATCATGACTGAAATCAATAGTCTCGGCGCACAAACCGCGATCAACAGCCAAGCCATGAAGGCGGGCATGACCGGCGATCTGCTGCGCCTGATCCAGGCGCAGCCTGGGTTGCTAGCTCCTGGCg                                                                                                               >  6:183171/1‑148 (MQ=255)                                gCCGGGCTCGACAGGTTGTATACAATCATGACTGAAATCAATAGTCTCGGCGCACAAACCGCGATCAACAGCCAAGCCATGAAGGCGGGCATGACCGGCGATCTGCTGCGCCTGATCCAGGCGCAGCCTGGGTTGCTAGCTCCTGGCg                                                                                                               >  7:2351/1‑148 (MQ=255)                                     gCTCGACAGGTTGTATACAATCATGACTGAAATCAATAGTCTCGGCGCACAAACCGCGATCAACAGCCAAGCCATGAAGGCGGGCATGACCGGCGATCTGCTGCGCCTGATCCAGGCGCAGCCTGGGTTGCTAGCTCCTGGCGAAACcg                                                                                                         <  8:2351/149‑1 (MQ=255)                                                  tataCAATCATGACTGAAATCAATAGTCTCGGCGCACAAACCGCGATCAACAGCCAAGCCATGAAGGCGGGCATGACCGGCGATCTGCTGCGCCTGATCCAGGCGCAGCCTGGGTTGCTAGCTCCTGGCGAAACCGCGCAAGCCGAgg                                                                                             <  5:183171/148‑1 (MQ=255)                                                    taCAATCATGACTGAAATCAATAGTCTCGGCGCACAAACCGCGATCAACAGCCAAGCCATGAAGGCGGGCATGACCGGCGATCTGCTGCGCCTGATCCAGGCGCAGCCTGGGTTGCTAGCTCCTGGCGAAACCGCGCAAGCCGAGGtt                                                                                           <  3:98008/148‑1 (MQ=255)                                                     aCAATCATGACTGAAATCAATAGTCTCGGCGCACAAACCGCGATCAACAGCCAAGCCATGAAGGCGGGCATGACCGGCGATCTGCTGCGCCTGATCCAGGCGCAGACTGGGTTGCTAGCTCCTGGCGAAACCGCGCAAGCCGAGGttt                                                                                          >  1:138023/1‑148 (MQ=255)                                                      taaTCATGACTGAAATCAATAGTCTCGGCGCACAAACCGCGATCAACAGCCAAGCCATGAAGGCAGGCATGACCGGCGATCTGCTGCGCCTGATCCAGGCGCAGCCTGGGTTGCTAGCTCCTGGCGAAACCGCGCAAGCCGAGGTTTTa                                                                                        <  5:115565/148‑1 (MQ=255)                                                        aTCATGACTGAAATCAATAGTCTCGGCGCACAAACCGCGATCAACAGCCAAGCCATGAAGGCAGGCATGACCGGCGATCTGCTGCGCCTGATCCAGGCGCAGCCTGGGTTGCTAGCTCCTGGCGAAACCGCGCAAGCCGAGGTTTTa                                                                                        >  6:115565/1‑147 (MQ=255)                                                         tCATGACTGAAATCAATAGTCTCGGCGCACAAACCGCGATCAACAGCCAAGCCATGAAGGCGGGCATGACCGGCGATCTGCTGCGCCTGATCCAGGCGCAGCCTGGGTTGCTAGCTCCTGGCGAAACCGCGCAAGCCGAGGTTTTATCg                                                                                     >  5:4045/1‑149 (MQ=255)                                                           aTGACTGAAATCAATAGTCTCGGCGCACAAACCGCGATCAACAGCCAAGCCATGAAGGCGGGCATGACCGGCGATCTGCTGCGCCTGATCCAGGCGCAGCCTGGGTTGCTAGCTCCTGGCGAAACCGCGCAAGCCGAGGTTTTATCGtt                                                                                   >  1:258050/1‑149 (MQ=255)                                                                         tAGTCTCGGCGCACAAACCGCGATCAACAGCCAAGCCATGAAGGCGGGCATGACCGGCGATCTGCTGCGCCTGATCCAGGCGCAGCCTGGGTTGCTAGCTCCTGGCGAAACCGCGCAAGCCGAGGTTTTATCGTTGCGTCAGACGGGc                                                                      >  5:180553/1‑148 (MQ=255)                                                                              tcGGCGCACAAACCGCGATCAACAGCCAAGCCATGAAGGCGGGCATGACCGGCGATCTGCTGCGCCTGGATCAGGCGCa                                                                                                                                      <  2:163890/79‑1 (MQ=39)                                                                              tcGGCGCACAAACCGCGATCAACAGCCAAGCCATGAAGGCGGGCATGACCGGCGATCTGCTGCGCCTGGATCAGGCGCa                                                                                                                                      >  1:163890/1‑79 (MQ=39)                                                                                ggCGCACAAACCGCGATCAACAGCCAAGCCATGAAGGCGGGCATGACCGGCGATCTGCTGCGCCTGATCCAGGCGCAGCCTGGGTTGCTAGCTCCTGGCGAAACCGCGCAAGCCGAGGTTTTATCGTTGCGTCAGACGGGCCAGGAc                                                                <  6:4045/147‑1 (MQ=255)                                                                                    cacaAACCGCGATCAACAGCCAAGCCATGAAGGCGGGCATGACCGGCGATCTGCTGCGCCTGATCCAGGCGCAGCCTGGGTTGCTAGCTCCTGGCGAAACCGCGCAAGCCGAGGTTTTATCGTTGCGTCAGACGGGCCAGGACttt                                                             >  2:72581/1‑146 (MQ=255)                                                                                    cacaAACCGCGATCAACAGCCAAGCCATGAAGGCGGGCATGACCGGCGATCTGCTGCGCCTGATCCAGGCGCAGCCTGGGTTGCTAGCTCCTGGCGAAACCGCGCAAGCCGAGGTTTTATCGTTGCGTCAGACGGGCCAGGACttt                                                             <  1:72581/146‑1 (MQ=255)                                                                                     acaAACCGCGATCAACAGCCAAGCCATGAAGGCGGGCATGACCGCCGATCTGCTGCGCCTGATACAGGCGCAGCCTGGGTTGCTAGCTCCTGGCGAAACCGCGCAAGCCGAGGTTTTATCGTTGCGTCAGACGCGCCAGGACTTTCAGc                                                         >  5:237367/1‑149 (MQ=255)                                                                                       aaaCCGCGATCAACAGCCAAGCCATGAAGGCGGGCATGACCGGCGATCTGCTGCGCCTGATCCAGGCGCAGCCTGGGTTGCTAGCTCCTGGCGAAACCGCGCAAGCCGAGGTTTTATCGTTGCGTCAGACGGGCCAGGACTTTCAGCTg                                                       >  5:14679/1‑149 (MQ=255)                                                                                         aCCGCGATCATCAGCCAAGCCATGAAGGCGGGCTTGACCGGGGATCTGCTGCGCCTGATGCAGGGGCAGCCTGGGTTGCTAGCTCCTGGCGAAACCGCGCAAGCCGAGGTTTTATCGTTGCGTCAGACGGGCCAGGACTTTCAGCTGGt                                                     <  6:237367/149‑1 (MQ=255)                                                                                         aCCGCGATCAACAGCCAAGCCATGAAGGCGGGCATGACCGGCGATCTGCTGCGCCTGATCCAGGCGCAGCCTGGGTTGCTAGCTCCTGGCGAAACCGCGCAAGCCGAGGTTTTATCGTTGCGTCAGACGGGCCAGGACTTTCAGCTGGt                                                     <  1:34285/149‑1 (MQ=255)                                                                                          ccGCGATCAACAGCCAAGCCATGAAGGCGGGCATGACCGGCGATCTGCTGCGCCTGATCCAGGCGCAGCCTGGGTTGCTAGCTCCTGGCGAAACCGCGCAAGCCGAGGTTTTATCGTTGCGTCAGACGGGCCAGGAc                                                                <  4:196791/137‑1 (MQ=255)                                                                                          ccGCGATCAACAGCCAAGCCATGAAGGCGGGCATGACCGGCGATCTGCTGCGCCTGATCCAGGCGCAGCCTGGGTTGCTAGCTCCTGGCGAAACCGCGCAAGCCGAGGTTTTATCGTTGCGTCAGACGGGCCAGGAc                                                                >  3:196791/1‑137 (MQ=255)                                                                                              gATCAACAGCCAAGCCATGAAGGCGGGCATGACCGGCGATCTGCTGCGCCTGATCCAGGCGCAGCCTGGGTTGCTAGCTCCTGGCGAAACCGCGCAAGCCGAGGTTTTATCGTTGCGTCAGACGGGCCAGGACTTTCAGCTGGTGCTAc                                                <  6:249777/149‑1 (MQ=255)                                                                                                  aaCAGCCAAGCCATGAAGGCGGGCATGACCGGCGATCTGCTGCGCCTGATCCAGGCGCAGCCTGGGTTGCTAGCTCCTGGCGAAACCGCGCAAGCCGAGGTTTTATCGTTGCGTCAGACGGGCCAGGACTTTCAGCTGGTGCTACGGc                                             <  2:258050/148‑1 (MQ=255)                                                                                                  aaCAGCCAAGCCATGAAGGCGGGCATGACCGGCGATCTGCTGCGCCTGATCCAGGCGCAGCCTGGGTTGCTAGCTCCTGGCGAAACCGCGCAAGCCGAGGTTTTATCGTTGCGTCAGACGGGCCAGGACTTTCAGCTGGTGCTACGGCt                                            >  1:131957/1‑149 (MQ=255)                                                                                                      gccaagccaTGAAGGCGGGCATGACCGGCGATCTGCTGCGCCTGATCCAGGCGCAGCCTGGGTTGCTAGCTCCTGGCGAAACCGCGCAAGCCGAGGTTTTATCGTTGCGTCAGACGGGCCAGGACTTTCAGCTGGTGCTACGGCTGATg                                        >  7:83268/1‑149 (MQ=255)                                                                                                      gccaagccaTGAAGGCGGGCATGACCGGCGATCTGCTGCGCCTGATCCAGGCGCAGCCTGGGTTGCTAGCTCCTGGCGAAACCGCGCAAGCCGAGGTTTTATCGTTGCGTCAGACGGGCCAGGACTTTCAGCTGGTGCTACGGCTGATg                                        >  3:34049/1‑149 (MQ=255)                                                                                                         aagccaTGAAGGCGGGCATGACCGGCGATCTGCTGCGCCTGATCCAGGCGCAGCCTGGGTTGCTAGCTCCTGGCGAAACCGCGCAAGCCGAGGTTTTATCGTTGCGTCAGACGGGCCAGGACTTTCAGCTGGTGCTACGGCTGATGCAg                                     <  8:83268/149‑1 (MQ=255)                                                                                                         aagccaTGAAGGCGGGCATGACCGGCGATCTGCTGCGCCTGATCCAGGCGCAGCCTGGGTTGCTAGCTCCTGGCGAAACCGCGCAAGCCGAGGTTTTATCGTTGCGTCAGACGGGCCAGGACTTTCAGCTGGTGCTACGGCTGATGCAg                                     <  4:34049/149‑1 (MQ=255)                                                                                                           gccaTGAAGGCGGGCATGACCGGCGATCTGCTGCGCCTGATCCAGGCGCAGCCTGGGTTGCTAGCTCCTGGCGAAACCGCGCAAGCCGAGGTTTTATCGTTGCGTCAGACGGGCCAGGACTTTCAGCTGGTGCTACGGCTGATGCAGGc                                   >  6:144221/1‑149 (MQ=255)                                                                                                                             aCCGGCGATCTGCTGCGCCTGATCCAGGCGCAGCCTGGGTTGCTAGCTCCTGGCGAAACCGCGCAAGCCGAGGTTTTATCGTTGCGTCAGACGGGCCAGGACTTTCAGCTGGTGCTACGGCTGATGCAGGCCAACGGCGTGCAGAccc                  <  2:138023/148‑1 (MQ=255)                                                                                                                                ggCGATCTGCTGCGCCTGATCCAGGCGCAGCCTGGGTTGCTAGCTCCTGGCGAAACCGCGCAAGCCGAGGTTTTATCGTTGCGTCAGACGGGCCAGGACTTTCAGCTGGTGCTACGGCTGATGCAGGCCAACGGCGTGCAGACCCAGCt              >  7:75594/1‑149 (MQ=255)                                                                                                                                        gctgcGCCTGGATCAGGCGCAGCCTGGGTTGCTAGCTCCTGGCGAAACCGCGCAAGCCGAGGTTTTATCGTTGCGTCAGACGGGCCAGGACTTTCAGCTGGTGCTACGGCTGATGCAGGCCAACGGCGTGCAGACCCAGCTTCAGGcc       >  6:129359/1‑148 (MQ=255)                                                                                                                                            cgcCTGATCCAGGCGCAGCCTGGGTTGCTAGCTCCTGGCGAAACCGCGCAAGCCGAGGTTTTATCGTTGCGTCAGACGGGCCAGGACTTTCAGCTGGTGCTACGGCTGATGCAGGCCAACGGCGTGCAGACCCAGCTTCAGGCCAACgc  >  4:100735/1‑149 (MQ=255)                                                                                                                                              ccTGATCCAGGCGCAGCCTGGGTTGCTAGCTCCTGGCGAAACCGCGCAAGCCGAGGTTTTATCGTTGCGTCAGACGGGCCAGGACTTTCAGCTGGTGCTACGGCTGATGCAGGCCAACGGCGTGCAGACCCAGCTTCAGGCCAACgc  >  4:249419/1‑147 (MQ=255)                                                                                                                                                  |                                                                                                                                               ATTATTACACGCGCCGCCGCGCTGCCAAGAGGCCGGGCTCGACAGGTTGTATACAATCATGACTGAAATCAATAGTCTCGGCGCACAAACCGCGATCAACAGCCAAGCCATGAAGGCGGGCATGACCGGCGATCTGCTGCGCCTGATCCAGGCGCAGCCTGGGTTGCTAGCTCCTGGCGAAACCGCGCAAGCCGAGGTTTTATCGTTGCGTCAGACGGGCCAGGACTTTCAGCTGGTGCTACGGCTGATGCAGGCCAACGGCGTGCAGACCCAGCTTCAGGCCAACGC  >  NC_002947/4917544‑4917831

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