Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,621,600	T→C	*240Q (TAA→CAA)	deoD →	purine nucleoside phosphorylase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,621,600	0	T	C	100.0%	180.4 / NA	50	*240Q (TAA→CAA)	deoD	purine nucleoside phosphorylase
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (18/32);  total (18/32)

GAAAGCCCTGACCATCTGCACCGTATCTGACCACATCCGCACTCACGAGCAGACCACTGCCGCTGAGCGTCAGACTACCTTCAACGACATGATCAAAATCGCACTGGAATCCGTTCTGCTGGGCGATAAAGAGTAATTGTGTTTCGCTGCAAGGCGATTGCCTTGTGAAGCCGGAGCGGGAGACTGCTCCGGCTTTTTAGTATCTATTCATTTTTCTCTCCAGCTTGAATATTTTCGCTATACTTTTCAGTGAAATGTGTTAATAAAT  >  NC_000913/4621467‑4621734                                                                                                                                      |                                                                                                                                       gAAAGCCCTGACCATCTGCACCGTATCTGACCACATCCGCACTCACGAGCAGACCACTGCCGCTGAGCGTCAGACTACCTTCAACGACATGATCAAAATCGCACTGGAATCCGTTCTGCTGGGCGATAAAGAGCAATtgt                                                                                                                                  <  1:372742/140‑1 (MQ=255)           aCCATCTGCACCGTATCTGACCACATCCGCACTCACGAGCAGACCACTGCCGCTGAGCGTCAGACTACCTTCAACGACATGATCAAAATCGCACTGGAATCCGTTCTGCTGGGCGATAAAGAGCAATTGTGTTTCGCTGc                                                                                                                        <  2:832907/140‑1 (MQ=255)           aCCATCTGCACCGTATCTGACCACATCCGCACTCACGAGCAGACCACTGCCGCTGAGCGTCAGACTACCTTCAACGACATGATCAAAATCGCACTGGAATCCGTTCTGCTGGGCGATAAAGAGCAATTGTGTTTCGCTGc                                                                                                                        <  1:823987/140‑1 (MQ=255)            ccATCTGCACCGTATCTGACCACATCCGCACTCACGAGCAGACCACTGCCGCTGAGCGTCAGACTACCTTCAACGACATGATCAAAATCGCACTGGAATCCGTTCTGCTGGGCGATAAAGAGCAATTGTGTTTCGCTGCa                                                                                                                       <  1:509243/140‑1 (MQ=255)               tCTGCACCGTATCTGACCACATCCGCACTCACGAGCAGACCACTGCCGCTGAGCGTCAGACTACCTTCAACGACATGATCAAAATCGCACTGGAATCCGTTCTGCTGGGCGATAAAGAGCAATTGTGTTTCGCTGCAAgg                                                                                                                    >  1:172459/1‑140 (MQ=255)               tCTGCACCGTATCTGACCACATCCACACTCACGAGCAGACCACTGCCGCTGAGCGTCAGACTACCTTCAACGACATGATCAAAATCGCACTGGAATCCGTTCTGCTGGGCGATAAAGAGCAATTGTGTTTCGCTGCAAgg                                                                                                                    >  2:345305/1‑140 (MQ=255)                  gCACCGTATCTGACCACATCCGCACTCACGAGCAGACCACTGCCGCTGAGCGTCAGACTACCTTTAACGACATGATCAAAATCGCACTGGAATCCGTTCTGCTGGGCGATAAAGAGCAATTGTGTTTCGCTGCAAGGCGa                                                                                                                 <  2:672782/140‑1 (MQ=255)                        tATCTGACCACATCCGCACTCACGAGCAGACCACTGCCGCTGAGCGTCAGACTACCTTCAACGACATGATCAAAATCGCACTGGAATCCGTTCTGCTGGGCGATAAAGAGCAATTGTGTTTCGCTGCAAGGCGATTGCCt                                                                                                           >  2:64076/1‑140 (MQ=255)                          tCTGACCACATCCGCACTCACGAGCAGACCACTGCCGCTGAGCGTCAGACTACCTTCAACGACATGATCAAAATCGCACTGGAATCCGTTCTGCTGGGCGATAAAGAGCAATTGTGTTTCGCTGCAAGGCGATTGCCTtg                                                                                                         >  1:111719/1‑140 (MQ=255)                                  caTCCGCACTCACGAGCAGACCACTGCCGCTGAGCGTCAGACTACCTTCAACGACATGATCAAAATCGCACTGGAATCCGTTCTGCTGGGCGATAAAGAGCAATTGTGTTTCGCTGCAAGGCg                                                                                                                  >  1:318320/1‑123 (MQ=255)                                  caTCCGCACTCACGAGCAGACCACTGCCGCTGAGCGTCAGACTACCTTCAACGACATGATCAAAATCGCACTGGAATCCGTTCTGCTGGGCGATAAAGAGCAATTGTGTTTCGCTGCAAGGCg                                                                                                                  <  2:318320/123‑1 (MQ=255)                                  caTCCGCACTCACGAGCAGACCACTGCCGCTGAGCGTCAGACTACCTTCAACGACATGATCAAAATCGCACTGGAATCCGTTCTGCTGGGCGATAAAGAGCAATTGTGTTTCGCTGCAAGGCGATTGCCTTGTGAAGCCg                                                                                                 <  2:785188/140‑1 (MQ=255)                                  caTCCGCACTCACGAGCAGACCACTGCCGCTGAGCGTCAGACTACCTTCAACGACATGATCAAAATCGCACTGGAATCCGTTCTGCTGGGCGATAAAGAGCAATTGTGTTTCGCTGCAAGGCGATTGCCTTGTGAAGCCg                                                                                                 >  2:148029/1‑140 (MQ=255)                                  caTCCGCACTCACGAGCAGACCACTGCCGCTGAGCGTCAGACTACCTTCAACGACATGATCAAAATCGCACTGGAATCCGTTCTGCTGGGCGATAAAGAGCAATTGTGTTTCGCTGCAAGGCGATTGCCTTGTGAAGCCg                                                                                                 >  2:1081817/1‑140 (MQ=255)                                   aTCCGCACTCACGAGCAGACCACTGCCGCTGAGCGTCAGACTACCTTCAACGACATGATCAAAATCGCACTGGAATCCGTTCTGCTGGGCGATAAAGAGCAATTGTGTTTCGCTGCAAGGCGATTGCCTTGTGAAGCCgg                                                                                                >  1:716329/1‑140 (MQ=255)                                       gCACTCACGAGCAGACCACTGCCGCTGAGCGTCAGACTACCTTCAACGACATGATCAAAATCGCACTGGAATCCGTTCTGCTGGGCGATAAAGAGCAATTGTGTTTCGCTGCAAGGCGATTGCCTTGTGAAGCCGGAGCg                                                                                            <  1:1076356/140‑1 (MQ=255)                                       gCACTCACGAGCAGACCACTGCCGCTGAGCGTCAGACTACCTTCAACGACATGATCAAAATCGCACTGGAATCCGTTCTGCTGGGCGATAAAGAGCAATTGTGTTTCGCTGCAAGGCGATTGCCTTGTGAAGCCGGAGCg                                                                                            <  1:148029/140‑1 (MQ=255)                                            cACGAGCAGACCACTGCCGCTGAGCGTCAGACTACCTTCAACGACATGATCAAAATCGCACTGGAATCCGTTCTGCTGGGCGATAAAGAGCAATTGTGTTTCGCTGCAAGGCGATTGCCTTGTGAAGCCGGAGCGGgaga                                                                                       <  1:467637/140‑1 (MQ=255)                                               gAGCAGACCACTGCCGCTGAGCGTCAGACTACCTTCAACGACATGATCAAAATCGCACTGGAATCCGTTCTGCTGGGCGATAAAGAGCAATTGTGTTTCGCTGCAAGGCGATTGCCTTGTGAAGCCGGAGCGGGAGACTg                                                                                    >  2:33063/1‑140 (MQ=255)                                                    gACCACTGCCGCTGAGCGTCAGACTACCTTCAACGACATGATCAAAATCGCACTGGAATCCGTTCTGCTGGGCGATAAAGAGCAATTGTGTTTCGCTGCAAGGCGATTGCCTTGTGAAGCCGGAGCGGGAGACTGCTCCg                                                                               <  1:1081817/140‑1 (MQ=255)                                                     aCCACTGCCGCTGAGCGTCAGACTACCTTCAACGACATGATCAAAATCGCACTGGAATCCGTTCTGCTGGGCGATAAAGAGCAATTGTGTTTCGCTGCAAGGCGATTGCCTTGTGAAGCCGGAGCGGGAGACTGCTCCgg                                                                              <  1:632746/140‑1 (MQ=255)                                                     aCCACTGCCGCTGAGCGTCAGACTACCTTCAACGACATGATCAAAATCGCACTGGAATCCGTTCTGCTGGGCGATAAAGAGCAATTGTGTTTCGCTGCAAGGCGATTGCCTTGTGAAGCCGGAGCGGGAGACTGCTCCgg                                                                              >  1:131866/1‑140 (MQ=255)                                                     aCCACTGCCGCTGAGCGTCAGACTACCTTCAACGACATGATCAAAATCGCACTGGAATCCGTTCTGCTGGGCGATAAAGAGCAATTGTGTTTCGCTGCAAGGCGATTGCCTTGTGAAGCCGGAGCGGGAGACTGCTCCgg                                                                              >  1:988843/1‑140 (MQ=255)                                                     aCCACTGCCGCTGAGCGTCAGACTACCTTCAACGACATGATCAAAATCGCACTGGAATCCGTTCTGCTGGGCGATAAAGAGCAATTGTGTTTCGCTGCAAGGCGATTGCCTTGTGAAGCCGGAGCGGGAGACTGCTCCgg                                                                              <  2:172459/140‑1 (MQ=255)                                                      ccACTGCCGCTGAGCGTCAGACTACCTTCAACGACATGATCAAAATCGCACTGGAATCCGTTCTGCTGGGCGATAAAGAGCAATTGTGTTTCGCTGCAAGGCGATTGCCTTGTGAAGCCGGAGCGGGAGACTGCTCCGGc                                                                             <  1:308165/140‑1 (MQ=255)                                                       cACTGCCGCTGAGCGTCAGACTACCTTCAACGACATGATCAAAATCGCACTGGAATCCGTTCTGCTGGGCGATAAAGAGCAATTGTGTTTCGCTGCAAGGCGATTGCCTTGTGAAGCCGGAGCGGGAGACTGCTCCGGCt                                                                            <  1:521667/140‑1 (MQ=255)                                                       cACTGCCGCTGAGCGTCAGACTACCTTCAACGACATGATCAAAATCGCACTGGAATCCGTTCTGCTGGGCGATAAAGAGCAATTGTGTTTCGCTGCAAGGCGATTGCCTTGTGAAGCCGGAGCGGGAGACTGCTCCGGCt                                                                            <  2:716329/140‑1 (MQ=255)                                                           gCCGCTGAGCGTCAGACTACCTTCAACGACATGATCAAAATCGCACTGGAATCCGTTCTGCTGGGCGATAAAGAGCAATTGTGTTTCGCTGCAAGGCGATTGCCTTGTGAAGCCGGAGCGGGAGACTGCTCCGGCttttt                                                                        <  1:183923/140‑1 (MQ=255)                                                            ccGCTGAGCGTCAGACTACCTTCAACGACATGATCAAAATCGCACTGGAATCCGTTCTGCTGGGCGATAAAGAGCAATTGTGTTTCGCTGCAAGGCGATTGCCTTGTGAAGCCGGAGCGGGAGACTGCTCCGGCTTTTTa                                                                       <  2:111719/140‑1 (MQ=255)                                                             cGCTGAGCGTCAGACTACCTTCAACGACATGATCAAAATCGCACTGGAATCCGTTCTGCTGGGCGATAAAGAGCAATTGTGTTTCGCTGCAAGGCGATTGCCTTGTGAAGCCGGAGCGGGAGACTGCTc                                                                                 <  1:997162/129‑1 (MQ=255)                                                             cGCTGAGCGTCAGACTACCTTCAACGACATGATCAAAATCGCACTGGAATCCGTTCTGCTGGGCGATAAAGAGCAATTGTGTTTCGCTGCAAGGCGATTGCCTTGTGAAGCCGGAGCGGGAGACTGCTc                                                                                 >  2:997162/1‑129 (MQ=255)                                                                tGAGCGTCAGACTACCTTCAACGACATGATCAAAATCGCACTGGAATCCGTTCTGCTGGGCGATAAAGAGCAATTGTGTTTCGCTGCAAGGCGATTGCCTTGTGAAGCCGGAGCGGGAGACTGCTCCGGCTTTTTAGTAt                                                                   <  1:294001/140‑1 (MQ=255)                                                                  aGCGTCAGACTACCTTCAACGACATGATCAAAATCGCACTGGAATCCGTTCTGCTGGGCGATAAAGAGCAATTGTGTTTCGCTGCAAGGCGATTGCCTTGTGAAGCCGGAGCGGGAGACTGCTCCGGCTTTTTAGTATCt                                                                 >  2:346411/1‑140 (MQ=255)                                                                  aGCGTCAGACTACCTTCAACGACATGATCAAAATCGCACTGGAATCCGTTCTGCTGGGCGATAAAGAGCAATTGTGTTTCGCTGCAAGGCGATTGCCTTGTGAAGCCGGAGCGGGAGACTGCTCCGGCTTTTTAGTATCt                                                                 <  2:208051/140‑1 (MQ=255)                                                                    cGTCAGACTACCTTCAACGACATGATCAAAATCGCACTGGAATCCGTTCTGCTGGGCGATAAAGAGCAATTGTGTTTCGCTGCAAGGCGATTGCCTTGTGAAGCCGGAGCGGGAGACTGCTCCGGCTTTTTAGTATCTAt                                                               >  1:406814/1‑140 (MQ=255)                                                                           cTACCTTCAACGACATGATCAAAATCGCACTGGAATCCGTTCTGCTGGGCGATAAAGAGCAATTGTGTTTCGCTGCAAGGCGATTGCCTTGTGAAGCCGGAGCGGGAGACTGCTCCGGCTTTTTAGTATCTATTCAtttt                                                        <  1:300766/140‑1 (MQ=255)                                                                                 tCAACGACATGATCAAAATCGCACTGGAATCCGTTCTGCTGGGCGATAAAGAGCAATTGTGTTTCGCTGCAAGGCGATTGCCTTGTGAAGCCGGAGCGGGAGACTGCTCCGGCTTTTTAGTATCTATTCATTTTtctctc                                                  <  1:632110/140‑1 (MQ=255)                                                                                  cAACGACATGATCAAAATCGCACTGGAATCCGTTCTGCTGGGCGATAAAGAGCAATTGTGTTTCGCTGCAAGGCGATTGCCTTGTGAAGCCGGAGCGGGAGACTGCTCCGGCTTTTTAGTATCTATTCATTTTTCTCTcc                                                 <  2:988843/140‑1 (MQ=255)                                                                                     cGACATGATCAAAATCGCACTGGAATCCGTTCTGCTGGGCGATAAAGAGCAATTGTGTTTCGCTGCAAGGCGATTGCCTTGTGAAGCCGGAGCGGGAGACTGCTCCGGCTTTTTAGTATCTATTCATTTTTCTCTCCAGc                                              <  1:33063/140‑1 (MQ=255)                                                                                        cATGATCAAAATCGCACTGGAATCCGTTCTGCTGGGCGATAAAGAGCAATTGTGTTTCGCTGCAAGGCGATTGCCTTGTGAAGCCGGAGCGGGAGACTGCTCCGGC‑TTTTAGTATCTATTCATTTTTCTCTCCAGCTTGa                                          >  2:298937/1‑140 (MQ=255)                                                                                             tCAAAATCGCACTGGAATCCGTTCTGCTGGGCGATAAAGAGCAATTGTGTTTCGCTGCAAGGCGATTGCCTTGTGAAGCCGGAGCGGGAGACTGCTCCGGCTTTTTAGTATCTATTCATTTTTCTCTCCAGCTTGAatat                                      <  1:671434/140‑1 (MQ=255)                                                                                             tCAAAATCGCACTGGAATCCGTTCTGCTGGGCGATAAAGAGCAATTGTGTTTCGCTGCAAGGCGATTGCCTTGTGAAGCCGGAGCGGGAGACTGCTCCGGCTTTTTAGTATCTATTCATTTTTCTCTCCAGCTTGAatat                                      <  1:64076/140‑1 (MQ=255)                                                                                                    cGCACTGGAATCCGTTCTGCTGGGCGATGAAGAGCAATTGTGTTTCGCTGCAAGGCGATTGCCTTGTGAAGCCGGAGCGGGAGACTGCTCCGGCTTTTTAGTATCTATTCATTTTTCTCTCCAGCTTGAATATTTTCGCt                               >  1:318776/1‑140 (MQ=255)                                                                                                                  ttCTGCTGGGCGATAAAGAGCAATTGTGTTTCGCTGCAAGGCGATTGCCTTGTGAAGCCGGAGCGGGAGACTGCTCCGGC‑TTTTAGTATCTATTCATTTTTCTCTCCAGCTTGAATATTTTCGCTATACTTTTCAGTGaa                <  1:298937/140‑1 (MQ=255)                                                                                                                     tgctgGGCGATAAAGAGCAATTGTGTTTCGCTGCAAGGCGATTGCCTTGTGAAGCCGGAGCGGGAGACTGCTCCGGCTTTTTAGTATCTATTCATTTTTCTCTCCAGCTTGAATATTTTCGCTATACTTTTCAGTGAAAt              >  1:482145/1‑140 (MQ=255)                                                                                                                          ggCGATAAAGAGCAATTGTGTTTCGCTGCAAGGCGATTGCCTTGTGAAGCCGGAGCGGGAGACTGCTCCGGCTTTTTAGTATCTATTCATTTTTCTCTCCAGCTTGAATATTTTCGCTATACTTTTCAGTGAAATGTGtt         <  2:757904/140‑1 (MQ=255)                                                                                                                               tAAAGAGCAATTGTGTTTCGCTGCAAGGCGATTGCCTTGTGAAGCCGGAGCGGGAGACTGCTCCGGCTTTTTAGTATCTATTCATTTTTCTCTCCAGCTTGAATATTTTCGCTATACTTTTCAGTGAAATGTGTtaataa    <  2:359666/140‑1 (MQ=255)                                                                                                                                 aaGAGCAATTGTGTTTCGCTGCAAGGCGATTGCCTTGTGAAGCCGGAGCGGGAGACTGCTCCGGCTTTTTAGTATCTATTCATTTTTCTCTCCAGCTTGAATATTTTCGCTATACTTTTCAGTGAAATGTGTTAATAAAt  <  2:522859/140‑1 (MQ=255)                                                                                                                                 aaGAGCAATTGTGTTTCGCTGCAAGGCGATTGCCTTGTGAAGCCGGAGCGGGAGACTGCTCCGGCTTTTTAGTATCTATTCATTTTTCTCTCCAGCTTGAATATTTTCGCTATACTTTTCAGTGAAATGTGTTAATAAAt  <  2:131866/140‑1 (MQ=255)                                                                                                                                 aaGAGCAATTGTGTTTCGCTGCAAGGCGATTGCCTTGTGAAGCCGGAGCGGGAGACTGCTCCGGCTTTTTAGTATCTATTCATTTTTCTCTCCAGCTTGAATATTTTCGCTATACTTTTCAGTGAAATGTGTTAATAAAt  >  1:522859/1‑140 (MQ=255)                                                                                                                                      |                                                                                                                                       GAAAGCCCTGACCATCTGCACCGTATCTGACCACATCCGCACTCACGAGCAGACCACTGCCGCTGAGCGTCAGACTACCTTCAACGACATGATCAAAATCGCACTGGAATCCGTTCTGCTGGGCGATAAAGAGTAATTGTGTTTCGCTGCAAGGCGATTGCCTTGTGAAGCCGGAGCGGGAGACTGCTCCGGCTTTTTAGTATCTATTCATTTTTCTCTCCAGCTTGAATATTTTCGCTATACTTTTCAGTGAAATGTGTTAATAAAT  >  NC_000913/4621467‑4621734

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