Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,005,368	C→T	Q201* (CAA→TAA)	pyrD →	dihydroorotate dehydrogenase, type 2

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,005,368	0	C	T	100.0%	201.1 / NA	57	Q201* (CAA→TAA)	pyrD	dihydroorotate dehydrogenase, type 2
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (30/27);  total (30/27)

TCTGATTTGTATGGAAAAAATCTATGCCTATGCGGGATATATCGCCATCAATATTTCATCGCCGAATACCCCAGGATTACGCACGCTGCAATATGGTGAAGCGCTGGATGATCTCTTAACCGCGATTAAAAATAAGCAAAATGATTTGCAAGCGATGCACCATAAATATGTGCCGATCGCAGTGAAGATCGCGCCGGATCTTTCTGAAGAAGAATTGATCCAGGTTGCCGATAGTTTAGTTCGCCATAATATTGATGGCGTTATTGCAACCA  >  NC_000913/1005232‑1005503                                                                                                                                         |                                                                                                                                        tCTGATTTGTATGGAAAAAATCTATGCCTATGCGGGATATATCGCCATCAATATTTCATCGCCGAATACCCCAGGATTACGCACGCTGCAATATGGTGAAGCGCTGGATGATCTCTTAACCGCGATTAAAAATAAGTaaa                                                                                                                                      >  2:300251/1‑140 (MQ=255)          tATGGAAAAAATCTATGCCTATGCGGGATATATCGCCATCAATATTTCATCGCCGAATACCCCAGGATTACGCACGCTGCAATATGGTGAAGCGCTGGATGATCTCTTAACCGCGATTAAAAATAAGTAAAATGATTTGc                                                                                                                             >  1:637911/1‑140 (MQ=255)           aTGGAAAAAATCTATGCCTATGCGGGATATATCGCCATCAATATTTCATCGCCGAATACCCCAGGATTACGCACGCTGCAATATGGTGAAGCGCTGGATGATCTCTTAACCGCGATTAAAAATAAGTAAAATGATTTGc                                                                                                                             <  1:101003/139‑1 (MQ=255)           aTGGAAAAAATCTATGCCTATGCGGGATATATCGCCATCAATATTTCATCGCCGAATACCCCAGGATTACGCACGCTGCAATATGGTGAAGCGCTGGATGATCTCTTAACCGCGATTAAAAATAAGTAAAATGATTTGc                                                                                                                             >  2:101003/1‑139 (MQ=255)           aTGGAAAAAATCTATGCCTATGCGGGATATATCGCCATCAATATTTCATCGCCGAATACCCCAGGATTACGCACGCTGCAATATGGTGAAGCGCTGGATGATCTCTTAACCGCGATTAAAAATAAGTAAAATGATTTGCa                                                                                                                            >  2:366279/1‑140 (MQ=255)           aTGGAAAAAATCTATGCCTATGCGGGATATATCGCCATCAATATTTCATCGCCGAATACCCCAGGATTACGCACGCTGCAATATGGTGAAGCGCTGGATGATCTCTTAACCGCGATTAAAAATAAGTAAAATGATTTGCa                                                                                                                            <  1:851519/140‑1 (MQ=255)             ggAAAAAATCTATGCCTATGCGGGATATATCGCCATCAATATTTCATCGCCGAATACCCCAGGATTACGCACGCTGCAATATGGTGAAGCGCTGGATGATCTCTTAACCGCGATTAAAAATAAGTAAAATGATTTGCAAg                                                                                                                          <  1:300251/140‑1 (MQ=255)                aaaaaTCTATGCCTATGCGGGATATATCGCCATCAATATTTCATCGCCGAATACCCCAGGATTACGCACGCTGCAATATGGTGAAGCGCTGGATGATCTCTTAACCGCGATTAAAAATAAGTAAAATGATTTGCAAGCGa                                                                                                                       <  1:943630/140‑1 (MQ=255)                  aaaTCTATGCCTATGCGGGATATATCGCCATCAATATTTCATCGCCGAATACCCCAGGATTACGCACGCTGCAATATGGTGAAGCGCTGGATGATCTCTTAACCGCGATTAAAAATAAGTAAAATGATTTGCAAGCGATg                                                                                                                     <  2:953221/140‑1 (MQ=255)                  aaaTCTATGCCTATGCGGGATATATCGCCATCAATATTTCATCGCCGAATACCCCAGGATTACGCACGCTGCAATATGGTGAAGCGCTGGATGATCTCTTAACCGCGATTAAAAATAAGTAAAATGATTTGCAAGCGATg                                                                                                                     <  1:687111/140‑1 (MQ=255)                        atgcctatgcGGGATATATCGCCATCAATATTTCATCGCCGAATACCCCAGGATTACGCACGCTGCAATATGGTGAAGCGCTGGATGATCTCTTAACCGCGATTAAAAATAAGTAAAATGATTTGCAAGCGATGCACCAt                                                                                                               >  2:812357/1‑140 (MQ=255)                                      tatatCGCCATCAATATTTCATCGCCGAATACCCCAGGATTACGCACGCTGCAATATGGTGAAGCGCTGGATGATCTCTTAACCGCGATTAAAAATAAGTAAAATGATTTGCAAGCGATGCACCATAAATATGTGCCGAt                                                                                                 >  1:652335/1‑140 (MQ=255)                                       atatatccATCAATATTTCATCGCCGAATACCCCAGGATTACGCACGCTGCAATATGGTGAAGCGCTGGATGATCTCTTAACCGCGATTAAAAATAAGTAAAATGATTTGCAAGCGATGCACCATAAATATGTGCCGATc                                                                                                <  1:397840/134‑1 (MQ=255)                                       atatCGCCATCAATATTTCATCGCCGAATACCCCAGGATTACGCACGCTGCAATATGGTGAAGCGCTGGATGATCTCTTAACCGCGATTAAAAATAAGTAAAATGATTTGCAAGCGATGCACCATAAATATGTGCCGATc                                                                                                <  1:812357/140‑1 (MQ=255)                                         atCGCCATCAATATTTCATCGCCGAATACCCCAGGATTACGCACGCTGCAATATGGTGAAGCGCTGGATGATCTCTTAACCGCGATTAAAAATAAGTAAAATGATTTGCAAGCGATGCACCATAAATATGTGCCGATCGc                                                                                              >  2:323626/1‑140 (MQ=255)                                         atCGCCATCAATATTTCATCGCCGAATACCCCAGGATTACGCACGCTGCAATATGGTGAAGCGCTGGATGATCTCTTAACCGCGATTAAAAATAAGTAAAATGATTTGCAAGCGATGCACCATAAATATGTGCCGATCGc                                                                                              >  2:405666/1‑140 (MQ=255)                                                tCAATATTTCATCGCCGAATACCCCAGGATTACGCACGCTGCAATATGGTGAAGCGCTGGATGATCTCTTAACCGCGATTAAAAATAAGTAAAATGATTTGCAAGCGATGCACCATAAATATGTGCCGATCGCAGTGAAg                                                                                       >  2:315276/1‑140 (MQ=255)                                                  aaTATTTCATCGCCGAATACCCCAGGATTACGCACGCTGCAATATGGTGAAGCGCTGGATGATCTCTTAACCGCGATTAAAAATAAGTAAAATGATTTGCAAGCGATGCACCATAAATATGTGCCGATCGCAGTGAAGAt                                                                                     >  1:434056/1‑140 (MQ=255)                                                  aaTATTTCATCGCCGAATACCCCAGGATTACGCACGCTGCAATATGGTGAAGCGCTGGATGATCTCTTAACCGCGATTAAAAATAAGTAAAATGATTTGCAAGCGATGCACCATAAATATGTGCCGATCGCAGTGAAGAt                                                                                     <  1:315276/140‑1 (MQ=255)                                                  aaTATTTCATCGCCGAATACCCCAAGATTACGCACGCTGCAATATGGTGAAGCGCTGGATGATCTCTTAACCGCGATTAAAAATAAGTAAAATGATTTGCAAGCGATGCACCATAAATATGTGCCGATCGCAGTGAAGAt                                                                                     >  1:492690/1‑140 (MQ=255)                                                      tttCATCGCCGAATACCCCAGGATTACGCACGCTGCAATATGGTGAAGCGCTGGATGATCTCTTAACCGCGATTAAAAATAAGTAAAATGATTTGCAAGCGATGCACCATAAATATGTGCCGATCGCAGTGAAGATcgcg                                                                                 <  1:818359/140‑1 (MQ=255)                                                              ccGAATACCCCAGGATTACGCACGCTGCAATATGGTGAAGCGCTGGATGATCTCTTAACCGCGATTAAAAATAAGTAAAATGATTTGCAAGCGATGCACCATAAATATGTGCCGATCGCAGTGAAGATcgcg                                                                                 <  1:423025/132‑1 (MQ=255)                                                              ccGAATACCCCAGGATTACGCACGCTGCAATATGGTGAAGCGCTGGATGATCTCTTAACCGCGATTAAAAATAAGTAAAATGATTTGCAAGCGATGCACCATAAATATGTGCCGATCGCAGTGAAGATcgcg                                                                                 >  2:423025/1‑132 (MQ=255)                                                              ccGAATACCCCAGGATTACGCACGCTGCAATATGGTGAAGCGCTGGATGATCTCTTAACCGCGATTAAAAATAAGTAAAATGATTTGCAAGCGATGCACCATAAATATGTGCCGATCGCAGTGAAGATCGCGCCGGATCt                                                                         >  2:70360/1‑140 (MQ=255)                                                              ccGAATACCCCAGGATTACGCACGCTGCAATATGGTGAAGCGCTGGATGATCTCTTAACCGCGATTAAAAATAAGTAAAATGATTTGCAAGCGATGCACCATAAATATGTGCCGATCGCAGTGAAGATCGCGCCGGATCt                                                                         >  2:596987/1‑140 (MQ=255)                                                              ccGAATACCCCAGGATTACGCACGCTGCAATATGGTGAAGCGCTGGATGATCTCTTAACCGCGATTAAAAATAAGTAAAATGATTTGCAAGCGATGCACCATAAATATGTGCCGATCGCAGTGAAGATCGCGCCGGATCt                                                                         >  2:381981/1‑140 (MQ=255)                                                                           gATTACGCACGCTGCAATATGGTGAAGCGCTGGATGATCTCTTAACCGCGATTAAAAATAAGTAAAATGATTTGCAAGCGATGCACCATAAATATGTGCCGATCGCAGTGAAGATCGCGCCGGATCTTTCTgaagaagaa                                                            <  2:492690/140‑1 (MQ=255)                                                                           gATTACGCACGCTGCAATATGGTGAAGCGCTGGATGATCTCTTAACCGCGATTAAAAATAAGTAAAATGATTTGCAAGCGATGCACCATAAATATGTGCCGATCGCAGTGAAGATCGCGCCGGATCTTTCTgaagaagaa                                                            >  2:379364/1‑140 (MQ=255)                                                                           gATTACGCACGCTGCAATATGGTGAAGCGCTGGATGATCTCTTAACCGCGATTAAAAATAAGTAAAATGATTTGCAAGCGATGCACCATAAATATGTGCCGATCGCAGTGAAGATCGCGCCGGATCTTTCTgaagaagaa                                                            <  2:690813/140‑1 (MQ=255)                                                                             ttACGCACGCTGCAATATGGTGAAGCGCTGGATGATCTCTTAACCGCGATTAAAAATAAGTAAAATGATTTGCAAGCGATGCACCATAAATATGTGCCGATCGCAGTGAAGATCGCGCCGGATCTTTCTGAAGAAGAAtt                                                          >  1:113671/1‑140 (MQ=255)                                                                              tACGCACGCTGCAATATGGTGAAGCGCTGGATGATCTCTTAACCGCGATTAAAAATAAGTAAAATGATTTGCAAGCGATGCACCATAAATATGTGCCGATCGCAGTGAAGATCGCGCCGGATCTTTCTGAAGAAGAATTg                                                         >  1:735240/1‑140 (MQ=255)                                                                               acgcacgcTGCAATATGGTGAAGCGCTGGATGATCTCTTAACCGCGATTAAAAATAAGTAAAATGATTTGCAAGCGATGCACCATAAATATGTGCCGATCGCAGTGAAGATCGCGCCGGATCTTTCTGAAGAAGAATTGa                                                        >  2:1018129/1‑140 (MQ=255)                                                                                cgcacgcTGCAATATGGTGAAGCGCTGGATGATCTCTTAACCGCGATTAAAAATAAGTAAAATGATTTGCAAGCGATGCACCATAAATATGTGCCGATCGCAGTGAAGATCGCGCCGGATCTTTCTgaagaagaa                                                            <  2:961785/135‑1 (MQ=255)                                                                                cgcacgcTGCAATATGGTGAAGCGCTGGATGATCTCTTAACCGCGATTAAAAATAAGTAAAATGATTTGCAAGCGATGCACCATAAATATGTGCCGATCGCAGTGAAGATCGCGCCGGATCTTTCTgaagaagaa                                                            >  1:961785/1‑135 (MQ=255)                                                                                cgcacgcTGCAATATGGTGAAGCGCTGGATGATCTCTTAACCGCGATTAAAAATAAGTAAAATGATTTGCAAGCGATGCACCATAAATATGTGCCGATCGCAGTGAAGATCGCGCCGGATCTTTCTGAAGAAGAATTGAt                                                       >  2:253287/1‑140 (MQ=255)                                                                                           atatGGTGAAGCGCTGGATGATCTCTTAACCGCGATTAAAAATAAGTAAAATGATTTGCAAGCGATGCACCATAAATATGTGCCGATCGCAGTGAAGATCGCGCCGGATCTTTCTGAAGAAGAATTGATCCAGGTTGCCg                                            <  1:405666/140‑1 (MQ=255)                                                                                           atatGGTGAAGCGCTGGATGATCTCTTAACCGCGATTAAAAATAAGTAAAATGATTTGCAAGCGATGCACCACAAATATGTGCCGATCGCAGTGAAGATCGCGCCGGATCTTTCTGAAGAAGAATTGATCCAGGTTGCCg                                            <  1:366279/140‑1 (MQ=255)                                                                                               ggTGAAGCGCTGGATGATCTCTTAACCGCGATTAAAAATAAGTAAAATGATTTGCAAGCGATGCACCATAAATATGTGCCGATCGCAGTGAAGATCGCGCCGGATCTTTCTGAAGAAGAATTGATCCAGGTTGCCGAtag                                        <  1:253287/140‑1 (MQ=255)                                                                                                 tGAAGCGCTGGATGATCTCTTAACCGCGATTAAAAATAAGTAAAATGATTTGCAAGCGATGCACCATAAATATGTGCCGATCGCAGTGAAGATCGCGCCGGATCTTTCTGAAGAAGAATTGATCCAGGTTGCCGAtagtt                                      >  2:411035/1‑140 (MQ=255)                                                                                                 tGAAGCGCTGGATGATCTCTTAACCGCGATTAAAAATAAGTAAAATGATTTGCAAGCGATGCACCATAAATATGTGCCGATCGCAGTGAAGATCGCGCCGGATCTTTCTGAAGAAGAATTGATCCAGGTTGCCGAtagtt                                      >  2:346910/1‑140 (MQ=255)                                                                                                   aaGCGCTGGATGATCTCTTAACCGCGATTAAAAATAAGTAAAATGATTTGCAAGCGATGCACCATAAATATGTGCCGATCGCAGTGAAGATCGCGCCGGATCTTTCTGAAGAAGAATTGATCCAGGTTGCCGAtagttta                                    >  1:200542/1‑140 (MQ=255)                                                                                                      cgcTGGATGATCTCTTAACCGCGATTAAAAATAAGTAAAATGATTTGCAAGCGATGCACCATAAATATGTGCCGATCGCAGTGAAGATCGCGCCGGATCTTTCTGAAGAAGAATTGATCCAGGTTGCCGAtagtttagtt                                 <  1:381981/140‑1 (MQ=255)                                                                                                                   ctTAACCGCGATTAAAAATAAGTAAAATGATTTGCAAGCGATGCACCATAAATATGTGCCGATCGCAGTGAAGATCGCGCCGGATCTTTCTGAAGAAGAATTGATCCAGGTTGCCGATAGTTTAGTTCGCCATAATATTg                    >  2:688625/1‑140 (MQ=255)                                                                                                                   ctTAACCGCGATTAAAAATAAGTAAAATGATTTGCAAGCGATGCACCATAAATATGTGCCGATCGCAGTGAAGATCGCGCCGGATCTTTCTGAAGAAGAATTGATCCAGGTTGCCGATAGTTTAGTTCGCCATAATATTg                    >  2:1058432/1‑140 (MQ=255)                                                                                                                           cgATTAAAAATAAGTAAAATGATTTGCAAGCGATGCACCATAAATATGTGCCGATCGCAGTGAAGATCGCGCCGGATCTTTCTGAAGAAGAATTGATCCAGGTTGCCGATAGTTTAGTTCGCCATAATATTGATGGCGtt            >  2:343460/1‑140 (MQ=255)                                                                                                                           cgATTAAAAATAAGTAAAATGATTTGCAAGCGATGCACCATAAATATGTGCCGATCGCAGTGAAGATCGCGCCGGATCTTTCTGAAGAAGAATTGATCCAGGTTGCCGATAGTTTAGTTCGCCATAATATTGATGGCGtt            >  2:411909/1‑140 (MQ=255)                                                                                                                            gATTAAAAATAAGTAAAATGATTTGCAAGCGATGCACCATAAATATGTGCCGATCGCAGTGAAGATCGCGCCGGATCTTTCTGAAGAAGAATTGATCCAGGTTGCCGATAGTTTAGTTCGCCATAATATTGATGGCg              >  2:967798/1‑137 (MQ=255)                                                                                                                            gATTAAAAATAAGTAAAATGATTTGCAAGCGATGCACCATAAATATGTGCCGATCGCAGTGAAGATCGCGCCGGATCTTTCTGAAGAAGAATTGATCCAGGTTGCCGATAGTTTAGTTCGCCATAATATTGATGGCg              <  1:967798/137‑1 (MQ=255)                                                                                                                              ttAAAAATAAGTAAAATGATTTGCAAGCGATGCACCATAAATATGTGCCGATCGCAGTGAAGATCGCGCCGGATCTTTCTGAAGAAGAATTGATCCAGGTTGCCGATAGTTTAGTTCGCCATAATATTGATGGCGTTAtt         >  1:553228/1‑140 (MQ=255)                                                                                                                               tAAAAATAAGTAAAATGATTTGCAAGCGATGCACCATAAATATGTGCCGATCGCAGTGAAGATCGCGCCGGATCTTTCTGAAGAAGAATTGATCCAGGTTGCCGATAGTTTAGTTCGCCATAATATTGATGGCGTTATTg        <  2:637911/140‑1 (MQ=255)                                                                                                                               tAAAAATAAGTAAAATGATTTGCAAGCGATGCACCATAAATATGTGCCGATCGCAGTGAAGATCGCGCCGGATCTTTCTGAAGAAGAATTGATCCAGGTTGCCGATAGTTTAGTTCGCCATAATATTGATGGCGTTATTg        <  1:926697/140‑1 (MQ=255)                                                                                                                               tAAAAATAAGTAAAATGATTTGCAAGCGATGCACCATAAATATGTGCCGATCGCAGTGAAGATCGCGCCGGATCTTTCTGAAGAAGAATTGATCCAGGTTGCCGATAGTTTAGTTCGCCATAATATTGATGGCGTTATTg        <  1:1018129/140‑1 (MQ=255)                                                                                                                                 aaaaTAAGTAAAATGATTTGCAAGCGATGCACCATAAATATGTGCCGATCGCAGTGAAGATCGCGCCGGATCTTTCTGAAGAAGAATTGATCCAGGTTGCCGATAGTTTAGTTCGCCATAATATTGATGGCGTTATTGCa      <  1:346910/140‑1 (MQ=255)                                                                                                                                   aaTAAGTAAAATGATTTGCAAGCGATGCACCATAAATATGTGCCGATCGCAGTGAAGATCGCGCCGGATCTTTCTGAAGAAGAATTGATCCAGGTTGCCGATAGTTTAGTTCGCCATAATATTGATGGCGTTATTGCAAc    <  2:289068/140‑1 (MQ=255)                                                                                                                                    aTAAGTAAAATGATTTGCAAGCGATGCACCATAAATATGTGCCGATCGCAGTGAAGATCGCGCCGGATCTTTCTGAAGAAGAATTGATCCAGGTTGCCGATAGTTTAGTTCGCCATAATATTGATGGCGTTATTGCAAcc   <  2:113671/140‑1 (MQ=255)                                                                                                                                     tAAGTAAAATGATTTGCAAGCGATGCACCATAAATATGTGCCGATCGCAGTGAAGATCGCGCCGGATCTTTCTGAAGAAGAATTGATCCAGGTTGCCGATAGTTTAGTTCGCCATAATATTGATGGCGTTATTGCAACCa  <  1:688625/140‑1 (MQ=255)                                                                                                                                     tAAGTAAAATGATTTGCAAGCGATGCACCATAAATATGTGCCGATCGCAGTGAAGATCGCGCCGGATCTTTCTGAAGAAGAATTGATCCAGGTTGCCGATAGTTTAGTTCGCCATAATATTGATGGCGTTATTGCAACCa  <  1:1058432/140‑1 (MQ=255)                                                                                                                                         |                                                                                                                                        TCTGATTTGTATGGAAAAAATCTATGCCTATGCGGGATATATCGCCATCAATATTTCATCGCCGAATACCCCAGGATTACGCACGCTGCAATATGGTGAAGCGCTGGATGATCTCTTAACCGCGATTAAAAATAAGCAAAATGATTTGCAAGCGATGCACCATAAATATGTGCCGATCGCAGTGAAGATCGCGCCGGATCTTTCTGAAGAAGAATTGATCCAGGTTGCCGATAGTTTAGTTCGCCATAATATTGATGGCGTTATTGCAACCA  >  NC_000913/1005232‑1005503

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