Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	2,937,810	G→T	17.6%	R430S (CGC→AGC)	PP_2570 ←	transposase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	2,937,810	0	G	T	17.6%	37.8 / 2.2	22	R430S (CGC→AGC)	PP_2570	transposase
Reads supporting (aligned to +/- strand):  ref base G (8/10);  new base T (2/2);  total (10/12)
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.81e-01

AAGATCGGGATATCACTCACACCTGAAGGTCGAGTTTGCCTGTGCCGGCCCCTTCGCGGGTGAACCCGCTCTCATCAAACAAAATGCAACCCATTGTTTTCTAAGCACTCAAGCGGCCACCAGTTGATACCCCAGTTTCCGGGCTTTGCGAAGCAAAGCCCGGACCTGCCGGTCTTGGCTTCTGGCTTCGAACTCCTCAAG  >  NC_002947/2937661‑2937861                                                                                                                                                      |                                                    aaGATCGGGATATCACTCACACCTGAAGGTCGAGTTTGCCTGTGCCGGCCCCTTCGCGGGTGAACCCGCTCTCATCCAACAAAATGCAACCCTTTGTTTTCTAAGCACTCAAGCGGCCACCAGTTGATACCCCAGTTTCCGGGCTTTGCg                                                     <  1:156844/150‑1 (MQ=255)          atatCACTCACACCTGAAGGTCGAGTTTGCCTGTGCCGGCCCCTTCGCTGGTGAACCCGCTCTCATCAAACAAAATGCAACCCATTGTTTTCTAAGCACTCAAGCGGCCACCAGTTGATACCCCAGTTTCCGGGCTTTGCGAAGCAAAGc                                            >  1:239520/1‑150 (MQ=255)            atCACTCACACCTGAAGGTCGAGTTTGCCTGTGCCGGCCCCTTCGCTGGTGAACCCGCTCTCATCAAACAAAATGCAACCCATTGTTTTCTAAGCACTCAAGCGGCCACCAGTTGATACCCCAGTTTCCGGGCTTTGCGAAGCAAAGccc                                          <  2:239520/150‑1 (MQ=255)            atCACTCACACCTGAAGGTCGAGTTTGCCTGTGCCGGCCCCTTCGCGGGTGAACCCGCTCTCATCAAACAAAATGCAACCCATTGTTTTCTAAGCACTCAAGCGGCCACCAGTTGATACCCCAGTTTCCGGGCTTTGCGAAGCAAAGccc                                          <  1:928681/150‑1 (MQ=255)                ctcacACCTGAAGGTCGAGTTTGCCTGTGCCGGCCCCTTCGCGGGTGAACCCGCTCTCATCAAACAAAATGCAACCCATTGTTTTCTAAGCACTCAAGCGGCCACCAGTTGATACCCCAGTTTCCGGGCTTTGCGAAGCAAAGCCCGGAc                                      <  1:1153540/150‑1 (MQ=255)                   acacCTGAAGGTCGAGTTTGCCTGTGCCGGCCCCTTCGCGGGTGAACCCGCTCTCATCAAACAAAATGCAACCCATTGTTTAAGAGTCACTCAAGCGGCCACCAGTTGATACCCCAGTTTCCGGGCTTTGCGAAGCAAAGTCCGGACCTg                                   >  1:720493/1‑150 (MQ=17)                    cacCTGAAGGTCGAGTTTGCCTGTGCCGGCCCCTTCGCGGGTGAACCCGCTCTCATCAAACAAAATGCAACCCATTGTTTTCTAAGCACTCAAGCGGCCACCAGTTGATACCCCAGTTTCCGGGCTTTGCGAAGCAAAGCCCGGACCTGc                                  <  2:220448/150‑1 (MQ=255)                      ccTGAAGGTCGAGTTTGCCTGTGCCGGCCCCTTCGCGGGTGAACCCGCTCTCATCAAACAAAATGCAACCCATTGTTTTCTAAGCACTCAAGCGGCCACCAGTTGATCCCCCAGTTTCCGGGCTTTGCGAAGCAAAGCCCGGACCTGCCg                                <  2:493993/150‑1 (MQ=255)                      ccTGAAGGTCGAGTTTGCCTGTGCCGGCCCCTTCGCGGGTGAACCCGCTCTCATCAAACAAAATGCAACCCATTGTTTTCTAAGCACTCAAGCGGCCAACAGTTGATACCCCAGTTTCCGGGCTTTGCGACGCAAAGCCCGGACCTGCCg                                <  2:688473/150‑1 (MQ=255)                       cTGAAGGTCGAGTTTGCCTGTGCCGGCCCCTTCGCGGGTGAACCCGCTCTCATCAAACAAAATGCAACCCATTGTTTTCTAAGCACTCAAGCGGCCACCAGTTGATACCCCAGTTTCCGGGCTTTGCGAAGCAAAGCCCGGACCTGCCgg                               >  1:80192/1‑150 (MQ=255)                       cTGAAGGTCGAGTTTGACTGTGCCGGCCCCTTCGCGGGTGAACCCGCTAATCTCAAACCAAATGAAACCCATTGTTTTCTAAGCACTCAAGCGGCCACCAGTTGATACCCCAGTTTCCGGGCTGTGCGAAGCAAAGCCCGGACCTGCCgg                               >  2:500142/1‑150 (MQ=17)                           aGGGCGAGTTTGCCTGTGACGGAAACTTCGCGGTTTTACTCGCTCTCATCAAACAAAATGCAACCCATTGTTTTCTAAGCACTCAAGCGGCCACCAGTTGATACCCCAGTTTCCGGGCTTTGCGAAGCAAAGCCCGGACCTGCCGGTCtt                           <  2:11703/150‑1 (MQ=11)                             gTCGAGTTTGCCTGTGCCGTCCCCTTCGCGGGTGAACCCGCTCTCATCAAACAAAATGCAACCCATTGTTTTCTAAGCACTCAAGCGGCCACCAGTTGATACCCCAGTTTCCGGGCTTTGCTTCGCAAAGCCCg                                         <  1:204581/134‑1 (MQ=17)                             gTCGAGTTTGCCTGTGCCGGCCCCTTCGCGGGTGAACCCGCTCTCATCAAACAAAATGCAACCCATTGTTTTCTAAGCACTCAAGCGGCCACCAGTTGATACCCCAGTTTCCGGGCTTTGCTTCGCAAAGCCCg                                         >  2:204581/1‑134 (MQ=255)                                   tttGCCTGTGCCGGCCCCTTCGCGGGTGAACCCGCTCTCATCAAACAAAATGCAACCCATTGTTTTCTAAGCACTCAAGCGGCCACCAGTTGATACCCCAGTTTCCGGGCTTTGCGAAGCAAAGCCCGGACCTGCCGGTCTtggcttctg                   <  1:949080/150‑1 (MQ=255)                                      gCCTGTGCCGGCCCCTTCGCGGGTGAACCCGCTCTCATCAAACAAAATGCAACCCATTGTTTTCTAAGCACTCAAGCGGCCACCAGTTGATACCCCAGTTTCCGGGCTTTGCGAAGCAAAGCCCGGACCCGCCGGTCTtggcttctggct                >  1:264552/1‑150 (MQ=255)                                        cTGTGCCGGCCCCTTCGCGGGTGAACCCGCTCTCATCAAACAAAATGCAACCCATTGTTTTCTAAGCACTCAAGCGGCCAGCAGTTGATACCCCAGTTTCCGGGCTTTGCGAAGCAAAGCCCGGACCTGCCGGTCtt                           <  2:832113/137‑1 (MQ=255)                                        cTGTGCCGGCCCCTTCGCGGGTGAACCCGCTCTCATCAAACAAAATGCAACCCATTGTTTTCTAAGCACTCAAGCGGCCACCAGTTGATACCCCAGTTTCCGGGCTTTGCTTCGCAAAGCCCgg                                        >  1:561039/1‑124 (MQ=255)                                        cTGTGCCGGCCCCTTCGCGGGTGAACCCGCTCTCATCAAACAAAATGCAACCCATTGTTTTCTAAGCACTCAAGCGGCCACCAGTTGATACCCCAGTTTCCGGGCTTTGCTTCGCAAAGCCCgg                                        <  2:561039/124‑1 (MQ=255)                                        cTGTGCCGGCCCCTTCGCGGGTGAACCCGCTCTCATCAAACAAAATGCAACCCATTGTTTTCTAAGCACTCAAGCGGCCACCAGTTGATACCCCAGTTTCCGGGCTTTGCGAAGCAAAGCCCGGACCTGCCGGTCtt                           >  1:832113/1‑137 (MQ=39)                                               ggCCCCTTCGCGGGTGAACCCGCTCTCATCAAACAAAATGCAACCCATTGTTTTCTAAGCACTCAAGCGGCCACCAGTTGATACCCCAGTTTCCGGGCTTTGCGAAGCAAAGCCCGGACCTGCCGGTCTTGGCTTCTGGCTTCGAActcc       <  1:601007/150‑1 (MQ=32)                                                   ccTTCGCGGGTGAACCCGCTCTCATCAAACAAAATGCAACCCATTGTTTTCTAAGCACTCAAGCGGCCACCAGTTGATACCCCAGTTTCCGGGCTTTGCGAAGCAAAGCCCGGACCTGCCGGTCTTGGCTTCTGGCTTCGAACTCCTCaa   >  2:854874/1‑150 (MQ=31)                                                    cTTCGCGGGTGAACCCGCTCTCATCAAACAAAATGCAACCCATTGTTTTCTAAGCACTCAAGCGGCCACCAGTTGATACCCCAGTTTCCGGGCTTTGCGAAGCAAAGCCCGGACCTGCCGGTCTTGGCTTCTGGCTTCGAACTCCTCAAg  >  2:521572/1‑150 (MQ=11)                                                                                                                                                      |                                                    AAGATCGGGATATCACTCACACCTGAAGGTCGAGTTTGCCTGTGCCGGCCCCTTCGCGGGTGAACCCGCTCTCATCAAACAAAATGCAACCCATTGTTTTCTAAGCACTCAAGCGGCCACCAGTTGATACCCCAGTTTCCGGGCTTTGCGAAGCAAAGCCCGGACCTGCCGGTCTTGGCTTCTGGCTTCGAACTCCTCAAG  >  NC_002947/2937661‑2937861

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