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%	150.6 / NA	40	I225T (ATC→ACC)	pyrC	dihydro‑orotase
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (16/24);  total (16/24)

AAGCAGCCCGCGCAGCCGCAACTGCTCTCTTTGCGATGACGTGCATGTGGCGCAGAATCCGTACCGAGGAATACTCGATTAAAACCGCTGGCGACCAGTTCACGCAATGCCTGTTGGTGAATATTACGTTTGAGGATGGGTAGACAATACAGGTGCGGACGCACGCCTCCAACCAGCATATGGTTGCGGTTAAACATCAGATGCTGCGGAGTGATGGTGGCAGCCAGCCGTTCATTTCCGTCACGGACATAGTCGGCAGCATCTTTGGTGG  >  NC_000913/1121799‑1122069                                                                                                                                        |                                                                                                                                        aaGCAGCCCGCGCAGCCGCAACTGCTCTCTTTGCGATGACGTGCATGTGGCGCAGAATCCGTACCGAGGAATACTCGATTAAAACCGCTGGCGACCAGTTCACGCAATGCCTGTTGGTGAATATTACGTTTGAGGGTgg                                                                                                                                      <  2:316554/139‑1 (MQ=255) aaGCAGCCCGCGCAGCCGCAACTGCTCTCTTTGCGATGACGTGCATGTGGCGCAGAATCCGTACCGAGGAATACTCGATTAAAACCGCTGGCGACCAGTTCACGCAATGCCTGTTGGTGAATATTACGTTTGAGGGTgg                                                                                                                                      <  1:512182/139‑1 (MQ=255) aaGCAGCCCGCGCAGCCGCAACTGCTCTCTTTGCGATGACGTGCATGTGGCGCAGAATCCGTACCGAGGAATACTCGATTAAAACCGCTGGCGACCAGTTCACGCAATGCCTGTTGGTGAATATTACGTTTGAGGGTgg                                                                                                                                      <  1:223584/139‑1 (MQ=255)    cagcCCGCGCAGCCGCAACTGCTCTCTTTGCGATGACGTGCATGTGGCGCAGAATCCGTACCGAGGAATACTCGATTAAAACCGCTGGCGACCAGTTCACGCAATGCCTGTTGGTGAATATTACGTTTGAGGGTGGGTa                                                                                                                                   <  1:458487/139‑1 (MQ=255)     agcCCGCGCAGCCGCAACTGCTCTCTTTGCGATGACGTGCATGTGGCGCAGAATCCGTACCGAGGAATACTCGATTAAAACCGCTGGCGACCAGTTCACGCAATGCCTGTTGGTGAATATTACGTTTGAGGGTGGGTAg                                                                                                                                  >  2:888883/1‑139 (MQ=255)      gcCCGCGCAGCCGCAACTGCTCTCTTTGCGATGACGTGCATGTGGCGCAGAATCCGTACCGAGGAATACTCGATTAAAACCGCTGGCGACCAGTTCACGCAATGCCTGTTGGTGAATATTACGTTTGAGGGTGGGTAGa                                                                                                                                 <  2:295319/139‑1 (MQ=255)      gcCCGCGCAGCCGCAACTGCTCTCTTTGCGATGACGTGCATGTGGCGCAGAATCCGTACCGAGGAATACTCGATTAAAACCGCTGGCGACCAGTTCACGCAATGCCTGTTGGTGAATATTACGTTTGAGGGTGGGTAGa                                                                                                                                 <  2:527837/139‑1 (MQ=255)                                tGCGATGACGTGCATGTGGCGCAGAATCCGTACCGAGGAATACTCGATTAAAACCGCTGGCGACCAGTTCACGCAATGCCTGTTGGTGAATATTACGTTTGAGGGTGGGTAGACAATACAGGTGCGGACGCACGCCTcc                                                                                                       <  1:22330/139‑1 (MQ=255)                                   gATGACGTGCATGTGGCGCAGAATCCGTACCGAGGAATACTCGATTAAAACCGCTGGCGACCAGTTCACGCAATGCCTGTTGGTGAATATTACGTTTGAGGGTGGGTAGACAATACAGGTGCGGACGCACGCCTCCAAc                                                                                                    <  2:798654/139‑1 (MQ=255)                                         gTGCATGTGGCGCAGAATCCGTACCGAGGAATACTCGATTAAAACCGCTGGCGACCAGTTCACGCAATGCCTGTTGGTGAATATTACGTTTGAGGGTGGGTAGACAATACAGGTGCGGACGCACGCCTCCAACCAGCat                                                                                              <  2:572264/139‑1 (MQ=255)                                                 ggCGCAGAATCCGTACCGAGGAATACTCGATTAAGACCGCTGGCGACCAGTTCACGCAATGCCTGTTGGTGAATATTACGTTTGAGGGTGGGTAGACAATACAGGTGCGGACGCACGCCTCCAACCAGCATATGGTTGc                                                                                      >  1:232049/1‑139 (MQ=255)                                                    gcAGAATCCGTACCGAGGAATACTCGATTAAAACCGCTGGCGACCAGTTCACGCAATGCCTGTTGGTGAATATTACGTTTGAGGGTGGGTAGACAATACAGGTGCGGACGCACGCCTCCAACCAGCATATGGTTGCGGt                                                                                   <  1:632638/139‑1 (MQ=255)                                                    gcAGAATCCGTACCGAGGAATACTCGATTAAAACCGCTGGCGACCAGTTCACGCAATGCCTGTTGGTGAATATTACGTTTGAGGGTGGGTAGACAATACAAGTGCGGACGCACGCCTCCAACCAGCATATGGTTGCGGt                                                                                   <  1:159787/139‑1 (MQ=255)                                                       gAATTCGTACCGAGGAATACTTGAATAAAACCGCTTGCGAACAATTCACGCAATGCCTGTTGGGTAATATTACGTTTGAGGGTTGGTAGACAATACAGGTGCGGACGCACGCCTCCAACCAGCATATGGTTGCGGTTaa                                                                                <  2:578198/139‑1 (MQ=255)                                                        aaTCCGTACCGAGGAATACTCGATTAAAACCGCTGGCGACCAGTTCACGCAATGCCTGTTGGTGAATATTACGTTTGAGGGTgg                                                                                                                                      >  2:368160/1‑84 (MQ=255)                                                        aaTCCGTACCGAGGAATACTCGATTAAAACCGCTGGCGACCAGTTCACGCAATGCCTGTTGGTGAATATTACGTTTGAGGGTgg                                                                                                                                      <  1:368160/84‑1 (MQ=255)                                                        aaTCCGTACCGAGGAATACTCGATTAAAACCGCTGGCGACCAGTTCACGCAATGCCTGTTGGTGAATATTACGTTTGAGGGTGGGTAGACAATACAGGTGCGGACGCACGCCTCCAACCAGCATATGGTTGCGGTTaaa                                                                               >  2:247439/1‑139 (MQ=255)                                                               aCCGAGGAATACTCGATTAAAACCGCTGGCGACCAGTTCACGCAATGCCTGTTGGTGAATATTACGTTTGAGGGTGGGTAGACAATACAGGTGCGGacgca                                                                                                              >  1:829971/1‑101 (MQ=255)                                                               aCCGAGGAATACTCGATTAAAACCGCTGGCGACCAGTTCACGCAATGCCTGTTGGTGAATATTACGTTTGAGGGTGGGTAGACAATACAGGTGCGGacgca                                                                                                              <  2:829971/101‑1 (MQ=255)                                                                 cGAGGAATACTCGATTAAAACCGCTGGCGACCAGTTCACGCAATGCCTGTTGGTGAATATTACGTTTGAGGGTGGGTAGACAATACAGGTGCGGACGCACGCCTCCAACCAGCATATGGTTGCGGTTAAACATCAGAtg                                                                      <  1:247439/139‑1 (MQ=255)                                                                    ggAATACTCGATTAAAACCGCTGGCGACCAGTTCACGCAATGCCTGTTGGTGAATATTACGTTTGAGGGTGGGTAGACAATACAGGTGCGGACGCACGCCTCCAAc                                                                                                    >  1:859823/1‑106 (MQ=255)                                                                    ggAATACTCGATTAAAACCGCTGGCGACCAGTTCACGCAATGCCTGTTGGTGAATATTACGTTTGAGGGTGGGTAGACAATACAGGTGCGGACGCACGCCTCCAAc                                                                                                    <  2:859823/106‑1 (MQ=255)                                                                        tACTCGATTAAAACCGCTGGCGACCAGTTCACGCAATGCCTGGTGGTGAATATTACGTTTGAGGGTGGGTAGACAATACAGGTGCGGACGCACGCCTCCAACCAGCATATGGTTGCGGTTAAACATCAGATGCTGCGGa                                                               <  1:1052560/139‑1 (MQ=255)                                                                          cTCGATTAAAACCGCTGGCGACCAGTTCACGCAATGCCTGTTGGTGAATATTACGTTTGAGGGTGGGTAGACAATACAGGTGCGGACGCACGCCTCCAACCAGCATATGGTTGCGGTTAAACATCAGATGCTGCGGAGt                                                             >  2:1016231/1‑139 (MQ=255)                                                                           tCGATTAAAACCGCTGGCGACCAGTTCACGCAATGCCTGTTGGTGAATATTACGTTTGAGGGTGGGTAGACAATACAGGTGCGGacgca                                                                                                              >  2:224072/1‑89 (MQ=255)                                                                           tCGATTAAAACCGCTGGCGACCAGTTCACGCAATGCCTGTTGGTGAATATTACGTTTGAGGGTGGGTAGACAATACAGGTGCGGacgca                                                                                                              <  1:224072/89‑1 (MQ=255)                                                                           tCGATTAAAACCGCTGGCGACCAGTTCACGCAATGCCTGTTGGTGAATATTACGTTTGAGGGTGGGTAGACAATACAGGTGCGGACGCACGCCTCCAACCAGCATATGGTTGCGGTTAAACATCAGATGCTGCGGAGTg                                                            >  2:514737/1‑139 (MQ=255)                                                                               ttAAAACCGCTGGCGACCAGTTCACGCAATGCCTGTTGGTGAATATTACGTTTGAGGGTGGGTa                                                                                                                                   <  2:17293/64‑1 (MQ=255)                                                                               ttAAAACCGCTGGCGACCAGTTCACGCAATGCCTGTTGGTGAATATTACGTTTGAGGGTGGGTa                                                                                                                                   >  1:17293/1‑64 (MQ=255)                                                                                                         cAATGCCTGTTGGTGAATATTACGTTTGAGGGTGGGTAGACAATACAGGTGCGGACGCACGCCTCCAACCAGCATATGGTTGCGGTTAAACATCAGATGCTGCGGAGTGATGGTGGCAGCCAGCCGTTCATTTCCGTCa                              <  1:1016231/139‑1 (MQ=255)                                                                                                             gCCTGTTGGTGAATATTACGTTTGAGGGTGGGTAGACAATACAGGTGCGGACGCACGCCTCCAAc                                                                                                    <  2:536173/65‑1 (MQ=255)                                                                                                             gCCTGTTGGTGAATATTACGTTTGAGGGTGGGTAGACAATACAGGTGCGGACGCACGCCTCCAAc                                                                                                    >  1:536173/1‑65 (MQ=255)                                                                                                                        aaTATTACGTTTGAGGGTGGGTAGACAATACAGGTGCGGACGCACGCCTCCAACCAGCATATGGTTGCGGTTAAACATCAGATGCTGCGGAGTGATGGTGGCAGCCAGCCGTTCATTTCCGTCACGGACATAGTCGgca               >  2:821297/1‑139 (MQ=255)                                                                                                                               cGTTTGAGGGTGGGTAGACAATACAGGTGCGGACGCACGCCTCCAACCAGCATATGGTTGc                                                                                      <  2:1029913/61‑1 (MQ=255)                                                                                                                               cGTTTGAGGGTGGGTAGACAATACAGGTGCGGACGCACGCCTCCAACCAGCATATGGTTGc                                                                                      >  1:1029913/1‑61 (MQ=255)                                                                                                                               cGTTTGAGGGTGGGTAGACAATACAGGTGCGGACGCACGCCTCCAACCAGCATATGGTTGCGGTTAAACATCAGATGCTGCGGAGTGATGGTGGca                                                   >  1:832287/1‑96 (MQ=255)                                                                                                                               cGTTTGAGGGTGGGTAGACAATACAGGTGCGGACGCACGCCTCCAACCAGCATATGGTTGCGGTTAAACATCAGATGCTGCGGAGTGATGGTGGca                                                   <  2:832287/96‑1 (MQ=255)                                                                                                                               cGTTTGAGGGTGGGTAGACAATACAGGTGCGGACGCACGCCTCCAACCAGCATATGGTTGCGGTTAAACATCAGATGCTGCGGAGTGATGGTGGCAGCCAGCCGTTCATTTCCGTCACGGACATAGTCGGCAGCATCtt        <  1:821297/139‑1 (MQ=255)                                                                                                                                    gAGGGTGGGTAGACAATACAGGTGCGGACGCACGCCTCCAACCAGCATATGGTTGCGGTTAAACATCAGATGCTGCGGAGTGATGGTGGCAGCCAGCCGTTCATTTCCGTCACGGACATAGTCGGCAGCATCTTtggtg   >  2:1091046/1‑139 (MQ=255)                                                                                                                                     aGGGTGGGTAGACAATACAGGTGCGGACGCACGCCTCCAACCAGCATATGGTTGCGGTTAAACATCAGATGCTGCGGAGTGATGGTGGCAGCCAGCCGTTCATTTCCGTCACGGACATAGTCGGCAGCATCTTtggtgg  >  1:76186/1‑139 (MQ=255)                                                                                                                                        |                                                                                                                                        AAGCAGCCCGCGCAGCCGCAACTGCTCTCTTTGCGATGACGTGCATGTGGCGCAGAATCCGTACCGAGGAATACTCGATTAAAACCGCTGGCGACCAGTTCACGCAATGCCTGTTGGTGAATATTACGTTTGAGGATGGGTAGACAATACAGGTGCGGACGCACGCCTCCAACCAGCATATGGTTGCGGTTAAACATCAGATGCTGCGGAGTGATGGTGGCAGCCAGCCGTTCATTTCCGTCACGGACATAGTCGGCAGCATCTTTGGTGG  >  NC_000913/1121799‑1122069

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