Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,121,934	A→G	I225T (ATC→ACC)	pyrC ←	dihydro‑orotase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,121,934	0	A	G	100.0%	168.0 / NA	47	I225T (ATC→ACC)	pyrC	dihydro‑orotase
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (17/30);  total (17/30)

AAGCAGCCCGCGCAGCCGCAACTGCTCTCTTTGCGATGACGTGCATGTGGCGCAGAATCCGTACCGAGGAATACTCGATTAAAACCGCTGGCGACCAGTTCACGCAATGCCTGTTGGTGAATATTACGTTTGAGGATGGGTAGACAATACAGGTGCGGACGCACGCCTCCAACCAGCATATGGTTGCGGTTAAACATCAGATGCTGCGGAGTGATGGTGGCAGCCAGCCGTTCATTTCCGTCACGGACATAGTCGGCAGCATCTT  >  NC_000913/1121799‑1122063                                                                                                                                        |                                                                                                                                  aaGCAGCCCGCGCAGCCGCAACTGCTCTCTTTGCGATGACGTGCATGTGGCGCAGAATCCGTACCGAGGAATACTCGATTAAAACCGCTGGCGACCCGTTCACGCAATGCCTGTTGGTGAATATTACGTTTGAGGGTgg                                                                                                                                <  1:526243/139‑1 (MQ=255) aaGCAGCCCGCGCAGCCGCAACTGCTCTCTTTGCGATGACGTGCATGTGGCGCAGAATCCGTACCGAGGAATACTCGATTAAAACCGCTGGCGACCAGTTCACGCAATGCCTGTTGGTGAATATTACGTTTGAGGGTgg                                                                                                                                <  2:330297/139‑1 (MQ=255) aaGCAGCCCGCGCAGCCGCAACTGCTCTCTTTGCGATGACGTGCATGTGGCGCAGAATCCGTACCGAGGAATACTCGATTAAAACCGCTGGCGACCAGTTCACGCAATGCCTGTTGGTGAATATTACGTTTGAGGGTgg                                                                                                                                <  1:669354/139‑1 (MQ=255)    cagcCCGCGCAGCCGCAACTGCTCTCTTTGCGATGACGTGCATGTGGCGCAGAATCCGTACCGAGGAATACTCGATTAAAACCGCTGGCGACCAGTTCACGCAATGCCTGTTGGTGAATATTACGTTTGAGGGTGGGTa                                                                                                                             <  2:714874/139‑1 (MQ=255)    cagcCCGCGCAGCCGCAACTGCTCTCTTTGCGATGACGTGCATGTGGCGCAGAATCCGTACCGAGGAATACTCGATTAAAACCGCTGGCGACCAGTTCACGCAATGCCTGTTGGTGAATATTACGTTTGAGGGTGGGTa                                                                                                                             <  1:455071/139‑1 (MQ=255)      gcCCGCGCAGCCGCAACTGCTCTCTTTGCGATGACGTGCATGTGGCGCAGAATCCGTACCGAGGAATACTCGATTAAAACCGCTGGCGACCAGTTCACGCAATGCCTGTTGGTGAATATTACGTTTGAGGGTGGGTAGa                                                                                                                           <  1:747722/139‑1 (MQ=255)      gcCCGCGCAGCCGCAACTGCTCTCTTTGCGATGACGTGCATGTGGCGCAGAATCCGTACCGAGGAATACTCGATTAAAACCGCTGGCGACCAGTTCACGCAATGCCTGTTGGTGAATATTACGTTTGAGGGTGGGTAGa                                                                                                                           <  1:447379/139‑1 (MQ=255)        ccGCGCAGCCGCAACTGCTCTCTTTGCGATGACGTGCATGTGGCGCAGAATCCGTACCGAGGAATACTCGATTAAAACCGCTGGCGACCAGTTCACGCAATGCCTGTTGGTGAATATTACGTTTGAGGGTGGGTAGACa                                                                                                                         >  2:560477/1‑139 (MQ=255)                            tctTTGCGATGACGTGCATGTGGCGCAGAATCCGTACCGAGGAATACTCGATTAAAACCGCTGGCGACCAGTTCACGCAATGCCTGTTGGTGAATATTACGTTTGAGGGTGGGTAGACAATACAGGTGCGGacgcacgc                                                                                                     <  2:39363/139‑1 (MQ=255)                                  cGATGACGTGCATGTGGCGCAGAATCCGTACCGAGGAATACTCGATTAAAACCGCTGGCGACCAGTTCACGCAATGCCGGTTGGTGAATATTACGTTTGAGGGGGGGTAGACAATACAGGTGCGGACGCACGCCTCCaa                                                                                               <  2:604298/139‑1 (MQ=255)                                        cGTGCATGTGGCGCAGAATCCGTACCGAGGAATACTCGATTAAAACCGCTGGCGACCAGTTCACGCAATGCCTGTTGGTGAATATTACGTTTGAGGGTGGGTAGACAATACAGGTGCGGACGCACGCCTcc                                                                                                 >  1:333031/1‑131 (MQ=255)                                        cGTGCATGTGGCGCAGAATCCGTACCGAGGAATACTCGATTAAAACCGCTGGCGACCAGTTCACGCAATGCCTGTTGGTGAATATTACGTTTGAGGGTGGGTAGACAATACAGGTGCGGACGCACGCCTcc                                                                                                 <  2:333031/131‑1 (MQ=255)                                                 ggCGCAGAATCCGTACCGAGGAATACTCGATTAAAACCGCTGGCGACCAGTTCACGCAATGCCTGTTGGTGAATATTACGTTTGAGGGTGGGTAGACAATACAGGTGCGGACGCACGCCTCCAACCAGCATATGGTTGc                                                                                >  1:289537/1‑139 (MQ=255)                                                 ggCGCAGAATCCGTACCGAGGAATACTCGATTAAAACCGCTGGCGACCAGTTCACGCAATGCCTGTTGGTGAATATTACGTTTGAGGGTGGGTAGACAATACAGGTGCGGACGCACGCCTCCAACCAGCATATGGTTGc                                                                                >  1:338535/1‑139 (MQ=255)                                                     cAGAATCCGTATCGAGGAATACTCGATTAAAACCGCTGGCGACCAGTTCACGCAATGCCTGTTGGTGAATATTACGTTTGAGGGTGGGTAGACAATACAGGTGCGGACGCACGCCTCCAACCAGCATATGGTTGCGGtt                                                                            <  1:522097/139‑1 (MQ=255)                                                                  gAGGAATACTCGATTAAAACCGCTGGCGACCAGTTCACGCAATGCCTGTTGGTGAATATTACGTTTGAGGGTGGGTAGACAATACAGGTGCGGACGCACGCCTCCAACCAGCATATGGTTGCGGTTAAACATCAGAtgc                                                               >  1:459707/1‑139 (MQ=255)                                                                    ggAATACTCGATTAAAACCGCTGGCGCCCAGTTCACGCAATGCCTGTTGGTGAATATTACGTTTGAGGGTGGGTAGACAATACAGGTGCGGACGCACGCCTCCAACCAGCATATGGTTGCGGTTAAACATCAGAtgctg                                                             <  1:300219/139‑1 (MQ=255)                                                                    ggAATACTCGATTAAAACCGCTGGCGACCAGTTCACGCAATGCCTGTTGGTGAATATTACGTTTGAGGGTGGGTAGACAATACAGGTGCGGACGCACGCCTCCAACCAGCATATGGTTGCGGTTAAACATCAGAtgctg                                                             <  2:459707/139‑1 (MQ=255)                                                                     gAATACTCGATTAAAACCGCTGGCGACCAGTTCACGCAATGCCTGTTGGTGAATATTACGTTTGAGGGTGGGTAGa                                                                                                                           >  1:393532/1‑76 (MQ=255)                                                                     gAATACTCGATTAAAACCGCTGGCGACCAGTTCACGCAATGCCTGTTGGTGAATATTACGTTTGAGGGTGGGTAGa                                                                                                                           <  2:393532/76‑1 (MQ=255)                                                                     gAATACTCGATTAAAACCGCTGGCGACCAGTTCACGCAATGCCTGTTGGTGAATATTACGTTTGAGGGTGGGTAGACAATACAGGTGCGGACGCACGCCTCCAACCAGCATATGGTTGCGGTTAAACATCAGAtgctgc                                                            <  1:515720/139‑1 (MQ=255)                                                                            cGATTAAAACCGCTGGCGACCAGTTCACGCAATGCCTGTTGGTGAATATTACGTTTGAGGGTGGGTAGACAATACAGGTGCGGACGCACGCCTCCAACCAGCATATGGTTGCGGTTAAACATCAGATGCTGCGGAGTGa                                                     <  2:702510/139‑1 (MQ=255)                                                                             gATTAAAACCGCTGGCGACCAGTTCACGCAATGCCTGTTGGTGAATATTACGTTTGAGGGTGGGTAGACAATACAGGTGCGGACGCACGCCTCCAACCAGCATATGGTTGCGGt                                                                             >  2:436367/1‑114 (MQ=255)                                                                             gATTAAAACCGCTGGCGACCAGTTCACGCAATGCCTGTTGGTGAATATTACGTTTGAGGGTGGGTAGACAATACAGGTGCGGACGCACGCCTCCAACCAGCATATGGTTGCGGt                                                                             <  1:436367/114‑1 (MQ=255)                                                                               ttAAAACCGCTGGCGACCAGTTCACGCAATGCCTGTTGGTGAATATTACGTTTGAGGGTGGGTAGACAATACAGGTGCGGACGCACGCCTCCAACCAGCATATGGTTGCGGTTAAACATCAGATGCTGCGGa                                                         >  1:437547/1‑132 (MQ=255)                                                                               ttAAAACCGCTGGCGACCAGTTCACGCAATGCCTGTTGGTGAATATTACGTTTGAGGGTGGGTAGACAATACAGGTGCGGACGCACGCCTCCAACCAGCATATGGTTGCGGTTAAACATCAGATGCTGCGGa                                                         <  2:437547/132‑1 (MQ=255)                                                                                       gCTGGCGACCAGTTCACGCAATGCCTGTTGGTGAATATTACGTTTGAGGGTgg                                                                                                                                >  1:117965/1‑53 (MQ=255)                                                                                       gCTGGCGACCAGTTCACGCAATGCCTGTTGGTGAATATTACGTTTGAGGGTgg                                                                                                                                <  2:117965/53‑1 (MQ=255)                                                                                       gCTGGCGACCAGTTCACGCAATGCCTGTTGGTGAATATTACGTTTGAGGGTGGGTAGACAATACAGGTGCGGACGCACGCCTCCAACCAGCATATGGTTGCGGTTAAACATCAg                                                                   <  1:584090/114‑1 (MQ=255)                                                                                       gCTGGCGACCAGTTCACGCAATGCCTGTTGGTGAATATTACGTTTGAGGGTGGGTAGACAATACAGGTGCGGACGCACGCCTCCAACCAGCATATGGTTGCGGTTAAACATCAg                                                                   >  2:584090/1‑114 (MQ=255)                                                                                       gCTGGCGACCAGTTCACGCAATGCCTGTTGGTGAATATTACGTTTGAGGGTGGGTAGACAATACAGGTGCGGACGCACGCCTCCAACCAGCATATGGTTGCGGTTAAACATCAGATGCTGCGGAGTGATGGTGGcagcc                                          <  2:338535/139‑1 (MQ=255)                                                                                       gCTGGCGACCAGTTCACGCAATGCCTGTTGGTGAATATTACGTTTGAGGGTGGGTAGACAATACAGGTGCGGACGCACGCCTCCAACCAGCATATGGTTGCGGTTAAACATCAGATGCTGCGGAGTGATGGTGGcagcc                                          <  1:63542/139‑1 (MQ=255)                                                                                       gCTGGCGACCAGTTCACGCAATGCCTGTTGGTGAATATTACGTTTGAGGGTGGGTAGACAATACAGGTGCGGACGCACGCCTCCAACCAGCATATGGTTGCGGTTAAACATCAGATGCTGCGGAGTGATGGTGGcagc                                           <  1:108133/138‑1 (MQ=255)                                                                                       gCTGGCGACCAGTTCACGCAATGCCTGTTGGTGAATATTACGTTTGAGGGTGGGTAGACAATACAGGTGCGGACGCACGCCTCCAACCAGCATATGGTTGCGGTTAAACATCAGATGCTGCGGAGTGATGATGGcagc                                           >  2:108133/1‑138 (MQ=255)                                                                                                   ttCACGCAATGCCTGTTGGTGAATATTACGTTTGAGGGTGGGTAGACAATACAGGTGCGGACGCACGCCTCCAACCAGCATATGGTTGCGGTTAAACATCAGATGCTGCGGAGTGATGGTGGCAGCCAGCCGTTCAttt                              <  1:620662/139‑1 (MQ=255)                                                                                                             gCCTGTTGGTGAATATTACGTTTGAGGGTGGGTAGACAATACAGGTGCGGACGCACGCCTCCAACCAGCATATGGTTGCGGTTAAACATCAGATGCTGCGGAGTGATGGTGGcagccag                                        <  2:297917/119‑1 (MQ=255)                                                                                                             gCCTGTTGGTGAATATTACGTTTGAGGGTGGGTAGACAATACAGGTGCGGACGCACGCCTCCAACCAGCATATGGTTGCGGTTAAACATCAGATGCTGCGGAGTGATGGTGGcagccag                                        >  1:297917/1‑119 (MQ=255)                                                                                                                  ttGGTGAATATTACGTTTGAGGGTGGGTAGACAATACAGGTGCGGACGCACGCCTCCAACCAGCATATGGTTGCGGTTAAACATCAGATGCTGCGGAGTGa                                                     >  2:658455/1‑101 (MQ=255)                                                                                                                  ttGGTGAATATTACGTTTGAGGGTGGGTAGACAATACAGGTGCGGACGCACGCCTCCAACCAGCATATGGTTGCGGTTAAACATCAGATGCTGCGGAGTGa                                                     <  1:658455/101‑1 (MQ=255)                                                                                                                    ggTGAATATTACGTTTGAGGGTGGGTAGACAATACAGGTGCGGACGCACGCCTCCAAc                                                                                              >  2:95278/1‑58 (MQ=255)                                                                                                                    ggTGAATATTACGTTTGAGGGTGGGTAGACAATACAGGTGCGGACGCACGCCTCCAAc                                                                                              <  1:95278/58‑1 (MQ=255)                                                                                                                      tGAATATTACGTTTGAGGGTGGGTAGACAATACAGGTGCGGACGCACGCCTCCAACCAGCATATGGTTGCGGTTAAACATCAGATGCTGCg                                                           <  1:555513/91‑1 (MQ=255)                                                                                                                      tGAATATTACGTTTGAGGGTGGGTAGACAATACAGGTGCGGACGCACGCCTCCAACCAGCATATGGTTGCGGTTAAACATCAGATGCTGCg                                                           >  2:555513/1‑91 (MQ=255)                                                                                                                            ttACGTTTGAGGGTGGGTCGACAATACAGGCCCGGACCCCCGCCTGCAACCAGCATCTGGTTACGGATAGACATCAGATGCTCCGGCGTGATGGTGGCAGCCACGCGTTCGTTTACGTTACGGACATAGTCGGCGGCAt     >  2:35473/1‑139 (MQ=255)                                                                                                                               cGTTTGAGGGTGGGTAGACAATACAGGTGCGGACGCACGCCTCCAACCAGCATATGGTTGCGGTTAAACATCAGATGCTGCGGAGTGATGGTGGCAGCCAGCCGTTCATTTCCGTCACGGACATAGTCGGCAGCATCtt  <  2:289537/139‑1 (MQ=255)                                                                                                                                    gAGGGTGGGTAGACAATACAGGTGCGGACGCACGCCTCCAACCAGCATATGGTTGCGGTTaa                                                                          <  1:139109/62‑1 (MQ=255)                                                                                                                                    gAGGGTGGGTAGACAATACAGGTGCGGACGCACGCCTCCAACCAGCATATGGTTGCGGTTaa                                                                          >  2:139109/1‑62 (MQ=255)                                                                                                                                        |                                                                                                                                  AAGCAGCCCGCGCAGCCGCAACTGCTCTCTTTGCGATGACGTGCATGTGGCGCAGAATCCGTACCGAGGAATACTCGATTAAAACCGCTGGCGACCAGTTCACGCAATGCCTGTTGGTGAATATTACGTTTGAGGATGGGTAGACAATACAGGTGCGGACGCACGCCTCCAACCAGCATATGGTTGCGGTTAAACATCAGATGCTGCGGAGTGATGGTGGCAGCCAGCCGTTCATTTCCGTCACGGACATAGTCGGCAGCATCTT  >  NC_000913/1121799‑1122063

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