Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	64,352	G→C	R477G (CGC→GGC)	polB ←	DNA polymerase II

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	64,352	0	G	C	97.6%	151.5 / NA	42	R477G (CGC→GGC)	polB	DNA polymerase II
Reads supporting (aligned to +/- strand):  ref base G (0/1);  new base C (16/25);  total (16/26)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ATGACCACGCATGGTGATCGACGATGCCAGCCGCGGATCGAAGAAGCGGCAGGCGGTGGTGCCGAGCACGCCATAAAAGGCATTCATGATGATTTTCAGCGCCTGCGACAGCGGTTTGTTACCCTGGCGTTTGGCTTCATCGCGCCCGTGCCAGATGTTAGTCACAATCTCCGGCAGGCAATGTTTTTCTCGCGAGAACCAGGCATCGAGAAAACCTTCGGTACTGTGCTCTGGATCAGGCTGCGCCATGCCTTCCACCAGCCCGACGGGATCAATCAGAAAGGTGCGGA  >  NC_000913/64209‑64498                                                                                                                                                |                                                                                                                                                   atgaCCACGCATGGTGATCGACGATGCCAGCCGCGGATCGAAGAAGCGGCAGGCGGTGGTGCCGAGCACGCCATAAAAGGCATTCATGATGATTTTCAGCGCCTGCGACAGCGGTTTGTTACCCTGGCGTTTGGCTTCATCGCCCCCGTg                                                                                                                                              <  2:480032/150‑1 (MQ=255)          cATGGTGATCGACGATGCCAGCCGCGGATCGAAGAAGCGGCAGGCGGTGGTGCCGAGCACGCCATAAAAGGCATTCATGATGATTTTCAGCGCCTGCGACAGCGGTTTGTTACCCTGGCGTTTGGCTTCATCGCCCCCGTGCCCGATGtt                                                                                                                                     >  1:599866/1‑150 (MQ=255)              gTGATCGACGATGCCAGCCGCGGATCGAAGAAGCGGCAGGCGGTGGTGCCGAGCACGCCATAAAAGGCATTCATGATGATTTTCAGCGCCTGCGACAGCGGTTTGTTACCCTGGCGTTTGGCTTCATCGCCCCCGTGCCAGATGTTAGTc                                                                                                                                 >  1:433083/1‑150 (MQ=255)                gATCGACGATGCCAGCCGCGGATCGAAGAAGCGGCAGGCGGTGGTGCCGAGCACGCCATAAAAGGCATTCATGATGATTTTCAGCGCCTGCGACAGCGGTTTGTTACCCTGGCGTTTGGCTTCATCGCCCCCGTGCCAGATGTTAGTcaca                                                                                                                              <  1:710222/151‑1 (MQ=255)                   cgacgaTGCCAGCCGCGGATCGAAGAAGCGGCAGGCGGTGGTGCCGAGCACGCCATAAAAGGCATTCATGATGATTTTCAGCGCCTGCGACAGCGGTTTGTTACCCTGGCGTTTGGCTTCATCGCCCCCGTGCCAGATGTTAGTCACAAtc                                                                                                                           >  2:323488/1‑151 (MQ=255)                      cgaTGCCAGCCGCGGATCGAAGAAGCGGCAGGCGGTGGTGCCGAGCACGCCATAAAAGGCATTCATGATGATTTTCAGCGCCTGCGACAGCGGTTTGTTACCCTGGCGTTTGGCTTCATCGCCCCCGTGCCAGATGTTAGTCACAAtct                                                                                                                          <  1:503640/149‑1 (MQ=255)                                   ggATCGAAGAAGCGGCAGGCGGTGGTGCCGAGCACGCCATAAAAGGCATTCATGATGATTTTCAGCGCCTGCGACAGCGGTTTGTTACCCTGGCGTTTGGCTTCATCGCCCCCGTGCCAGATGTTAGTCACAATCTCCGGCAGGCAATGt                                                                                                            <  2:387323/150‑1 (MQ=255)                                        gaagaaGCGGCAGGCGGTGGTGCCGAGCACGCCATAAAAGGCATTCATGATGATTTTCAGCGCCTGCGACAGCGGTTTGTTACCCTGGCGTTTGGCTTCATCGCCCCCGTGCCAGATGTTAGTCACAATCTCCGGCAGGCAATGTTTTtct                                                                                                      >  1:81807/1‑151 (MQ=255)                                        gaagaaGCGGCAGGCGGTGGTGCCGAGCACGCCATAAAAGGCATTCATGATGATTTTCAGCGCCTGCGACAGCGGTTTGTTACCCTGGCGTTTGGCTTCATCGCCCCCGTGCCAGATGTTAGTCACAATCTCCGGCAGGCAATGTTTTtct                                                                                                      >  2:189750/1‑151 (MQ=255)                                         aagaagCGGCAGGCGGTGGTGCCGAGCACGCCATAAAAGGCATTCATGATGATTTTCAGCGCCTGCGACAGCGGTTTGTTACCCTGGCGTTTGGCTTCATCGCCCCCGTGCCAGATGTTAGTCACAATCTCCGGCAGGCAATGTTTTtct                                                                                                      <  1:441439/150‑1 (MQ=255)                                          agaagCGGCAGGCGGTGGTGCCGAGCACGCCATAAAAGGCATTCATGATGATTTTCAGCGCCTGCGACAGCGGTTTGTTACCCTGGCGTTTGGCTTCATCGCCCCCGTGCCAGATGTTAGTCACAATCTCCGGCAGGCAATGTTTTTCTcg                                                                                                    <  2:370913/151‑1 (MQ=255)                                            aagCGGCAGGCGGTGGTGCCGAGCACGCCATAAAAGGCATTCATGATGATTTTCAGCGCCTGCGACAGCGGTTTGTTACCCTGGCGTTTGGCTTCATCGCCCCCGTGCCAGATGTTAGTCACAATCTCCGGCAGGCAATGTTTTTCTcgcg                                                                                                  >  1:691636/1‑151 (MQ=255)                                                ggCAGGCGGTGGTGCCGAGCACGCCATAAAAGGCATTCATGATGATTTTCAGCGCCTGCGACAGCGGTTTGTTACCCTGGCGTTTGGCTTCATCGCCCCCGTGCCAGATGTTAGTCACAATCTCCGGCAGGCAATGTTTTTCTCGCgaga                                                                                               <  1:425465/150‑1 (MQ=255)                                                  cAGGCGGTGGTGCCGAGCACGCCATAAAAGGCATTCATGATGATTTTCAGCGCCTGCGACAGCGGTTTGTTACCCTGGCGTTTGGCTTCATCGCCCCCGTGCCAGATGTTAGTCACAATCTCCGGCAGGCAATGTTTTTCTCGCGAGAAcc                                                                                            <  2:45000/151‑1 (MQ=255)                                                    ggCGGTGGTGCCGAGCACGCCATAAAAGGCATTCATGATGATTTTCAGCGCCTGCGACAGCGGTTTGTTACCCTGGCGTTTGGCTTCATCGCCCCCGTGCCAGATGTTAGTCACAATCTCCGGCAGGCAATGTTTTTCTCGCGAGAACCa                                                                                           <  2:199782/150‑1 (MQ=255)                                                    ggCGGTGGTGCCGAGCACGCCATAAAAGGCATTCATGATGATTTTCAGCGCCTGCGACAGCGGTTTGTTACCCTGGCGTTTGGCTTCATCGCCCCCGTGCCAGATGTTAGTCACAATCTCCGGCAGGCAATGTTTTTCTCGCGAGAACCa                                                                                           <  2:617235/150‑1 (MQ=255)                                                    ggCGGTGGTGCCGAGCACGCCATAAAAGGCATTCATGATGATTTTCAGCGCCTGCGACAGCGGTTTGTTACCCTGGCGTTTGGCTTCATCGCCCCCGTGCCAGATGTTAGTCACAATCTCCGGCAGGCAATGTTTTTCTCGCGAGAACCAg                                                                                          <  2:295959/151‑1 (MQ=255)                                                          ggtgCCGAGCACGCCATAAAAGGCATTCATGATGATTTTCAGCGCCTGCGACAGCGGTTTGTTACCCTGGCGTTTGGCTTCATCGCCCCCGTGCCAGATGTTAGTCACAATCTCCGGCAGGCAATGTTTTTCTCGCGAGAACCAGGCATCg                                                                                    <  1:596548/151‑1 (MQ=255)                                                            tgCCGAGCACGCCATAAAAGGCATTCATGATGATTTTCAGCGCCTGCGACAGCGGTTTGTTACCCTGGCGTTTGGCTTCATCGCCCCCGTGCCAGATGTTAGTCACAATCTCCGGCAGGCAATGTTTTTCTCGCGAGAACCAGGCATCg                                                                                    >  1:32013/1‑149 (MQ=255)                                                                    aCGCCATAAAAGGCATTCATGATGATTTTCAGCGCCTGCGACAGCGGTTTGTTACCCTGGCGTTTGGCTTCATCGCCCCCGTGCCAGATGTTAGTCACAATCTCCGGCAGGCAATGTTTTTCTCGCGAGAACCAGGCATCGAGAAAACCtt                                                                          >  1:17153/1‑151 (MQ=255)                                                                           aaaaGGCATTCATGATGATTTTCAGCGCCTGCGACAGCGGTTTGTTACCCTGGCGTTTGGCTTCATCGCCCCCGTGCCAGATGTTAGTCACAATCTCCGGCAGGCAATGTTTTTCTCGCGAGAACCAGGCATCGAGAAAACCTTCGGTACt                                                                   <  2:285530/151‑1 (MQ=255)                                                                           aaaaGGCATTCATGATGATTTTCAGCGCCTGCGACAGCGGTTTGTTACCCTGGCGTTTGGCTTCATCGCCCCCGTGCCAGATGTTAGTCACAATCTCCGGCAGGCAATGTTTTTCTCGCGAGAACCAGGCATCGAGAAAACCTTCGGTACt                                                                   <  2:360855/151‑1 (MQ=255)                                                                                     cATGATGATTTTCAGCGCCTGCGACAGCGGTTTGTTACCCTGGCGTTTGGCTTCATCGCCCCCGTGCCAGATGTTAGTCACAATCTCCGGCAGGCAATGTTTTTCTCGCGAGAACCAGGCATCGAGAAAACCTTCGGTACTGTGCTCTGGa                                                         <  1:729227/151‑1 (MQ=255)                                                                                          tgatTTTCAGCGCATGCGACAGCGGTTTGTTACCCTGGCGTTTGGCTTCATCGCCCCCGTGCCAGATGTTAGTCACAATCTCCGGCAGGCAATGTTTTTCTCGCGAGAACCAGGCATCGAGAAAACCTTCGGTACTGTGCTCTGGATCAgg                                                    >  1:104399/1‑151 (MQ=255)                                                                                            atTTTCAGCGCCTGCGACAGCGGTTTGTTACCCTGGCGTTTGGCTTCATCGCCCCCGTGCCAGATGTTAGTCACAATCTCCGGCAGGCAATGTTTTTCTCGCGAGAACCAGGCATCGAGAAAACCTTCGGTACTGTGCTCTGGATCAGGCt                                                  >  1:54282/1‑151 (MQ=255)                                                                                             ttttCAGCGCCTGCGACAGCGGTTTGTTACCCTGGCGTTTGGCTTCATCGCCCCCGTGCCAGATGTTAGTCACAATCTCCGGCAGGCAATGTTTTTCTCGCGAGAACCAGGCATCGAGAAAACCTTCGGTACTGTGCTCTGGATCAGGCTg                                                 >  2:659109/1‑151 (MQ=255)                                                                                              tttCAGCGCCTGCGACAGCGGTTTGTTACCCTGGCGTTTGGCTTCATCGCCCCCGTGCCAGATGTTAGTCACAATCTCCGGCAGGCAATGTTTTTCTCGCGAGAACCAGGCATCGAGAAAACCTTCGGTACTGTGCTCTGGATCAGGCTg                                                 <  1:506367/150‑1 (MQ=255)                                                                                               ttCAGCGCCTGCGACAGCGGTTTGTTACCCTGGCGTTTGGCTTCATCGCCCCCGTGCCAGATGTTAGTCACAATCTCCGGCAGGCAATGTTTTTCTCGCGAGAACCAGGCATCGAGAAAACCTTCGGTACTGTGCTCTGGATCAGGCTgc                                                >  2:602566/1‑150 (MQ=255)                                                                                                   gcgcCTGCGACAGCGGTTTGTTACCCTGGCGTTTGGCTTCATCGCCCCCGTGCCAGATGTTAGTCACAATCTCCGGCAGGCAATGTTTTTCTCGCGAGAACCAGGCATCGAGAAAACCTTCGGTACTGTGCTCTGGATCAGGCTGCGCCAt                                           <  2:65606/151‑1 (MQ=255)                                                                                                      ccTGCGACAGCGGTTTGTTACCCTGGCGTTTGGCTTCATCGCGCCCGTGCCAGATGTTAGTCACAATCTCCGGCAGGCAATGTTTTTCTCGCGAGAACCAGGCATCGAGAAAACCTTCGGTACTGTGCTCTGGATCAGGCTGCGCCATGcc                                        <  1:81426/151‑1 (MQ=255)                                                                                                                 ggTTTGTTACCCTGGCGTTTGGCTTCATCGCCCCCGTGCCAGATGTTAGTCACAATCTCCGGCAGGCAATGTTTTTCTCGCGAGAACCAGGCATCGAGAAAACCTTCGGTACTGTGCTCTGGATCAGGCTGCGCCATGCCTTCCACCAGcc                             >  2:427526/1‑151 (MQ=255)                                                                                                                   tttGTTACCCTGGCGTTTGGCTTCATCGCCCCCGTGCCAGATGTTAGTCACAATCTCCGGCAGGCAATGTTTTTCTCGCGAGAACCAGGCATCGAGAAAACCTTCGGTACTGTGCTCTGGATCAGGCTGCGCCATGCCTTCCACCAGCCCg                           >  1:216410/1‑151 (MQ=255)                                                                                                                    ttGTTACCCTGGCGTTTGGCTTCATCGCCCCCGTGCCAGATGTTAGTCACAATCTCCGGCAGGCAATGTTTTTCTCGCGAGAACCAGGCATCGAGAAAACCTTCGGTACTGTGCTCTGGATCAGGCTGCGCCATGCCTTCCACCAGCCCg                           <  1:592816/150‑1 (MQ=255)                                                                                                                      gTTACCCTGGCGTTTGGCTTCATCGCCCCCGTGCCAGATGTTAGTCACAATCTCCGGCAGGCAATGTTTTTCTCGCGAGAACCAGGCATCGAGAAAACCTTCGGTACTGTGCTCTGGATCAGGCTGCGCCATGCCTTCCACCAGCCCGAc                         <  2:584087/150‑1 (MQ=255)                                                                                                                      gTTACCCTGGCGTTTGGCTTCATCGCCCCCGTGCCAGATGTTAGTCACAATCTCCGGCAGGCAATGTTTTTCTCGCGAGAACCAGGCATCGAGAAAACCTTCGGTACTGTGCTCTGGATCAGGCTGCGCCATGCCTTCCACCAGCCCGACg                        <  2:545042/151‑1 (MQ=255)                                                                                                                        tACCCTGGCGTTTGGCTTCATCGCCCCCGTGCCAGATGTTAGTCACAATCTCCGGCAGGCAATGTTTTTCTCGCGAGAACCAGGCATCGAGAAAACCTTCGGTACTGTGCTCTGGATCAGGCTGCGCCATGCCTTCCACCAGCCCGACgg                       <  1:323420/150‑1 (MQ=255)                                                                                                                              ggCGTTTGGCTTCATCGCCCCCGTGCCAGATGTTAGTCACAATCTCCGGCAGGCAATGTTTTTCTCGCGAGAACCAGGCATCGAGAAAACCTTCGGTACTGTGCTCTGGATCAGGCTGCGCCATGCCTTCCACCAGCCCGACGGGatcaa                 <  1:72476/150‑1 (MQ=255)                                                                                                                              ggCGTTTGGCTTCATCGCCCCCGTGCCAGATGTTAGTCACAATCTCCGGCAGGCAATGTTTTTCTCGCGAGAACCAGGCATCGAGAAAACCTTCGGTACTGTGCTCTGGATCAGGCTGCGCCATGCCTTCCACCAGCCCGACGGGatca                  <  2:552577/149‑1 (MQ=255)                                                                                                                                cttttGGCTTCATCGCCCCCGTGCCAGATGTTAGTCACAATCTCCGGCAGGCAATGTTTTTCTCGCGAGAACCAGGCATCGAGAAAACCTTCGGTACTGTGCTCTGGATCAGGCTGCGCCATGCCTTCCACCAGCCCGACGGGatcaatc               >  2:301586/3‑150 (MQ=255)                                                                                                                                 gTTTGGCTTCATCGCCCCCGTGCCAGATGTTAGTCACAATCTCCGGCAGGCAATGTTTTTCTCGCGAGAACCAGGCATCGAGAAAACCTTCGGTACTGTGCTCTGGATCAGGCTGCGCCATGCCTTCCACCAGCCCGACGGGatcaatc               >  2:73888/1‑149 (MQ=255)                                                                                                                                         tCATCGCCCCCGTGCCAGATGTTAGTCACAATCTCCGGCAGGCAATGTTTTTCTCGCGAGAACCAGGCATCGAGAAAACCTTCGGTACTGTGCTCTGGATCAGGCTGCGCCATGCCTTCCACCAGCCCGACGGGATCAATCAGAAAGGTg      <  2:223588/150‑1 (MQ=255)                                                                                                                                         tCATCGCCCCCGTGCCAGATGTTAGTCACAATCTCCGGCAGGCAATGTTTTTCTCGCGAGAACCAGGCATCGAGAAAACCTTCGGTACTGTGCTCTGGATCAGGCTGCGCCATGCCTTCCACCAGCCCGACGGGATCAATCAGAAAGGTg      <  1:6786/150‑1 (MQ=255)                                                                                                                                             cgccCCCGTGCCAGATGTTAGTCACAATCTCCGGCAGGCAATGTTTTTCTCGCGAGAACCAGGCATCGAGAAAACCTTCGGTACTGTGCTCTGGATCAGGCTGCGCCATGCCTTCCACCAGCCCGACGGGATCAATCAGAAAGGTGCGga  <  1:61341/150‑1 (MQ=255)                                                                                                                                                |                                                                                                                                                   ATGACCACGCATGGTGATCGACGATGCCAGCCGCGGATCGAAGAAGCGGCAGGCGGTGGTGCCGAGCACGCCATAAAAGGCATTCATGATGATTTTCAGCGCCTGCGACAGCGGTTTGTTACCCTGGCGTTTGGCTTCATCGCGCCCGTGCCAGATGTTAGTCACAATCTCCGGCAGGCAATGTTTTTCTCGCGAGAACCAGGCATCGAGAAAACCTTCGGTACTGTGCTCTGGATCAGGCTGCGCCATGCCTTCCACCAGCCCGACGGGATCAATCAGAAAGGTGCGGA  >  NC_000913/64209‑64498

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