Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	3,996,941	G→T	19.4%	P89H (CCT→CAT)	PP_3525 ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	3,996,941	0	G	T	19.4%	68.9 / 9.1	31	P89H (CCT→CAT)	PP_3525	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (13/12);  new base T (3/3);  total (16/15)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.79e-01

CACGCCCTGCTGCAGTTGCAACCACTGCACCTGGCCCTCCTGCCCGGGGTACCGGCTGGCCACTCGCCCATCACGGTCGAGCAGCACCCGGTAGCTGTAATCGCGCATGGCCGGGACCGCGAACAGTTTGCTGATCAGCGCAGGCATGCGCTGAATGTCGGCCACGAATACAGCGGCGCGTGCTTCCAGGTAGCCCTTGGGCTGCTCGGCCAGTGCCGCCTTGACCAGCTTGGCACCGTCCATATCGCGCGCCACCAGTAAAATCTGCGTGTCGGCACCAAGGCTGT  >  NC_002947/3996799‑3997085                                                                                                                                               |                                                                                                                                                 cACGCCCTGCTGCAGTTGCAACCACTGCACCTGGCCCTCCTGCCCGGGGTACCGGCTGGCCACTCGCCCATCACGGTCGAGCAGCACCCGGTAGCTGTAATCGCGCATGCCGGGGACCGCGAACAGTTTGCTGATCAGCGCAGGCATgcg                                                                                                                                           >  1:383864/1‑150 (MQ=255)        tgctgcAGTTGCAACCACTGCACCTGGCCCTCCTGCCCGGGGTACCGGCTGGCCACTCGCCCATCACGGTCGAGCAGCACCCGGTAGCTGTAATCGCGCATGGCCGGGACCGCGAACAGTTTGCTGATCAGCGCAGGCATGCGCTGAATg                                                                                                                                    >  1:338867/1‑150 (MQ=255)               gTTGCAACCACTGCACCTGGCCCTCCTGCCCGGGGTACCGGCTGGCCACTCGCCCATCACGGTCGAGCAGCACCCGGTAGCTGTAATCGCGCATGGCCGGGACCGCGAACAGTTTGCTGATCAGCGCAGGCATGCGCTGAATGTCg                                                                                                                                 >  1:735890/1‑146 (MQ=255)               gTTGCAACCACTGCACCTGGCCCTCCTGCCCGGGGTACCGGCTGGCCACTCGCCCATCACGGTCGAGCAGCACCCGGTAGCTGTAATCGCGCATGGCCGGGACCGCGAACAGTTTGCTGATCAGCGCAGGCATGCGCTGAATGTCg                                                                                                                                 <  2:735890/146‑1 (MQ=255)                   cAACCACTGCACCTGGCCCTCCTGCCCGGGGTACCGGCTGGCCACTCGCCCATCACGGTCGAGCAGCACCCGGTAGCTGTAATCGCGCATGGCCGGGACCGCGAACAGTTTGCTGATCAGCGCAGGCATGCGCTGAATGTCGGCCACGaa                                                                                                                         >  1:29148/1‑150 (MQ=255)                        aCTGCACCTGGCCCTCCTGCCCGGGGTACCGGCTGGCCACTCGCCCATAACGGTCGAGCAGCACCCGGTAGCTGTAATCGCGCATGGCCGGGACCGCGAACAGTTTGCTGATCAGCGCAGGCATGCGCTGAATGTCGGCCACGAATACAg                                                                                                                    >  1:429515/1‑150 (MQ=255)                                   cccTCCTGCCCGGGGTACCGGCTGGCCACTCGCCCATCACGGTCGAGCAGCACCCGGTAGCTGTAATCGCGCATGGCCGGGACCGCGAACAGTTTGCTGATCAGCGCAGGCATGCGCTGAATGTCGGCCACGAATACAGCGGCGCGTGCt                                                                                                         <  1:257359/150‑1 (MQ=255)                                    cctcctGCCCGGGGTACCGGCTGGCCACTCGCCCATCACGGTCGAGCAGCACCCGGTAGCTGTAATCGCGCATGGCCGGGACCGCGAACAGTTTGCTGATCAGCGCAGGCATGCGCTGAATGTCGGCCACGAATACAgcggc                                                                                                                >  2:255541/1‑142 (MQ=255)                                    cctcctGCCCGGGGTACCGGCTGGCCACTCGCCCATCACGGTCGAGCAGCACCCGGTAGCTGTAATCGCGCATGGCCGGGACCGCGAACAGTTTGCTGATCAGCGCAGGCATGCGCTGAATGTCGGCCACGAATACAgcggc                                                                                                                <  1:255541/142‑1 (MQ=255)                                    cctcctGCCCGGGGTACCGGCTGGCCACTCGCCCATCACGGTCGAGCAGCACCCGGTAGCTGTAATCGCGCATGGCCGGGACCGCGAACAGTTTGCTGATCAGCGCAGGCATGCGCTGAATGTCGGCCACGAATACAGCGGCGCGTGCtt                                                                                                        >  2:633608/1‑150 (MQ=255)                                           cccGGGGTACCGGCTGGCCACTCGCCCATCACGGTCGAGCAGCACCCGGTAGCTGTAATCGCGCATGGCCGGGACCGCGAACAGTTTGCTGATCAGCGCATGCCTGCGCTGat                                                                                                                                      <  1:289051/113‑2 (MQ=255)                                           cccGGGGTACCGGCTGGCCACTCGCCCATCACGGTCGAGCAGCACCCGGTAGCTGTAATCGCGCATGGCCGGGACCGCGAACAGTTTGCTGATCAGCGCATGCCTGCGCTGat                                                                                                                                      >  2:289051/1‑112 (MQ=255)                                                             cACTCGCCCATCACGGTCGAGCAGCACCCGGTAGCTGTAATCGCGCATGGCCGGGACCGCGAACAGTTTGCTGATCAGCGCAGGCATGCGCTGAATGTCGGCCACGAATACAGCGGCGCGTGCTTCCa                                                                                                     >  2:79619/1‑128 (MQ=255)                                                             cACTCGCCCATCACGGTCGAGCAGCACCCGGTAGCTGTAATCGCGCATGGCCGGGACCGCGAACAGTTTGCTGATCAGCGCAGGCATGCGCTGAATGTCGGCCACGAATACAGCGGCGCGTGCTTCCa                                                                                                     <  1:79619/128‑1 (MQ=255)                                                                 cGCCCATCACGGTCGAGCAGCACCCGGTAGCTGTAATCGCGCATGGCCGGGACCGCGAACAGTTTGCTGATCAGCGCAGGCATGCGCTGAATGTCGGCCACGAATACAGCGGCGCGTGCTTCCAGGTAGCCCTTGGGCTGCTCGGCCAGt                                                                           <  2:460478/150‑1 (MQ=255)                                                                  gCCCATCACGGTCGAGCAGCACCCGGTAGCTGTAATCGCGCATGGCCGGGACCGCGAACAGTTTGCTGATCAGCGCAGGCATGCGCTGAATGTCGGCCACGAATACAGCGGCGCGTGCTTCCAGGTAg                                                                                                <  1:592709/128‑1 (MQ=255)                                                                  gCCCATCACGGTCGAGCAGCACCCGGTAGCTGTAATCGCGCATGGCCGGGACCGCGAACAGTTTGCTGATCAGCGCAGGCATGCGCTGAATGTCGGCCACGAATACAGCGGCGCGTGCTTCCAGGTAg                                                                                                >  2:592709/1‑128 (MQ=255)                                                                               gAGCAGCACCCGGTAGCTGTAATCGCGCATGGCCGGGACCGCGAACAGTTTGTTGATCAGCGCATGCCTGCGCTGatc                                                                                                                                     >  2:142300/1‑76 (MQ=38)                                                                               gAGCAGCACCCGGTAGCTGTAATCGCGCATGGCCGGGACCGCGAACAGTTTGTTGATCAGCGCATGCCTGCGCTGatc                                                                                                                                     <  1:142300/78‑3 (MQ=38)                                                                                                            gggCCGGTACCGCGAACAGTTTGATGCTCGGCGCAGGCGTGCGTGGAATGTCGTCCACGAATACAGCGGCGCGTGCTTCCAGGTAGCCCTTGGGCTGCTCGGCCAGTGCCGCCTTGACCAGCTTGGCACCGTCCATATCGCGCGccacca                                <  2:383864/149‑1 (MQ=255)                                                                                                              gCCGGGACCGCGAACAGTTTGCTGATCAGCGCAGGCATGCGCTGAATGTCGGCCACGAATACAGCGGCGCGTGCTTCCAGGTAGCCCTTGGGCTGCTCGGCCAGTGCCGCCTTGACCAGCTTGGCACCGTCCATATCGCGCGCCACCAGt                              >  2:191150/1‑150 (MQ=255)                                                                                                               ccGGGACCGCGAACAGTTTGCTGATCAGCGCAGGCATGCGCTGAATGTCGGCCACGAATACAGCGGCGCGTGCTTCCAGGTAGCCCTTGGGCTGCTCGGCCAGTGCCGCCTTGACCAGCTTGGCACCGTCCATATCGCGCGCCACCAGTa                             <  1:191150/150‑1 (MQ=255)                                                                                                                      cgcgAACAGTTTGCTGATCAGCGCAGGCATGCGCTGAATGTCGGCCACGAATACAGCGGCGCGTGCTTCCAGGTAGCCCTTGGGCTGCTCGGCCAGTGCCGCCTTGACCAGCTTGGCACCGTCCATATCGCGCGCCACCAGTAAAATCTg                      >  2:479176/1‑150 (MQ=255)                                                                                                                         gAACAGTTTGCTGATCAGCGCAGGCATGCGCTGAATGTCGGCAACGAATACAGCGGCGCGTGCTTCCAGGTAGCCCTTGGGCTGCTCGGCCAGTGCCGCCTTGACCAGCTTGGCACCGTCCATATCGCGCGCCACCAGTAAAATCTGCgt                   >  1:207719/1‑150 (MQ=255)                                                                                                                                 tGCTGATCAGCGCAGGCATGCGCTGAATGTCGGCCACGAATACAGCGGCGCGTGCTTCCAGGTAGCCCTTGGGCTGCTCGGCCAGTGCCGCCTTGACCAGCTTGGCACCGTCCATATCGCGCGCCACCAGTAAAATCTGCGTGTCGGCAc           <  2:429515/150‑1 (MQ=255)                                                                                                                                    tGATCAGCGCAGGCATGCGCTGAATGTCGGCCACGAATACAGCGGCGCGTGCTTCCAGGTAGCCCTTGGGCTGCTCGGCCAGTGCCGCCTTGACCAGCTTGGCACCGTCCATATCGCGCGCCACCAGTAAAATCTGCGTGTCGGCACCaa        >  2:651383/1‑150 (MQ=255)                                                                                                                                     gATCAGCGCAGGCATGCGCTGAATGTCGGCCACGAATACAGCGGCGCGTGCTTCCAGGTAGCCCTTGGGCTGCTCGGCCAGTGCCGCCTTGACCAGCTTGGCACCGTCCATATCGCGCGCCACCAGTAAAATCTGCGTGTCGGCACCAAg       <  2:612785/150‑1 (MQ=255)                                                                                                                                       tCAGCGCATGCCTGCGCTGAATGTCGGCCACGAATACAGCGGCGCGTGCTTCCAGGTAg                                                                                                >  2:268198/1‑59 (MQ=38)                                                                                                                                       tCAGCGCATGCCTGCGCTGAATGTCGGCCACGAATACAGCGGCGCGTGCTTCCAGGTAg                                                                                                <  1:268198/59‑1 (MQ=38)                                                                                                                                        cAGCGCAGGCATGCGCTGAATGTCGGCCACGAATACAGCGGCGCGTGCTTCCAGGTAGCCCTTGGGCTGCTCGGCCAGTGCCGCCTTGACCAGCTTGGCACCGTCCATATCGCGCGCCACCAGTATAATCTGCGTGTAGGCACCAAGGCt    <  2:29148/150‑1 (MQ=255)                                                                                                                                          gcgcAGGCATGCGCTGAATGTCGGCCACGAATACAGCGGCGCGTGCTTCCAGGTAGCCCTTGGGCTGCTCGGCCAGTGCCGCCTTGACCAGCTTGGCACCGTCCATATCGCGCGCCACCAGTAAAATCTGCGTGTCGGCACCAAGGCTGt  <  1:633608/150‑1 (MQ=255)                                                                                                                                               |                                                                                                                                                 CACGCCCTGCTGCAGTTGCAACCACTGCACCTGGCCCTCCTGCCCGGGGTACCGGCTGGCCACTCGCCCATCACGGTCGAGCAGCACCCGGTAGCTGTAATCGCGCATGGCCGGGACCGCGAACAGTTTGCTGATCAGCGCAGGCATGCGCTGAATGTCGGCCACGAATACAGCGGCGCGTGCTTCCAGGTAGCCCTTGGGCTGCTCGGCCAGTGCCGCCTTGACCAGCTTGGCACCGTCCATATCGCGCGCCACCAGTAAAATCTGCGTGTCGGCACCAAGGCTGT  >  NC_002947/3996799‑3997085

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