Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	4,917,691	C→T	11.3%	I30I (ATC→ATT) ‡	PP_4328 →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	4,917,691	0	C	T	11.3%	117.2 / 5.5	46	I30I (ATC→ATT) ‡	PP_4328	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (19/20);  new base T (3/2);  total (23/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.94e-01

ATTATTACACGCGCCGCCGCGCTGCCAAGAGGCCGGGCTCGACAGGTTGTATACAATCATGACTGAAATCAATAGTCTCGGCGCACAAACCGCGATCAACAGCCAAGCCATGAAGGCGGGCATGACCGGCGATCTGCTGCGCCTGATCCAGGCGCAGCCTGGGTTGCTAGCTCCTGGCGAAACCGCGCAAGCCGAGGTTTTATCGTTGCGTCAGACGGGCCAGGACTTTCAGCTGGTGCTACGGCTGATGCAGGCCAACGGCGTGCAGACCCAGCTTCAGGCCAACGCCAGCCAGC  >  NC_002947/4917544‑4917839                                                                                                                                                    |                                                                                                                                                     attattACACGCGCCGCCGCGCTGCCAAGAGGCCGGGCTCGACAGGTTGTATACAATCATGACTGAAATCAATAGTCTCGGCGCACAAACCGCGATCAACAGCCAAGCCATGAAGGCGGGCATGACCGGCGATCTGCTGCGCCTGATc                                                                                                                                                      >  6:111305/1‑148 (MQ=255)              ccgccgcGCTGCCAAGAGGCCGGGCTCGACAGGTTGTATACAATCATGACTGAAATCAATAGTCTCGGCGCACAAACCGCGATCAACAGCCAAGCCATGAAGGCGGGCATGACCGGCGATCTGCTGCGCCTGATCCAGGCGCAGCCTgg                                                                                                                                        <  4:263671/149‑1 (MQ=255)                 ccgcGCTGCCAAGAGGCCGGGCTCGACAGGTTGTATACAATCATGACTGAAATCAATAGTCTCGGCGCACAAACCGCGATCAACAGCCAAGCCATGAAGGCGGGCATGACCGGCGATCTGCTGCGCCTGATCCAGGCGCAGCCTGGGtt                                                                                                                                     <  2:149161/149‑1 (MQ=255)                   gcgcTGCCAAGAGGCCGGGCTCGACAGGTTGTATACAATCATGACTGAAATCAATAGTCTCGGCGCACAAACCGCGATCAACAGCCAAGCCATGAAGGCGGGCATGACCGGCGATCTGCTGCGCCTGGATCAGGCGCAg                                                                                                                                             >  1:252735/1‑139 (MQ=255)                   gcgcTGCCAAGAGGCCGGGCTCGACAGGTTGTATACAATCATGACTGAAATCAATAGTCTCGGCGCACAAACCGCGATCAACAGCCAAGCCATGAAGGCGGGCATGACCGGCGATCTGCTGCGCCTGGATCAGGCGCAg                                                                                                                                             <  2:252735/139‑1 (MQ=255)                    cgcTGCCAAGAGGCCGGGCTCGACAGGTTGTATACAATCATGACTGAAATCAATAGTCTCGGCGCACAAACCGCGATCAACAGCCAAGCCATGAAGGCGGGCATGACCGGCGATCTGCTGCGCCTGATCCAGGCGCAGCCTGGGTTGCt                                                                                                                                  <  5:111305/149‑1 (MQ=255)                    cgcTGCCAAGAGGCCGGGCTCGACAGGTTGTATACAATCATGACTGAAATCAATAGTCTCGGCGCACAAACCGCGATCAACAGCCAAGCCATGAAGGCGGGCATGACCGGCGATCTGCTGCGCCTGATCCAGGCGCAGCCTGGGTTGCt                                                                                                                                  <  2:239264/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                                                                                                                                              >  1:163890/1‑79 (MQ=39)                                                                              tcGGCGCACAAACCGCGATCAACAGCCAAGCCATGAAGGCGGGCATGACCGGCGATCTGCTGCGCCTGGATCAGGCGCa                                                                                                                                              <  2:163890/79‑1 (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)                                                                                                                                                    ccAGGCGCAGCCTGGGTTGCTAGCTCCTGGCGAAACCGCGCAAGCCGAGGTTTTATCGTTGCGTCAGACGGGCCAGGACTTTCAGCTGGTGCTACGGCTGATGCAGGCCAACGGCGTGCAGACCCAGCTTCAGGCCAACGCcagccagc  >  3:192728/1‑149 (MQ=255)                                                                                                                                                    |                                                                                                                                                     ATTATTACACGCGCCGCCGCGCTGCCAAGAGGCCGGGCTCGACAGGTTGTATACAATCATGACTGAAATCAATAGTCTCGGCGCACAAACCGCGATCAACAGCCAAGCCATGAAGGCGGGCATGACCGGCGATCTGCTGCGCCTGATCCAGGCGCAGCCTGGGTTGCTAGCTCCTGGCGAAACCGCGCAAGCCGAGGTTTTATCGTTGCGTCAGACGGGCCAGGACTTTCAGCTGGTGCTACGGCTGATGCAGGCCAACGGCGTGCAGACCCAGCTTCAGGCCAACGCCAGCCAGC  >  NC_002947/4917544‑4917839

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