Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	4,342,644	C→A	13.8%	V97V (GTC→GTA)	PP_3812 →	lipase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	4,342,644	0	C	A	13.8%	86.0 / 3.3	29	V97V (GTC→GTA)	PP_3812	lipase
Reads supporting (aligned to +/- strand):  ref base C (14/11);  new base A (2/2);  total (16/13)
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.05e-01

AACGACGCCCGCGCGGTACTGACCGGTGCCCAGGCTTCGGTAAAGGTCGACCTATCGGGGGTCGAGGTGAAGGAACGCACGATTCAGGCCAATGGCCAGTCGATCAAGCTGCAAGTGGTGCGCCCGACTAATGTCAAAGGTGATCTGCCGGTGTTCATGTTCTTCCACGGCGGCGGTTGGGTGCTGGGCGACTTCCCGACCCACCAGCGACTGATCCGCGACCTGGTGGTCGGCTCGGGTGCGGTGGCGGTCTATGTGGACTAC  >  NC_002947/4342510‑4342773                                                                                                                                       |                                                                                                                                  aaCGACGCCCGCGCGGTACTGACCGGTGCCCAGGCTTCGGTAAAGGTCGACCTATCGGGGGTCGAGGTGAAGGAACGCACGATTCGGGCCAATGGCCAGTCGATCAAGCTGCAAGTGGTGCGCCCGACTAATGTCa                                                                                                                                  <  2:159875/136‑1 (MQ=255)          cgcgcgGTACTGACCGGTGCCCAGGCTTCGGTAAAGGTCGACCTATCGGGGGTCGAGGTGAAGGAACGCACGATTCAGGCCAATGGCCAGTCGATCAAGCTGCAAGTGGTGCGCCCGACTAATGTCAAAGGTGATc                                                                                                                         <  1:141439/136‑1 (MQ=255)          cgcgcgGTACTGACCGGTGCCCAGGCTTCGGTAAAGGTCGACCTATCGGGGGTCGAGGTGAAGGAACGCACGATTCAGGCCAATGGCCAGTCGATCAAGCTGCAAGTGGTGCGCCCGACTAATGTCAAAGGTGATc                                                                                                                         >  1:351848/1‑136 (MQ=255)          cgcgcgGTACTGACCGGTGCCCAGGCTTCGGTAAAGGTCGACCTATCGGGGGTCGAGGTGAAGGAACGCACGATTCAGGCCAATGGCCAGTCGATCAAGCTGCAAGTGGTGCGCCCGACTAATGTCAAAGGAGATc                                                                                                                         >  1:51559/1‑136 (MQ=255)                              ccAGGCTTCGGTAAAGGTCGACCTATCGGGGGTCGAGGTGAAGGAACGCACGATTCAGGCCAATGGCCAGTCGATCAAGCTGCAAGTGGTGCGCCCGACTAATGTCAAATGTGATCTGCCGGTGTTCATGttcttc                                                                                                     >  1:83017/1‑136 (MQ=255)                              ccAGGCTTCGGTAAAGGTCGACCTATCGGGGGTCGAGGTGAAGGAACGCACGATTCAGGCCAATGGCCAGTCGATCAAGCTGCAAGTGGTGCGCCCGACTAATGTCAAAGGTGATCTGCCGGTGTTCATGttcttc                                                                                                     >  2:296006/1‑136 (MQ=255)                              ccAGGCTTCGGTAAAGGTCGACCTATCGGGGGTCGAGGTGAAGGAACGCACGATTCAGGCCAATGGCCAGTCGATCAAGCTGCAAGTGGTGCGCCCGACTAATGTCAAAGGTGATCTGCCGGTGTTCATGttcttc                                                                                                     >  2:213628/1‑136 (MQ=255)                              ccAGGCTTCGGTAAAGGTCGACCTATCGGGGGTCGAGGTGAAGGAACGCACGATTCAGGCCAATGGCCAGTCGATCAAGCTGCAAGTGGTGCGCCCGACTAATGTAAAAGGTGATCTGCCGGTGTTCATGttcttc                                                                                                     >  2:352535/1‑136 (MQ=255)                              ccAGGCTTCGGTAAAGGTCGACCTATCGGGGGTCGAGGTGAAGGAACGCACGATTCAGGCCAATGGCCAGTCGATCAAGCTGCAAGTGGTGCGCCCGACTAATGTAAAAGGTGATCTGCCGGTGTTCATGttcttc                                                                                                     >  2:288984/1‑136 (MQ=255)                                   cTTCGGTAAAGGTCGACCTATCGGGGGTCGAGGTGAAGGAACGCACGATTCAGGCCAATGGCCAGTCGATCAAGCTGCAAGTGGTGCGCCCGACTAATGTCAAAGGTGATCTGCCGGTGTTCATGTTCTTCCAcgg                                                                                                <  2:351848/136‑1 (MQ=255)                                                     tATCGGGGGTCGAGGTGAAGGAACGCACGATTCAGGCCAATGGCCAGTCGATCAAGCTGCAAGTGGTGCGCCCGACTAATGTCAAAGGTGATCTGCCGGTGTTCATGTTCTTCAACGGCGGCGGTTGGGTGCTggg                                                                              >  1:320361/1‑136 (MQ=255)                                                        cGGGGGTCGAGGTGAAGGAACGCACGATTCAGGCCAATGGCCAGTCGATCAAGCTGCAAGTGGTGCGCCCGACTAATGTCAAAGGTGATCTGCCGGTGTTCATGTTCTTCCACGGCGGCGGTTGGGTGCTGGGCGa                                                                           <  1:250967/136‑1 (MQ=255)                                                          ggggTCGAGGTGAAGGAACGCACGATTCAGGCCAATGGCCAGTCGATCAAGCTGCAAGTGGTGCGCCCGACTAATGTCAAAGGTGATCTGCCGGTGTTCATGTTCTTCCACGGCGGCGGTTGGGTGCTGGGCGACt                                                                         >  2:229788/1‑136 (MQ=255)                                                          ggggTCGAGGTGAAGGAACGCACGATTCAGGCCAATGGCCAGTCGATCAAGCTGCAAGTGGTGCGCCCGACTAATGTCAAAGGTGATCTGCCGGTGTTCATGTTCTTCCACGGCGGCGGTTGGGTGCTGGGCGACt                                                                         <  2:160944/136‑1 (MQ=255)                                                                 aGGTGAAGGAACGCACGATTCAGGCCAATGGCCAGTCGATCAAGCTGCAAGTGGTGCGCCCGACTAATGTCAAAGGTGATCTGCCGGTGTTCATGTTCTTCCACGGCGGCGGTTGGGTGCTGGGCGACTTCCCGAc                                                                  <  1:132236/136‑1 (MQ=255)                                                                               aCGATTCAGGCCAATGGCCAGTCGATCAAGCTGCAAGTGGTGCGCCCGACTAATGTCAAAGGTGATCTGCCGGTGTTCATGTTCTTCCACGGCGGCGGTTGGGTGCTGGGCGACTTCCCGACCCACCAGCGACTGa                                                    >  2:387463/1‑136 (MQ=255)                                                                                 gATTCAGGCCAATGGCCAGTCGATCAAGCTGCAAGTGGTGCGCCCGACTAATGTCAAAGGTGATCTGCCGGTGTTCATGTTCTTCCACGGCGGCGGTTGGGTGCTGGGCGACTTCCCGACCCACCAGCGACTGATc                                                  >  1:390593/1‑136 (MQ=255)                                                                                            ttGGCCAGTCGATCAAGCTGCAAGTGGTGCGCCCGACTAATGTCAAAGGTGATCTGCCGGTGTTCATGTTCTTCCACGGCGGCGGTTGGGTGCTGGGCGACTTCCCGACCCACCAGCGACTGATCCGCGACCtggt                                       <  1:260587/135‑1 (MQ=255)                                                                                                ccAGTCGATCAAGCTGCAAGTGGTGCGCCCGACTAATGTCAAAGGTGATCTGCCGGTGTTCATGTTCTTCCACGGCGGCGGTTGGGTGCTGGGCGACTTCCCGACCCACCAGCGACTGATCCGCGACCTGGTGGTc                                   >  1:184961/1‑136 (MQ=255)                                                                                                     cGATCAAGCTGCAAGTGGTGCGCCCGACTAATGTCAAAGGTGATCTCCCGGTCTTCATGTTCTTCCACGGCGGCGGTTGGGTGCTGGGCGACTTCCCGACCCACCAGCGACTGATCCGCGACCTGGAGGACGGCTc                              >  1:36608/1‑136 (MQ=255)                                                                                                            gCTGCAAGTGGTGCGCCCGACTAATGTCAAAGGTGATCTGCCGGTGTTCATGTTCTTCCACGGCGGCGGTTGGGTGCTGGGCGACTTCCCGACCCACCAGCGACTGATCCGCGACCTGGTGGTCGGCTCGGGTGCg                       <  1:296006/136‑1 (MQ=255)                                                                                                              ggCAAGTGGTGCGCCCGACTAATGTCAAAGGTGATCTGCCGGTGTTCATGTTCTTCCACGGCGGCGGTTGGGTGCTTGGCGACTTCCCGACCCACCAGCGACTGATCCGCGACCTGGTGGTCGGCTCGGGTGCGGt                     <  1:362633/135‑1 (MQ=255)                                                                                                                         cgcCCGACTAATGTCAAAGGTGATCTGCCGGTGTTCATGTTCTTCCACGGCGGCGGTTGGGTGCTGGGCGACTTCCCGACCCACCAGCGACTGATCCGCGACCTGGTGGTCGGCTCGGGTGCGGTGGCGGTCTAtg          >  1:286485/1‑136 (MQ=255)                                                                                                                          gcCCGACTAATGTCAAAGGTGATCTGCCGGTGTTCATGTTCTTCCACGGCGGCGGTTGGGTGCTGGGCGACTTCCCGACCCACCAGCGACTGATCCGCGACCTGGTGGTCGGCTCGGGTGCGGTGGCGGTCTAtgt         >  1:243535/1‑136 (MQ=255)                                                                                                                          gcCCGACTAATGTCAAAGGTGATCTGCCGGTGTTCATGTTCTTCCACGGCGGCGGTTGGGTGCTGGGCGACTTCCCGACCCACCAGCGACTGATCCGCGACCTGGTGGTCGGCTCGGGTGCGGTGGCGGTCTAtgt         >  1:199095/1‑136 (MQ=255)                                                                                                                           cccGACTAATGTCAAAGGTGATCTGCCGGTGTTCATGTTCTTCCACGGCGGCGGTTGGGTGCTGGGCGACTTCCCGACCCACCAGCGACTGATCCGCGACCTGGTGGTCGGCTCGGGTGCGGTGGCGGTCTAtgtg        <  1:222613/136‑1 (MQ=255)                                                                                                                               aCTAATGTAAAAGGTGATCTGCCGGTGTTCATGTTCTTCCACGGCGGCGGTTGGGTGCTGGGCGACTTCCCGACCCACCAGCGACTGATCCGCGACCTGGTGGTCGGCTCGGGTGCGGTGGCGGTCTATGTGGACt    <  1:352535/136‑1 (MQ=255)                                                                                                                               aCTAATGTAAAAGGTGATCTGCCGGTGTTCATGTTCTTCCACGGCGGCGGTTGGGTGCTGGGCGACTTCCCGACCCACCAGCGACTGATCCGCGACCTGGTGGTCGGCTCGGGTGCGGTGGCGGTCTATGTGGACt    <  1:288984/136‑1 (MQ=255)                                                                                                                                 tAATGTCAAAGGTGATCTGCCGGTGTTCATGTTCTTCCACGGCGGCGGTTGGGTGCTGGGCGACTTCCCGACCCACCAGCGACTGATCCGCGACCTGGTGGTCGGCTCGGGTGCGGTGGCGGTCTATGTGGACTac  <  1:213628/136‑1 (MQ=255)                                                                                                                                       |                                                                                                                                  AACGACGCCCGCGCGGTACTGACCGGTGCCCAGGCTTCGGTAAAGGTCGACCTATCGGGGGTCGAGGTGAAGGAACGCACGATTCAGGCCAATGGCCAGTCGATCAAGCTGCAAGTGGTGCGCCCGACTAATGTCAAAGGTGATCTGCCGGTGTTCATGTTCTTCCACGGCGGCGGTTGGGTGCTGGGCGACTTCCCGACCCACCAGCGACTGATCCGCGACCTGGTGGTCGGCTCGGGTGCGGTGGCGGTCTATGTGGACTAC  >  NC_002947/4342510‑4342773

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