Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	1,180,137	C→A	18.7%	P86H (CCC→CAC)	PP_1030 →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	1,180,137	0	C	A	18.7%	72.2 / 2.9	26	P86H (CCC→CAC)	PP_1030	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (8/13);  new base A (3/2);  total (11/15)
Fisher's exact test for biased strand distribution p-value = 6.20e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TGGGCCATGCCCTGGCGGCAGGGGCGCGCAAAGTGGCCCTGCACCCGGAAAAACTGCCCTTCGGCGAACCGGTCAACGGCCTGGAAGTGATGCTCAAACGCTGCATCTATACCCCCACCGAAGGCTTCCTCGAAGAAGCCGGTTGCCCGGAGTGCCGGCGCGAGGTAGGTGAGCCGTTGTTCGAAAGCCTGGAAGAGTGGATGCCGGCAGTGAGCGACAACTTTACCTGCCCACTGTGCGGCTTCGAAGA  >  NC_002947/1180023‑1180272                                                                                                                   |                                                                                                                                        agggCCATGCCCTGGCGGCAGGGGCGCACAAAGTTGCCCTGCAACCCGAAAAACTGCCCATCGGCGAACCGGTCAACGGCCAGCAAGTGATTCTAAAACTCTGCAACAATACCCCCACCCAAGCCTTCCTCgaagc                                                                                                                    >  2:219856/2‑135 (MQ=255)          cccTGGCGGCAGGGGCGCGCAAAGTGGCCCTGCACCCGGAAAAACTGCCCTTCGGCGAACCGGTCAACGGCCTGGAAGTGATGCTCAAACGCTGCATCTATACCCCCACCGAAGGCTTCCTCGAAGAAGCCGGTTg                                                                                                           >  1:188636/1‑136 (MQ=255)          cccTGGCGGCAGGGGCGCGCAAAGTGGCCATGCACCCGGAAAAACTGCCCTTCGGCGAACCGGTCAACGGCCTGGAAGTGATGCTCAAACGCTGCATCTATACCCCCACCGAAGGCTTCCTCGAAGAAGCCGGTTg                                                                                                           >  2:333178/1‑136 (MQ=255)                      gggCGCGCATAGTGTCCCTGCTAACGTAAAAACTGACCTGTGGAGAACCGTTCAACGGCCGGGATGGGATGCTCAAACGCTGCATCTATACCCCCACCTAAGGCTTCCTCTAATCATCCGGTTGCCCGTAGTGACg                                                                                               <  2:321071/136‑1 (MQ=255)                      gggCGCGCAAAGTGGCCCTGCACCCGGAAAAACTGCCCTTCGGCGAACCGGTCAACGGCCTGGAAGTGATGCTCAAACGCTGCATCTATACCCCCACCGAAGGCTTCCTCGAAGAAGCCGGTTGCCCGGAGTGCCg                                                                                               <  2:2008/136‑1 (MQ=255)                        gcgcgcAAAGTGGCCCTGCACCCGGAAAAACTGCCCTTCGGCGAACCGGTCAACGGCCTGGAAGTGATGCTCAAACGCTGCATCTATACCCCCACCGAAGGCTTCCTCGAAGAAGCCGGTTGCCCGGAGTGCCg                                                                                               >  2:515074/1‑134 (MQ=255)                        gcgcgcAAAGTGGCCCTGCACCCGGAAAAACTGCCCTTCGGCGAACCGGTCAACGGCCTGGAAGTGATGCTCAAACGCTGCATCTATACCCCCACCGAAGGCTTCCTCGAAGAAGCCGGTTGCCCGGAGTGCCg                                                                                               <  1:515074/134‑1 (MQ=255)                              aaaGTGGCCCTGCACCCGGAAAAACTGCCCTTCGGCGAACCGGTCAACGGCCTGGAAGTGATGCTCAAACGCTGCATCTATACCCCCCCCGAAGGCTTCCTCGAAGAAGCCGGTTGCCCGGAGTGCCGGCGCGagg                                                                                       >  1:47291/1‑136 (MQ=255)                              aaaGTGGCCCTGCACCCGGAAAAACTGCCCTTCGGCGAACCGGTCAACGGCCTGGAAGTGATGCTCAAACGCTGCATCTATACCCCCACCGAAGGCTTCCTCGAAGAAGCCGGTTGCCCGGAGTGCCGGCGCGagg                                                                                       >  1:68282/1‑136 (MQ=255)                              aaaGTGGCCCTGCACCCGGAAAAACTGCCCTTCGGCGAACCGGTCAACGGCCTGGAAGTGATGCTCAAACGCTGCATCTATACCCACACCGAAGGCTTCCTCGAAGAAGCCGGTTGCCCGGAGTGCCGGCGCGagg                                                                                       >  2:527212/1‑136 (MQ=255)                              aaaGTGGCCCTGCACCCGGAAAAACTGCCCTTCGGCGAACCGGTCAACGGCCTGGAAGTGATGCTCAAACGCTGCATCTATACCCACACCGAAGGCTTCCTCGAAGAAGCCGGTTGCCCGGAGTGCCGGCGCGagg                                                                                       >  2:527214/1‑136 (MQ=255)                                              cGGAAAAACTGCCCTTCGGCGAACCGGTCAACGGCCTGGAAGTGATGCTCAAACGCTGCATCTATACCCCCACCGAAGGCTTCCTCGAAGAAGCCGGTTGCGCGGAGTGCCGGCGCGAGGTAGGTGAGCCgttgtt                                                                       <  1:219856/136‑1 (MQ=255)                                                            ttCGGCGAACCGGTCAACGGCCTGGAAGTGATGCTCAAACGCTGCATCTATACCCCCACCGAAGGCTTCCTCGAAGAAGCCGGTTGCCCGGAGTGCCGGCGCGAGGTAGGTGAGCCGTTGTTCGAAAGCCTGGAag                                                         >  2:539246/1‑136 (MQ=255)                                                             tCGGCGAACCGGTCAACGGCCTGGAAGTGATGCTCAAACTCTGCATCTATACCGACACCGAAGTCTTCCTCTAAGAAGCCGGTTGCCCGGAGTGCCGGCTCGAGGTAGGTGAGCCGTTGTTCGAAAGCCTGGAaga                                                        <  1:527214/136‑1 (MQ=255)                                                             tCGGCGAACCGGTCAACGGCCTGGAAGTGATGCTCAAACGCTGCATCTATACCCCCACCGAAGGCTTCCTCGAAGAAGCCGGTTGCCCGGAGTGCCGGCGCGAGGTAGGTGAGCCGTTGTTCGAAAGCCTGGAaga                                                        <  2:527090/136‑1 (MQ=255)                                                             tCGGCGAACCGGTCAACGGCCTGGAAGTGATGCTCAAACGCTGCATCTATACCCCCACCGAAGGCTTCCTCGAAGAAGCCGGTTGCCCGGAGTGCCGGCGCGAGGTAGGTGAGCCGTTGTTCGAAAGCCTGGAaga                                                        <  2:84418/136‑1 (MQ=255)                                                             tCGGCGAACCGGTCAACGGCCTGGAAGTGATGCTCAAACGCTGCATCTATACCCCCACCGAAGGCTTCCTCGAAGAAGCCGGTTGCCCGGAGTGCCGGCGCGAGGTAGGTGAGCCGTTGTTCGAAAGCCTGGAaga                                                        <  2:458602/136‑1 (MQ=255)                                                             tCGGCGAACCGGTCAACGGCCTGGAAGTGATGCTCAAACGCTGCATCTATACCCCCACCGAAGGCTTCCTCGAAGAAGCCGGTTGCCCGGAGTGCCGGCGCGAGGTAGGTGAGCCGTTGTTCGAAAGCCTGGAaga                                                        <  1:333178/136‑1 (MQ=255)                                                             tCGGCGAACCGGTCAACGGCCTGGAAGTGATGCTCAAACGCTGCATCTATACCCCCACCGAAGGCTTCCTCGAAGAAGCCGGTTGCCCGGAGTGCCGGCGCGAGGTAGGTGAGCCGTTGTTCGAAAGCCTGGAaga                                                        <  1:231909/136‑1 (MQ=255)                                                             tCGGCGAACCGGTCAACGGCCTGGAAGTGATGCTCAAACGCTGCATCTATACCCACACCGAAGGCTTCCTCGAAGAAGCCGGTTGCCCGGAGTGCCGGCGCGAGGTAGGTGAGCCGTTGTTCGAAAGCCTGGAaga                                                        <  1:527212/136‑1 (MQ=255)                                                                              ggCCTGGAAGTGATGCTCAAACGCTGCATCTATACCCCCACCGATGGCTTCCTCGAAGAAGCGGGTTGCCCGGAGTGCCGGCGCGAGGTAGGTGAGCCGTTGTTCGAAAGCCTGGAAGAGTGGATGCCGGCAGTGa                                       <  1:379924/136‑1 (MQ=255)                                                                                                  aCGCTGCATCTATAACCACACCGAAGGCTTCCTCGAAGAAGCCGCTTGCCCGGAGTGCAGGCGCGAGGAAGGTGAGCCGATGTTCGAAACCCAGGAAGAGTGGATCCCGGCAGTGAGCGACAACTTAACCTGCCCa                   >  1:514323/1‑136 (MQ=255)                                                                                                         aTCTATACCCCCACCGAAGGCTTCCTCGAAGAAGCCGGTTGCCCGGAGTGCCGGCGCGAGGTAGGTGAGCCgttgtt                                                                       >  2:344144/1‑77 (MQ=255)                                                                                                         aTCTATACCCCCACCGAAGGCTTCCTCGAAGAAGCCGGTTGCCCGGAGTGCCGGCGCGAGGTAGGTGAGCCgttgtt                                                                       <  1:344144/77‑1 (MQ=255)                                                                                                            tataCCCCCACCGAAGGCTTCCTCGAAGAAGCCGGTTGCCCGGAGTGCCGGCGCGAGGTAGGTGAGCCGTTGTTCGAAAGCCTGGAAGAGTGGATGCCGGCAGTGAGCGACAACTTTACCTGCCCACTGTGCGGCt         <  1:525479/136‑1 (MQ=255)                                                                                                            aataCCCCCACCGAAGGCTTCCTCGAAGAAGCCGGTTGCCCGGAGTGCCGGCGCGAGGTAGGTGAGCCGTTGTTCGAAAGCCTGGAAGAGTGGATGCCGGCAGTGAGCGACAACTTTACCTGCCCACTGTGCGGCt         <  1:48096/135‑1 (MQ=255)                                                                                                                   ccACCGAAGGCTTCCTCGAAGAAGCCGGTTGCCCGGAGTGCCGGCGCGAGGTAGGTGAGCCGTTGTTCGAAAGCCTGGAAGAGTGGATGCCGGCAGTGAGCGACAACTTTACCTGCCCACTGTGCGGCTTCGAAGa  >  2:108902/1‑136 (MQ=255)                                                                                                                   ccACCGAAGGCTTCCTCGAAGAAGCCGGTTGCCCGGAGTGCCGGCGCGAGGTAGGTGAGCCGTTGTTCGAAAGCCTGGAAGAGTGGATGCCGGCAGTGAGCGACAACTTTACCTGCCCACTGTGCGGCTTCGAAGa  >  2:359687/1‑136 (MQ=255)                                                                                                                   |                                                                                                                                        TGGGCCATGCCCTGGCGGCAGGGGCGCGCAAAGTGGCCCTGCACCCGGAAAAACTGCCCTTCGGCGAACCGGTCAACGGCCTGGAAGTGATGCTCAAACGCTGCATCTATACCCCCACCGAAGGCTTCCTCGAAGAAGCCGGTTGCCCGGAGTGCCGGCGCGAGGTAGGTGAGCCGTTGTTCGAAAGCCTGGAAGAGTGGATGCCGGCAGTGAGCGACAACTTTACCTGCCCACTGTGCGGCTTCGAAGA  >  NC_002947/1180023‑1180272

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