Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	3,344,389	C→G	19.2%	Q261H (CAG→CAC)	PP_2941 ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	3,344,389	0	C	G	19.2%	59.8 / 9.9	26	Q261H (CAG→CAC)	PP_2941	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (13/8);  new base G (3/2);  total (16/10)
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.82e-01

AACAGTACGTGGAAAGTCAGAAAAGCGGAGCTTGTGCTGCCAGGTGAGTTCGCCGAGGTCTGATCGGACCCACCGCAACCACGGGCAGGTCGACCACCTGCGAATCTGTGCAGCCGCTTCACTGAAGCGGCTCTATCAGCTCAGCGCCTTGATTGCGAACATTGGCAATGGCTGTGCTCACTCGATACCACTCGAAGGCCTCGGCTGGCTCGCCCAGGTTGAGCGCCATCTGCTCGGCTTGTTCGCCCGGCGTCGCGGGTTGCA  >  NC_002947/3344268‑3344531                                                                                                                          |                                                                                                                                               aaCAGTACGTGGAAAGTCAGAAAAGCGGAGCTTGTGCTGCCAGGTGAGTTCGCCGAGGTCTGATCGGACCCACCGCAACCACGGGCAGGTCGACCACCTGCGAATCTGTGCAGCCGCTTCACTGAAGCGGCTCTATCAGCTCAGCGCCtt                                                                                                                    <  1:364324/150‑1 (MQ=255)      tACGTGGAAAGTCAGAAAAGCGGAGCTTGTGCTGCCAGGTGAGTTCGCCGAGGTCTGATCGGACCCACCGCAACCACGGGCAGGTCGACCACCTGCGAATCTGTGCAGCCGCTTCACTGAAGCGGCTCTATCAGCTCAGCGCCTTGATTg                                                                                                               >  1:241697/1‑150 (MQ=255)        cGTGGAAAGTCAGAAAAGCGGAGCTTGTGCTGCCAGGTGAGTTCGCCGAGGTCTGATCGGACCCACCGCAACCACGGGCAGGTCGACCACCTGCGAATCTGTGCAGCCGCTTCACTGAAGCGGCTCTATCAGCTCAGCGCCTTGATTGCg                                                                                                             <  1:196561/150‑1 (MQ=255)              aaGTCAGAAAAGCGGAGCTTGTGCTGCCAGGTGAGTTCGCCGAGGTCTGATCGGACCCACCGCAACCACGGGCAGGTCGACCACCTGCGAATCTGTGCAGCCGCTTCAGTGAAGCGGct                                                                                                                                      >  2:41280/1‑119 (MQ=255)              aaGTCAGAAAAGCGGAGCTTGTGCTGCCAGGTGAGTTCGCCGAGGTCTGATCGGACCCACCGCAACCACGGGCAGGTCGACCACCTGCGAATCTGTGCAGCCGCTTCAGTGAAGCGGct                                                                                                                                      >  2:41345/1‑119 (MQ=255)              aaGTCAGAAAAGCGGAGCTTGTGCTGCCAGGTGAGTTCGCCGAGGTCTGATCGGACCCACCGCAACCACGGGCAGGTCGACCACCTGCGAATCTGTGCAGCCGCTTCAGTGAAGCGGct                                                                                                                                      <  1:41280/119‑1 (MQ=255)              aaGTCAGAAAAGCGGAGCTTGTGCTGCCAGGTGAGTTCGCCGAGGTCTGATCGGACCCACCGCAACCACGGGCAGGTCGACCACCTGCGAATCTGTGCAGCCGCTTCAGTGAAGCGGct                                                                                                                                      <  1:41345/119‑1 (MQ=255)                           ggAGCTTGTGCTGCCAGGTGAGTTCGCCGAGGTCTGATCGGACCCACCGCAACCACGGGCAGGTCGACCACCTGCGAATCTGTGCAGCCGCTTCACTGAAGCGGCTCTATCAGCTCAGACCCTTGATTACTTACATTGTCAATGGCtgtg                                                                                          >  1:137084/1‑150 (MQ=255)                                            gTGAGTTCGCCGAGGTCTGATCGGACCCACCGCAACCACGGGCAGGTCGACCACCTGCGAATCTGTGCAGCCGCTTCACTGAAGCGGCTCTATCAGCTCAGCGCCTTGATTGCGAACATTGGCAATGGCTGTGCTCACTCGATACCACTc                                                                         >  2:413218/1‑150 (MQ=255)                                                 ttCGCCGAGGTCTGATCGGACCCACCGCAACCACGGGCAGGTCGACCACCTGCGAATCTGTGCAGCCGCTTCACTGAAGCGGCTCTATCAGCTCAGCGCCTTGATTGCGAACATTGGCAATGGCTGTGCTCACTCGATACCACTCGAAgg                                                                    >  1:572422/1‑150 (MQ=255)                                                           tCTGATCGGACCCACCGCAACCACGGGCAGGTCGACCACCTGCGAATCTGTGCAGCCGCTTCACTGAAGCGGCTCTATCAGCTCAGCGCCTTGATTGCGAACATTGGCAATGGCTGTGCTCACTCGATACCACTCGAAGGCCTCggctgg                                                          >  2:431353/1‑150 (MQ=255)                                                                        aCCGCAACCACGGGCAGGTCGACCACCTGCGAATCTGTGCAGCCGCTTCACTGAAGCGGCTCTATCAGCTCAGCGCCTTGATTGCGAACATTGGCAATGGCTGTGCTCACTCGATACCACTCGAAGGCCTCGGCTGGCTCGCCCAGGTTg                                             >  1:511098/1‑150 (MQ=255)                                                                               ccACGGGCAGGTCGACCACCTGCGAATCTGTGCAGCCGCTTCACTGAAGCGGctc                                                                                                                                     >  1:470071/1‑55 (MQ=255)                                                                               ccACGGGCAGGTCGACCACCTGCGAATCTGTGCAGCCGCTTCACTGAAGCGGctc                                                                                                                                     <  2:470071/55‑1 (MQ=255)                                                                                           cGACCACCTGCGAATCTGTGCAGCCGCTTCACTGAAGCGGCTCTATCAGCTCAGCGCCTTGATTGCGAACATTGGCAATGGCTGTGCTCACTCGATACCACTCGAAGGCCTCGGCTGGCTCGCCCAGGTTGAGCGCCATCTGCTCGGCtt                          <  2:470888/150‑1 (MQ=255)                                                                                               caccTGCGAATCTGTGCAGCCGCTTCACTGAAGCGGCTCTATCAGCTCAGCGCCTTGATTGCGAACATTGGCAATGGCTGTGCTCACTCGATACCACTCGAAGGCCTCGGCTGGCTCGCCCAGGTTGAGCGCCATCTGCTCGGCTTGTTc                      >  1:91084/1‑150 (MQ=255)                                                                                                accTGCGAATCTGTGCAGCCGCTTCACTGAAGCGGCTCTATCAGCTCAGCGCCTTGATTGCGAACATTGGCAATGGCTGTGCTCACTCGATACCACTCGAAGGCCTCGGCTGGCTCGCCCAGGTTGAGCGCCATCTGCTCGGCTTGTTCg                     >  1:117097/1‑150 (MQ=255)                                                                                                        aTCTGTGCAGCCGCTTCACTGAAGCGGCTCTATCAGCTCAGCGCCTTGATTGCGAACATTGGCAATGGCTGTGCTCACTCGATACCACTCGAAGGCCTCGGCTGGCTCGCCCAGGTTGAGCGCCATCTGCTCGGCTTGTTCGCCCGGc               >  2:239170/1‑148 (MQ=255)                                                                                                        aTCTGTGCAGCCGCTTCACTGAAGCGGCTCTATCAGCTCAGCGCCTTGATTGCGAACATTGGCAATGGCTGTGCTCACTCGATACCACTCGAAGGCCTCGGCTGGCTCGCCCAGGTTGAGCGCCATCTGCTCGGCTTGTTCGCCCGGc               <  1:239170/148‑1 (MQ=255)                                                                                                           tgtgCAGCCGCTTCACTGAAGCGGCTCTATCAGCTCAGCGCCTTGATTGCGAACATTGGCAATGGCTGTGCTCACTCGATACCACTCGAAGGCCTCGGCTGGCTCGCCCAGGTTGAGCGCCATCTGCTCGGCTTGTTCGCCCGGCGTcgc          >  1:221937/1‑150 (MQ=255)                                                                                                           tgtgCAGCCGCTTCACTGAAGCGGCTCTATCAGCTCAGCGCCTTGATTGCGAACATTGGCAATGGCTGTGCTCACTCGATACCACTCGAAGGCCTCGGCTGGCGCGCCCAGGTTGAGCGCCATCTGCTCGGCTTGTTCGCCCGGCGTcgc          <  2:241697/150‑1 (MQ=255)                                                                                                              gCAGCCGCTTCACTGAAGCGGCTCTATCAGCTCAGCGCCTTGATTGCGAACATTGGCAATGGCTGTGCTCACTCGATACCACTCGAAGGCCTCGGCTGGCTCGCCCAGGTTGAGCGCCATCTGCTCGGCTTGTTCGCCCGGCGTCGCggg       >  1:118193/1‑150 (MQ=255)                                                                                                              gCAGCCGCTTCACTGAAGCGGCTCTATCAGCTCAGCGCCTTGATTGCGAACATTGGCAATGGCTGTGCTCACTCGATACCACTCGAAGGCCTCGGCTGGCTCGCCCAGGTTGAGCGCCATCTGCTCGGCTTGTTCGCCCGGCGTCGCggg       <  2:118193/150‑1 (MQ=255)                                                                                                               cAGCCGCTTCACTGAAGCGGCTCTATCAGCTCAGCGCCTTGATTGCGAACATTGGCAATGGCTGTGCTCACTCGATACCACTCGAAGGCCTCGGCTGGCTCGCCCAGGTTGAGCGCCATCTGCTCGGCTTGTTCGCCCGGCGTCGCGGGt      >  1:26632/1‑150 (MQ=255)                                                                                                                 gCCGCTTCAGTGAAGCGGCTCTATCAGCTCAGCGCCTTGATTGCGAACATTGGCAATGGCTGTGCTCACTCGATACCACTCGAAGGCCTCGGCTGGCTCGCCCAGGTTGAGCGCCATCTGCTCGGCTTGTTCGCCCGGCGTCGCGGGTTg    >  1:529244/1‑150 (MQ=255)                                                                                                                   cGCTTCACTGAAGCGGCTCTATCAGCTCAGCGCCTTGATTGCGAACATTGGCAATGGCTGTGCTCACTCGATACCACTCGAAGGCCTCGGCTGGCTCGCCCAGGTTGAGCGCCATCTGCTCGGCTTGTTCGCCCGGCGTCGCGGGTTGCa  <  2:341592/150‑1 (MQ=255)                                                                                                                          |                                                                                                                                               AACAGTACGTGGAAAGTCAGAAAAGCGGAGCTTGTGCTGCCAGGTGAGTTCGCCGAGGTCTGATCGGACCCACCGCAACCACGGGCAGGTCGACCACCTGCGAATCTGTGCAGCCGCTTCACTGAAGCGGCTCTATCAGCTCAGCGCCTTGATTGCGAACATTGGCAATGGCTGTGCTCACTCGATACCACTCGAAGGCCTCGGCTGGCTCGCCCAGGTTGAGCGCCATCTGCTCGGCTTGTTCGCCCGGCGTCGCGGGTTGCA  >  NC_002947/3344268‑3344531

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