Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	1,712,542	G→T	7.6%	Q92H (CAG→CAT)	adk →	adenylate kinase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	1,712,542	0	G	T	7.6%	162.3 / 4.4	53	Q92H (CAG→CAT)	adk	adenylate kinase
Reads supporting (aligned to +/- strand):  ref base G (24/25);  new base T (2/2);  total (26/27)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.60e-01

CCTGGAGCTGAAGAAAGTCATGGATGCCGGCCAGCTGGTCTCCGACGAGTTGATCATCAGCCTGGTGAAGGAGCGCATCGCTCAGCCTGATTGCGCCAACGGCTGCCTGTTCGACGGCTTCCCACGCACCATCCCACAGGCCGAAGCCATGGTCGCTGCCGGTGTCGACATCGATGCCGTGGTCGAAATTGCCGTGGACGACGAGGAAATCGTCGGCCG‑CATGGCTGGCCGCCGCGTGCACCTGGCCTCGGGCCGTACCTACCACATTCAGTACAACCCGCCGAAA  >  NC_002947/1712404‑1712689                                                                                                                                           |                                                                                                                                                     ccTGGAGCTGAAGAAAGTCATGGATGCCGGCCAGCTGGTCTCCGACGAGTTGATCATCAGCCTGGTGAAGGAGCGCATCGCTCAGCCTGATTGCGCCAACGGCTGCCTGTTCGACGGCTTCCCACGCACCATCCCACAGGCCGAAGCCa                                                                                                                                            <  8:95949/149‑1 (MQ=255)  cTGGAGCTGAAGAAAGTCATGGATGCCGGCCAGCTGGTCTCCGACGAGTTGATCATCAGCCTGGTGAAGGAGCGCATCGCTCAGCCTGATTGCGCCAACGGCTGCCTGTTCGACGGCTTCCCACGCACCATCCCACAGGCCGAAGCCAt                                                                                                                                           <  8:306702/149‑1 (MQ=255)       gCTGAAGAAAGTCATGGATGCCGGCCAGCTGGTCTCCGACGAGTTGATCATCAGCCTGGTGAAGGAGCGCATCGCTCAGCCTGATTGCGCCAACGGCTGCCTGTTCGACGGCTTCCCACGCACCATCCCACAGGCCGAAGCCATGGTc                                                                                                                                       <  4:4080/148‑1 (MQ=255)                  tCATGGATGCCGGCCAGCTGGTCTCCGACGAGTTGATCATCAGCCTGGTGAAGGAGCGCATCGCTCAGCCTGATTGCGCCAACGGCTGCCTGTTCGACGGCTTCCCACGCACCATCCCACAGGCCGAAGCCATGGTCGCTGCCGGTGTc                                                                                                                           >  5:315235/1‑149 (MQ=255)                      ggATGCCGGCCAGCTGGTCTCCGACGAGTTGATCATCAGCCTGGTGAAGGAGCGCATCGCTCAGCCTGATTGCGCCAACGGCTGCCTGTTCGACGGCTTCCCACGCACCATCCCACAGGCCGAAGCCATGGTCGCTGCCGgtgt                                                                                                                            >  8:79240/1‑144 (MQ=255)                      ggATGCCGGCCAGCTGGTCTCCGACGAGTTGATCATCAGCCTGGTGAAGGAGCGCATCGCTCAGCCTGATTGCGCCAACGGCTGCCTGTTCGACGGCTTCCCACGCACCATCCCACAGGCCGAAGCCATGGTCGCTGCCGgtgt                                                                                                                            <  7:79240/144‑1 (MQ=255)                          gCCGGCCAGCTGGTCTCCGACGAGTTGATCATCAGCCTGGTGAAGGAGCGCATCGCTCAGCCTGATTGCGCCAACGGCTGCCTGTTCGACGGCTTCCCACGCACCATCCCACAGGCCGAAGCCATGGTCGCTGCCGGTGTCGACATCGa                                                                                                                   >  1:268814/1‑149 (MQ=255)                          gCCGGCCAGCTGGTCTCCGACGAGTTGATCATCAGCCTGGTGAAGGAGCGCATCGCTCAGCCTGATTGCGCCAACGGCTGCCTGTTCGACGGCTTCCCACGCACCATCCCACAGGCCGAAGCCATGGTCGCTGCCGGTGTCGACATCGa                                                                                                                   >  5:76644/1‑149 (MQ=255)                           ccGGCCAGCTGGTCTCCGACGAGTTGATCATCAGCCTGGTGAAGGAGCGCATCGCTCAGCCTGATTGCGCCAACGGCTGCCTGTTCGACGGCTTCCCACGCACCATCCCACAGGCCGAAGCCATGGTCGCTGCCGGTGTCGACATCGa                                                                                                                   <  6:1156/148‑1 (MQ=255)                            cGGCCAGCTGGTCTCCGACGAGTTGATCATCAGCCTGGTGAAGGAGCGCATCGCTCAGCCTGATTGCGCCAACGGCTGCCTGTTCGACGGCTTCCCACGCACCATCCCACAGGCCGAAGCCATGGTCGCTGCCGGTGTCGACATCGATg                                                                                                                 >  7:46734/1‑149 (MQ=255)                                           cgacgaGTTGATCATCAGCCTGGTGAAGGAGCGCATCGCTCAGCCTGATTGCGCCAACGGCTGCCTGTTCGACGGCTTCCCACGCACCATCCCACAGGCCGAAGCCATGGTCGCTGCCGGTGTCGACATCGATGCCGTg                                                                                                            <  4:147151/139‑1 (MQ=255)                                           cgacgaGTTGATCATCAGCCTGGTGAAGGAGCGCATCGCTCAGCCTGATTGCGCCAACGGCTGCCTGTTCGACGGCTTCCCACGCACCATCCCACAGGCCGAAGCCATGGTCGCTGCCGGTGTCGACATCGATGCCGTg                                                                                                            >  3:147151/1‑139 (MQ=255)                                            gacgaGTTGATCATCAGCCTGGTGAAGGAGCGCATCGCTCAGCCTGATTGCGCCAACGGCTGCCTGTTCGACGGCTTCCCACGCACCATCCCACAGGCCGAAGCCATGGTCGCTGCCGGTGTCGACATCGATGCCGTGGTCGAAATTGc                                                                                                 >  7:282660/1‑149 (MQ=255)                                                  ttGATCATCAGCCTGGTGAAGGAGCGCATCGCTCAGCCTGATTGCGCCAACGGCTGCCTGTTCGACGGCTTCCCACGCACCATCCCACAGGCCGAAGCCATGGTCGCTGCCGgtgt                                                                                                                            >  4:164145/1‑116 (MQ=255)                                                  ttGATCATCAGCCTGGTGAAGGAGCGCATCGCTCAGCCTGATTGCGCCAACGGCTGCCTGTTCGACGGCTTCCCACGCACCATCCCACAGGCCGAAGCCATGGTCGCTGCCGgtgt                                                                                                                            <  3:164145/116‑1 (MQ=255)                                                        atcaGCCTGGTGAAGGAGCGCATCGCTCAGCCTGATTGCGCCAACGGCTGCCTGTTCGACGGCTTCCCACGCACCATCCCACAGGCCGAAGCCATGGTCGCTGCCGGTGTCGACATCGATGCCGTGGTCGAAATTGCCGTGGAcgacga                                                                                     <  4:44073/149‑1 (MQ=255)                                                         tcaGCCTGGTGAAGGAGCGCATCGCTCAGCCTGATTGCGCCAACGGCTGCCTGTTCGACGGCTTCCCACGCACCATCCCACAGGCCGAAGCCATGGTCGCTGCCGGTGTCGACATCGATGCCGTGGTCGAAATTGCCGTGGACGACGAg                                                                                    <  6:76644/149‑1 (MQ=255)                                                                 gTGAAGGAGCGCATCGCTCAGCCTGATTGCGCCAACGGCTGCCTGTTCGACGGCTTCCCACGCACCATCCCACAGGCCGAAGCCATGGTCGCTGCCGGTGTCGACATCGATGCCGTGGTCGAAATTGCCGTGGACGACGAGGAAAtc                                                                              <  7:23433/147‑1 (MQ=255)                                                                  tGAAGGAGCGCATCGCTCAGCCTGATTGCGCCAACGGCTGCCTGTTCGACGGCTTCCCACGCACCATCCCACAGGCCGAAGCCATGGTCGCTGCCGGTGTCGACATCGATGCCGTGGTCGAAATTGCCGTGGACGACGAGGAAAtcgtc                                                                           <  5:79492/149‑1 (MQ=255)                                                                   gAAGGAGCGCATCGCTCAGCCTGATTGCGCCAACGGCTGCCTGTTCGACGGCTTCCCACGCACCATCCCACAGGCCGAAGCCATGGTCGCTGCCGGTGTCGACATCGATGCCGTGGTCGAAATTGCCGTGGACGACGAGGAAAtcgtcg                                                                          >  2:58435/1‑149 (MQ=255)                                                                            cATCGCTCAGCCTGATTGCGCCAACGGCTGCCTGTTCGACGGCTTCCCACGCACCATCCCACAGGCCGAAGCCATGGTCGCTGCCGGTGTCGACATCGATGCCGTGGTCGAAATTGCCGTGGACGACGAGGAAATCGTCGGCCG‑CAtgg                                                                <  8:46734/149‑1 (MQ=255)                                                                             aTCGCTCAGCCTGATTGCGCCAACGGCTGCCTGTTCGACGGCTTCCCACGCACCATCCCACAGGCCGAAGCCATGGTCGCTGACGGTGTCGACATCGATGCCGTGGTCGAAATTGCCGTGGACGACGAGGAAATCGTCGGCCG‑CAtggc                                                               >  6:171170/1‑149 (MQ=255)                                                                                      ccTGATTGCGCCAACGGCTGCCTGTTCGACGGCTTCCCACGCACCATCCCACAGGCCGAAGCCATGGTCGCTGCCGGTGTCGACATCGATGCCGTGGTCGAAATTGCCGTGGACGACGAGGAAATCGTCGGCCG‑CATGGCTGGccgccg                                                      >  8:12453/1‑149 (MQ=255)                                                                                             gcgcCAACGGCTGCCTGTTCGACGGCTTCCCACGCACCATCACACAGGCCGAAGCCATGGTCGCTGCAGGTGTCGACATCGATGCCGTGGTCGAAATTGCCGTGGACGACGAGGAAATCGTCGGCCG‑CATGGCTCGCCGCCGCGTGCAc                                               >  6:268803/1‑149 (MQ=255)                                                                                              cgcCAACGGCTGCCTGTTCGACGGCTTCCCACGCACCATCCCACAGGCCGAAGCCATGGTCGCTGCCGGTGTCGACATCGATGCCGTGGTCGAAATTGCCGTGGACGACGAGGAAATCGTCGGCCG‑CATGGCTGGCCGCCGCGTGCAc                                               <  6:315235/148‑1 (MQ=255)                                                                                                ccAACGGCTGCCTGTTCGACGGCTTCCCACTCACCATCCCACATGCCGAAGCCATGGTCGCTGCCGGTGTCGACATCGATGCCGTGGTCGAAATTGCCGTGGACGACGAGGAAATCGTCGGCCG‑CATGGCTGGCcgccgc                                                     >  3:33048/1‑140 (MQ=255)                                                                                                ccAACGGCTGCCTGTTCGACGGCTTCCCACGCACCATCCCACATTCCTAAGCCATGGTCGCTGCCGGTGTCGACATCGATGCCGTGGTCGAAATTGCCTTGGACGACGAGGAAATCGTCGGCCG‑CATGGCTTGCcgccgc                                                     <  4:33048/140‑1 (MQ=255)                                                                                                ccAACGGCTGCCTGTTCGACGGCTTCCCACGCACCATCCCACATGCCGAAGCCATGGTCGCTGCCGGTGTCGACATCGATGCCGTGGTCGAAATTGCCGTGGACGACGAGGAAATCGTCGGCCG‑CATGGCTGGCcgccgc                                                     >  3:37097/1‑140 (MQ=255)                                                                                                ccAACGGCTGCCTGTTCGACGGCTTCCCACGCACCATCCCACATGCCGAAGCCATGGTCGCTGCCGGTGTCGACATCGATGCCGTGGTCGAAATTGCCGTGGACGACGAGGAAATCGTCGGCCG‑CATGGCTGGCcgccgc                                                     <  4:37097/140‑1 (MQ=255)                                                                                                     ggCTGCCTGTTCGACGGCTTCCCACGCACCATCCCACAGGCCGAAGCCATGGTCGCTGCCGGTGTCGa                                                                                                                         >  8:62001/1‑68 (MQ=255)                                                                                                     ggCTGCCTGTTCGACGGCTTCCCACGCACCATCCCACAGGCCGAAGCCATGGTCGCTGCCGGTGTCGa                                                                                                                         <  7:62001/68‑1 (MQ=255)                                                                                                     ggCTGCCTGTTCGACGGCTTCCCACGCACCATCCCACAGGCCGAAGCCATGGTCGCTGCCGGTGTCGACATCGATGCCGTGGTCGAAATTGCCGTGGACGACGAGGAAATCGTCGGCCG‑CATGGCTGGCCGCCGCGTGCACCTGGCCTc                                       >  2:196401/1‑149 (MQ=255)                                                                                                      gCTGCCTGTTCGACGGCTTCCCACGCACCATCCCACAGGCCGAAGCCATGGTCGCTGCCGGTGTCGACATCGATGCCGTGGTCGAAATTGCCGTGGACGACGAGGAAATCGTCGGCCG‑CATGGCTGGCCGCCGCGTGCACCTGGCCTCg                                      <  2:201857/149‑1 (MQ=255)                                                                                                             gTTCGACGGCTTCCCACGCACCATCCCACAGGCCGAAGCCATGGTCGCTGCCGGTGTCGACATCGATGCCGTGGTCGAAATTGCCGTGGACGACGAGGAAATCGTCGGCCG‑CATGGCTGGCCGCCGCGTGCACCTGGCCTCGGGCCGta                               >  8:141689/1‑149 (MQ=255)                                                                                                                 gACGGCTTCCCACGCACCATCCCACAGGCCGAAGCCATGGTCGCTGCCGGTGTCGACATCGATGCCGTGGTCGAAATTGCCGTGGACGACGAGGAAATCGTCGGCCG‑CATGGCTGGCCGCCGCGTGCACCTGGCCTCGGGCCGtaccta                           >  5:302333/1‑149 (MQ=255)                                                                                                                 gACGGCTTCCCACGCACCATCCCACAGGCCGAAGCCATGGTCGCTGCCGGTGTCGACATCGATGCCGTGGTCGAAATTGCCGTGGACGACGAGGAAATCGTCGGCCG‑CATGGCTGGCCGCCGCGTGCACCTGGCCTCGGGCCGtaccta                           <  1:58435/149‑1 (MQ=255)                                                                                                                     gCTTCCCACGCACCATCCCACAGGCCGAAGCCATGGTCGCTGGCGGTGTCGACATCGATGCCGTGGTCGAAATTGCCGTGGACGACGAGGAAATCGTCGGCCGCCATGGCTGGCCGCCGCGTGCACCTGGCCTCgg                                     <  1:40626/136‑1 (MQ=255)                                                                                                                     gCTTCCCACGCACCATCCCACAGGCCGAAGCCATGGTCGCTGCCGGTGTCGACATCGATGCCGTGGTCGAAATTGCCGTGGACGACGAGGAAATCGTCGGCCG‑CATGGCTGGCCGCCGCGTGCACCTGGCCTCgg                                     >  2:40626/1‑135 (MQ=255)                                                                                                                       ttCCCACGCACCATCCCACAGGCCGAAGCCATGGTCGCTGCCGGTGTCGACATCGATGCCGTGGTCGAAATTGCCGTGGACGACGAGGAAATCGTCGGCCG‑CATGGCTGGCCGCCGCGTGCACCTGGCCTCGGGCCGTACCTACCACAt                     >  2:249560/1‑149 (MQ=255)                                                                                                                          ccACGCACCATCCCACAGGCCGAAGCCATGGTCGCTGCCGGTGTCGACATCGATGCCGTGGTCGAAATTGCCGTGGACGACGAGGAAATCGTCGGCCG‑CATGGCTGGCCGCCGCGTGCACCTGGCCTCGGGCCGTACCTACcaca                      <  3:189799/145‑1 (MQ=255)                                                                                                                          ccACGCACCATCCCACAGGCCGAAGCCATGGTCGCTGCCGGTGTCGACATCGATGCCGTGGTCGAAATTGCCGTGGACGACGAGGAAATCGTCGGCCG‑CATGGCTGGCCGCCGCGTGCACCTGGCCTCGGGCCGTACCTACcaca                      >  4:189799/1‑145 (MQ=255)                                                                                                                           cACGCACCATCCCACAGGCCGAAGCCATGGTCGCTGCCGGTGTCGACATCGATGCCGTGGTCGAAATTGCCGTGGACGACGAGGAAATCGTCGGCCG‑CATGGCTGGCCGCCGCGTGCACCTGGCCTCGGGCCGTACCTACCACATTCAg                 <  1:231490/149‑1 (MQ=255)                                                                                                                               cACCATCCCACAGGCCGAAGCCATGGTCGCTGCCGGTGTCGACATTGATGCCGTGGTCGAAATTGCCGTGGACGACGAGGAAATCGTCGGCCG‑CATGGCTGGCCGCCGCGTGCACCTGGCCTCGGGCCGtacct                            <  2:234797/134‑1 (MQ=255)                                                                                                                               cACCATCCCACAGGCCGAAGCCATGGTCGCTGCCGGTGTCGACATTGATGCCGTGGTCGAAATTGCCGTGGACGACGAGGAAATCGTCGGCCG‑CATGGCTGGCCGCCGCGTGCACCTGGCCTCGGGCCGtacct                            >  1:234797/1‑134 (MQ=255)                                                                                                                                  cATCCCACAGGCCGAAGCCATGGTCGCTGCCGGTGTCGACATCGATGCCGTGGTCGAAATTGCCGTGGACGACGAGGAAATCGTCGGCCG‑CATGGCTGGCCGCCGCGt                                                   >  8:277731/1‑108 (MQ=255)                                                                                                                                  cATCCCACAGGCCGAAGCCATGGTCGCTGCCGGTGTCGACATCGATGCCGTGGTCGAAATTGCCGTGGACGACGAGGAAATCGTCGGCCG‑CATGGCTGGCCGCCGCGt                                                   <  7:277731/108‑1 (MQ=255)                                                                                                                                  cATCCCACAGGCCGAAGCCATGGTCGCTGCCGGTGTCGACATCGATGCCGTGGTCGAAATTGCCGTGGACGACGAGGAAATCGTCGGCCG‑CATGGCTGGCCGCCGCGTGCACCTGGCCTCGGGCCGTACCTACca                        <  2:207513/135‑1 (MQ=255)                                                                                                                                  cATCCCACAGGCCGAAGCCATGGTCGCTGCCGGTGTCGACATCGATGCCGTGGTCGAAATTGCCGTGGACGACGAGGAAATCGTCGGCCG‑CATGGCTGGCCGCCGCGTGCACCTGGCCTCGGGCCGTACCTACca                        >  1:207513/1‑135 (MQ=255)                                                                                                                                   aTCCCACAGGCCGAAGCCATGGTCGCTGCCGGTGTCGACATCGATGCCGTGGTCGAAATTGCCGTGGACGACGAGGAAATCGTCGGCCG‑CATGGCTGGCCGCCGCGTGCACCTGGCCTCGGGCCGTACCTACCACATTCAGt                <  8:194823/142‑1 (MQ=255)                                                                                                                                   aTCCCACAGGCCGAAGCCATGGTCGCTGCCGGTGTCGACATCGATGCCGTGGTCGAAATTGCCGTGGACGACGAGGAAATCGTCGGCCG‑CATGGCTGGCCGCCGCGTGCACCTGGCCTCGGGCCGTACCTACCACATTCAGt                >  7:194823/1‑142 (MQ=255)                                                                                                                                       cacaGGCCGAAGCCATGGTCGCTGCCGGTGTCGACATCGATGCCGTGGTCGAAATTGCCGTGGACGACGAGGAAATCGTCGGCCG‑CATGGCTGGCCGCCGCGTGCACCTGGCCTCGGGCCGTACCTAACACATTCAGTACAACccgcc      >  7:58441/1‑148 (MQ=255)                                                                                                                                        acaGGCCGAAGCCATGGTCGCTGCCGGTGTCGACATCGATGCCGTGGTCGAAATTGCCGTGGACGACGAGGAAATCGTCGGCCG‑CATGGCTGGCCGCCGCGTGCACCTGGCCTCGGGCCGTACCTACCACATTCAGTACAACCCGCCGa    <  2:268814/149‑1 (MQ=255)                                                                                                                                          aGGCCGAAGCCATGGTCGCTGCCGGTGTCGACATCGATGCCGTGGTCGAAATTGCCGTGGACGACGAGGAAATCGTCGGCCG‑CATGGCTGGCCGCCGCGTGCACCTGGCCTCGGGCCGTACCTACCACATTCAGTACAACCCGCCGaaa  <  4:240845/149‑1 (MQ=255)                                                                                                                                           |                                                                                                                                                     CCTGGAGCTGAAGAAAGTCATGGATGCCGGCCAGCTGGTCTCCGACGAGTTGATCATCAGCCTGGTGAAGGAGCGCATCGCTCAGCCTGATTGCGCCAACGGCTGCCTGTTCGACGGCTTCCCACGCACCATCCCACAGGCCGAAGCCATGGTCGCTGCCGGTGTCGACATCGATGCCGTGGTCGAAATTGCCGTGGACGACGAGGAAATCGTCGGCCG‑CATGGCTGGCCGCCGCGTGCACCTGGCCTCGGGCCGTACCTACCACATTCAGTACAACCCGCCGAAA  >  NC_002947/1712404‑1712689

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