Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	2,964,941	A→G	9.1%	Y71Y (TAT→TAC)	thyA ←	thymidylate synthetase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,964,941	0	A	G	9.1%	118.1 / 2.4	44	Y71Y (TAT→TAC)	thyA	thymidylate synthetase
Reads supporting (aligned to +/- strand):  ref base A (20/20);  new base G (2/2);  total (22/22)
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.60e-01

TTTTCAGCTGGTTCAGTACCGTAGTGATCTGGTCAATATGACGACCATCTGGCGTTGGCCAGGCGCGCCACTGTTTACCATACACTGGCCCGAGGTCGCCGTTTTCATCGGCCCATTCGTCCCAGATGGTGACATTGTTTTCGTGTAGATAAGCAATGTTAGTGTCGCCCTGCAGAAACCACAGCAGTT‑CATGGATGATGGAACGCAGGTGGCAACGTTTAGTTGTCACCAGCGGGAATCCATCTTGCAGGTTAAAACGCATCTGATGACCAAAAATGGAAAGCGTTCCGG  >  NC_000913/2964793‑2965083                                                                                                                                                     |                                                                                                                                                ttttCAGCTGGTTCAGTACCGTAGTGATCTGGTCAATATGACGACCATCTGGCGTTGGCCAGGCGCGCCACTGTTTACCATACACTGGCCCGAGGTCGCCGTTTTCATCGGCCCATTCGTCCCAGATGGTGACATTGTTTTCGTGTAGa                                                                                                                                                 >  2:731613/1‑149 (MQ=255)    tCAGCTGGTTCAGTACCGTAGTGATCTGGTCAATATGACGACCATCTGGCGTTGGCCAGGCGCGCCACTGTTTACCATACACTGGCCCGAGGTCGCCGTTTTCATCGGCCCATTCGTCCCAGATGGTGACATTGTTTTCGTGTAGATaa                                                                                                                                              <  1:351897/149‑1 (MQ=255)     cAGCTGGTTCAGTACCGTAGTGATCTGGTCAATATGACGACCATCTGGCGTTGGCCAGGCGCGCCACTGTTTACCATACACTGGCCCGAGGTCGCCGTTTTCATCGGCCCATTCGTCCCAGATGGTGACATTGTTTTCGTGTAGATAAg                                                                                                                                             <  2:630405/149‑1 (MQ=255)      aGCTGGTTCAGTACCGTAGTGATCTGGTCAATATGACGACCATCTGGCGTTGGCCAGGCGCGCCACTGTTTACCATACACTGGCCCGAGGTCGCCGTTTTCATCGGCCCATTCGTCCCAGATGGTGACATTGTTTTCGTGTAGATAAGc                                                                                                                                            >  2:734543/1‑149 (MQ=255)               aGTACCGTAGTGATCTGGTCAATATGACGACCATCTGGCGTTGGCCAGGCGCGCCACTGTTTACCATACACTGGCCCGAGGTCGCCGTTTTCATCGGCCCATTCGTCCCAGATGGTGACATTGTTTTCGTGTAGATAAGCAATGTTAgt                                                                                                                                   <  1:287154/149‑1 (MQ=255)               aGTACCGTAGTGATCTGGTCAATATGACGACCATCTGGCGTTGGCCAGGCGCGCCACTGTTTACCATACACTGGCCCGAGGTCGCCGTTTTCATCGGCCCATTCGTCCCAGATGGTGACATTGTTTTCGTGTAGATAAGCAATGTTAgt                                                                                                                                   <  1:735719/149‑1 (MQ=255)                 tACCGTAGTGATCTGGTCAATATGACGACCATCTGGCGTTGGCCAGGCGCGCCACTGTTTACCATACACTGGCCCGAGGTCGCCGTTTTCATCGGCCCATTCGTCCCAGATGGTGACATTGTTTTCGTGTAGGTAAGCAATGTTAgtgt                                                                                                                                 >  2:305826/1‑149 (MQ=255)                 tACCGTAGTGATCTGGTCAATATGACGACCATCTGGCGTTGGCCAGGCGCGCCACTGTTTACCATACACTGGCCCGAGGTCGCCGTTTTCATCGGCCCATTACTCCCAGAAGGTGACATCGTTTTCGTGTAGGTAAGCAATGTTAgtgt                                                                                                                                 >  2:86943/1‑149 (MQ=255)                   ccGGAGTGATCTGGTCAATATGACGACCATCTGGCGTTGGCCAGGCGCGCCACTGTTTACCATACACTGGCCCGAGGTCGCCGTTTTCATCGGCCCATTCGTCCCAGATGGTGACATTGTTTTCGTGTAGATAAGCAATGTTAGTGTCg                                                                                                                               <  1:165554/149‑1 (MQ=255)                     gaaGTGATGCGGCCAAGATGACGACCATCTGGCGTTGGCCAGGCGCGCCACTGTTTACCATACACTGGCCCGAGGTCGCCGTTTTCATCGGCCCATTCGTCCCAGATGGTGACATTGTTTTCGTGTAGATAAGCAATGTTAGTGTCGcc                                                                                                                             <  1:14759/147‑1 (MQ=255)                        gtgtTCTGGTCAATCGGACGACCAGCTGGCGTTGGCCAGGCGCGCCACTGTTTACCATACACTGGCCCGAGGTCGCCGTTTTCATCGGCACATTCGTCCCAGATGGTGACATTGTTTTCGTGTAGATAAGCAATGTTAGTGTCGCCCTg                                                                                                                          <  2:56007/149‑1 (MQ=255)                         tGATCTGGTCAATATGACGACCATCTGGCGTTGGCCAGGCGCGCCACTGTTTACCATACACTGGCCCGAGGTCGCCGTTTTCATCGGCCCATTCGTCCCAGATGGTGCCCTTGTTTTCGTGTAGATAAGCAACGTTAGTGTCGCCCTGc                                                                                                                         >  2:667932/1‑149 (MQ=255)                         tGATCTGGTCAATATGACGACCATCTGGCGTTGGCCAGGCGCGCCACTGTTTACCATACACTGGCCCGAGGTCGCCGTTTTCATCGGCCCATTCGTCCCAGATGGTGACATTGTTTTCGTGTAGATAAGCAATGTTAGTGTCGCCCTGc                                                                                                                         >  2:269522/1‑149 (MQ=255)                          gATCTGGTCAATATGACGACCATCTGGCGTTGGCCAGGCGCGCCACTGTTTACCATACACTGGCCCGAGGTCGCCGTTTTCATCGGCCCTTTCGTCCCAGATGGTGACATTGTTTTCGTGTAGATAAGCAATGTTAGTGTCGCCCTGCa                                                                                                                        <  1:519652/149‑1 (MQ=255)                          gATCTGGTCAATATGACGACCATCTGGCGTTGGCCAGGCGCGCCACTGTTTACCATACACTGGCCCGAGGTCGCCGTTTTCATCGGCCCATTCGTCCCAGATGGTGACATTGTTTTCGTGTAGATAAGCAATGTTAGTGTCGCCCTGCa                                                                                                                        <  1:86720/149‑1 (MQ=255)                               ggTCAATATGACGACCATCTGGCGTTGGCCAGGCGCGCCACTGTTTACCATACACTGGCCCGAGGTCGCCGTTTTCATCGGCCCATTCGTCCCAGATGGTGACATTGTTTTCGTGTAGATAAGCAATGTTAGTGTCGCCCTGCAGAAAc                                                                                                                   <  1:280421/149‑1 (MQ=255)                                 tCAATATGACGGCCATCTGGCGTTGGCCAGGCGCGCCACTGTTTACCATACACTGGCCCGAGGTCGCCGTTTTCATCGGCCCATTCGTCCCAGATGGTGACATTGTTTTCGTGTAGATAAGCAATGTTAGTGTCGCCCTGCAGAAACca                                                                                                                 <  2:46017/149‑1 (MQ=255)                                           gacCATCTGGCGTTGGCCGGGCGGGCCACTGTTTACCATACACTGGCCCGAGGTCGCCGTTTTCATCGGCCCATTCGTCCCAGATGGTGACATTGTTTTCGTGTAGATAAGCAATGTTAGTGTCGCCCTGCAGAAACCACAGCAGTT‑Ca                                                                                                      <  2:202276/149‑1 (MQ=255)                                            acCATCTGGCGTTGGCCAGGCGCGCCACTGTTTACCATACACTGGCCCGAGGTCGCCGTTTTCATCGGCCCATTCGTCCCAGATGGTGACATTGTTTTCGTGTAGATAAGCAATGTTAGTGTCGCCCTGCAGCAACCACATCAGTT‑CAt                                                                                                     <  2:678201/149‑1 (MQ=255)                                                 cTGGCGTTGGCCAGGCGCGCCACTGTTTACCATACACTGGCCCGAGGTCGCCGTTTTCATCGGCCCATTCGTCCCAGATGGTGACATTCTTTTCGTGTAGATAAGCAATGTTAGTGTCGCCCTGCAGAAACCACAGCAGTT‑CATGgatg                                                                                                <  2:557981/149‑1 (MQ=255)                                                    gCGTTGCCCAGGCGCTCCACCGTTTACCATACACTGGCCCGAGGTCGCCGTTTTCTTCGGCCCATTCGTCCCAGATGGTGACATTGTTTTCGTGTAGATAAGCAATGTTAGTGTCGCCCTGCAGAAACCACAGCAGTT‑CATGGatgatg                                                                                             <  2:698196/149‑1 (MQ=255)                                                          gCCAGGCGCGCCACTGTTTACCATACACTGGCCCGAGGTCGCCGTTTTCATCGGCCCATTCGTCCCAGATGGTGACATTGTTTTCGTGTAGATAAGCAATGTTAGTGTCGCCCTGCAGAAACCACAGCAGTT‑CATGGATGATGGAACGc                                                                                       >  1:348102/1‑149 (MQ=255)                                                          gCCAGGCGCGCCACTGTTTACCATACACTGGCCCGAGGTCGCCGTTTTCATCGGCCCATTCGTCCCAGATGGTGACATTGTTTTCGTGTAGATAAGCAATGTTAGTGTCGCCCTGCAGAAACCACAGCAGTT‑CATGGATGATGGAACGc                                                                                       >  1:348149/1‑149 (MQ=255)                                                           ccAGGCGCGCCACTGTTTACCATACACTGGCCCGAGGTCGCCGTTTTCATCGGCCCATTCGTCCCAGATGGTGACATTGTTTTCGTGTAGATAAGCAATGTTAGTGTCGCCCTGCAGAAACCACAGCAGTT‑CATGGATGATGGAACGCa                                                                                      >  2:749496/1‑149 (MQ=255)                                                            cAGGCGCGCCACTGTTTACCATACACTGGCCCGAGGTCGCCGTTTTCATCGGCCCATTCGTCCCAGATGGTGACATTGTTTTCGTGTAGATAAGCAATGTTAGTGTCGCCCTGCAGAAACCACAGCAGTT‑CATGGATGATGGAACGCAg                                                                                     >  2:488463/1‑149 (MQ=255)                                                               gcgcgcCACTGTTTACCATACACTGGCCCGAGGTCGCCGTTTTCATCGGCCCATTCGTCCCAGATGGTGACATTGTTTTCGTGTAGATAAGCAATGTTAGTGTCGCCCTGGAGAAACCACAGCAGTT‑CATGGATGATGGAACGCAGGTg                                                                                  >  1:264284/1‑149 (MQ=255)                                                                   gcCACTGTTTACCATACACTGGCCCGAGGTCGCCGTTTTCATCGGCCCATTCGTCCCAGATGGTGACATTGTTTTCGTGTAGATAAGCAATGTTAGTGTCGCCCTGCAGAAACCACAGCAGTT‑CATGGATGATGGCACGCAGGTGGCaa                                                                              >  1:624996/1‑149 (MQ=255)                                                                   gcCACTGTTTACCATACACTGGCCCGAGGTCGCCGTTTTCATCGGCCCATTCGTCCCAGATGGTGACATTGTTTTCGTGTAGATAAGCAATGTTAGTGTCGCCCTGCAGAAACCACAGCAGTT‑CATGCATGATGGACCGCAGGCGGCaa                                                                              >  2:782656/1‑149 (MQ=255)                                                                   gacACTGTTTACCATACACTGGCCAGAGGTCGCCGTTTTCATCGGCCCATTCGTCCCAGATGGTGACATTGTTTTCGTGTAGATAAGCAATGTTAGTGTCGCCCTACAGAAA‑CACAACAATTGCATGGATGATGGAACACACGAGGCaa                                                                              >  2:748219/3‑149 (MQ=255)                                                                    ccACTGTTTACCATACACTGGCCCGAGGTCGCCGTTTTCATCGGCCCATTCGTCCCAGATGGTGACATTGTTTTCGTGTAGATAAGCAATGTTAGTGTCGCCCTGCAGAAACCACAGCAGTT‑CATGGATGATGGAACGCAGGTGGCAAc                                                                             <  1:761045/149‑1 (MQ=255)                                                                     cACTGTTTACCATACACTGGCCCGAGGTCGCCGTTTTCATCGGCCCATTCGTCCCAGATGGTGACATTGTTTTCGTGTAGATAAGCAATGTTAGTGTCGCCCTGCAGAAACCACAGCAGTG‑CATGGATGATGGAACGCAGGTGGCAACg                                                                            >  1:229393/1‑149 (MQ=255)                                                                      aCTGTTTACCATACACTGGCCCGAGGTCGCCGTTTTCATCGGCCCATTCGTCCCAGATGGTGACATTGTTTTCGTGTAGATAAGCAATGTTAGTGTCGCCCTGCAGAAACCACAGCAGTT‑CATGGATGATGGAACGCAGGTGGCAACGt                                                                           >  2:168797/1‑149 (MQ=255)                                                                          tttACCATACACTGGCCCGAGGTCGCCGTTTTCATCGGCCCATTCGTCCCAGATGGTGACATTGTTTTCGTGTAGATAAGCAATGTTAGTGTCGCCCTGCAGAAACCACAGCAGTT‑CATGGATGATGGAACGCAGGTGGCAACGTTTAg                                                                       <  1:605420/149‑1 (MQ=255)                                                                                    acTGGCCCGAGGTCGCCGTTTTCATCGGCCCATTCGTCCCAGATGGTGACATTGTTTTCGTGTAGATAAGCAATGTTAGTGTCGCCCTGCAGAAACCACAGCAGTT‑CATGGATGATGGAACGCAGGTGGCAACGTTTAGTTGTCACCAg                                                             <  1:470408/149‑1 (MQ=255)                                                                                          ccGAGGTCGCCGTTTTCATCGGCCCATTCGTCCCAGATGGTGACATTGTTTTCGTGTAGATAAGCAATGTTAGTGTCGCCCTGCAGAAACCACAGCAGTT‑CATGGATGATGGAACGCAGGTGGCAACGTTTAGTTGTCACCAGCGGGaa                                                       >  2:237205/1‑149 (MQ=255)                                                                                                        ttCATCGGCCCATTCGTCCCAGATGGTGACATTGTTTTCGTGTAGATAAGCAATGTTAGTGTCGCCCTGCAGAAACCACAGCAGTT‑CATGGATGATGGAACGCAGGTGGCAACGTTTAGTTGTCACCAGCGGGAATCCATCTTGCAGGt                                         >  1:502452/1‑149 (MQ=255)                                                                                                        ttCATCGGCCCATTCGTCCCAGATGGTGACATTGTTTTCGTGTAGATAAGCAATGTTAGTGTCGCCCTGCAGAAACCACAGCAGTT‑CATGGATGATGGAACGCAGGTGGCAACGTTTAGTTGTCACCAGCGGGAATCCAACTTGCAGGt                                         >  1:561449/1‑149 (MQ=255)                                                                                                         tCATCGGCCCATTCGTCCCAGATGGTGACATTGTTTTCGTGTAGATAAGCAATGTTAGTGTCGCCCTGCAGAAACCACAGCAGTT‑CATGGATGATGGAACGCAGGTGGCAACGTTTAGTTGTCACCAGCGGGAATCCATCTTGCAGGtt                                        >  1:46734/1‑149 (MQ=255)                                                                                                           aTCGGCCCAGTCGTCCCAGATGGTGACATTGTTTTCGTGTAGATAAGCAATGTTAGTGTCGCCCTGCAGAAACCACAGCAGTT‑CATGGATGATGGAACGCAGGTGGCAACGTTTAGTTGTCACCAGCGGGAATCCATCTTGCAGGTTaa                                      <  2:120956/149‑1 (MQ=255)                                                                                                                 ccATTCGTCCCAGATGGTGACATTGTTTTCGTGTAGATAAGCAATGTTAGTGTCGCCCTGCAGAAACCACAGCAGTT‑CATGGATGATGGAACGCAGGTGGCAACGTTTAGTTGTCATCAGCGGGAATCCATCTTGCAGGTTAAAACGCa                                >  2:439723/1‑149 (MQ=255)                                                                                                                         tccAGGTGGGGACATTGTTTTCGTGTAGGTGAGCAATGTTAGTGGCGCCCTGCAGAAACCACTGCAGGT‑CATGGATGATGGAACGCAGGTGGCAACGTTTAGTTGTCACCAGCGGGAATCCATCTTGCAGGTTAAAACGCATCtgatga                        <  1:86943/148‑1 (MQ=255)                                                                                                                         cccAGATGGTGACATTGTTTTCGTGTAGGTAAGCAATGTTAGTGTCGCCCTGCAGAAACCACAGCAGTT‑CATGGATGATGGAACGCAGGTGGCAACGTTTAGTTGTCACCAGCGGGAATCCATCTTGCAGGTTAAAACGCATCtgatga                        <  1:305826/149‑1 (MQ=255)                                                                                                                                      aTTGTTTTCGTGTAGATAAGCAATGTTAGTGTCGCCCTGCAGAAACCCCAGCCGTT‑CATGGATGATGGAACGCAGGTGGCAACGTTTAGTTGTCACCAGCGGGAATCCATCTTGCAGGTTAAAACGCATCCGATGACCAAAAATGGaaa           >  1:667593/1‑149 (MQ=255)                                                                                                                                       ttGTTTTCGTGTAGATAAGCAATGTTAGTGTCGCCCTGCAGAAACCACAGCAGTT‑CATGGATGATGGAACGCAGGTGGCAACGTTTAGTTGTCACCAGCGGGAATCCATCTTGCAGGTTAAAACGCATCTGATGACCAAAAATGGAAAg          <  1:237205/149‑1 (MQ=255)                                                                                                                                               gtgtAGATAAGCAATGTTAGTGTCGCCCTGCAGAAACCACAGCAGTT‑CATGGATGATGGAACGCAGGTGGCAACGTTTAGTTGTCACCAGCGGGAATCCATCTTGCAGGTTAAAAAGCATATGATGACCAACAATGGAAAGGgttccgg  >  2:648790/1‑149 (MQ=255)                                                                                                                                                     |                                                                                                                                                TTTTCAGCTGGTTCAGTACCGTAGTGATCTGGTCAATATGACGACCATCTGGCGTTGGCCAGGCGCGCCACTGTTTACCATACACTGGCCCGAGGTCGCCGTTTTCATCGGCCCATTCGTCCCAGATGGTGACATTGTTTTCGTGTAGATAAGCAATGTTAGTGTCGCCCTGCAGAAACCACAGCAGTT‑CATGGATGATGGAACGCAGGTGGCAACGTTTAGTTGTCACCAGCGGGAATCCATCTTGCAGGTTAAAACGCATCTGATGACCAAAAATGGAAAGCGTTCCGG  >  NC_000913/2964793‑2965083

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