Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	63,556	T→C	Y742C (TAC→TGC)	polB ←	DNA polymerase II

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	63,556	0	T	C	98.2%	190.2 / NA	55	Y742C (TAC→TGC)	polB	DNA polymerase II
Reads supporting (aligned to +/- strand):  ref base T (0/1);  new base C (28/26);  total (28/27)
Fisher's exact test for biased strand distribution p-value = 4.91e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TGCTCAAAATAGCCCAAGTTGCCCGGTCATAAGTGTAGCAAAATTATCCTCAATAAAAGGGAGTATTCCCTCCGCCACGGGTTGTAGCTGGCGGGTCAGATAGTGTTCGTAATCCAGTGGTGAACGTTGGTAGTCCAGCGGCTCCGGGCCGTTGGTGGTCCATACGTACTTAATGGTGCCGCGATTCTGATATTGCAAGGGGCGACCACGCTTTTGGTTTTCTTCATCGGCAAGGCGAGCGGCGCGTACATGAGGCGGCACATTACG  >  NC_000913/63426‑63692                                                                                                                                   |                                                                                                                                         cgCTCAAAATAGCCCAAGTTGCCCGGTCATAAGTGTAGCAAAATTATCCTCAATAAAAGGGAGTATTCCCTCCGCCACGGGTTGTAGCTGGCGGGTCAGATAGTGTTCGTAATCCAGTGGTGAACGTTGGCAGTCCAGCg                                                                                                                                 >  2:510396/2‑140 (MQ=255)    tCAAAATAGCCCAAGTTGCCCGGTCATAAGTGTAGCAAAATTATCCTCAATAAAAGGGAGTATTCCCTCCGCCACGGGTTGTAGCTGGCGGGTCAGATAGTGTTCGTAATCCAGTGGTGAACGTTGGCAGTCCAGCGGCt                                                                                                                              >  1:85988/1‑140 (MQ=255)      aaaaTAGCCCAAGTTGTCGGGTCATAAGTGTAGCAAAATTATCCTCAATAAAATGGAGTATTCCCTCCGCCCCGGGTTGTTGCTGGCGGGTCAGATAGTGTTCGTAATCCAGTGGTGAACGTTGGCAGTCCAGCGGCTcc                                                                                                                            <  2:96224/140‑1 (MQ=255)      aaaaTAGCCCAAGTTGCCCGGTCATAAGTGTAGCAAAATTATCCTCAATAAAAGGGAGTATTCCCTCCGCCACGGGTTGTAGCTGGCGGGTCAGATAGTGTTCGTAATCCAGTGGTGAACGTTGGCAGTCCAGCGGCTcc                                                                                                                            >  1:391932/1‑140 (MQ=255)                aaGTTGCCCGGTCATAAGTGTAGCAAAATTATCCTCAATAAAAGGGAGTATTCCCTCCGCCACGGGTTGTAGCTGGCGGGTCAGATAGTGTTCGTAATCCAGTGGTGAACGTTGGCAGTCCAGCGGCTCCGGGCCGTtgg                                                                                                                  <  1:351371/140‑1 (MQ=255)                        cGGTCATAAGTGTAGCAAAATTATCCTCAATAAAAGGGAGTATTCCCTCCGCCACGGGTTGTAGCTGGCGGGTCAGATAGTGTTCGTAATCCAGTGGTGAACGTTGGCAGTCCAGCGGCTCCGGGCCGTTGGTGGTCCAt                                                                                                          <  2:618785/140‑1 (MQ=255)                        cGGTCATAAGTGTAGCAAAATTATCCTCAATAAAAGGGAGTATTCCCTCCGCCACGGGTTGTAGCTGGCGGGTCAGATAGTGTTCGTAATCCAGTGGTGAACGTTGGCAGTCCAGCGGCTCCGGGCCGTTGGTGGTCCAt                                                                                                          <  1:733394/140‑1 (MQ=255)                        cGGTCATAAGTGTAGCAAAATTATCATCAATAAAAGGGAGTATTCCCTCCGCCACGGGTTGTAGCTGGCGGGTCAGATAGTGTTCGTAATCCAGTGGTGAACGTTGGCAGTCCAGCGGCTCCGGGCCGTTGGTGGTCCAt                                                                                                          <  2:574841/140‑1 (MQ=255)                         ggTCATAAGTGTAGCAAGATTACCCTCAATAAAAGTGAGTATTCCCTCCGCCACGGGTTGTAGCTGGCGGGTCAGATAGTGTTCGTAATCCAGTGGTGAACGTTGGCAGTCCAGCGGCTCCGGGCCGTTGGTGGTCCATa                                                                                                         <  2:241334/140‑1 (MQ=255)                          gTCATAAGTGTAGCAAACTTATTCTCAATAACAGGTAGTCTTCCCTCCCCCACGGGTTGTAGCTGGCGGGTCAGATAGTGTTCGTAATCCAGTGGTGAACGTTGGCAGTCCAGCGGCTCCGGGCCGTTGGTGGTCCATAc                                                                                                        <  2:191647/140‑1 (MQ=255)                          gTCATAAGTGTAGCAAAATTATCCTCAATCAAAGGGAATATTCCCTCCGCCACGGGTTGTAGCTGGCGGGTCAGATAGTGTTCGTAATCCAGTGGTGAACGTTGGCAGTCCAGCGGCTCCGGGCCGTTGGTGGTCCATAc                                                                                                        <  1:273241/140‑1 (MQ=255)                          gTCATAAGTGTAGCAAAATTATCCTCAATAAAAGGGAGTATTCCCTCCGCCACGGGTTGTAGCTGGCGGGTCAGATAGTGTTCGTAATCCAGTGGTGAACGTTGGCAGTCCAGCGGCTCCGGGCCGTTGGTGGTCCATAc                                                                                                        <  2:140160/140‑1 (MQ=255)                           tCATAAGTGTAGCAAAATTATCCTCAATAAAAGGGAGTATTCCCTCCGCCACGGGTTGTAGCTGGCGGGTCAGATAGTGTTCGTAATCCAGTGGTGAACGTTGGCAGTCCAGCGGCTCCGGGCCGTTGGTGGTCCATACg                                                                                                       <  1:469556/140‑1 (MQ=255)                             aTAAGTGTAGCAAAATTATCCTCAATAAAAGGGAGTATTCCCTCCGCCACGGGTTGTAGCTGGCGGGTCAGATAGTGTTCGTAATCCAGTGGTGAACGTTGGCAGTCCAGCGGCTCCGGGCCGTTGGTGGTCCATACGTa                                                                                                     <  1:494750/140‑1 (MQ=255)                                    tAGCAAAATTATCCTCAATAAAAGGGAGTATTCCCTCCGCCACGGGTTGTAGCTGGCGGGTCAGATAGTGTTCGTAATCCAGTGGTGAACGTTGGCAGTCCAGCGGCTCCGGGCCGTTGGTGGTCCATACGTACTTAATg                                                                                              >  1:62308/1‑140 (MQ=255)                                     aGCAAAATTATCCTCAATAAAAGGGAGTATTCCCTCCGCCACGGGTTGTAGCTGGCGGGTCAGATAGTGTTCGTAATCCAGTGGTGAACGTAGGCAGTCCAGCGGCTCTGGGCCGTTGGTGGGCCATACGTACTTAATgg                                                                                             >  2:293677/1‑140 (MQ=255)                                       cAAAATTATCCTCAATAAAAGGGAGTATTCCCTCCGCCACGGGTTGTAGCTGGCGGGTCAGATAGTGTTCGTAATCCAGTGGTGAACGTTGGCAGTCCAGCGGCTCCGGGCCGTTGGTGGTCCATACGTACTTAATgg                                                                                             >  1:33492/1‑138 (MQ=255)                                       cAAAATTATCCTCAATAAAAGGGAGTATTCCCTCCGCCACGGGTTGTAGCTGGCGGGTCAGATAGTGTTCGTAATCCAGTGGTGAACGTTGGCAGTCCAGCGGCTCCGGGCCGTTGGTGGTCCATACGTACTTAATgg                                                                                             <  2:33492/138‑1 (MQ=255)                                          aaTTATCCTCAATAAAAGGGAGTATTCCCTCCGCCACGGGTTGTAGCTGGCGGGTCAGATAGTGTTCGTAATCCAGTGGTGAACGTTGGCATTCCAGCGGGTCCGGGCCGTTGGTGGACCATTAGTACTTACTGGTGCcg                                                                                        >  2:330113/1‑140 (MQ=255)                                          aaTTATCCTCAATAAAAGGGAGTATTCCCTCCGCCACGGGTTGTAGCTGGCGGGTCAGATAGTGTTCGTAATCCAGTGGTGAACGTTGGCAGTCCAGCGGCTCCGGGCCGTTGGTGGTCCATACGTACTTAATGGTGCcg                                                                                        <  2:85988/140‑1 (MQ=255)                                            ttATCCTCAATAAAAGGGAGTATTCCCTCCGCCACGGGTTGTAGCTGGCGGGTCAGATAGTGTTCGTAATCCAGTGGTGAACGTTGGCAGTCCAGCGGCTCCGGGCCGTTGGTGGTCCATACGTACTTAATGGTGCcgcg                                                                                      >  2:695396/1‑140 (MQ=255)                                                    aaTAAAAGGGAGTATTCCCTCCGCCACGGGTTGTAGCTGGCGGGTCAGATAGTGTTCGTAATCCAGTGGTGAACGTTGGCAGTCCAGCGGCTCCGGGCCGTTGGTGGTCCATACGTACTTAATGGTGCCGCGATTCTGat                                                                              >  2:721064/1‑140 (MQ=255)                                                     aTAAAAGGGAGTATTCCCTCCGCCACGGGTTGTAGCTGGCGGGTCAGATAGTGTTCGTAATCCAGTGGTGAACGTTGGCAGTCCAGCGGCTCCGGGCCGTTGGTGGTCCATACGTACTTAATGGTGCCGCGATTc                                                                                  >  2:493833/1‑135 (MQ=255)                                                     aTAAAAGGGAGTATTCCCTCCGCCACGGGTTGTAGCTGGCGGGTCAGATAGTGTTCGTAATCCAGTGGTGAACGTTGGCAGTCCAGCGGCTCCGGGCCGTTGGTGGTCCATACGTACTTAATGGTGCCGCGATTc                                                                                  <  1:493833/135‑1 (MQ=255)                                                           gggAGTATTCCCTCCGCCACGGGTTGTAGCTGGCGGGTCAGATAGTGTTCGTAATCCAGTGGTGAACGTTGGCAGTCCAGCGGCTCCGGGCCGTTGGTGGTCCATACGTACTTAATGGTGCCGCGATTCTGATATTGCaa                                                                       <  2:375410/140‑1 (MQ=255)                                                              aGTATTCCCTCCGCCACGGGTTGTAGCTGGCGGGTCAGATAGTGTTCGTAATCCAGTGGTGAACGTTGGCAGTCCAGCGGCTCCGGGCCGTCGGTGGTCCATACCTACATAATGGTGCCGCGGTTCTGATATTGCAAggg                                                                    >  2:256420/1‑140 (MQ=255)                                                                tATTCCCTCCGCCACGGGTTGTAGCTGGCGGGTCAGATAGTGTTCGTAATCCAGTGGTGAACGTTGGCAGTCCAGCGGCTCCGGGCCGTTGGTGGTCCATACGTACTTAATGGTGCCGCGATTCTGATATTGCa                                                                        >  2:48/1‑134 (MQ=255)                                                                tATTCCCTCCGCCACGGGTTGTAGCTGGCGGGTCAGATAGTGTTCGTAATCCAGTGGTGAACGTTGGCAGTCCAGCGGCTCCGGGCCGTTGGTGGTCCATACGTACTTAATGGTGCCGCGATTCTGATATTGCa                                                                        <  1:48/134‑1 (MQ=255)                                                                        ccGCCACGGGTTGTAGCTGGCGGGTCAGATAGTGTTCGTAATCCAGTGGTGAACGTTGGCAGTCCAGCGGCTCCGGGCCGTTGGTGGTCCATACGTACTTAATGGTGCCGCGATTCTGATATTGCAAGGGGCGACCACGc                                                          >  1:582827/1‑140 (MQ=255)                                                                         cGCCACGGGTTGTAGCTGGCGGGTCAGATAGTGTTCGTAATCCAGTGGTGAACGTTGGCAGTCCAGCGGCTCCGGGCCGTTGGTGGTCCATACGTACTTAATGGTGCCGCGATTCTGATATTGCAAGGGGCGACCACGCt                                                         >  2:654604/1‑140 (MQ=255)                                                                          gCCACGGGTTGTAGCTGGCGGGTCAGATAGTGTTCGTAATCCAGTGGTGAACGTTGGCAGTCCAGCGGCTCCGGGCCGTTGGTGGTCCATACGTACTTAATGGTGCCGCGATTCTGATATTGCAAGGGGCGACCACGCtt                                                        >  1:445701/1‑140 (MQ=255)                                                                              cGGGTTGTAGCTGGCGGGTCAGATAGTGTTCGTAATCCAGTGGTGAACGTTGGCAGTCCAGCGGCTCCGGGCCGTTGGTGGTCCATACGTACTTAATGGTGCCGCGATTCTGATATTGCAAGGGGCGACCACGCTTTTgg                                                    >  1:565692/1‑140 (MQ=255)                                                                              cGGGTTGTAGCTGGCGGGTCAGATAGTGTTCGTAATCCAGTGGTGAACGTTGGCAGTCCAGCGGCTCCGGGCCGTTGGTGGTCCATACGTACTTAATGGTGCCGCGATTCTGATATTGCAAGGGGCGACCACGCTTTTgg                                                    >  2:63689/1‑140 (MQ=255)                                                                                        cTGGCGGGTCAGATAGTGTTCGTAATCCAGTGGTGAACGTTGGCAGTCCAGCGGCTCCGGGCCGTTGGTGGTCCATACGTACTTAATGGTGCCGCGATTCTGATATTGCAAGGGGCGACCACGCTTTTGGTTTTCTTCAt                                          >  1:571651/1‑140 (MQ=255)                                                                                          ggCGGGTCAGATAGTGTTCGTAATCCAGTGGTGAACGTTGGCAGTCCAGCGGCTCCGGGCCGTTGGTGGTCCATACGTACTTAATGGTGCCGCGATTCTGATATTGCAAGGGGCGACCACGCTTTTGGTTTTCTTCATCg                                        >  2:504044/1‑140 (MQ=255)                                                                                           gCGGGTCAGATAGTGTTCGTAATCCAGTGGTGAACGTTGGCAGTCCAGCGGCTCCGGGCCGTTGGTGGTCCATACGTACTTAATGGTGCCGCGATTCTGATATTGCAAGGGGCGACCACGCTTTTGGTTTTCTTCATCgg                                       >  2:56618/1‑140 (MQ=255)                                                                                           gCGGGTCAGATAGTGTTCGTAATCCAGTGGTGAACGTTGGCAGTCCAGCGGCTCCGGGCCGTTGGTGGTCCATACGTACTTAATGGTGCCGCGATTCTGATATTGCAAGGGGCGACCACGCTTTTGGTTTTCTTCATCgg                                       >  2:557466/1‑140 (MQ=255)                                                                                               gTCAGATAGTGTTCGTAATCCAGTGGTGAACGTTGGCAGTCCAGCGGCTCCGGGCCGTTGGTGGTCCATACGTACTTTATGGTGCCGCGATTCTGATATTGCTAGGGGCGACCAGGCGTTTGGTATTCTTCATCGGCAAg                                   >  1:263069/1‑140 (MQ=255)                                                                                                   gATAGTGTTCGTAATCCAGTGGTGAACGTTGGCAGTCCAGCGGCTCCGGGCCGTTGGTGGTCCATACGTACTTAATGGTGCCGCGATTCTGATATTGCAAGGGGCGACCACGCTTTTGGTTTTCTTCATCGGCAAGGCGa                               <  2:391932/140‑1 (MQ=255)                                                                                                    aTAGTGTTCGTAATCCAGTGGTGAACGTTGGCAGTCCAGCGGCTCCGGGCCGTTGGTGGTCCATACGTACTTAATGGTGCCGCGATTCTGATATTGCAAGGGGCGACCACGCTTTTGGTTTTCTTCATCGGCAAGGCGAg                              <  1:504044/140‑1 (MQ=255)                                                                                                    aTAGTGTTCGTAATCCAATGGTTAACGTTGGTAGTTCAGCGGCTCCGGCCCGTTGGTGGTCCATACGTACTTAATGGTGCCGCGATTCTGATATTGCAAGGGGCGACCACGCTTTTGGTTTTCTTCATCGGCAAGGCGAg                              <  2:233973/140‑1 (MQ=255)                                                                                                        tgtTCGTAATCCAGTGGTGAACGTTGGCAGTCCAGCGGCTCCGGGCCGTTGGTGGTCCATACGTACTTAATGGTGCCGCGATTCTGATATTGCAAGGGGCGACCACGCTTTTGGTTTTCTTCATCGGCAAGGCGAgcggc                          <  2:445701/140‑1 (MQ=255)                                                                                                        tgtTCGTAATCCAGTGGTGAACGTTGGCAGTCCAGCGGCTCCGGGCCGTTGGTGGTCCATACGTACTTAATGGTGCCGCGATTCTGATATTGCAAGGGGCGACCACGCTTTTGGTTTTCTTCATCGGCAAGGCGAgcggc                          <  1:56618/140‑1 (MQ=255)                                                                                                        tgtTCGTAATCCAGTGGTGAACGTTGGCAGTCCAGCGGCTCCGGGCCGTTGGTGGTCCATACGTACTTAATGGTGCCGCGATTCTGATATTGCAAGGGGCGACCACGCTTTTGGTTTTCTTCATCGGCAAGGCGAgcggc                          >  1:37835/1‑140 (MQ=255)                                                                                                         gtTCGTAATCCAGTGGTGAACGTTGGCAGTCCAGCGGCTCCGGGCCGTTGGTGGTCCATACGTACTTAATGGTGCCGCGATTCTGATATTGCAAGGGGCGACGACGCTTTGGGTGTTGTTCATCTGCAAGGCGAgcggcg                         >  1:233973/1‑140 (MQ=255)                                                                                                              tAATCCAGTGGTGAACGTTGGCAGTCCAGCGGCTCCGGGCCGTTGGTGGTCCATACGTACTTAATGGTGCCGCGATTCTGATATTGCAAGGGGCGACCACGCTTTTGGTTTTCTTCATCGGCAAGGCGAgcggcg                         <  1:466796/135‑1 (MQ=255)                                                                                                              tAATCCAGTGGTGAACGTTGGCAGTCCAGCGGCTCCGGGCCGGTGGTGGTCCATACGTACTTAATGGTCCCGCGATTCTGATATTGCAAGGGGCGACCACGCTTTTGGTTTTCTTCATCGGCAAGGCGAgcggcg                         >  2:466796/1‑135 (MQ=255)                                                                                                                 tCCAGTGGTGAACGTTGGCAGTCCAGCGGATCCGGGCTGTTGGTGGTCCATACGTACTTAATGGTGCCGCGATTCTGATATTGCAAGGGGCGACCACGCTTTTGGTTTTCTTCATCGGCAAGGCGAGTGGCGCGTACATg                 <  1:256420/140‑1 (MQ=255)                                                                                                                    aGTGGTGAACGTTGGCAGTCCAGCGGCTCCGGGCCGTTGGTGGTCCATACGTACTTAATGGTGCCGCGATTCTGATATTGCAAGGGGCGACCACGCTTTTGGTTTTCTTCATCGGCAAGGCGAGCGGCGCGTACATGAgg              >  1:561015/1‑140 (MQ=255)                                                                                                                      tggtgAACGTTGGCAGTCCAGCGGCTCCGGGCCGTTGGTGGTCCATACGTACTTAATGGTGCCGCGATTCTGATATTGCAAGGGGCGACCACGCTTTTGGTTTTCTTCATCGGCAAGGCGAGCGGCGCGTACATGAggcg            <  1:695396/140‑1 (MQ=255)                                                                                                                      tggtgAACGTTGGCAGTCCAGCGGCTCCGGGCCGTTGGTGGTCCATACGTACTTAATGGTGCCGCGATTCTGATATTGCAAGGGGCGACCACGCTTATGGTTTTCTTCATCGGCAAGGCGAGCGGCGCGTACATGAggcg            <  1:596374/140‑1 (MQ=255)                                                                                                                        gtgAACGTTGGCAGTCCAGCGGCTCCGGGCCGTTGGTGGTCCATACGTACTTAATGGTGCCGCGATTCTGATATTGCAAGGGGCGACCACGCTTTTGGTTTTCTTCATCGGCAAGGCGAGCGGCGCGTACATGAggcggc          <  1:654604/140‑1 (MQ=255)                                                                                                                           aaCGTTGGCAGTCCAGCGGCTCCGGGCCGTTGGTGGTCCATACGTACTTAATGGTGCCGCGATTCCGATATTGCAAGGGGCGACCACGCTTTTGGTTTTCTTCATCGGCAAGGCGAGCGGCGCGTACATGAGGCGGcaca       >  2:433878/1‑140 (MQ=255)                                                                                                                             cGTTGGCAGTCCAGCGGCTCCGGGCCGTTGGTGGTCCATACGTACTTAATGGTGCCGCGATTCTGATATTGCAAGGGGCGACCACGCTTTTGGTTTTCTTCATCGGCAAGGCGAGCGGCGCGTACATGAGGCGGCACAtt     <  1:510396/140‑1 (MQ=255)                                                                                                                                tGGCAGTCCAGCGGCTCCGGGCCGTTGGTGGTCCATACGTACTTAATGGTGCCGCGATTCTGATATTGCAAGGGGCGACCACGCTTTTGGTTTTCTTCATCGGCAAGGCGAGCGGCGCGTACATGAGGCGGCACATTACg  >  2:669577/1‑140 (MQ=255)                                                                                                                                   |                                                                                                                                         TGCTCAAAATAGCCCAAGTTGCCCGGTCATAAGTGTAGCAAAATTATCCTCAATAAAAGGGAGTATTCCCTCCGCCACGGGTTGTAGCTGGCGGGTCAGATAGTGTTCGTAATCCAGTGGTGAACGTTGGTAGTCCAGCGGCTCCGGGCCGTTGGTGGTCCATACGTACTTAATGGTGCCGCGATTCTGATATTGCAAGGGGCGACCACGCTTTTGGTTTTCTTCATCGGCAAGGCGAGCGGCGCGTACATGAGGCGGCACATTACG  >  NC_000913/63426‑63692

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