Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,005,691	T→A	I308I (ATT→ATA)	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,691	0	T	A	100.0%	46.8 / NA	16	I308I (ATT→ATA)	pyrD	dihydroorotate dehydrogenase, type 2
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base A (6/10);  total (6/10)

TCCAGGTTGCCGATAGTTTAGTTCGCCATAATATTGATGGCGTTATTGCAACCAATACCACACTCGATCGTTCTCTTGTTCAGGGAATGAAAAATTGCGATCAAACCGGTGGCTTAAGTGGTCGTCCGCTTCAGTTAAAAAGCACCGAAATTATTCGCCGCTTGTCACTGGAATTAAACGGTCGCTTACCGATCATCGGTGTTGGCGGCATCGACTCGGTTATCGCTGCGCGTGAAAAGATTGCTGCGGGTGCCTCACTGGTGCAAATTTATTCTGGTTTTATTTTTAAAGGTCCGCCGCTGATTAAAGAAATCGTTACCCATATCTAATTATTTCTTCGAATTCCTCTTTATAAGATAATCAGGGCTTTATTTTCAGCCCTGGTTGTTTTATATTCATCGCTGTTGCTTATTTAGACATTTTGTACTTTTATTGTTGAGGTTATTA  >  NC_000913/1005450‑1005896                                                                                                                                                                                                                                                  |                                                                                                                                                                                                              tCCAGGTTGCCGATAGTTTAGTTCGCCATAATATTGATGGCGTTATTGCAACCAATACCACACTCGATCGTTCTCTTGTTCAGGGAATGAAAAATTGCGATCAAACCGGTGGCTTAAGTGGTCGTCCGCTTCAGTTAAAAAGCACCGAAATTATTCGCCGCTTGTCACTGGAATTAAACGGTCGCTTACCGATCATCGGTGTTGGCGGCATCGACTCGGTTATCGCTGCGCGTGAAAAGATAGCTGCgggg                                                                                                                                                                                                      >  1:111170/1‑250 (MQ=255)                    agttCGCCATAATATTGATGGCGTTATTGCAACCAATACCACACTCGATCGTTCTCTTGTTCAGGGAATGAAAAATTGCGATCAAACCGGTGGCTTAAGTGGTCGTCCGCTTCAGTTAAAAAGCACCGAAATTATTCGCCGCTTGTCACTGGAATTAAACGGTCGCTTACCGATCATCGGTGTTGGCGGCATCGACTCGGTTATCGCTGCGCGTGAAAAGATAGCTGCGGGTGCTTCACTGGTGCAAAtt                                                                                                                                                                                    <  1:201033/250‑1 (MQ=255)                                     aTGGCGTTATTGCAACCAATACCACACTCGATCGTTCTCTTGTTCAGGGAATGAAAAATTGCGATCAAACCGGTGGCTTAAGTGGTCGTCCGCTTCAGTTAAAAAGCACCGAAATTATTCGCCGCTTGTCACTGGAATTAAACGGTCGCTTACCGATCATCGGTGTTGGCGGCATCGACTCGGTTATCGCTGCGCGTGAAAAGATAGCTGCGGGTGCTTCCCTGGTGCCATTTTATTCTGGtttttttttt                                                                                                                                                                  >  2:290743/1‑251 (MQ=255)                                                          ccACACTCGATCGTTCTCTTGTTCAGGGAATGAAAAATTGCGATCAAACCGGTGGCTTAAGTGGTCGTCCGCTTCAGTTAAAAAGCACCGAAATTATTCGCCGCTTGTCACTGGAATTAAACGGTCGCTTACCGATCATCGGTGTTGGCGGCATCGACTCGGTTATCGCTGCGCGTGAAAAGATAGCTGCGGGTGCTTCACTGGTGCAAATTTATTCTGGTTTTATTTTTAAAGGTCCGCCGCTGATTaaa                                                                                                                                             <  1:290743/251‑1 (MQ=255)                                                                  gATCGTTCTCTTGTTCAGGGAATGAAAAATTGCGATCAAACCGGTGGCTTAAGTGGTCGTCCGCTTCAGTTAAAAAGCACCGAAATTATTCGCCGCTTGTCACTGGAATTAAACGGTCGCTTACCGATCATCGGTGTTGGCGGCATCGACTCGGTTATCGCTGCGCGTGAAAAGATAGCTGCGGGTGCTTCACTGGTGCAAATTTATTCTGGTTTTATTTTTAAAGGTCCGCCGCTGATTAAAGAAATCGt                                                                                                                                     >  1:71684/1‑251 (MQ=255)                                                                       ttCTCTTGTTCAGGGAATGAAAAATTGCGATCAAACCGGTGGCTTAAGTGGTCGTCCGCTTCAGTTAAAAAGCACCGAAATTATTCGCCGCTTGTCACTGGAATTAAACGGTCGCTTACCGATCATCGGTGTTGGCGGCATCGACTCGGTTATCGCTGCGCGTGAAAAGATAGCTGCGGGTGCTTCACTGGGGCAAATTTATTCTGGTTTTATTTTTAAAGGTCCGCCGCTGATTAAAGAAATCGTTAccc                                                                                                                                >  1:326156/1‑251 (MQ=255)                                                                                      aaTGAAAAATTGCGATCAAACCTGTGGCTTAAGTGGTCGTCCGCTTCAGTTAAAAAGCACCGAAATTATTCGCCGCTTGTCACTGGAATTAAACGGTCGCTTACCGATCATCGGTGTTGGCGGCATCGACTCGGTTATCGCTGCGCGTGAAAAGATAGCTGCGGGTGCTTCACTGGTGCAAATTTATTCTGGTTTTATTTTTAAAGGTCCGCCGCTGATTAAAGAAATCGTTACCCATATCTAATTATttc                                                                                                                 <  1:134258/251‑1 (MQ=255)                                                                                                   gATCAAACCGGTGGCTTAAGTGGTCGTCCGCTTCAGTTAAAAAGCACCGAAATTATTCGCCGCTTGTCACTGGAATTAAACGGTCGCTTACCGATCATCGGTGTTGGCGGCATCGACTCGGTTATCGCTGCGCGTGAAAAGATAGCTGCGGGTGCTTCACTGGTGCAAATTTATTCTGGTTTTATTTTTAAAGGTCCGCCGCTGATTAAAGAAATCGTTACCCATATCTAATTATttcttc                                                                                                              <  1:116296/241‑1 (MQ=255)                                                                                                   gATCAAACCGGTGGCTTAAGTGGTCGTCCGCTTCAGTTAAAAAGCACCGAAATTATTCGCCGCTTGTCACTGGAATTAAACGGTCGCTTACCGATCATCGGTGTTGGCGGCATCGACTCGGTTATCGCTGCGCGTGAAAAGATAGCTGCGGGTGCTTCACTGGTGCAAATTTATTCTGGTTTTATTTTTAAAGGTCCGCCGCTGATTAAAGAAATCGTTACCCATATCTAATTATttcttc                                                                                                              >  2:116296/1‑241 (MQ=255)                                                                                                                   tAAGTGGTCGTCCGCTTCAGTTAAAAAGCACCGAAATTATTCGCCGCTTGTCACTGGAATTAAACGGTCGCTTACCGATCATCGGTGTTGGCGGCATCGACTCGGTTATCGCTGCGCGTGAAAAGATAGCTGCGGGTGCTTCACTGGTGCAAATTTATTCTGGTTTTATTTTTAAAGGTCCGCCGCTGATTAAAGAAATCGTTACCCATATCTAATTATTTCTTCGAATTCCTCTTTATAAGATAATCa                                                                                      <  1:339712/249‑1 (MQ=255)                                                                                                                                                            cgccgcTTGTCACTGGAATTAAACGGTCGCTTACCGATCATCGGTGTTGGCGGCATCGACTCGGTTATCGCTGCGCGTGAAAAGATAGCTGCGGGTGCTTCACTGGTGCAAATTTATTCTGGTTTTATTTTTAAAGGTCCGCCGCTGATTAAAGAAATCGTTACCCATATCTAATTATTTCTTCGAATTCCTCTTTATAAGATAATCAGGGCTTTATTTTCAGCCCTGGTTGTTTTATATTCATCGCTGtt                                           <  1:344556/251‑1 (MQ=255)                                                                                                                                                                      cACTGGAATTAAACGGTCGCTTACCGATCATCGGTGTTGGCGGCATCGACTCGGTTATCGCTGCGCGTGAAAAGATAGCTGCGGGTGCTTCACTGGTGCAAATTTATTCTGGTTTTATTTTTAAAGGTCCGCCGCTGATTAAAGAAATCGTTACCCATATCTAATTATTTCTTCGAATTCCTCTTTATAAGATAATCAGGGCTTTATTTTCAGCCCTGGTTGTTTTATATTCATCGCTGTTGCTTATTTAg                                 <  1:283707/251‑1 (MQ=255)                                                                                                                                                                                 aaCGGTCGCTTACCGATCATCGGTGTTGGCGGCATCGACTCGGTTATCGCTGCGCGTGAAAAGATAGCTGCGGGTGCTTCACTGGTGCAAATTTATTCTGGTTTTATTTTTAAAGGTCCGCCGCTGATTAAAGAAATCGTTACCCATATCTAATTATTTCTTCGAATTCCTCTTTATAAGATAATCAGGGCTTTATTTTCAGCCCTGGTTGTTTTATATTCATCGCTGTTGCTTATTTAGACATTTTGTAc                      <  1:233970/251‑1 (MQ=255)                                                                                                                                                                                                   accgggtttGGCGGCATCGACTCGGTTATCGCTGCGCGGGAAAAGATAGCTGCGTGTGCTTCACTGGTGCAAATTTATTCTGTTTTTATTTTTAAAGGTCCGCCGCTGATTAAAGAAATCGTTACCCATATCTAATTATTTCTTCGAATTCCTCTTTATAAGATAATCAGGGCTTTATTTTCAGCCCTGGTTGTTTTATATTCATCGCTGTTGCTTATTTAGACATTTTGTACTTTTATTGTTGAGGttat    <  2:71684/244‑1 (MQ=255)                                                                                                                                                                                                     cGGTGTTGGCGGCATCGACTCGGTTATCGCTGCGCGTGAAAAGATAGCTGCGGGTGCTTCACTGGTGCAAATTTATTCTGGTTTTATTTTTAAAGGTCCGCCGCTGATTAAAGAAATCGTTACCCATATCTAATTATTTCTTCGAATTCCTCTTTATAAGATAATCAGGGCTTTATTTTCAGCCCTGGTTGTTTTATATTCATCGCTGTTGCTTATTTAGACATTTTGTACTTTTATTGTTGAGGttatta  <  1:213221/251‑1 (MQ=255)                                                                                                                                                                                                      ggTGTTGGCGGCATCGACTCGGTTATCGCTGCGCGTGAAAAGATAGCTGCGGGTGCTTCACTGGTGCAAATTTATTCTGGTTTTATTTTTAAAGGTCCGCCGCTGATTAAAGAAATCGTTACCCATATCTAATTATTTCTTCGAATTCCTCTTTATAAGATAATCAGGGCTTTATTTTCAGCCCTGGTTGTTTTATATTCATCGCTGTTGCTTATTTAGACATTTTGTACTTTTATTGTTGAGGttatta  >  1:92548/1‑250 (MQ=255)                                                                                                                                                                                                                                                  |                                                                                                                                                                                                              TCCAGGTTGCCGATAGTTTAGTTCGCCATAATATTGATGGCGTTATTGCAACCAATACCACACTCGATCGTTCTCTTGTTCAGGGAATGAAAAATTGCGATCAAACCGGTGGCTTAAGTGGTCGTCCGCTTCAGTTAAAAAGCACCGAAATTATTCGCCGCTTGTCACTGGAATTAAACGGTCGCTTACCGATCATCGGTGTTGGCGGCATCGACTCGGTTATCGCTGCGCGTGAAAAGATTGCTGCGGGTGCCTCACTGGTGCAAATTTATTCTGGTTTTATTTTTAAAGGTCCGCCGCTGATTAAAGAAATCGTTACCCATATCTAATTATTTCTTCGAATTCCTCTTTATAAGATAATCAGGGCTTTATTTTCAGCCCTGGTTGTTTTATATTCATCGCTGTTGCTTATTTAGACATTTTGTACTTTTATTGTTGAGGTTATTA  >  NC_000913/1005450‑1005896

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