Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	2,308,961	T→G	10.6%	K51T (AAG→ACG)	PP_2031 ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	2,308,961	0	T	G	10.6%	84.8 / 5.7	38	K51T (AAG→ACG)	PP_2031	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (19/15);  new base G (2/2);  total (21/17)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.48e-01

GAAGATCGGCCGCTCGCGTTCTTCGTGGAACAGCTTGGTGTGCGGCTCCTGGGTACGCGCAGTGACACGCCAGTCGATGTTGCGCACATCATCGCCGGCCTGGTAAACGCGCACCTG‑GTCGAAGTCCACGCCGCGCCCGCGCAGCTTGGAATGGTGCAGGCCCACCAGTGGGCTGCGCTGGCCGGGCTTGGAAAACAGCTGGATTTCGCGCACGCGGTGGCGCATGTCGATCAGCTCGGCAAGGCCAATGCGGATGCCGTGTTCGGCCAGCTGCGCAGTGGGCATGGGTTCAG  >  NC_002947/2308816‑2309108                                                                                                                                                   |                                                                                                                                                    gAAGATCGGCCGCTCGCGTTCTTCGTGGAACAGCTTGGTGTGCGGCTCCTGGGTACGCGCAGTGACACGCCAGTCGATGTTGCGCACATCATCGCCGGCCTGGTAAACGCGCACCTG‑GTCGAAGTCCACGCCGCGCCCGCGCAGCTTGGa                                                                                                                                                 >  2:318535/1‑150 (MQ=255)  aaGATCGGCCGCTCGCGTTCTTCGTGGAACAGCTTGGTGTGCGGCTCCTGGGTACGCGCAGTGACACGCCAGTCGATGTTGCGCACATCATCGCCGGCCTGGTAAACGCGCACCTG‑GTCGAAGTCCACGCCGCGCCCGCGCAGCTTGGaa                                                                                                                                                <  1:516238/150‑1 (MQ=255)    gATCGGCCGCTCGCGTTCTTCGTGGAACAGCTTGGTGTGCGGCTCCTGGGTACGCGCAGTGACACGCCAGTCGATGTTGCGCACATCATCGCCGGCCTGGTAAACGCGCACCTG‑GTCGAAGTCCACGCCGCGCCCGCGCAGCTTGGAATg                                                                                                                                              >  1:30353/1‑150 (MQ=255)           cGCTCGCGTTCTTCGTGGAACAGCTTGGTGTGCGGCTCCTGGGTACGCGCAGTGACACGCCAGTCGATGTTGCGCACATCATCGCCGGCCTGGTAAACGCGCACCTG‑GTCGAAGTCCACGCCGCGCCCGCGCAGCTTGGAATGGTGCAgg                                                                                                                                       <  2:288205/150‑1 (MQ=255)                  gTTCTTCGTGGACAAGCTTGGTGTGCGGCTCCTGGGTACGCGCAGTGACACGCCAGTCGATGTTGCGCACATCATCGCCGGCCTGGTAAACGCGCACCTG‑GTCGAAGTCCACGCCGCGCCCGCGCAGCTTGGAATGGTGCAGGCccacca                                                                                                                                >  2:184454/1‑150 (MQ=255)                  gTTCTTCGTGGAACAGCTTGGTGTGCGGCTCCTGGGTACGCGCAGTGACACGCCAGTCGATGTTGCGCACATCATCGCCGGCCTGGTAAACGCGCACCTG‑GTCGAAGTCCACGCCGCGCCCGCGCAGCTTGGAATGGTGCAGGCccacca                                                                                                                                >  2:376186/1‑150 (MQ=255)                          tGGAACAGCTTGGTGTGCGGCTCCTGGGTACGGGCATTGACACGCCATTCGATGTTGCGCCCATCATCGCCGGCCCGCCAGTCACAC‑GCTGCACCTAAGTCCACGCCGCGCCCGCGCAGCTTGGAATGGTGACGGCCCAGTAGTGGGCTg                                                                                                                        >  1:734162/1‑150 (MQ=255)                                                           gcAGTGACACGCCAGTCGATGTTGCGCACATCATCGCCGGCCTGGTAAACGCGCACCTG‑GTCGAAGTCCACGCCGCGCCCGCGCAGCTTGGAATGGTGCAGGCCCACCAGTGGGCTGCGCTGGCCGGGCTTGGAAAACAGCTGGATTTcg                                                                                       >  1:286301/1‑150 (MQ=255)                                                                  cacGCCAGTCGATGTTGCGCACATCATCGCCGGCCTGGTAAACGCGCACCTG‑GTCGAAGTCCACGCCGCGCCCGCGCAGCTTGGAATGGTGCAGGCCCACCAGTGGGCTGCGCTGGCCGGGCTTGGAAAACAGCTGGATTTCGCGCAcgc                                                                                <  2:286301/150‑1 (MQ=255)                                                                   acGCCAGTCGATGTTGCGCACATCATCGCCGGCCTGGTAAACGCGCACCTG‑GTCGAAGTCCACGCCGCGCCCGCGCAGCTTGGAATGGTGCAGGCCCACCAGTGGGCTGCGCTGGCCGGGCTTGGAAAACAGCTGGATTTCGCGCAcgcg                                                                               >  1:554678/1‑150 (MQ=255)                                                                     gCCAGTCGATGTTGCGCACATCATCGCCGGCCTGGTAAACGCGCACCTG‑GTCGAAGTCCACGCCGCGCCCGCGCAGCTTGGAATGGTGCAGGCCCACCAGTGGGCTGCGCTGGCCGGGCTTGGAAAACAGCTGGATTTCGCGCACGCGGt                                                                             >  1:252544/1‑150 (MQ=255)                                                                           cGATGTTGCGCACATCATCGCCGGCCTGGTAAACGCGCACCTG‑GTCGAAGTCCACGCCGCGCCCGCGCAGCTTGGAATGGTGCAGGCCCACCAGTGGGCTGCGCTGGCCGGGCTTGGAAAACAGCTGGATTTCGCGCACGCGGTGGCGCa                                                                       <  2:30353/150‑1 (MQ=255)                                                                               gTTGCGCACATCATCGCCGGCCTGGTAAACGCGCACCTG‑GTCGAAGTCCACGCCGCGCCCGCGCAGCTTGGAATGGTGCAGGCCCACCAGTGGGCTGCGCTGGCCGGGCTTGGAAAACAGCTGGATTTCGCGCACGCGGTGGCGCATGTc                                                                   <  1:184454/150‑1 (MQ=255)                                                                                  gcgcACATCATCGCCGGCCTGGTAAACGCGCACCTG‑GTCGAAGTCCACGCCGCGCCCGCGCAGCTTGGAATGGTGCAGGCCCACCAGTGGGCTGCGCTGGCCGGGCTTGGAAAACAGCTGGATTTCGCGCACGCGGTGGCGCATGTCGAt                                                                <  2:734162/150‑1 (MQ=255)                                                                                    gcACATCATCGCCGGCCTGGTAAACGCGCACCTG‑GTCGAAGTCCACGCCGCGCCCGCGCAGCTTGGAATGGTGCAGGCCCACCAGTGGGCTGCGCTGGCCGGGCTTGGAAAACAGCTGGATTTCGCGCACGCGGTGGCGCATGTCGATCa                                                              >  2:634334/1‑150 (MQ=255)                                                                                        atcatcGCCGGCCTGGTAAACGCGCACCTG‑GTCGAAGTCCACGCCGCGCCCGCGCAGCTTGGAATGGTGCAGGCCCACCAGTGGGCTGCGCTGGCCGGGCtt                                                                                                          <  1:51825/102‑1 (MQ=255)                                                                                        atcatcGCCGGCCTGGTAAACGCGCACCTG‑GTCGAAGTCCACGCCGCGCCCGCGCAGCTTGGAATGGTGCAGGCCCACCAGTGGGCTGCGCTGGCCGGGCtt                                                                                                          >  2:51825/1‑102 (MQ=255)                                                                                         tcatcGCCGGCCTGGTAAACGCGCACCTG‑GTCGAAGTCCACGCCGCGCCCGCGCAGCTTGGAATGGTGCAGGCCCACCAGTGGGCTGCGCTGGCCGGGCTTGGAAAACAGCTGGATTTCGCGCACGCGGTGGCGCATGTCGATCAGCTCg                                                         >  1:40969/1‑150 (MQ=255)                                                                                            tcGCCGGCCTGGTAAACGCGCACCTG‑GTCGAAGTCCACGCCGCGCCCGCGCAGCTTGGAATGGTGCAGGCCCACCAGTGGGCTGCGCTGGCCGGGCTTGGAAAACAGCTGGATTTCGCGCACGCGGTGGCGCATGTCGATCAGCTCGGCa                                                      >  1:58486/1‑150 (MQ=255)                                                                                                 agcgcagcccactggtgggcctgcaccatTCCAAGCTGCGCGGGCGCGGCGTGGAATGGTGCAGGCCCACCAGTGGGCTGCGCTGGCCGGGCTTGGAAAACAGCTGGATTTCGCGCa                                                                                   <  1:217633‑M1/88‑1 (MQ=255)                                                                                                 agcgcagcccactggtgggcctgcaccatTCCAAGCTGCGCGGGCGCGGCGTGGAATGGTGCAGGCCCACCAGTGGGCTGCGCTGGCCGGGCTTGGAAAACAGCTGGATTTCGCGCa                                                                                   >  2:217633‑M1/30‑117 (MQ=255)                                                                                                  gCCTGGTAAACGCGCACCTG‑GTCGAAGTCCACGCCGCGCCCGCGCAGCTTGGAATGGTGCAGGCCCACCAGTGGGCTGCGCTGGCCGGGCTTGGAAAACAGCTGGATTTCGCGCACGCGGTGGCGCATGTCGATCAGCTCGGCaa                                                     <  1:321628/145‑1 (MQ=255)                                                                                                  gCCTGGTAAACGCGCACCTG‑GTCGAAGTCCACGCCGCGCCCGCGCAGCTTGGAATGGTGCAGGCCCACCAGTGGGCTGCGCTGGCCGGGCTTGGAAAACAGCTGGATTTCGCGCACGCGGTGGCGCATGTCGATCAGCTCGGCaa                                                     >  2:321628/1‑145 (MQ=255)                                                                                                       gTAAACGCGCACCTG‑GTCGAAGTCCACGCCGCGCCCGCGCAGCTTGGAATGGTGCAGGCCCACCAGTGGGCGGCGCTGGCCGGGCTTGGAAAACAGCTGGATTTCGCGCACGCGGTGGCGCATGTCGATCAGCTCGGCAAGGCCAATGCg                                           >  2:705434/1‑150 (MQ=255)                                                                                                                 gggcctgcaccatTCCAAGCTGCGCGGGCGCGGCGTGGAATGGTGCAGGCCCACCAGTGGGCTGCGCTGGCCGGGCTTGGAAAACAGCTGGATTTCGCGCACGCGGTGGCGCATGTCGATCAGCTCGGCAAGGCCAATGCGGATGCCgtg                                  <  1:688731‑M1/137‑1 (MQ=255)                                                                                                                       gcaccatTCCAAGCTGCGCGGGCGCGGCGTGGAATGGTGCAGGCCCACCAGTGGGCTGCGCTGGCCGGGCTTGGAAAACAGCTGGATTTCGCGCACGCGGTGGCGCATGTCGATCAGCTCGGCAAGGCCAATGCGGATGCCGTGTTCGGc                            >  2:688731‑M1/8‑150 (MQ=255)                                                                                                                       gTCGAAGTCCACGCCGCGCCCGCGCAGCTTGGAATGGTGCAGGCCCACCAGTGGGCTGCGCTGGCCTGGCTTGGAAAACAGCTGGATTTCGCGCACGCGGTGGCGCATGTCGATCAGCTCGGCAAGGCCAATGCGGATGCCGTGTTCGGc                            >  1:266765/1‑150 (MQ=255)                                                                                                                       gTCGAAGTCCACGCCGCGCCCGCGCAGCTTGGAATGGTGCAGGCCCACCAGTGGGCTGCGCTGGCCGGGCTTGGAAAACAGCTGGATTTCGCGCACGCGGTGGCGCATGTCGATCAGCTCGGCAAGGCCAATGCGGATGCCGTGTTCGGc                            >  1:573277/1‑150 (MQ=255)                                                                                                                            aGTCCACGCCGCGCCCGCGCAGCTTGGAATGGTGCAGGCCCACCAGTGGGCTGCGCTGGCCGGGCTTGGAAAACAGCTGGATTTCGCGCACGCGGTGGCGCATGTCGATCAGCTCGGCAAGGCCAATGCGGATGCCGTGTTCGGCCAGCt                       <  1:576546/150‑1 (MQ=255)                                                                                                                              TCCACGCCGCGCCCGCGCAGCTTGGAATGGTGCAGGCCCACCAGTGGGCTGCGCTGGCCGGGCTTGGAAAACAGCTGGATTTCGCGCAGGCGGCGGCGCATGTCGATCAGCTCGGCAAGGCCAATGCGGATGCCGTGTTCGGCCAGCTgc                     <  1:357790‑M1/150‑1 (MQ=255)                                                                                                                                 ACGCCGCGCCCGCGCAGCTTGGAATGGTGCAGGCCCACCAGTGGGCTGCGCTGGCCGGGCTTGGAAAACAGCTGGATTTCGCGCACGCGGTGGCGCATGTCGATCAGCTCGGCAAGGCCAATGCGGATGCCGTGTTCGGCCAGCTGCGCa                  >  1:376186‑M1/1‑150 (MQ=255)                                                                                                                                    ccgcgcccgcgcAGCTTGGAATGGTGCAGGCCCACCAGTGGGCTGCGCTGGCCGGGCTTGGAAAACAGCTTGATTTCGCGCACGCGGTGGCGCATGTCGATCAGCTCGGCAAGGCCAATGCGGATGCCGTGTTCGGCCAGCTGCGCAGTg               <  1:288784/150‑1 (MQ=255)                                                                                                                                           cgcgcAGCTTGGAATGGTGCAGGGCTGCCAGTGGGCTGCGCTGGCCGGGCTTTGAAAACAGCTGGATTTCGCGCACGCGGTGGCGCATGTCGATCAGCTCGGCAAGGCCAATGCGGATGCCGTGTTCGGCCAGCTGCGCAGTGGGCATgg        <  2:266765/150‑1 (MQ=255)                                                                                                                                            gcgcAGCTGGGAATGGTGCAGGCCCACCAGTGGGCTGCGCTGGCCGGGCTTGGAAAACAGCTGGATTTCGCGCACGCGGTGGCGCATGTCGATCAGCTCGGCACGGCCACTGCTGATGCCGTGTTCGGCCAGCTGCGCAGTGGGCATggg       >  2:552040/1‑150 (MQ=255)                                                                                                                                              gcAGCTTGGAATGGTGCAGGCCCACCAGTGGGCTGCGCTGGCCGGGCTTGGAAAACAGCTGGATTTCGCGCACGCGGTGGCGCATGTCGATCAGCTCGGCAAGGCCAATGCGGATGCCGTTTTCGGCCAGCTGCGCAGTGGGCATGGGtt     <  2:554678/150‑1 (MQ=255)                                                                                                                                                 gCTTGGAATGGTGCAGGCCCCCCAGTGGGCTGCGCTGGCCGGGCTTGGAGAACAGCTGGATTTCGCGCACGCGGTGGCGCATGTCGATCAGCTCGGCAAGGCCAATGCGGATGCCGTGTTCGGCCAGCTGCGCAGTGGGCATGGGTTCAg  <  2:6135/150‑1 (MQ=255)                                                                                                                                                 gCTTGGAATGGTGCAGGCCCACCAGTGGGCTGCGCTGGCCGGGCTTGGAAAACAGCTGGATTTCGCGCACGCGGTGGCGCATGTCGATCAGCTCGGCAAGGCCAATGCGGATGCCGTGTTCGGCCAGCTGCGCAGTGGGCATGGGTTCAg  <  1:318535/150‑1 (MQ=255)                                                                                                                                                 gCTTGGAATGGTGCAGGCCCACCAGTGGGCTGCGCTGGCCGGGCTTGGAAAACAGCTGGATTTCGCGCACGCGGTGGCGCATGTCGATCAGCTCGGCAAGGCCAATGCGGATGCCGTGTTCGGCCAGCTGCGCAGTGGGCATGGGTTCAg  >  1:71815/1‑150 (MQ=255)                                                                                                                                                   |                                                                                                                                                    GAAGATCGGCCGCTCGCGTTCTTCGTGGAACAGCTTGGTGTGCGGCTCCTGGGTACGCGCAGTGACACGCCAGTCGATGTTGCGCACATCATCGCCGGCCTGGTAAACGCGCACCTG‑GTCGAAGTCCACGCCGCGCCCGCGCAGCTTGGAATGGTGCAGGCCCACCAGTGGGCTGCGCTGGCCGGGCTTGGAAAACAGCTGGATTTCGCGCACGCGGTGGCGCATGTCGATCAGCTCGGCAAGGCCAATGCGGATGCCGTGTTCGGCCAGCTGCGCAGTGGGCATGGGTTCAG  >  NC_002947/2308816‑2309108

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