Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,284,507	C→G	14.3%	S288R (AGC→AGG)	PP_4658 →	methyl‑accepting chemotaxis protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	5,284,507	0	C	G	14.3%	74.8 / 6.0	28	S288R (AGC→AGG)	PP_4658	methyl‑accepting chemotaxis protein
Reads supporting (aligned to +/- strand):  ref base C (13/11);  new base G (2/2);  total (15/13)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.98e-01

GTACGCCAGTTCCGCGACGCGGTTGCCAACATCACCACCTCGCGCGCCGAAATGACCGTGCAGGGCGCCGACATCGTCAAGCGCAGTGACGCGTTGTACCAGATCCAGCTGGAGCGCCGTGACATCGAGAGCACCCAGGCCCGCAGCCTGCAGGCCATCGCCACACTGTTGGCGCTGCTGGTCGGCGTGCTGGCGGCGGTGCTGATCACCCGCCAGATCACGCGCCCGTTGCAAGACACCCTGGTGGCCGTGGAGAAAATCGCCAGTGGCGACCTTACCCAACACATGCGC  >  NC_002947/5284361‑5284651                                                                                                                                                   |                                                                                                                                                 gTACGCCAGTTCCGCGACGCGGTTGCCAACATCACCACCTCGCGCGCCGAAATGACCGTGCAGGGCGCCGACATCGTCAAGCGCAGTGACGCGTTGTACCAGATCCAGCTGGAGCGCCGTGACATCGAGAGCACCCAGGCCTGCAGGCTg                                                                                                                                               >  1:510171/1‑150 (MQ=255)       cAGTTCCGCGACGCGGTTGCCAACATCACCACCTCGCGCGCCGAAATGACCGTGCAGGGCGCCGACATCGTCAAGCGCAGTGACGCGTTGTACCAGATCCAGCTGGAGCGCCGTGACATCGAGAGCACCCAGGCCCGCAGCCTGCAGGcc                                                                                                                                         >  1:174137/1‑150 (MQ=255)        aGTTCCGCGACGCGGTTGCCAACATCACCACCTCGCGCGCCGAAATGACCGTGCAGGGCGCCGACCTCTTCAAGCGCAGTGACGCGTTGTACCAGATCCAGCTGGAGCGCCGTGACATCGAGAGCACCCAGGCCTGCAGGCTGCGGGcct                                                                                                                                        <  1:422111/150‑2 (MQ=255)          ttCCGCGACGCGGTTGCCAACATCACCACCTCGCGCGCCGAAATGACCGTGCAGGGCGCCGACATCGTCAAGCGCAGTGACGCGTTGTACCAGATCCAGCTGGAGCGCCGTGACATCGAGAGCACCCAGGCCTGCAGGCTGCGGGcctgg                                                                                                                                      <  2:510171/150‑4 (MQ=255)            ccGCGACGCGGTTGCCAACATCACCACCTCGCGCGCCGAAATGACCGTGCAGGGCGCCGACATCGTCAAGCGCAGTGACGCGTTGTACCAGATCCAGCTGGAGCGCCGTGACATCGAGAGCACCCAGGCCCGCAGCCTGCAGGCCATCGc                                                                                                                                    <  1:615501/150‑1 (MQ=255)                               aTCACCACCTCGCGCGCCGAAATGACCGTGCAGGGCGCCGACATCGTCAAGCGCAGTGACGCGTTGTACCAGATCCAGCTGGAGCGCCGTGACATCGAGAGCACCCAGGCCCGCAGCCTGCAGGCCATCGCCACACTGTTGGCGctgctg                                                                                                                 >  2:382339/1‑150 (MQ=255)                                 caccacCTCGCGCGCCGAAATGACCGTGCAGGGCGCCGACATCGTCAAGCGCAGTGACGCGTTGTACCAGATCCAGCTGGAGCGCCGTGACATCGAGAGCACCCAGGCCCGCAGCCTGCAGGCCATCGCCACACTGTTGGCGCTGCTGGt                                                                                                               >  2:284885/1‑150 (MQ=255)                                             cgcCGAAATGACCGTGCAGGGCGCCGACATCGTCAAGCGCAGTGACGCGTTGTACCAGATCCAGCTGGAGCGCCGTGACATCGAGAGCACCCAGGCCCGCAGCCTGCAGGCCATCGCCACACTGTTGGCGCTGCTGGTCGGCGTGCTggc                                                                                                   >  2:703315/1‑150 (MQ=255)                                               ccGAAATGACCGTGCAGGGCGCCGACATCGTCAAGCGCAGTGACGCGTTGTACCAGATCCAGCTGGAGCGCCGTGACATCGAGAGCACCCAGGCCCGCAGCCTGCAGGCCATCGCCACACTGTTGGCGCTGCTGGTCGGCGTGCTggcgg                                                                                                 <  2:553024/150‑1 (MQ=255)                                                     tGACCGTGCAGGGCGCCGACATCGTCAAGCGCAGTGACGCGTTGTACCAGATCCAGCTGGAGCGCCGTGACATCGAGAGCACCCAGGCCCGCAGCCTGCAGGCCATCGCCACACTGTTGGCGCTGCTGGTCGGCGTGCTGGCGGCGGTGc                                                                                           <  1:284885/150‑1 (MQ=255)                                                                             tCAAGCGCAGTGACGCGTTGTACCAGATCCAGCTGGAGCGCCGTGACATCGAGAGCACCCAGAGCCGCAGCCCGCAGGCCATCGCCACACTGTTGACGCTGCTGGTCGGCGTGCTGGCGGCGGTGCTGATCACCCGCCCACTCACGCGcc                                                                   >  2:24221/1‑150 (MQ=255)                                                                                 gcgcAGTGACGCGTTGTACCAGATCCAGCTGGAGCGCCGTGACATCGAGAGCACCCAGGCCCGCAGCCTGCAGGCCATAGCCACACTGTTGGCGCTGCTGGTCGGCGTGCTGGCGGCGGTGCTGATCACCCGCCAGATCACGCGCCCGtt                                                               >  2:587682/1‑150 (MQ=255)                                                                                     aGTGACGCGTTGTACCAGATCCAGCTGGAGCGCCGTGACATCGAGAGCACCCAGGCCCGCGGCCTGCAGGCCATCGCCACACTGTTGGCGCTGCTGGTCGGCGTGCTGGCGGCGGTGCTGATCACCCGCCAGATCACGCGCCCGTTGCaa                                                           <  1:464647/150‑1 (MQ=255)                                                                                                      gATCCAGCTGGAGCGCCGTGACATCGAGAGCACCCAGGCCCGCAGCCTGCAGGCCATCGCCACACTGTTGGCGCTGCTGGTCGGCGTGCTGGCGGCGGTGCTGATCACCCGCCAGATCACGCGCCCGTTGCAAGACACCCTGGTGGCCGt                                          <  2:127252/150‑1 (MQ=255)                                                                                                              tGGAGCGCCGTGACATCGAGAGCACCCAGGCCCGCAGCCTGCAGGCCATCGCCACACTGTTGGCGCTGCTGGTCGGCGTGCTGGCGGCGGTGCTGATCACCCGCCAGATCACGCGCCCGTTGCAAGACCCCCTGGTGGCCGTGCAGcaaa                                  >  2:189730/1‑150 (MQ=255)                                                                                                                       gTGACATCGAGAGCACCCAGGCCCGCAGCCTGCAGGCCATCGCCACACTGTTGGCGCTGCTGGTCGGCGTGCTGGCGGCGGTGCTGATCACCCGCCAGATCACGCGCCCGTTGCAAGACACCCTGGTGGCCGTGGAGAAAATCGc                              <  2:448992/145‑1 (MQ=255)                                                                                                                       gTGACATCGAGAGCACCCAGGCCCGCAGCCTGCAGGCCATCGCCACACTGTTGGCGCTGCTGGTCGGCGTGCTGGCGGCGGTGCTGATCACCCGCCAGATCACGCGCCCGTTGCAAGACACCCTGGTGGCCGTGGAGAAAATCGc                              >  1:448992/1‑145 (MQ=255)                                                                                                                          aCATCGAGAGCACCCAGGCCCGCAGCCTGCAGGCCATCGCCACACTGTTGGCGCTGCTGGTCGGCGTGCTGGCGGCGGTATTGATCACACGCCAGATCACGCGCCCGTTGCAAGACACCCTGGTGGCCGTGGAGAAAATCGCCAGTgg                        >  1:660627/1‑148 (MQ=255)                                                                                                                          aCATCGAGAGCACCCAGGCCCGCAGCCTGCAGGCCATCGCCACACTGTTGGCGCTGCTGGTCGGCGTGCTGGCGGCGGTATTGATCACACGCCAGATCACGCGCCCGTTGCAAGACACCCTGGTGGCCGTGGAGAAAATCGCCAGTgg                        <  2:660627/148‑1 (MQ=255)                                                                                                                             tCGAGAGCACCCAGGCCCGCAGCCTGCAGGCCATCGCCACACTGTTGGCGCTGCTGGTCGGCGTGCTGGCGGCGGTGCTGATCACCCGCCAGATCACGCGCCCGTTGCAAGACACCCTGGTGGCCGTGGAGAAAATCGCCAGTGGCGAcc                   >  2:267196/1‑150 (MQ=255)                                                                                                                             tCGAGAGCACCCAGGCCCGCAGCCTGCAGGCCATCGCCACACTGTTGGCGCTGCTGGTCGGCGTGCTGGCGGCGGTGCTGATCACCCGCCAGATCACGCGCCCGTTGCAAGACACCCTGGTGGCCGTGGAGAAAATCGCCAGTGGCGAcc                   >  2:598370/1‑150 (MQ=255)                                                                                                                                     aCCCAGGCCCGCAGCCTGCAGGCCATCGCCACACTGTTGGCGCTGCTGGTCGGCGTGCTGGCGGCGGTGCTGATCACCCGCCAGATCACGCGCCCGTTGCAAGACACCCTGGTGGCCGTGGAGAAAATCGCCAGTGGCGACCTTACCCaa           <  1:703315/150‑1 (MQ=255)                                                                                                                                       tttgtcccGCAGCCTGCAGGCCATCGCCACACTGTTGGCGCTGCTGGTCGGCGTGCTGGCGGCGGTGCTGATCACCCGCCAGATCACGCGCCCGTTGCAAGACACCCTGGTGGCCGTGGAGAAAATCGCCAGTGGCGACCTTACCCAaca         >  1:215735/6‑150 (MQ=255)                                                                                                                                       ccAGGCCCGCAGCCTGCAGGCCATCGCCACACTGTTGGCGCTGCTGGTCGGCGTGCTGGCGGCGGTGCTGATCACCCGCCAGATCACGCGCCCGtt                                                               >  2:303110/1‑96 (MQ=255)                                                                                                                                       ccAGGCCCGCAGCCTGCAGGCCATCGCCACACTGTTGGCGCTGCTGGTCGGCGTGCTGGCGGCGGTGCTGATCACCCGCCAGATCACGCGCCCGtt                                                               <  1:303110/96‑1 (MQ=255)                                                                                                                                            cctgCAGGCTGCGGGCCATCGCCACACTGTTGGCGCTGCTGGTCGGCGTGCTGGCGGCGGTGCTGATCACCCGCCAGATCACGCGCCCGTTGCAAGACACCCTGGTGGCCGTGGAGAAAATCGCCAGTGGCGACCTTACCCAACACATgc    >  2:604271/4‑150 (MQ=255)                                                                                                                                            cccGCAGCCTGCAGGCCATCGCCACACTGTTGGCGCTGCTGGTCGGCGTGCTGGCGGCGGTGCTGATCACCCGCCAGATCACGCGCCCGTTGCAAGACACCCTGGTGGCCGTGGAGAAAATCGCCAGTGGCGACCTTACCCAACACATgc    <  1:598370/150‑1 (MQ=255)                                                                                                                                              cGCAGCCTGCAGGCCATCGCCACACTGTTGGCGCTGCTGGTCGGCGTGCTGGCGGCGGTGCTGATCACCCGCCAGATCACGCGCCCGTTGCAAGACACCCTGGTGGCCGTGGAGAAAATCGCCAGTGGCGACCTTACCCAACACATgcgc  <  1:382339/150‑1 (MQ=255)                                                                                                                                                   |                                                                                                                                                 GTACGCCAGTTCCGCGACGCGGTTGCCAACATCACCACCTCGCGCGCCGAAATGACCGTGCAGGGCGCCGACATCGTCAAGCGCAGTGACGCGTTGTACCAGATCCAGCTGGAGCGCCGTGACATCGAGAGCACCCAGGCCCGCAGCCTGCAGGCCATCGCCACACTGTTGGCGCTGCTGGTCGGCGTGCTGGCGGCGGTGCTGATCACCCGCCAGATCACGCGCCCGTTGCAAGACACCCTGGTGGCCGTGGAGAAAATCGCCAGTGGCGACCTTACCCAACACATGCGC  >  NC_002947/5284361‑5284651

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