Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	4,250,135	C→T	13.8%	R357H (CGC→CAC)	PP_3725 ←	acyl‑CoA dehydrogenase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	4,250,135	0	C	T	13.8%	67.7 / 3.2	29	R357H (CGC→CAC)	PP_3725	acyl‑CoA dehydrogenase
Reads supporting (aligned to +/- strand):  ref base C (11/14);  new base T (2/2);  total (13/16)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.55e-01

GCAGGCGCTTGAGGTTGTCACGAATCGATTGCGACATGAGTGTTCTCCTGGCGCTCAGCGCAGCAGCGGGCGCAGCAGTTCGCGGGCGATGATATGCCGCTGGATTTCCGAAGTGCCCTCCCAGATGCGCTCGATACGCGCGTTGCGCCAGATGCGCTCAACCGGCCCTTCATCCATCAGGCCCATGCCACCGAAGATCTGCACCGCTTCGTCGGCGACCTTGCCCAGCACCTCGCTGGCAAACAGCTTGGCCATGCCGGCCTCACCGTCGGTCATGTT  >  NC_002947/4249990‑4250268                                                                                                                                                  |                                                                                                                                      gcagGCGCTTGAGGTTGTCACGAATCGATTGCGACATGAGTGTTCTCCTGGCGCTCAGCGCAGCAGCGGGCGCAGCAGTTCGCGGGCGATGATATGCCGCTGGATTTCCGAAGTGCCCTCCCAGATGCGCTCGATACGCGCGTTGCGcc                                                                                                                                    <  6:315732/149‑1 (MQ=255)           gAGGTTGTCACGAATCGATTGCGACATGAGTGTTCTCCTGGCGCTCAGCGCAGCAGCGGGCGCAGCAGTTCGCGGGCGATGATATGCCGCTGGATTTCCGAAGTGCCCTCCCAGATGCGCTCGATACGCGCGTTGCGCCAGATgcgc                                                                                                                            <  1:135701/147‑1 (MQ=255)           gAGGTTGTCACGAATCGATTGCGACATGAGTGTTCTCCTGGCGCTCAGCGCAGCAGCGGGCGCAGCAGTTCGCGGGCGATGATATGCCGCTGGATTTCCGAAGTGCCCTCCCAGATGCGCTCGATACGCGCGTTGCGCCAGATGCGCt                                                                                                                           >  6:81376/1‑148 (MQ=255)              gTTGTCACGAATCGATTGCGACATGAGTGTTCTCCTGGCGCTCAGCGCAGCAGCGGGCGCAGCAGTTCGCGGGCGATGATATGCCGCTGGATTTCCGAAGTGCCCTCCCAGATGCGCTCGATACGCGCGTTGCGCCAGATGCGCTCAAc                                                                                                                       <  6:344095/149‑1 (MQ=255)                  tCACGAATCGATTGCGACATGAGTGTTCTCCTGGCGCTCAGCGCAGCAGCGGGCGCAGCAGTTCGCGGGCGATGATATGCCGCTGGATTTCCGAAGTGCCCTCCCAGATGCGCTCGATACGCGCGTTGCGCCAGATGCGCTCAACCgg                                                                                                                    <  6:189869/148‑1 (MQ=255)                    aCGAATCGATTGCGACATGAGTGTTCTCCTGGCGCTCAGCGCAGCAGCGGGCGCAGCAGTTCGCGGGCGATGATATGCCGCTGGATATCCGAAGTCACCTAACAGATGCGCTCGATACGCGCGTTGCGCCAGATGCGCTCAACCGGccc                                                                                                                 >  5:174388/1‑149 (MQ=255)                           gATTGCGACATGAGTGTTCTCCTGGCGCTCAGCGCAGCAGCGGGCGCAGCAGTTCGCGGGCGATGATATGCCGCTGGATTTCCGAAGTGCCCTCCCAGATGCGCTCGATACGCGCGTTGCGCCAGATGCGCTCAACCGGCCCTTcatcc                                                                                                          <  5:81376/149‑1 (MQ=255)                            aTTGCGACATGAGTGTTCTCCTGGCGCTCAGCGCAGCAGCGGGCGCAGCAGTTCGCGGGCGATGATATGCCGCTGGATTTCCGAAGTGCCCTCCCAGATGCGCTCGATACGCGCGTTGCGCCAGATGCGCTCAACCGGCCCTTcatcca                                                                                                         >  4:171866/1‑149 (MQ=255)                               gCGACATGAGTGTTCTCCTGGCGCTCAGCGCAGCAGCGGGCGCAGCAGTTCGCGGGCGATGATATGCCGCTGGATTTCCGAAGTGCCCTCCCAGATGCGCTCGATACGCGCGTTGCGCCAGATGCGCTCAACCGGCCCTTcatccatc                                                                                                       <  3:171866/148‑1 (MQ=255)                                 gACATGAGTGTTCTCCTGGCGCTCAGCGCAGCAGCGGGCGCAGCAGTTCGCGGGCGATGATATGCCGCTGGATTTCCGAAGTGCCCTCCCAGATGCGCTCGATACGCGCGTTGCGCCAGATGCGCTCAACCGGCCCTTCATCCATCAgg                                                                                                    >  1:55905/1‑149 (MQ=255)                                 gACATGAGTGTTCTCCTGGCGCTCAGCGCAGCAGCGGGCGCAGCAGTTCGCGGGCGATGATATGCCGCTGGATTTCCGAAGTGCCCTCCCAGATGCGCTCGATACGCGCGTTGCGCCAGATGCGCTCAACCGGCCCTTCATCCATCAgg                                                                                                    <  2:55905/149‑1 (MQ=255)                                                   gcgcTCAGCGCAGCAGCGGGCTCAGCAGTTCGCGGGCGATGATATGCCGCTGCTTTTCCGAAGTGCCCTCCCAGATGCGCTCGATACGCGCGTTGCGCCAGATGCGCTCAACCGGCCCTTCATCCATCAGGCCCATGCCACCGAAGAt                                                                                   <  1:294881/148‑1 (MQ=255)                                                             cagcagCGGGCGCAGCAGTTCGCGGGCGATGATATGCCGCTGGATTTCCGAAGTGCCCTCCCAGATGCGCTCGATACGCGCGTTGCGCCAGCTGCGCTCAACCGGCCCTTCATCCATCAGGCCCATGCAACCGAAGATCTGCACCGCtt                                                                        >  4:310956/1‑149 (MQ=255)                                                                          agcagTTCGCGGGCGATGATATGCCGCTGGATTTCCGAAGTGCCCTCCCAGATGCGCTCGATACGCGCGTTGCGCCAGATGCGCTCAACCGGCCCTTCATCCATCAGGCCCATGCCACCGAAGATCTGCACCGCTTCGTCGGCGACCtt                                                           >  8:65461/1‑149 (MQ=255)                                                                             agTGCGCGGGCGATGATATGCCGCTGGATTTCCGAAGTGCCCTCCCAGATGCGCTCGATACGCGCGTTGCGCCAGATGCGCTAAATCGTTTTTTTTTTTAATTTTACGTTGACATCGTAGATCTTCACCGCTTCGTCGGCGACCTTGcc                                                        <  6:174388/149‑1 (MQ=255)                                                                                 cgcgGGCGATGATATGCCGCTGGATTTCCGAAGTGCCCTCCCAGATGCGCTCGATACGCGCGTTGCGCCAGATGCGCTCAACCGGCCCTTCATCCATCAGGCCCATGCCACCGAAGATCTGc                                                                               <  4:298493/122‑1 (MQ=255)                                                                                 cgcgGGCGATGATATGCCGCTGGATTTCCGAAGTGCCCTCCCAGATGCGCTCGATACGCGCGTTGCGCCAGATGCGCTCAACCGGCCCTTCATCCATCAGGCCCATGCCACCGAAGATCTGc                                                                               >  3:298493/1‑122 (MQ=255)                                                                                          tgacATGCCCCTGGATTTCCGACGGGCCCTCCCACATGCGCTCGATACGCGCGTTGCGCCAGATGCGCTCAACCGTTTTTTTTTTTAATTATACGTTGACAACCTAGATCTTCACCGCTTCGTCGGCGACCTTGCCCAGCACCTCGCTg                                           <  6:294295/149‑1 (MQ=255)                                                                                             tatGCCGCTGGATTTCCGAAGTGCCCTCCCAGATGCGCTCGATACGCGCGTTGTGCCAGATGCGCTCAACCGGCCCTTCATCCATCAGGCCCATGCCACCGAAGATCTGCACCGCTTCGTCGGCGAcc                                                             >  4:132942/1‑128 (MQ=255)                                                                                             tatGCCGCTGGATTTCCGAAGTGCCCTCCCAGATGCGCTCGATACGCGCGTTGTGCCAGATGCGCTCAACCGGCCCTTCATCCATCAGGCCCATGCCACCGAAGATCTGCACCGCTTCGTCGGCGAcc                                                             >  4:127836/1‑128 (MQ=255)                                                                                             tatGCCGCTGGATTTCCGAAGTGCCCTCCCAGATGCGCTCGATACGCGCGTTGTGCCAGATGCGCTCAACCGGCCCTTCATCCATCAGGCCCATGCCACCGAAGATCTGCACCGCTTCGTCGGCGAcc                                                             <  3:132942/128‑1 (MQ=255)                                                                                             tatGCCGCTGGATTTCCGAAGTGCCCTCCCAGATGCGCTCGATACGCGCGTTGTGCCAGATGCGCTCAACCGGCCCTTCATCCATCAGGCCCATGCCACCGAAGATCTGCACCGCTTCGTCGGCGAcc                                                             <  3:127836/128‑1 (MQ=255)                                                                                                    cTGGATTTCCGAAGTGCCCTCCCAGATGCGCTCGATACGCGCGTTGCGCCAGATGCGCTCAACCGGCCCTTCATCCATCAGGCCCATGCCACCGAAGATCTGCACCGCTTCGTCGGCGACCTTGCCCAGCACCTCGCTGGCAAACAGCt                                 >  6:7660/1‑149 (MQ=255)                                                                                                       gATTTCCGAAGTGCCCTCCCAGATGCGCTCGATACGCGCGTTGCGCCAGATGCGCTCAACCGGCCCTTCATCCATCAGGCCCATGCCACCGAAGATCTGCACCGCTTCGTCGGCGACCTTGCCCAGCACCTCGCTGGCAAACAGCTTgg                              <  5:211499/149‑1 (MQ=255)                                                                                                       gATTTCCGAAGTGCCCTCCCAGATGCGCTCGATACGCGCGTTGCGCCAGATGCGCTCAACCGGCCCTTCATCCATCAGGCCCATGCCACCGAAGATCTGCACCGCTTCGTCGGCGACCTTGCCCAGCACCTCGCTGGCAAACAGCTTgg                              >  6:211499/1‑149 (MQ=255)                                                                                                          ttCCGAAGTGCCCTCCCAGATGCGCTCGATACGCGCGTTGCGCCAGATGCGCTCAACCGGCCCTTCATCCATCAGGCCCATGCCACCGAAGATCTGCACCGCTTCGTCGGCGACCTTGCCCAGCACCTCGCTGGCAAACAGCTTGGCCa                           >  1:308970/1‑149 (MQ=255)                                                                                                                       tCCCAGATGCGCTCGATACGCGCGTTGCGCCAGATGCGCTCAACCGGCCCTTCATCCATCAGGCCCATGCCACCGAAGATCTGCACCGCTTCGTCGGCGACCTTGCCCAGCACCTCGCTGGCAAACAGCTTGGCCATGCCGGCCTCAcc              >  4:172225/1‑149 (MQ=255)                                                                                                                                 gcTCGATACGCGCGTTGCGCCAGATGCGCTCAACCGGCCCTTCATCCATCAGGCCCATGCCACCGAAGATCTGCACCGCTTCGTCGGCGACCTTGCCCAGCACCTCGCTGGCAAACAGCTTGGCCATGCCGGCCTCACCGTCGGTCATg    <  5:7660/149‑1 (MQ=255)                                                                                                                                     gATACGCGCGTTGCGCCAGATGCGCTCAACCGGCCCTTCATCCATCAGGCCCATGCCACCGAAGATCTGCACCGCTTCGTCGGCGACCTTGCCCAGCACCTCGCTGGCAAACAGCTTGGCCATGCCGGCCTCACCGTCGGTCATGtt  <  2:308970/147‑1 (MQ=255)                                                                                                                                                  |                                                                                                                                      GCAGGCGCTTGAGGTTGTCACGAATCGATTGCGACATGAGTGTTCTCCTGGCGCTCAGCGCAGCAGCGGGCGCAGCAGTTCGCGGGCGATGATATGCCGCTGGATTTCCGAAGTGCCCTCCCAGATGCGCTCGATACGCGCGTTGCGCCAGATGCGCTCAACCGGCCCTTCATCCATCAGGCCCATGCCACCGAAGATCTGCACCGCTTCGTCGGCGACCTTGCCCAGCACCTCGCTGGCAAACAGCTTGGCCATGCCGGCCTCACCGTCGGTCATGTT  >  NC_002947/4249990‑4250268

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