Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	6,165,280	T→A	23.3%	D20V (GAT→GTT)	PP_5405 ←	transposase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	6,165,280	0	T	A	23.3%	29.9 / 4.2	17	D20V (GAT→GTT)	PP_5405	transposase
Reads supporting (aligned to +/- strand):  ref base T (8/5);  new base A (2/2);  total (10/7)
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.95e-01

GACTGCCTCCTCAAAATCAAACTGGCTCAGTTGGATCAGGTGGGTGCTGATGTGTTTCCCATCGTCTTTGAATTCGACCAAGCGAACACGGTTTTGCCGGCCATCGCTGACCAGTAGGTGATATACATCGCCTTGTGGAAGGTTACGAAACCCCTCGGGGGCTATGATTCGCGAGCCTGCCCTCATGGCGCGGCCTCTTCAAAAAGCGAGCCGAACTCCTGGAGCAGATCGCATCCACCGGTAATCAGC  >  NC_002947/6165153‑6165401                                                                                                                                |                                                                                                                          gACTGCCTCCTCAAAATCAAACTGGCTCAGTTGGATCAGGAAGGTGCTGATGTGTTTCCCATCGTCTTTGAATTCGACCAAGCGAACACGGTTTTGCCGGCCATCGCTGACCAGTAGGTGATATACATCGCCTTGTGGAAGGTTACGaaa                                                                                                     >  2:426498/1‑150 (MQ=255)    tGCCTCCTCAAAATCAAACTGGCTCAGTTGGATCAGGTGGGTGCTGATGTGTTTCCCATCGTCTTTGAATTCGACCAAGCGAACACGGTTTTGCCGGCCATCGCTGACCAGTAGGTGATATACATCGCCTTGTGGAAGGTTACGAAAccc                                                                                                  <  2:277479/150‑1 (MQ=255)      cctcctCAAAATCAAACTGGCTCAGTTGGATCAGGTGGGTGCTGATGTGTTTCCCATCGTCTTTGAATTCGACCAAGCGAACACGGTTTTGCCGGCCATCGCTGACCAGTAGGTGATATACAACGCCTTGTGGAAGGTTACGAAACCCCt                                                                                                >  1:574389/1‑150 (MQ=255)       ctcctcAAAATCAAACTGGCTCAGTTGGATCAGGTGGGTGCTGATGTGTTTCCCATCGTCTTTGAATTCGACCAAGCGAACACGGTTTTGCCGGCCATCGCTGACCAGTAGGTGATATACATCGCCTTGTGGAAGGTTACGAAACCCCTc                                                                                               >  2:532204/1‑150 (MQ=255)          ctcAAAATCAAACTGGCTCAGTTGGATCAGGTGGGTGCTGATGTGTTTCCCATCGTCTTTGAATTCGACCAAGAGAACACGGTTTTGCCGGCCATCGCTGACCAGTAGGTGATATACATCGCCTTGTGGAAGGTTACGAAACCCCTCggg                                                                                            >  2:91101/1‑150 (MQ=255)            cAAAATCAAACTGGCTCAGTTGGATCAGGTGGGTGCTGATGTGTTTCCCATCGTCTTTGAATTCGACCAAGCGAACACGGTTTTGCCGGCCATCGCTGACCAGTAGGTGATATACATCGCCTTGTGGAAGGTTACGAAACCCCTCggggg                                                                                          >  1:550546/1‑150 (MQ=255)                             aGTTGGATCAGGTGGGTGCTGATGTGTTTCCCATCGTCTTTGAATTCGACCAAGCGAACACGGTTTTGCCGGCCATCGCTGACCAGTAGGTGATATACATCGCCTTGTGGAAGGTTACGAAACCCCTCGGGGGCTATGATTCGCGAGCCt                                                                         <  2:550546/150‑1 (MQ=255)                                                     tgtTTCCCATCGTCTTTGAATTCGACCAAGCGAACACGGTTTTGCCGGCCATCGCTGACCAGTAGGTGATATACATCGCCTTGTGGAAGGTTACGAAACCCCTCGGGGGCTATGATTCGCGAGCCTGCCCTCATGGCGCGGCCTCTTCaa                                                 >  1:555580/1‑150 (MQ=255)                                                            cATCGTCTTTGAATTCGACCAAGCGAACACGGTTTTGCCGGCCATCGCTGACCAGTAGGTGATATACATCGCCTTGTGGAAGGTTACGAAACCCCTCGGGGGCTATGATTCGCGAGCCTGCCCTCATGGCGCGGCCTCTTCAAAAAGCGa                                          >  2:26249/1‑150 (MQ=255)                                                                            gACCAAGCGAACACGGTTTTGCCGGCCATCGCTGACCAGTAGGTGATATACAACGCCTTGTGGAAGGTTACGAAACCCCTCGGGGGCTATGATTCGCGAGCCTGCCCTCATGGCGCGGCCTCTTCAAAAAGCGAGCCGAACTCCTGGAGc                          <  2:574389/150‑1 (MQ=255)                                                                             aCCAAGCGAACACGGTTTTGCCGGCCATCGCTGACCAGTAGGTGATATACAACGCCTTGTGGAAGGTTACGAAACCCCTCGGGGGCTATGATTCGCGAGCCTGCCCTCATGGCGCGGCCTCTTCAAAAAGCGAGCCGAACTc                                 >  1:408157/1‑142 (MQ=255)                                                                             aCCAAGCGAACACGGTTTTGCCGGCCATCGCTGACCAGTAGGTGATATACAACGCCTTGTGGAAGGTTACGAAACCCCTCGGGGGCTATGATTCGCGAGCCTGCCCTCATGGCGCGGCCTCTTCAAAAAGCGAGCCGAACTc                                 <  2:408157/142‑1 (MQ=255)                                                                                 aGCGAACACGGTTTTGCCGGCCATCGCTGACCAGTAGGTGATATACATCGCCTTGTGGAAGGTTACGAAACCCCTCGGGGGCTATGATTCGCGAGCCTGCCCTCATGGCGCGGCCTCTTCaaaaa                                              >  2:26235/1‑125 (MQ=255)                                                                                 aGCGAACACGGTTTTGCCGGCCATCGCTGACCAGTAGGTGATATACATCGCCTTGTGGAAGGTTACGAAACCCCTCGGGGGCTATGATTCGCGAGCCTGCCCTCATGGCGCGGCCTCTTCaaaaa                                              <  1:26235/125‑1 (MQ=255)                                                                                                    gCCATCGGTTACCAGTAGGTGATATACATCGCCTTGTGGAAGGTTACGAAACCCCTCGGGGGCTATGATTCGCGAGCCTGCCCTCATGGCGCGGCCTCTTCAAAAAGCGAGCCGAACTCCTGGAGCAGATCGCATCCACCGGTAATCAGc  <  1:26249/150‑1 (MQ=255)                                                                                                                             aCATCGCCTTGTGGAAGGTTACGAAACCCCTCGGGGGCTATGATTCGCGAGCCTGCCCTCATGGCGCGGCCTCTTCAAAAAGCGAGcc                                       >  1:186232/1‑88 (MQ=255)                                                                                                                             aCATCGCCTTGTGGAAGGTTACGAAACCCCTCGGGGGCTATGATTCGCGAGCCTGCCCTCATGGCGCGGCCTCTTCAAAAAGCGAGcc                                       <  2:186232/88‑1 (MQ=255)                                                                                                                                |                                                                                                                          GACTGCCTCCTCAAAATCAAACTGGCTCAGTTGGATCAGGTGGGTGCTGATGTGTTTCCCATCGTCTTTGAATTCGACCAAGCGAACACGGTTTTGCCGGCCATCGCTGACCAGTAGGTGATATACATCGCCTTGTGGAAGGTTACGAAACCCCTCGGGGGCTATGATTCGCGAGCCTGCCCTCATGGCGCGGCCTCTTCAAAAAGCGAGCCGAACTCCTGGAGCAGATCGCATCCACCGGTAATCAGC  >  NC_002947/6165153‑6165401

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