Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_002947	2,613,426	A→G	T134A (ACC→GCC)	PP_2287 →	internal core protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	2,613,426	0	A	G	91.7%	71.8 / ‑1.5	25	T134A (ACC→GCC)	PP_2287	internal core protein
Reads supporting (aligned to +/- strand):  ref base A (1/1);  new base G (11/11);  total (13/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.92e-01

CAACATGGATATGGCCCTGGCCGACTACAACGGCGGTCCTCGGGCTGCCAAGGCGCTGGCCGCAGGGAGTCCCTGGCCGGAAACACGGGACTACCTGAGCAAGTTCAACGGCACCACCTCCGCAGCAGGGCAGACCCCCTCGGCAAAGCCTGACCCCCTGTCCTCCGCTTTCACCTCCGGCGATTCGGTCAACCCGGATCAGAGCGCTTCCGCCTCGGACCTGGTGCAGATGCAGGCCCAGCAGCGGGCAGAGAACGGG  >  NC_002947/2613293‑2613551                                                                                                                                      |                                                                                                                              cAACATGGATATGGCCCTGGCCGACTACAACGGCGGTCCTCGGGCTGCCAAGGCGCTGGCCGCAGGGAGTCCCTGGCCGGAAACACGGGACTACCTGAGCAAGTTCAACGGCACCACCTCCGCAGCAGGGCAGGcc                                                                                                                             >  1:926993/1‑136 (MQ=255) cAACATGGATATGGCCCTGGACGACTACAACTGCTGTCCTCGGGCTGCCAAGGCGCAGGCCGCAGGGATTCCCAGGCCGGAAACACGGCAATACCTCAGCAAGTTCAACTGCACCACCTCAGCAGCAGGGCAGTcc                                                                                                                             >  2:797094/1‑136 (MQ=255)       ggATATGGCCCTGGCCGACTACAACGGCGGTCCTCGGGCTGCCAAGGCGCTGGCCGCAGGGAGTCCCTGGCCGGAAACACGGGACTACCTGAGCAAGTTCAACGGCACCACCTCCGCAGCAGGGCAGGCCCCCTCg                                                                                                                       <  1:933119/136‑1 (MQ=255)                      cGACTACAACGGCGGTCCTCGGGCTGCCAAGGCGCTGGCCGCAGGGAGTCCCTGGCCGGAAACACGGGACTACCTGAGCAAGTTCAACGGCACCACCTCCGCAGCAGGGCAGGCCCCCTCGGCAAAGCCTGAcccc                                                                                                        >  2:815293/1‑136 (MQ=255)                      cGACTACAACGGCGGTCCTCGGGCTGCCAAGGCGCTGGCCGCAGGGAGTCCCTGGCCGGAAACACGGGACTACCTGAGCAAGTTCAACGGCACCACCTCCGCAGCAGGGCAGGCCCCCTCGGCAAAGCCTGAcccc                                                                                                        >  2:518369/1‑136 (MQ=255)                       gACTACAACGGCGGTCCTCGGGCTGCCAAGGCGCTGGCCGCAGGGAGTCCCTGGCCGGAAACACGGGACTACCTGAGCAAGTTCAACGGCACCACCTCCGCAGCAGGGCAGGCCCCCTCGGCAAAGCCTGAccccc                                                                                                       >  1:586999/1‑136 (MQ=255)                           aCAACGGCGGTCCTCGGGCTGCCAAGGCGCTGGCCGCAGGGAGTCCCTGGCCGGAAACACGGGACTACCTGAGCAAGTTCAACGGCACCACCTCCGCAGCAGGGCAGGCCCCCTCGGCAAAGCCTGACCCCCTGtc                                                                                                   >  2:965527/1‑136 (MQ=255)                                  cggTCCTCGGGCTGCCAAGGCGCTGGCCGCAGGGAGTCCCTGGCCGGAAACACGGGACTACCTGAGCAAGTTCAACGGCACCACCTCCGCAGCAGGGCAGGCCCCCTCGGCAAAGCCTGCCCCCCTGTCCTCCGCt                                                                                            >  2:653445/1‑136 (MQ=255)                                  cggTCCTCGGGCTGCCAAGGCGCTGGCCGCAGGGAGTCCCTGGCCGGAAACACGGGACTACCTGAGCAAGTTCAACGGCACCACCTCCGCAGCAGGGCAGGCCCCCTCGGCAAAGCCTGACCCCCTGTCCTCCGCt                                                                                            >  2:660001/1‑136 (MQ=255)                                   ggTCCTCGGGCTGCCAAGGCGCTGGCCGCAGGGAGTCCCTGGCCGGAAACACGGGACTACCTGAGCAAGTTCAACGGCACCACCTCCGCAGCAGGGCAGGCCCCCTCGGCAAAGCCTGACCCCCTGTCCTCCGCtt                                                                                           <  1:219009/136‑1 (MQ=255)                                       cTCGGGCTGCCAAGGCGCTGGCCGCAGGGAGTCCCTGGCCGGAAACACGGGACTACCTGAGCAAGTTCAACGGCACCACCTCCGCAGCAGGGCAGGCCCCCTCGGCAAAGCCTGACCCCCTGTCCTCCGCTTTCAc                                                                                       >  1:341455/1‑136 (MQ=255)                                              tGCCAAGGCGCTTGCCTCAGGGAGTCACTGGCCGGAAACACGGGACTACCTGAGCAAGTTCAACGGCAACACCTCCGCAGCAGGGCAGGCCCCCTCGGCAAAGCCTGACCCCCTGTCCTCCGCTTTCACCTCCGGc                                                                                <  2:661405/136‑1 (MQ=255)                                              tGCCAAGGCGCTGGCCGCAGGGAGTCCCTGGCCGGAAACACGGGACTACCTGAGCAAGTTCAACGGCACCACCTCCGCAGCAGGGCAGGCCCCCTCGGCAAAGCCTGACCTCCTGTCCTCCGCTTTCACCTCCGGc                                                                                <  1:182495/136‑1 (MQ=255)                                               gCCAAGGCGCTGGCCGCAGGGAGTCCCTGGCCGGAAACACGGGACTACCTGAGCAAGTTCAACGGCACCACCTCCGCAGCAGGGCAGACCCCCTCGGCAAAGCCTGACCCCCTGTCCTCCGCTTTCAc                                                                                       >  2:1186152/1‑128 (MQ=255)                                               gCCAAGGCGCTGGCCGCAGGGAGTCCCTGGCCGGAAACACGGGACTACCTGAGCAAGTTCAACGGCACCACCTCCGCAGCAGGGCAGACCCCCTCGGCAAAGCCTGACCCCCTGTCCTCCGCTTTCAc                                                                                       <  1:1186152/128‑1 (MQ=255)                                                             cGCAGGGAGTCCCTGGCCGGAAACACGGGACTACCTGAGCAAGTTCAACGGCACCACCTCCGCAGCAGGGCAGGCCCCCTCGGCAAAGCCTGACCCCCTGTCCTCCGCTTTCACCTCCGGCGATTCGGTCAACCCg                                                                 >  2:228709/1‑136 (MQ=255)                                                                             ccGGAAACACGGGACTACCTGAGCAAGTTCAACGGCACCACCTCCGCAGCAGGGCAGGCCCCCTCGGCAAAGCCTGACCCCCTGTCCTCCGCTTTCACCTCCGGCGATTCGGTCAACCCGGATCAGAGCGCTTCCg                                                 <  2:341455/136‑1 (MQ=255)                                                                             ccGGAAACACGCGACTACCTGAGCAAGTTCAACGGCACCACCTCCGCAGCAGGGCAGGCCCCCTCGGCAAAGCCTGACCCCCTGTCCTCCGCTTTCACCTCCGGCGATTCGGTCAACCCGGATCAGAGCGCTTCCg                                                 <  1:473634/136‑1 (MQ=255)                                                                              cGGAAACACGGGACTACCTGAGCAAGTTCAACGGCACCACCTCCGCAGCAGGGCAGGCCCCCTCGGCAAAGCCTGACCCCCTGTCCTCCGCTTTCACCTCCGGCGATTCGGTCAACCCGGATCAGAGCGCATCCGc                                                >  2:875894/1‑136 (MQ=255)                                                                                 aaaCACGGGACTACCTGAGCAAGTTCAACGGCACCACCTCCGCAGCAGGGCAGGCCCCCTCGGCAAAGCCTGACCCCCTGTCCTCCGCTTTCACCTCCGGCGATTCGGTCAACCCGGATCAGAGCGCTTCCGCCTc                                             <  2:970646/136‑1 (MQ=255)                                                                                              ccTGAGCAAGTTCAACGGCACCACCTCCGCAGCAGGGCAGGCCCCCTCGGCAAAGCCTGACCCCCTGTCCTCCGCTTTCACCTCCGGCGATTCGGTCAACCCGGATCAGAGCGCTTCCGCCTCGGACCTGGTGCAg                                >  2:675489/1‑136 (MQ=255)                                                                                                        ttCTACGGCACCACCTCCGCAGCAGGGCAGGCCCCCTCGGCAAAGCCTGACCCCCTGTCCTCCGCTTTCACCTCCGGCGATTCGGTCAACCCGGATCAGAGCGCTTCCGCCTCGGACCTGGTGCAGATGCAGGccc                      <  2:314832/136‑1 (MQ=255)                                                                                                                          gcagcaGGGCAGGCCCCCTCGGCAAAGCCTGACCCCCTGTCCTCCGCTTTCACCTCCGGCGATTCGGTCAACCCGGATCAGAGCGCTTCCGCCTCGGACCTGGTGCAGATGCAGGCCCAGCAGCGGGCAGAGAACg    <  1:518369/136‑1 (MQ=255)                                                                                                                            agcagGGCAGGCCCCCTCGGCAAAGCCTGACCCCCTGTCCTCCGCTTTCACCTCCGGCGATTCGGTCAACCCGGATCAGAGCGCTTCCGCCTCGGACCTGGTGCAGATGCAGGCCCAGCAGCGGGCAGAGAACggg  <  2:926993/136‑1 (MQ=255)                                                                                                                            agcagGGCAGGCCCCCTCGGCAAAGCCTGACCCCCTGTCCTCCGCTTTCACCTCCGGCGATTCGGTCAACCCGGATCAGAGCGCTTCCGCCTCGGACCTGGTGCAGATGCAGGCCCAGCAGCGGGCAGAGAACggg  <  1:188775/136‑1 (MQ=255)                                                                                                                                      |                                                                                                                              CAACATGGATATGGCCCTGGCCGACTACAACGGCGGTCCTCGGGCTGCCAAGGCGCTGGCCGCAGGGAGTCCCTGGCCGGAAACACGGGACTACCTGAGCAAGTTCAACGGCACCACCTCCGCAGCAGGGCAGACCCCCTCGGCAAAGCCTGACCCCCTGTCCTCCGCTTTCACCTCCGGCGATTCGGTCAACCCGGATCAGAGCGCTTCCGCCTCGGACCTGGTGCAGATGCAGGCCCAGCAGCGGGCAGAGAACGGG  >  NC_002947/2613293‑2613551

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