Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_002947	4,618,760	G→A	G1055S (GGC→AGC)	PP_4085 →	rhs family protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	4,618,760	0	G	A	96.2%	82.1 / ‑3.8	26	G1055S (GGC→AGC)	PP_4085	rhs family protein
Reads supporting (aligned to +/- strand):  ref base G (1/0);  new base A (8/17);  total (9/17)
Fisher's exact test for biased strand distribution p-value = 3.46e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.26e-01

CTGCACCAGATCAACCTGAACGGCCGGGTGATCAGCGATTTCGAGCGTGATGCGGTGCACGACGAGGTCCTGCGGACCCAAGGCAAACTGGTCACCCGCACGCGCTACGACACCAGCGGTCGCTTGGCCGGCAAGGCTATTCACTACCGCGATGCCCCTGCTGAAGTGCTGCCGCTGCTGGACAAGGCCTATCGCTACGATGCCAGCGACAACCTGATTGCCGAGGTACTGACCCAGACACAGCGCCGAGGG  >  NC_002947/4618631‑4618882                                                                                                                                  |                                                                                                                           cTGCACCATATCAACCTCCACGGCCGGGTGATCAGCGATTTCGAGCGTGATGCGGTGCACGACGAGGTCCTGCGGACCCAAGGCAAACTGGTCACCCGCACGCGCTACGACACCAGCGGTCGCTTGGCCAGCAAgg                                                                                                                      <  1:1063126/136‑1 (MQ=255) cTGCACCAGATCAACCTGAACGGCCGGGTGATCAGCGATTTCGAGCGTGATGCGGTGCACGACGAGGTCCTGCGGACCCAAGGCAAACTGGTCACCCGCACGCGCTACGACACCAGCGGTCGCTTGGCCAGCAAgg                                                                                                                      <  2:791738/136‑1 (MQ=255)      ccAGATCAACCTGAACGGCCGGGTGATCAGCGATTTCGAGCGTGATGCGGTGCACGACGAGGTCCTGCGGACCCAAGGCAAACTGGTCACCCGCACGCGCTACGACACCAGCGGTCGCTTGGCCAGCAAGGCTAtt                                                                                                                 <  1:659460/136‑1 (MQ=255)      ccAGATCAACCTGAACGGCCGGGTGATCAGCGATTTCGAGCGTGATGCGGTGCACGACGAGGTCCTGCGGACCCAAGGCAAACTGGTCACCCGCACGCGCTACGACACCAGCGGTCGCTTGGCCAGCAAGGCTAtt                                                                                                                 <  1:600232/136‑1 (MQ=255)            cAACCTGAACGGCCGGGTGATCAGCGATTTCGAGCGTGATGCGGTGCACGACGAGGTCCTGCGGACCCAAGGCAAACTGGTCACCCGCACGCGCTACGACACCAGCGGTCGCTTGGCCAGCAAGGCTATTCACTAc                                                                                                           >  2:348184/1‑136 (MQ=255)                   aaCGGGCGGTTGATCAGCGATTTCGAGCGTGTTGCGGTGCACGACGAGGTCCTGCGGACCCAAGGCAAACTGGTCACCCGCACGCGCTACGACACCAGCGGTCGCTTGGCCAGCAAGGCTATTCACTACCGCGATg                                                                                                    <  2:964922/136‑1 (MQ=255)                           ggTGATCAGCGATTTCGAGCGTGATGCGGTGCACGACGAGGTCCTGCGGACCCAAGGCAAACTGGTCACCCGCACGCGCTACGACACCAGCGGTCGCTTGGCCAGCAAGGCTATTCACTACCGCGa                                                                                                      <  2:332860/126‑1 (MQ=255)                           ggTGATCAGCGATTTCGAGCGTGATGCGGTGCACGACGAGGTCCTGCGGACCCAAGGCAAACTGGTCACCCGCACGCGCTACGACACCAGCGGTCGCTTGGCCAGCAAGGCTATTCACTACCGCGa                                                                                                      >  1:332860/1‑126 (MQ=255)                                  aGCGATTTCGCGCGTGATGCGGTGCACGACGAGGGCCTGCGGACCCAAGGCAAACTGGTCACCCGCACGCGCTACGACACCAGCGGTCGCTTGGCCAGCAAGGCTATTCACTACCGCGATGCCCCTGCTGAAGtgc                                                                                     <  1:695839/136‑1 (MQ=255)                                              ggTGATGCGGTGCACGACGAGGTCCTGCGGACCCAAGGCAAACTGGTCACCCGCACGCGCTACGACACCAGCGGTCGCTTGGCCAGCAAGGCTATTCACTACCGCGATGCCCCTGCTGAAGTGCTGCCGCTGCTgg                                                                         <  1:191391/135‑1 (MQ=255)                                                ggATGCGGTGCACGACGAGGTCCTGCGGACCCAAGGCAAACTGGTCACCCGCACGCGCTACGACACCAGCGGTCGCTTGGCCAGCAAGGCTATTCACTACCGCGATGCCCCTGCTGAAGTGCTGCCGCTGCTGGAc                                                                       <  1:603606/135‑1 (MQ=255)                                                             gacgaGGTCCTGCGGACCCAAGGCAAACTGGTCACCCGCACGCGCTACGACACCAGCGGTCGCTTGGCCAGCAAGGCTATTCACTACCGCGATGCCCCTGCTGAAGTGCTGCCGCTGCTGGACAAGGCCTATCGCt                                                          >  1:101902/1‑136 (MQ=255)                                                                 aGGTCCTGCGGACCCAAGGCAAACTGGTCACCCGCACGCGCTACGACACCAGCGGTCGCTTGGCCAGCAAGGCTATTCACTACCGCGATGCCCCTGCTGAAGTGCTGCCGCTGCTGGACAAGGCCTATCGCTACGa                                                      <  2:1091041/136‑1 (MQ=255)                                                                 aGGTCCTGCGGACCCAAGGCAAACTGGTCACCCGCACGCGCTACGACACCAGCGGTCGCTTGGCCAGCAAGGCTATTCACTACCGCGATGCCCCTGCTGAAGTGCTGCCGCTGCTGGACAAGGCCTATCGCTACGa                                                      <  2:146817/136‑1 (MQ=255)                                                                   gTCCTGCGGACCCAAGGCAAACTGGTCACCCGCACGCGCTACGACACCAGCGGTCGCTTGGCCAGCAAGGCTATTCACTACCGCGATGCCCCTGCTGAAGTGCTGCCGCTGCTGGACAAGGCCTATCGCTACGATg                                                    <  1:1199494/136‑1 (MQ=255)                                                                               cAAGGCAAACTGGTCACCCGCACGCGCTACGACACCAGCGGTCGCTTGGCCAGCAAGGCTATTCACTACCGCGATGCCCCTGCTGAAGTGCTGCCGCTGCTGGACAAGGCCTATCGCTACGATGCCAGCGACAAcc                                        >  2:695224/1‑136 (MQ=37)                                                                                 aGGCAAACTGGTCACCCGCACGCGCTACGACACCAGCGGTCGCTTGGCCAGCAAGGCTATTCACTACCGCGATGCCCCTGCTGAAGTGCTGCCGCTGCTGGACAAGGCCTATCGCTACGATGCCAGCGACAACCTg                                      >  1:466384/1‑136 (MQ=37)                                                                                          ggTCACCCGCACGCGCTACGACACCAGCGGTCGCTTGGCCAGCAAGGCTATTCACTACCGCGATGCCCCTGCTGAAGTGCTGCCGCTGCTGGACAAGGCCTATCGCTACGATGCCAGCGACAACCTGATTGCCGAg                             <  1:589776/136‑1 (MQ=37)                                                                                               cccGCACGCGCTACGACACCAGCGGTCGCTTGGCCAGCAAGGCTATTCACTACCGCGATGCCCCTGCTGAAGTGCTGCCGCTGCTGGACAAGGCCTATCGCt                                                          >  2:137017/1‑102 (MQ=255)                                                                                               cccGCACGCGCTACGACACCAGCGGTCGCTTGGCCAGCAAGGCTATTCACTACCGCGATGCCCCTGCTGAAGTGCTGCCGCTGCTGGACAAGGCCTATCGCt                                                          <  1:137017/102‑1 (MQ=255)                                                                                                      gcgcTACGACACCAGCGGTCGCTTGGCCGGCAAGGCTATTCACTACCGCGATGCCCCTGCTGAAGTGCTGCCGCTGCTGGACAAGGCCTATCGCTACGATGCCAGCGACAACCTGATTGCCGAGGTACTGACCCAg                 >  1:601890/1‑136 (MQ=39)                                                                                                      gcgcTACGACACCAGCGGTCGCTTGGCCAGCAAGGCTATTCACTACCGCGATGCCCCTGCTGAAGTGCTGCCGCTGCTGGACAAGGCCTATCGCTACGATGCCAGCGACAACCTGATTGCCGAGGTACTGACCCAg                 >  2:1019444/1‑136 (MQ=255)                                                                                                        gcTACGACACCAGCGGTCGCTTGGCCAGCAAGGCTATTCACTACCGCGATGCCCCTGCTGAAGTGCTGCCGCTGCTGGACAAGGCCTATCGCTACGATGCCAGCGACAACCTGATTGCCGAGGTACTGACCCAGac               <  1:266923/136‑1 (MQ=255)                                                                                                             gACACCAGCGGTCGCTTGGCCAGCAAGGCTATTCACTACCGCGATGCCCCTGCTGAAGTGCTGCCGCTGCTGGACAAGGCCTATCGCTACGATGCCAGCGACAACCTGATTGCCGAGGTACTGACCCAGACACAgc          <  2:1158772/136‑1 (MQ=255)                                                                                                                  cAGCGGTCGCTTGGCCAGCAAGGCTATTCACTACCGCGATGCCCCTGCTGAAGTGCTGCCGCTGCTGGACAAGGCCTATCGCTACGATGCCAGCGACAACCTGATTGCCGAGGTACTGACCCAGACACAGCGCCGa     <  2:101902/136‑1 (MQ=255)                                                                                                                     cGGTCGCTTGGCCAGCAAGGCTATTCACTACCGCGATGCCCCTGCTGAAGTGCTGCCGCTGCTGGACAAGGCCTATCGCTACGATGCCAGCGACAACCTGATTGCCGAGGTACTGACCCAGACACAGCGCCGAggg  >  2:855565/1‑136 (MQ=255)                                                                                                                                  |                                                                                                                           CTGCACCAGATCAACCTGAACGGCCGGGTGATCAGCGATTTCGAGCGTGATGCGGTGCACGACGAGGTCCTGCGGACCCAAGGCAAACTGGTCACCCGCACGCGCTACGACACCAGCGGTCGCTTGGCCGGCAAGGCTATTCACTACCGCGATGCCCCTGCTGAAGTGCTGCCGCTGCTGGACAAGGCCTATCGCTACGATGCCAGCGACAACCTGATTGCCGAGGTACTGACCCAGACACAGCGCCGAGGG  >  NC_002947/4618631‑4618882

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