Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	1,155,411	G→A	6.9%	V96V (GTG→GTA)	tmk →	thymidylate kinase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,155,411	0	G	A	6.9%	195.7 / 4.3	58	V96V (GTG→GTA)	tmk	thymidylate kinase
Reads supporting (aligned to +/- strand):  ref base G (19/35);  new base A (2/2);  total (21/37)
Fisher's exact test for biased strand distribution p-value = 6.15e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.82e-01

GATATCAAATCGGTAGGCGATGAAGTCATTACCGATAAAGCCGAAGTTCTGATGTTTTATGCCGCGCGCGTTCAACTGGTAGAAACGGTCATCAAACCAGCTCTGGCTAACGGCACCTGGGTGATTGGCGATCGCCACGATCTCTCCACTCAGGCGTATCAGGGCGGCGGACGTGGTATTGACCAACATATGCTGGCAACACTGCGTGATGCTGTTCTCGGGGATTTTCGCCCCGACTTAACGCTCTATCTCGATGTT  >  NC_000913/1155289‑1155546                                                                                                                           |                                                                                                                                        gATATCAAATCGGTAGGCGATGAAGTCATTACCGATAAAGCCGAAGTTCTGATGTTTTATGCCGCGCGCGTTCAACTGGTAGAAACGGTCATCAAACCAGCTCTGGCTAACGGCACCTGGGTGATTGGCGATCGCCACGa                                                                                                                        >  1:351089/1‑140 (MQ=255)      cAAATCGGTAGGCGATGAAGTCATTACCGATAAAGCCGAAGTTCTGATGTTTTATGCCGCGCGCGTTCAACTGGTAGAAACGGTCATCAAACCAGCTCTGGCTAACGGCACCTGGGTGATTGGCGATCGCCACGAtctct                                                                                                                   <  1:1086782/140‑1 (MQ=255)      cAAATCGGTAGGCGATGAAGTCATTACCGATAAAGCCGAAGTTCTGATGTTTTATGCCGCGCGCGTTCAACTGGTAGAAACGGTCATCAAACCAGCTCTGGCTAACGGCACCTGGGTGATTGGCGATCGCCACGAtctct                                                                                                                   >  1:1114043/1‑140 (MQ=255)      cAAATCGGTAGGCGATGAAGTCATTACCGATAAAGCCGAAGTTCTGATGTTTTATGCCGCGCGCGTTCAACTGGTAGAAACGGTCATCAAACCAGCTCTGGCTAACGGCACCTGGGTGATTGGCGATCGCCACGAtctct                                                                                                                   >  1:251781/1‑140 (MQ=255)           cGGTAGGCGATGAAGTCATTACCGATAAAGCCGAAGTTCTGATGTTTTATGCCGCGCGCGTTCAACTGGTAGAAACGGTCATCAAACCAGCTCTGGCTAACGGCACCTGGGTGATTGGCGATCGCCACGATCTCTCCACt                                                                                                              >  1:468088/1‑140 (MQ=255)               aGGCGATGAAGTCATTACCGATAAAGCCGAAGTTCTGATGTTTTATGCCGCGCGCGTTCAACTGGTAGAAACGGTCATCAAACCAGCTCTGGCTAACGGCACCTGGGTGATTTGCGATCGCCACGATCTCTCCACTCAgg                                                                                                          <  2:255646/140‑1 (MQ=255)               aGGCGATGAAGTCATTACCGATAAAGCCGAAGTTCTGATGTTTTATGCCGCGCGCGTTCAACTGGTAGAAACGGTCATCAAACCAGCTCTGGCTAACGGCACCTGGGTGATTGGCGATCGCCACGATCTCTCCACTCAgg                                                                                                          <  2:887262/140‑1 (MQ=255)               aGGCGATGAAGTCATTACAGATAAAGCCGAAGTTCTGATGTTTTATGCCGCGCGCGTTCAACTGGTAGAAACGGTCATCAAACCAGCTCTGGCTAACGGCACCTGGGTGATTGGCGATCGCCACGATCTCTCCACTCAgg                                                                                                          <  2:255811/140‑1 (MQ=255)                            aTTACCGATAAAGCCGAAGTTCTGATGTTTTATGCCGCGCGCGTTCAACTGGTAGAAACGGTCATCAAACCAGCTCTGGCTAACGGCACCTGGGTGATTGGCGATCGCCACGATCTCTCCACTCAGGCGTATCAGggcgg                                                                                             >  2:953410/1‑140 (MQ=255)                            aTTACCGATAAAGCCGAAGTTCTGATGTTTTATGCCGCGCGCGTTCAACTGGTAGAAACGGTCATCAAACCAGCTCTGGCTAACGGCACCTGGGTGATTGGCGATCGCCACGATCTCTCCACTCAGGCGTATCAGggcgg                                                                                             >  2:1083039/1‑140 (MQ=255)                               aCCGATAAAGCCGAAGTTCTGATGTTTTATGCCGCGCGCGTTCAACTGGTAGAAACGGTCATCAAACCAGCTCTGGCTAACGGCACCTGGGTGATTGGCGATCGCCACGATCTCTCCACTCAGGCGTATCAGGGcggcgg                                                                                          >  2:512223/1‑140 (MQ=255)                                 cGATAAAGCCGAAGTTCTGATGTTTTATGCCGCGCGCGTTCAACTGGTAGAAACGGTCATCAAACCAGCTCTGGCTAACGGCACCTGGGTGATTGGCGATCGCCACGATCTCTCCACTCAGGCGTATCAGGGCGGCGGAc                                                                                        <  1:668682/140‑1 (MQ=255)                                       aGCCGAAGTTCTGATGTTTTATGCCGCGCGCGTTCAACTGGTAGAAACGGTCATCAAACCAGCTCTGGCTAACGGCACCTGGGTGATTGGCGATCGCCACGATCTCTCCACTCAGGCGTATCAGGGCGGCGGACGTGGTa                                                                                  <  1:1086070/140‑1 (MQ=255)                                           gcaGTTCTGATGTTTTATGCCGCGCGCGTTCAACTGGTAGAAACGGTCATCAAACCAGCTCTGGCTAACGGCACCTGGGTGATTGGCGATCGCCACGATCTCTCCACTCAGGCGTATCAGGGCGGCGGACGTGGTATTGa                                                                              <  1:317948/138‑1 (MQ=255)                                           gAAGTTCTGATGTTTTATGCCGCGCGCGTTCAACTGGTAGAAACGGTCATCAAACCAGCTCTGGCTAACGGCACCTGGGTGATTGGCGATCGCCACGATCTCTCCACTCAGGCGTATCAGGGCGGCGGACGTGGTATTGa                                                                              <  1:800590/140‑1 (MQ=255)                                           gAAGTTCTGATGTTTTATGCCGCGCGCGTTCAACTGGTAGAAACGGTCATCAAACCAGCTCTGGCTAACGGCACCTGGGTGATTGGCGATCGCCACGATCTCTCCACTCAGGCGTATCAGGGCGGCGGACGTGGTATTGa                                                                              <  1:800847/140‑1 (MQ=255)                                               ttCTGATGTTTTATGCCGCGCGCGTTCAACTGGTAGAAACGGTCATCAAACCAGCTCTGGCTAACGGCACCTGGGTGATTGGCGATCGCCACGATCTCTCCACTCAGGCGTATCAGGGCGGCGGACGTGGTATTGACCaa                                                                          <  2:225793/140‑1 (MQ=255)                                                    aTGTTTTATGCCGCGCGCGTTCAACTGGTAGATACGGTCATCAAACCAGCTCTGGCTAACGGCACCTGGTTGATTGGCGATCGGCAAGATCTCTCCACTCAGGCTTATCAGGGCGGCGGACGTGGTATTGACCAACatat                                                                     <  2:491916/140‑1 (MQ=255)                                                    aTGTTTTATGCCGCGCGCGTTCAACTGGTAGAAACGGTCATCAAACCAGCTCTGGCTAACGGCACCTGGGTGATTGGCGATCGCCACGATCTCTCCACTCAGGCGTATCAGGGCGGCGGACGTGGTATTGACCAACatat                                                                     <  1:667160/140‑1 (MQ=255)                                                     tGTTTTATGCCGCGCGCGTTCAACTGGTAGAAACGGTCATCAAACCAGCTCTGGCTAACGGCACCTGGGTGATTGGCGATCGCCACGATCTCTCCACTCAGGCGTATCAGGGCGGCGGACGTGGTATTGACCAACATATg                                                                    <  1:614085/140‑1 (MQ=255)                                                     tGTTTTATGCCGCGCGCGTTCAACTGGTAGAAACGGTCATCAAACCAGCTCTGGCTAACGGCACCTGGGTGATTGGCGATCGCCACGATCTCTCCACTCAGGCGTATCAGGGCGGCGGACGTGGTATTGACCAACATATg                                                                    <  1:614186/140‑1 (MQ=255)                                                     tGTTTTATGCCGCGCGCGTTCAACTGGTAGAAACGGTCATCAAACCAGCTCTGGCTAACGGCACCTGGGTGATTGGCGATCGCCACGATCTCTCCACTCAGGCGTATCAGGGCGGCGGACGTGGTATTGACCAACATATg                                                                    <  1:216465/140‑1 (MQ=255)                                                              ccGCGCGCGTTCAACTGGTCGAAACGGTCATCAAACCAGCTCGGGCTAACGGCACCTGGGTGATTGGCGATCGCCACGATCTCTCCACTCAGGCGTATCAGGGCGGCGGACGTGGTATTGACCAACATATGCTGGCAaca                                                           <  2:1017655/140‑1 (MQ=255)                                                              ccGCGCGCGTTCAACTGGTAGAAACGGTCATCAAACCAGCTCTGGCTAACGGCACCTGGGTGATTGGCGATCGCCACGATCTCTCCACTCAGGCGTATCAGGGCGGCGGACGTGGTATTGACCAACATATGCTGGCAaca                                                           <  2:826282/140‑1 (MQ=255)                                                              ccGCGCGCGTTCAACTGGTAGAAACGGTCATCAAACCAGCTCTGGCTAACGGCACCTGGGTGATTGGCGATCGCCACGATCTCTCCACTCAGGCGTATCAGGGCGGCGGACGTGGTATTGACCAACATATGCTGGCAaca                                                           <  2:825690/140‑1 (MQ=255)                                                                 cgcgcgTTCAACTGGTAGAAACGGTCATCAAACCAGCTCTGGCTAACGGCACCTGGGTGATTGGCGATCGCCACGATCTCTCCACTCAGGCGTATCAGGGCGGCGGACGTGGTATTGACCAACATATGCTGGCAACACTg                                                        >  1:231357/1‑140 (MQ=255)                                                                   cgcgTTCAACTGGTAGAAACGGTCATCAAACCAGCTCTGGCTAACGGCACCTGGGTGATTGGCGATCGCCACGATCTCTCCTCTCAGGCGTATCAGGGCGGCGGACGTGGTATTGACCAACATATGATGGCAACACTGCg                                                      <  2:897131/140‑1 (MQ=255)                                                                   cgcgTTCAACTGGTAGAAACGGTCATCAAACCAGCTCTGGCTAACGGCACCTGGGTGATTGGCGATCGCCACGATCTCTCCACTCAGGCGTATCAGGGCGGCGGACGTGGTATTGACCAACATATGCTGGCAACACTGCg                                                      <  2:113215/140‑1 (MQ=255)                                                                   cgcgTTCAACTGGTAGAAACGGTCATCAAACCAGCTCTGGCTAAAGGCACCTGGGTGATTGGCGATCGCCACGATCTCTCCACTCAGGCGTATCAGGGCGGCGGACGTGGTATTGACCAACATATGCTGGCAACACTGCg                                                      <  2:464781/140‑1 (MQ=255)                                                                        tCAACTGGTAGCAACGGTCATCAAACCAGCTCTGGCTAACGGCACCTGGGTGATTGGCGATCGCCACGATCTCTCCACTCAGGCGTATCAGGGCGGCGGACGTGGTTTTGACCAACATATGCTGGCAACACTGCGTGATg                                                 <  2:1058765/140‑1 (MQ=255)                                                                        tCAACTGGTAGAAACGGTGATCAATCCAGCTCTGGCTAACGGCACCTGGGTGATTGGCGATCGCCACGATCTCTCCACTCAGGCGTATCAGGGCGGCGGTCGTGGTATTGACCAACATATGCTGGCAACACTGCGTGATg                                                 <  2:159311/140‑1 (MQ=255)                                                                        acAACTGGTAGAAACGGTCATCAAACCAGCTCTGGCTAACGGCACCTGGGTGAATGGCGATCGCCACGATCTCTCCACTCAGGCGTATCAGGGCGGCGGACGTGGTATTGACCAACATATGCTGGCAACACTGCGTGATg                                                 <  2:1058828/139‑1 (MQ=255)                                                                                tAGAAACGGTCATCAAACCAGCTCTGGCTAACGGCACCTGGGTGATTGGCGATCGCCACGATCTCTCCACTCAGGCGTATCAGGGCGGCGGACGTGGTATTGACCAACATATGCTGGCAACACTGCGTGATGCTGTtctc                                         >  2:1161520/1‑140 (MQ=255)                                                                                         tcatcaAACCAGCTCTGGCTAACGGCACCTGGGTGATTGGCGATCGCCACGATCTCTCCACTCAGGCGTATCAGGGCGGCGGACGTGGTATTGACCAACATATGCTGGCAACACTGCGTGATGCTGTTCTCGGGGAtttt                                >  1:1083249/1‑140 (MQ=255)                                                                                         tcatcaAACCAGCTCTGGCTAACGGCACCTGGGTGATTGGCGATCGCCACGATCTCTCCACTCAGGCGTATCAGGGCGGCGGACGTGGTATTGACCAACATATGCTGGCAACACTGCGTGATGCTGTTCTCGGGGAtttt                                >  1:1083283/1‑140 (MQ=255)                                                                                          catcaAACCAGCTCTGGCTAACGGCACCTGGGTGATTGGCGATCGCCACGATCTCTCCACTCAGGCGTATCAGGGCGGCGGACGTGGTATTGACCAACATATGCTGGCAACACTGCGTGATGCTGTTCTCGGGGATTTTc                               >  2:417509/1‑140 (MQ=255)                                                                                          catcaAACCAGCTCTGGCTAACGGCACCTGGGTGATTGGCGATCGCCACGATCTCTCCACTCAGGCGTATCAGGGCGGCGGACGTGGTATTGACCAACATATGCTGGCAACACTGCGTGATGCTGTTCTCGGGGATTTTc                               >  2:864194/1‑140 (MQ=255)                                                                                          catcaAACCAGCTCTGGCTAACGGCACCTGGGTGATTGGCGATCGCCACGATCTCTCCACTCAGGCGTATCAGGGCGGCGGACGTGGTATTGACCAACATATGCTGGCAACACTGCGTGATGCTGTTCTCGGGGATTTTc                               >  2:447499/1‑140 (MQ=255)                                                                                                        tGGCTAACGGCACCTGGGTGATTGGCGATCGCCACGATCTCTCCACTCAGGCGTATCAGGGCGGCGGACGTGGTATTGACCAACATATGCTGGCAACACTGCGTGATGCTGTTCTCGGGGATTTTCGCCCCGACTTAACg                 <  2:231357/140‑1 (MQ=255)                                                                                                          gCTAACGGCACCTGGGTGATTGGCGTTCGCCACGATCTCTCCACTCTGGCGTTTCAGGGCGGCGGACGTGGTATTGACCTACATATGCTGGCAACTCTGAGTGTTGCTGTTGTCGGGGATTTTCGCCGCGAGTTAACGct               <  2:1105663/140‑1 (MQ=255)                                                                                                          gCTAACGGCACCTGGGTGATTGGCGATCGCCACGATCTCTCCACTCAGGCGTATCAGGGCGGCGGACGTGGTATTGACCAACATATGCTGGCAACACTGCGTGATGCTGTTCTCGGGGATTTTCGCCCCGACTTAACGct               <  1:759151/140‑1 (MQ=255)                                                                                                          gCTAACGGCACCTGGGTGATTGGCGATCGCCACGATCTCTCCACTCAGGCGTATCAGGGCGGCGGACGTGGTATTGACCAACATATGCTGGCAACACTGCGTGATGCTGTTCTCGGGGATTTTCGCCCCGACTTAACGct               <  2:1099360/140‑1 (MQ=255)                                                                                                          gCTAACGGCACCTGGGTGATTGGCGATCGCCACGATCTCTCCACTCAGGCGTATCAGGGCGGCGGACGTGGTATTGACCAACATATGCTGGCAACACTGCGTGATGCTGTTCTCGGGGATTTTCGCCCCGACTTAACGct               >  2:540046/1‑140 (MQ=255)                                                                                                          gCTAACGGCACCTGGGTGATTGGCGATCGCCACGATCTCTCCACTCAGGCGTATCAGGGCGGCGGACGTGGTATTGACCAACATATGCTGGCAACACTGCGTGATGCTGTTCTCGGGGATTTTCGCCCCGACTTAACGct               <  2:292731/140‑1 (MQ=255)                                                                                                          gCTAACGGCACCTGGGTGATTGGCGATCGCCACGATCTCTCCACTCAGGCGTATCAGGGCGGCGGACGTGGTATTGACCAACATATGCTGGCAACACTGCGTGATGCTGTTCTCGGGGATTTTCGCCCCGACTTAACGct               <  1:924812/140‑1 (MQ=255)                                                                                                          gCTAACGGCACCTGGGTGATTGGCGATCGCCACGATCTCTCCACTCAGGCGTATCAGGGCGGCGGACGTGGTATTGACCAACATATGCTGGCAACACTGCGTGATGCTGTTATCGGGGATTTTCGCCCCGACTTAACGct               <  1:1160605/140‑1 (MQ=255)                                                                                                               cGGCACCTGGGTGATTGGCGATCGCCACGATCTCTCCACTCAGGCGTATCAGGGCGGCGGACGTGGTATTGACCAACATATGCTGGCAACACTGCGTGATGCTGTTCTCGGGGATTTTCGccc                           <  1:108430/123‑1 (MQ=255)                                                                                                               cGGCACCTGGGTGATTGGCGATCGCCACGATCTCTCCACTCAGGCGTATCAGGGCGGCGGACGTGGTATTGACCAACATATGCTGGCAACACTGCGTGATGCTGTTCTCGGGGATTTTCGccc                           >  2:108430/1‑123 (MQ=255)                                                                                                               cGGCACCTGGGTGATTGGCGATCGCCACGATCTCTCCACTCAGGCGTATCAGGGCGGCGGACGTGGTATTGACCAACATATGCTGGCAACACTGCGTGATGCTGTTCTCGGGGATTTTCGCCCCGACTTAACGCTCTAtc          >  2:393407/1‑140 (MQ=255)                                                                                                               cGGCACCTGGGTAATTGGCGATCGCCACGATCTCTCCACTCAGGCGTATCAGGGCGGCGGACGTGGTATTGACCAACATATGCTGGCAACACTGCGTGATGCTGTTCTCGGGGATTTTCGCCCCGACTTa                    <  1:940738/130‑1 (MQ=255)                                                                                                               cGGCACCTGGGTAATTGGCGATCGCCACGATCTCTCCACTCAGGCGTATCAGGGCGGCGGACGTGGTATTGACCAACATATGCTGGCAACACTGCGTGATGCTGTTCTCGGGGATTTTCGCCCCGACTTa                    >  2:940738/1‑130 (MQ=255)                                                                                                               cGGCACCTGGGTAATTGGCGATCGCCACGATCTCTCCACTCAGGCGTATCAGGGCGGCGGACGTGGTATTGACCAACATATGCTGGCAACACTGCGTGATGCTGTTCTCGGGGATTTTCGCCCCGACTTa                    >  2:940706/1‑130 (MQ=255)                                                                                                               cGGCACCTGGGTAATTGGCGATCGCCACGATCTCTCCACTCAGGCGTATCAGGGCGGCGGACGTGGTATTGACCAACATATGCTGGCAACACTGCGTGATGCTGTTCTCGGGGATTTTCGCCCCGACTTa                    <  1:940706/130‑1 (MQ=255)                                                                                                                    ccTGGGTGATTGGCGATCGCCACGATCTCTCCACTCAGGCGTATCAGGGCGGCGGACGTGGTATTGACCAACATATGCTGGCAACACTGCGTGATGCTGTTCTCGGGGATTTTCGCCCCGACTTAACGCTCTATCTCGAt     >  1:560749/1‑140 (MQ=255)                                                                                                                      tGGGTGATTGGCGATCGCCACGATCTCTCCACTCAGGCGTATCAGGGCGGCGGACGTGGTATTGACCAACATATGCTGGCAACACTGCGTGATGCTGTTCTCGGGGATTTTCGCCCCGACTTAACGCTCTATCTCGATGt   >  1:105830/1‑140 (MQ=255)                                                                                                                      tGGGTGATTGGCGATCGCCACGATCTCTCCACTCAGGCGTATCAGGGCGGCGGACGTGGTATTGACCAACATATGCTGGCAACACTGCGTGATGCTGTTCTCGGGGATTTTCGCCCCGACTTAACGCTCTATCTCGATGt   <  1:478106/140‑1 (MQ=255)                                                                                                                      tGGGTGATTGGCGATCGCCACGATCTCTCCACTCAGGCGTATCAGGGCGGCGGACGTGGTATTGACCAACATATGCTGGCAACACTGCGTGATGCTGTTCTCGGGGATTTTCGCCCCGACTTAACGCTCTATCTCGATGt   <  1:973608/140‑1 (MQ=255)                                                                                                                       gggTGATTGGCGATCGCCACGATCTCTCCACTCAGGCGTATCAGGGCGGCGGACGTGGTATTGACCAACATATGCTGGCAACACTGCGTGATGCTGTTCTCGGGGATTTTCGCCCCGACTTAACGCTCTATCTCGATGtt  <  1:412373/140‑1 (MQ=255)                                                                                                                           |                                                                                                                                        GATATCAAATCGGTAGGCGATGAAGTCATTACCGATAAAGCCGAAGTTCTGATGTTTTATGCCGCGCGCGTTCAACTGGTAGAAACGGTCATCAAACCAGCTCTGGCTAACGGCACCTGGGTGATTGGCGATCGCCACGATCTCTCCACTCAGGCGTATCAGGGCGGCGGACGTGGTATTGACCAACATATGCTGGCAACACTGCGTGATGCTGTTCTCGGGGATTTTCGCCCCGACTTAACGCTCTATCTCGATGTT  >  NC_000913/1155289‑1155546

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