Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	2,531,947	T→A	21.0%	A53A (GCT→GCA)	PP_2222 →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	2,531,947	0	T	A	21.0%	40.0 / 5.8	19	A53A (GCT→GCA)	PP_2222	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (9/6);  new base A (2/2);  total (11/8)
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.72e-01

TCGATGTGGCTGTGATGATAGGGTCGGGTGTGCCGGCCACCTTGCAGGCGCTGGGCAAGCGTATTTGCTGGGTAGTGCTGGTCAATGGCGAGCGTCGCGGCACCGCCTTCGCCACGCGCGAAGAAGCGCTTGAGTGCCAGGCTGCCTGGCAGCAGCAGCTGAAGAAAGGGCTGGCTGCCTGAAAATGTTTATTGCCTCTCGACGGGTTTCGGTGTTTTCGGGCATATTTGTCCTACAGTAATGCGCACATGGATGCCAGGGCTTGTCTTCTAGCAGGCTCACGCC  >  NC_002947/2531805‑2532089                                                                                                                                               |                                                                                                                                               tCGATGTGGCTGTGATGATAGGGTCGGGTGTGCCGGCCACCTTGCAGGCGCTGGGCAAGCGTATTTGCTGGGTAGTGCTGGTCAATGGCGAGCGTCGCGGCACCGCCTTCGCCACGCGCGAAGAAGCGCTTGAGTGCCAGGCAGc                                                                                                                                              <  1:180436/145‑1 (MQ=255)  cGATGTGGCTGTGATGATAGGGTCGGGTGTGCCGGCCACCTTGCAGGCGCTGGGCAAGCGTATTTGCTGGGTAGTGCTGGTCAATGGCGAGCGTCGCGGCACCGCCTTCGCCACGCGCGAAGAAGCGCTTGAGTGCCAGGCTGCCTGgca                                                                                                                                        <  1:133707/150‑1 (MQ=255)   gATGTGGCTGTGATGATAGGGTCGGGTGTGCCGGCCACCTTGCAGGCGCTGGGCAAGCGTATTTGCTGGGTAGTGCTGGTCAATGGCGAGCGTCGCGGCACCGCCTTCGCCACGCGCGAAGAAGCGCTTGAGTGCCAGGCTGCCTGgcag                                                                                                                                       >  2:54905/1‑150 (MQ=255)           tgtgATGATAGGGTCGGGTGTGCCGGCCACCTTGCAGGCGCTGGGCAAGCGTATTTGCTGGGTAGTGCTGGTCAATGGCGAGCGTCGCGGCACCGCCTTCGCCGCGCGCGAAGAAGCGCTTGAGTGCCAGGCTGCCTGGCAGCAGCAGCt                                                                                                                               >  1:40972/1‑150 (MQ=255)             tgatgaTAGGGTCGGGTGTGCCGGCCACCTTGCAGGCGCTGGGCAAGCGTATTTGCTGGGTAGTGCTGGTCAATGGCGAGCGTCGCGGCACCGCCTTCGCCACGCGCGAAGAAGCGCTTGAGTGCCAGGCTGCCTGGCAGCAGCAGCTga                                                                                                                             >  1:275074/1‑150 (MQ=255)                                gCCGGCCACCTTGCAGGCGCTGGGCAAGCGTATTTGCTGGGTAGTGCTGGTCAATGGCGAGCGTCGCGGCACCGCCTTCGGCACGCGCGAAGAAGCGCTTGAGTGCCAGGCAGCCTGgcactcaag                                                                                                                                  >  1:155818/1‑120 (MQ=255)                                gCCGGCCACCTTGCAGGCGCTGGGCAAGCGTATTTGCTGGGTAGTGCTGGTCAATGGCGAGCGTCGCGGCACCGCCTTCGCCACGCGCGAAGAAGCGCTTGAGTGCCAGGCAGCCTGgcactcaag                                                                                                                                  <  2:155818/126‑7 (MQ=255)                                                       gCAAGCGTATTTGCTGGGTAGTGCTGGTCAATGGCGAGCGTCGCGGCACCGCCTTCGCCACGCGCGAAGAAGCGCTTGAGTGCCAGGCTGCCTGGCAGCAGCAGCTGAAGAAAGGGCTGGCTGCCTGAAAATGTTTATTGCCTCTCGACg                                                                                   >  1:56983/1‑150 (MQ=255)                                                                  tGCTGGGTAGTGCTGGTCAATGGCGAGCGTCGCGGCACCGCCTTCGCCGCGCGCGAAGAAGCGCTTGAGTGCCAGGCTGCCTGGCAGCAGCAGCTGAAGAAAGGGCTGGCTGCCTGAAAATGTTTATTGCCTCTCGACGGGTTTCGgtgt                                                                        <  2:40972/150‑1 (MQ=255)                                                                    cTGGGTAGTGCTGGTCAATGGCGAGCGTCGCGGCACCGCCTTCGCCACGCGCGAAGAAGCGCTTGAGTGCCAGGCTGCCTGGCAGCAGCAGCTGAAGAAAGGGCTGGCTGCCTGAAAATGTTTATTGCCTCTCGACGGGTTTCGGTGttt                                                                      <  2:275074/150‑1 (MQ=255)                                                                     tGGGTAGTGCTGGTCAATGGCGAGCGTCGCGGCACCGCCTTCGCCACGCGCGAAGAAGCGCTTGAGTGCCAGGCTGCCTGGCAGCAGCAGCTGAAGAAAGGGCTGGCTGCCTGAAAATGTTTATTGCCTCTCGACGGGTTTCGGTGtttt                                                                     <  1:54905/150‑1 (MQ=255)                                                                       ggTAGTGCTGGTCAATGGCGAGCGTCGCGGCACCGCCTTCGCCACGCGCGAAGAAGCGCTTGAGTGCCAGGCTGCCTGGCAGCAGCAGCTGAAGAAAGGGCTGGCTGCCTGAAAATGttt                                                                                                 >  1:389268/1‑120 (MQ=255)                                                                       ggTAGTGCTGGTCAATGGCGAGCGTCGCGGCACCGCCTTCGCCACGCGCGAAGAAGCGCTTGAGTGCCAGGCTGCCTGGCAGCAGCAGCTGAAGAAAGGGCTGGCTGCCTGAAAATGttt                                                                                                 <  2:389268/120‑1 (MQ=255)                                                                               tGGTCAATGGCGAGCGTCGCGGCAGCGCCTTCGCCACGCGCGAAGAAGCGCTTGAGTGCCAGGATACTTGACAACAGAAGCTGAAGAAAAGTCTGGCTGCCTGAAAATGTTTATTCCCTCTCGACGGGGTTGGGTGTTTTCggttgtatt                                                           >  2:211754/1‑143 (MQ=255)                                                                                 gTCAATGGCGAGCGTCGCGGCACCGCCTTCGCCACGCGCGAAGAAGCGCTTGAGTGCCAGGCTGCCTGGCAGCAGCAGCTGAAGAAAGGGCTGGCTGCCTGAAAATGTTTATTGCCTCTCGACGGGTTTCGGTGTTTTCGGGCATATTTg                                                         >  2:104049/1‑150 (MQ=255)                                                                                  tCAATGGCGAGCGTCGCGGCACCGCCTTCGCCACGCGCGAAGAAGCGCTTGAGTGCCAGGCTGCCTGGCAGCAGCAGCTGAAGAAAGGGCTGGCTGCCTGAAAATGTTTATTGCCTCTCGACGGGTTTCGGTGTTTTCGGGCATATTTGt                                                        >  2:216616/1‑150 (MQ=255)                                                                                                            ttCGCCACGCGCGAAGAAGCGCTTGAGTGCCAGGCTGCCTGGCAGCAGCAGCTGAAGAAAGGGCTGGCTGCCTGAAAATGTTTATTGCCTCTCGACGGGTTTCGGTGTTTTCGGGCATATTTGTc                                                       <  2:344870/125‑1 (MQ=255)                                                                                                            ttCGCCACGCGCGAAGAAGCGCTTGAGTGCCAGGCTGCCTGGCAGCAGCAGCTGAAGAAAGGGCTGGCTGCCTGAAAATGTTTATTGCCTCTCGACGGGTTTCGGTGTTTTCGGGCATATTTGTc                                                       >  1:344870/1‑125 (MQ=255)                                                                                                                                        gCCAGGCAGCCTGGCAGCAGCAGCTGAAGAAAGGGCTGGCTGCCTGAAAATGTTTATTGCCTCTCGACGGGTTTCGGTGTTTTCGGGCATATTTGTCCTACAGTAATGCGCACATGGATGCCAGGGCTTGTCTTCTAGCAGGCTCACGcc  >  1:416340/1‑150 (MQ=255)                                                                                                                                               |                                                                                                                                               TCGATGTGGCTGTGATGATAGGGTCGGGTGTGCCGGCCACCTTGCAGGCGCTGGGCAAGCGTATTTGCTGGGTAGTGCTGGTCAATGGCGAGCGTCGCGGCACCGCCTTCGCCACGCGCGAAGAAGCGCTTGAGTGCCAGGCTGCCTGGCAGCAGCAGCTGAAGAAAGGGCTGGCTGCCTGAAAATGTTTATTGCCTCTCGACGGGTTTCGGTGTTTTCGGGCATATTTGTCCTACAGTAATGCGCACATGGATGCCAGGGCTTGTCTTCTAGCAGGCTCACGCC  >  NC_002947/2531805‑2532089

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