Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,638,997	A→T	18.2%	H13L (CAC→CTC)	PP_4950 →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	5,638,997	0	A	T	18.2%	50.0 / 5.5	22	H13L (CAC→CTC)	PP_4950	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (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 = 9.95e-01

GCCCAACTGTGGCAAGGGCGAACCCAATCAGGTGGTGCGGGTCGGCCACGCCTCGCCGGCCTGTGTGTTCCGCCAGATCGATGTATTCGGAGGAGACGCCTGATGAAGCACGCTTTCGAAACCCTGGTCGGGGAGGTGCACGCCGCCCTTCAGGCGGGCGAGCAGTTCACCCTCGGTTACAGTGCCGAGCAATCGCAGTTCGTGCGTTTCAACCATGCCAAGGTGCGCCAGGCCGGCGAAGTCAGCCAGGCCAGCGCGCAACTGCGGCTGGTGCGCGATG  >  NC_002947/5638858‑5639137                                                                                                                                            |                                                                                                                                             gCCCAACTGTGGCAAGGGCGAACCCAATCAGGTGGTGCGGGTCGGCCACGCCTCGCCGGCCTGTGTGTTCCGCCAGATCGATGTATTCGGAGGAGACGCCTGATGAAGCACGCTTTCGAAACCCTGGTCGGGGAGGTGCACGCCGCCCtt                                                                                                                                    >  2:609542/1‑150 (MQ=255)             cAAGGGCGAACCCAATCAGGTGGTGCGGGTCGGCCACGCCTCGCCGGCCTGTGTGTTCCGCCAGATCGATGTATTCGGAGGAGACGCCTGATGAAGCACGCTTTCGAAACCCTGGTCGGGGAGGTGCACGCCGCCCTTCAGGCGGGCCAg                                                                                                                        <  1:609542/150‑1 (MQ=255)              aaGGGCGAACCCAATCAGGTGGTGCGGGTCGGCCACGCCTCGCCGGCCTGTGTGTTCCGCCAGATCGATGTATTCGGAGGAGACGCCTGATGAAGCACGCTTTCGAAACCCTGGTCGGGGAGGTGCACGCCGCCCTTCAGGCGGGCGAgc                                                                                                                       <  1:303595‑M2/150‑3 (MQ=255)                        ccAATCAGGTGGTGCGGGTCGGCCACGCCTCGCCGGCCTGTGTGTTCCGCCAGATCGATGTATTCGGAGGAGACGCCTGATGAAGCACGCTTTCGAAACCCTGGTCGGGGAGGTGCACGCCGCCCTTCAGGCGGGCGAGCAGTTCACCCt                                                                                                             >  1:215409/1‑150 (MQ=255)                           aTCAGGTGGTGCGGGTCGGCCACGCCTCGCCGGCCTGTGTGTTCCGCCAGATCGATGTATTCGGAGGAGACGCCTGATGAAGCACGCTTTCGAAACCCTGGTCGGGGAGGTGCACGCCGCCCTTCAGGCGGGCGAGCAGTTCACCCTCgg                                                                                                          <  1:220382/150‑1 (MQ=255)                                gtggtgCGGGTCGGCCACGCCTCGCCGGCCTGTGTGTTCCGCCAGATCGATGTATTCGGAGGAGACGCCTGATGAAGCACGCTTTCGAAACCCTGGTCGGGGAGGTGCACGCCGCCCTTCAGGCGGGCGAGCAGTTCACCCTCGGTTACa                                                                                                     <  1:270409/150‑1 (MQ=255)                                   gtgCGGGTCGGCCACGCCTCGCCGGCCTGTGTGTTCCGCCAGATCGATGTATTCGGAGGAGACGCCTGATGAAGCACGCTTTCGAAACCCTGGTCGGGGAGGTGCACGCCGCCCTTCAGGCGGGCGAGCAGTTCACCCTCGGTTACAGTg                                                                                                  >  2:280669/1‑150 (MQ=255)                                          tCGGCCACGCCTCGCCGGCCTGTGTGTTCCGCCAGATCGATGTATTCGGAGGAGACGCCTGATGAAGCACGCTTTCGAAACCCTGGTCGGGGAGGTGCTCGCCCGCCTGAAGGGCGGCGTgcacctccccgaccagggttt                                                                                                    >  1:125690‑M2/1‑120 (MQ=255)                                          tCGGCCACGCCTCGCCGGCCTGTGTGTTCCGCCAGATCGATGTATTCGGAGGAGACGCCTGATGAAGCACGCTTTCGAAACCCTGGTCGGGGAGGTGCTCGCCCGCCTGAAGGGCGGCGTgcacctccccgaccagggttt                                                                                                    <  2:125690‑M2/141‑22 (MQ=255)                                                            ccTGTGTGTTCCGACAGATCGATGTATTCGGAGGAGACGCCTGATGAAGCACGCTTTCGAAACCCTGGTCGGGGAGGTGCACGCCGCCCTTCAGGCGGGCGAGCAGTTCACCCTCGGTTACAGTGCCGAGCAATCGCAGTTCGTGCGttt                                                                         <  1:280669/150‑1 (MQ=255)                                                                  tgtTCCGCCAGATCGATGTATTCGGAGGAGACGCCTGATGAAGCACGCTTTCGAAACCCTGGTCGGGGAGGTGCACGCCGCCCTTCAGGCGGGCGAGCAGTTCACCCTCGGTTACAGTGCCGAGCAATCGCAGTTCGTGCGTTTCAACCa                                                                   >  2:600448/1‑150 (MQ=255)                                                                    ttCCGCCAGATCGATGTATTCGGAGGAGACGCCTGATGAAGCACGCTTTCGAAACCCTGGTCGGGGAGGTGCTCGCCCGCCTGAAGGGCGGCGTgcacctccccgaccagggtttcgaa                                                                                                <  1:277548‑M2/119‑26 (MQ=255)                                                                    ttCCGCCAGATCGATGTATTCGGAGGAGACGCCTGATGAAGCACGCTTTCGAAACCCTGGTCGGGGAGGTGCTCGCCCGCCTGAAGGGCGGCGTgcacctccccgaccagggtttcgaa                                                                                                >  2:277548‑M2/1‑94 (MQ=255)                                                                                      ttCGGAGGAGACGCCTGATGAAGCACGCTTTCGAAACCCTGGTCGGGGAGGTGCACGCCGCCCTTCAGGCGGGCGAGCAGTTCACCCTCGGTTACAGTGCCGAGCAATCGCAGTTCGTGCGTTTCAACCATGCCAAGGTGCGCCAg                                                   <  2:468323/146‑1 (MQ=255)                                                                                      ttCGGAGGAGACGCCTGATGAAGCACGCTTTCGAAACCCTGGTCGGGGAGGTGCACGCCGCCCTTCAGGCGGGCGAGCAGTTCACCCTCGGTTACAGTGCCGAGCAATCGCAGTTCGTGCGTTTCAACCATGCCAAGGTGCGCCAg                                                   >  1:468323/1‑146 (MQ=255)                                                                                               gaCGCCTGATGAAGCACGCTTTCGAAACCCTGGTCGGGGAGGTGCACGCCGCCCTTCAGGCGGGCGAGCAGTTCACCCTCGGTTACAGTGCCGAGCAATCGCAGTTCGTGCGTTTCAACCATGCCAAGGTGCGCCAGGCCGGCGAAGTCa                                      <  1:600448/150‑1 (MQ=255)                                                                                                        tgaAGCACGCTTTCGAAACCCTGGTCGGGGAGGTGCACGCCGCCCTTCAGGCGGGCGAGCAGTTCACCCTCGGTTACAGTGCCGAGCAATCGCAGTTCGTGCGTTTCAACCATGCCAAGGTGCGCCAGGCCGGCGAAGTCAgccaggcca                             <  2:215409/150‑1 (MQ=255)                                                                                                         gaAGCACGCTTTCGAAACCCTGGTCGGGGAGGTGCACGCCGCCCTTCAGGCGGGCGAGCAGTTCACCCTCGGTTACAGTGCCGAGCAATCGCAGTTCGTGCGTTTCAACCATGCCAAGGTGCGCCAGGCCGGCGAAGTCAgccaggccag                            >  1:63414/1‑150 (MQ=255)                                                                                                            gCACGCTTTCGAAACCCTGGTCGGGGAGGTGCACGCCGCCCTTCAGGCGGGCGAGCAGTTCACCCTCGGTTACAGTGCCGAGCAATCGCAGTTCGTGCGTTTCAACCATGCCAAGGTGCGCCAGGCCGGCGAAGTCAGCCAGGCCAgcgc                         >  2:473797/1‑150 (MQ=255)                                                                                                                          cccTGGTCGGGGAGGTGCACGCCGCCCTTCAGGCGGGCGAGCAGTTCACCCTCGGTTACAGTGCCGAGCAATCGCAGTTCGTGCGTTTCAACCATGCCAAGGTGCGCCAGGCCGGCGAAGTCAGCCAGGCCAGCGCGCAACTGCGGCTgg           <  2:63414/150‑1 (MQ=255)                                                                                                                             tGGTCGGGGAGGTGCACGCCGCCCTTCAGGCGGGCGAGCAGTTCACCCTCGGTTACAGTGCCGAGCAATCGCAGTTCGTGCGTTTCAACCATGCCAAGGTGCGCCAGGCCGGCGAAGTCAGCCAGGCCAGCGCGCAACTGCGGCTGGTgc        >  2:42946/1‑150 (MQ=255)                                                                                                                                   gggAGGTGCACGCCGCCCTTCAGGCGGGCGAGCAGTTCACCCTCGGTTACAGTGCCGAGCAATCGCAGTTCGTGCGTTTCAACCATGCCAAGGTGCGCCAGGCCGGCGAAGTCAGCCAGGCCAGCGCGCAACTGCGGCTGGTGCGCGATg  <  1:413654/150‑1 (MQ=255)                                                                                                                                            |                                                                                                                                             GCCCAACTGTGGCAAGGGCGAACCCAATCAGGTGGTGCGGGTCGGCCACGCCTCGCCGGCCTGTGTGTTCCGCCAGATCGATGTATTCGGAGGAGACGCCTGATGAAGCACGCTTTCGAAACCCTGGTCGGGGAGGTGCACGCCGCCCTTCAGGCGGGCGAGCAGTTCACCCTCGGTTACAGTGCCGAGCAATCGCAGTTCGTGCGTTTCAACCATGCCAAGGTGCGCCAGGCCGGCGAAGTCAGCCAGGCCAGCGCGCAACTGCGGCTGGTGCGCGATG  >  NC_002947/5638858‑5639137

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