Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,472,524	A→G	C4R (TGT→CGT)	pyrL ←	pyrBI operon leader peptide

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,472,524	0	A	G	100.0%	165.3 / NA	47	C4R (TGT→CGT)	pyrL	pyrBI operon leader peptide
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (24/23);  total (24/23)

TTAGCCATCTTTTATCTCCTGACGCCTGGGCAAAAAAAAGCCCCTCGATTGAGGGGCTGGGAATGGGTGATCAACGGGAAGAAAAACGGCAGGCCAGCGTCTTTTTTCAGACGCGGTAAGACAAAATGTCGAACACACTGAACCATACATCCTCCCGGCAAATTGTCCGGCATTATACTCATCGTCAGAAGCGGCGCAAGCATTTGATGCAATATTTTGTCAGCGCAAACGGTTTATTTGAATTAAAAGTCAAGGTATATGCATTTTTAT  >  NC_000913/4472388‑4472657                                                                                                                                         |                                                                                                                                      ttggCCATTTTTTATCCCCTGCCGCCTGGGCAAAAAAAAGCCCCTCGATTGAGGGGCTGGGAATGGGTGATCAACGGGAAGAAAAACGGCAGGCCAGCGTCTTTTTTCAGACGCGGTAAGACAAAATGTCGAACACGCTg                                                                                                                                    <  1:56034/137‑1 (MQ=255)                  ccTGACGCCTGGGCAAAAAAAAGCCCCTCGATTGAGGGGCTGGGAATGGGTGATCAACGGGAAGAAAAACGGCAGGCCAGCGTCTTTTTTCAGACGCGGTAAGACAAAATGTCGAACACGCTGAACCATACATCCTCCCg                                                                                                                   >  2:125439/1‑140 (MQ=255)                   cTGACGCCTGGGCAAAAAAAAGCCCCTCGATTTAGGGGCTGGGAATGGGTGATCAACGGGAAGAAAAACGGCAGGCCAGCGTCTTTTTTCAGACGCGGTAAGACAAAATGTCGAACACGCTGAACCATACAACCTCCCgg                                                                                                                  >  1:74332/1‑140 (MQ=255)                    tGACGCCTGGGCAAAAAAAAGCCCCTCGATTGAGGGGCTGGGAATGGGTGATCAACGGGAAGAAAAACGGCAGGCCAGCGTCTTTTTTCAGACGCGGTAAGACAAAATGTCGAACACGCTGAACCATACATCCTCCCGGc                                                                                                                 >  1:573500/1‑140 (MQ=255)                      aCGCCTGGGCAAAAAAAAGCCCCTCGATTGAGGGGCTGGGAATGGGTGATCAACGGGAAGAAAAACGGCAGGCCAGCGTCTTTTTTCAGACGCGGTAAGACAAAATGTCGAACACGCTGAACCATACATCCTCCCGGCaa                                                                                                               >  1:609664/1‑140 (MQ=255)                            gggCAAAAAAAAGCCCCTCGATTGAGGGGCTGGGAATGGGTGATCAACGGGAAGAAAAACGGCAGGCCAGCGTCTTTTTTCAGACGCGGTAAGACAAAATGTCGAACACGCTGAACCATACATCCTCCCGGCAAATTGTc                                                                                                         <  1:125439/140‑1 (MQ=255)                            gggCAAAAAAAAGCCCCTCGATTGAGGGGCTGGGAATGGGTGATCAACGGGAAGAAAAACGGCAGGCCAGCGTCTTTTTTCAGACGCGGTAAGACAAAATGTCGAACACGCTGAACCATACATCCTCCCGGCAAATTGTc                                                                                                         <  1:516386/140‑1 (MQ=255)                            gggCAAAAAAAAGCCCCTCGATTGAGGGGCTGGGAATGGGTGATCAACGGGAAGAAAAACGGCAGGCCAGCGTCTTTTTTCAGACGCGGTAAGACAAAATGTCGAACACGCTGAACCATACATCCTCCCGGCAAATTGTc                                                                                                         <  2:54509/140‑1 (MQ=255)                            gggCAAAAAAAAGCCCCCCGATTGAGGGGCTGGGATTGGGTGATCAACTGGAAGAAAAACGGCAGGCCAGCGTCTTTTTTCAGACGCGGTAAGACAAAATGTCGAACACGCTGAACCATACATCCTCCCGGCAAATTGTc                                                                                                         <  2:277226/140‑1 (MQ=255)                            gggCAAAAAAAAGCCCCCCGATTGAGGGGCTCGGATTGGCTGATCAACTGGAAGAAAAACGGCAGGCCAGCGTCTTTTTTCAGACGCGGTAAGACAAAATGTCGAACACGCTGAACCATACATCCTCCCGGCAAATTGTc                                                                                                         <  2:168586/140‑1 (MQ=255)                             ggCAAAAAAAAGGCCCTCGATTGAGGGGCTGGGAATGGGCGATCAACGGGAAGAAAAACGGCAGGCCAGCGTCTTTTTTCAGACGCGGTAAAACAAAATGACGAACACGCTGAACCACACATACTCCCGGCTAATTGTcc                                                                                                        >  2:323492/1‑140 (MQ=255)                             ggCAAAAAAAAGCCCCTCGATTGAGGGGCTGGGAATGGGTGATCAACGGGAAGAAAAACGGCAGGCCAGCGTCTTTTTTCAGACGCGGTAAGACAAAATGTCGAACACGCTGAAGCAGACATCCACCCGGCAAATTGTcc                                                                                                        >  1:195805/1‑140 (MQ=255)                             ggCAAAAAAAAGCCCCTCGATTGAGGGGCTGGGAATGGGTGATCAACGGGAAGAAAAACGGCAGGCCAGCGTCTTTTTTCAGACGCGGTAAGACAAAATGTCGAACACGCTGAACCATACATCCTCCCGGCAAATTGTcc                                                                                                        >  1:48251/1‑140 (MQ=255)                             ggCAAAAAAAAGCCCCTCGATTGAGGGGCTGGGAATGGGTGATCAACGGGAAGAAAAACGGCAGGCCAGCGTCTTTTTTCAGACGCGGTAAGACAAAATGTCGAACACGCTGAACCATACATCCTCCCGGCAAATTGTcc                                                                                                        >  2:148796/1‑140 (MQ=255)                             ggCAAAAAAAAGCCCCTCGATTGAGGGGCTGGGAATGGGTGATCAACGGGAAGAAAAACGGCAGGCCAGCGTCTTTTTTCAGACGCGGTAAGACAAAATGTCGAACACGCGGAAGCAGACATCCACCCGGCAAATTGTcc                                                                                                        >  1:168586/1‑140 (MQ=255)                                  aaaaaaGCCCCTCGATTGAGGGGCTGGGAATGGGTGATCAACGGGAAGAAAAACGGCAGGCCAGCGTCTTTTTTCAGACGCGGTAAGACAAAATGTCGAACACGCTGAACCATACATCCTCCCGGCAAATTGTCCGGCAt                                                                                                   >  1:658125/1‑140 (MQ=255)                                          cccTCGATTGAGGGGCTGGGAATGGGTGATCAACGGGAAGAAAAACGGCAGGCCAGCGTCTTTTTTCAGACGCGGTAAGACAAAATGTCGAACACGCTGAACCATACATCCTCCCGGCAAATTGTCCGGCATTATACTCa                                                                                           <  2:658125/140‑1 (MQ=255)                                            cTCGATTGAGGGGCTGGGAATGGGTGATCAACGGGAAGAAAAACGGCAGGCCAGCGTCTTTTTTCAGACGCGGTAAGACAAAATGTCGAACACGCTGAACCATACATCCTCCCGGCAAATTGTCCGGCATTATACTCATc                                                                                         >  1:499279/1‑140 (MQ=255)                                                    gggggCTGGGAATGGGTGATCAACGGGAAGAAAAACGGCAGGCCAGCGTCTTTTTTCAGACGCGGTAAGACAAAATGTCGAACACGCTGAACCATACATCCTCCCGGCAAATTGTCCGGCATTATACTCATCGTCAGAAg                                                                                 <  2:52955/139‑1 (MQ=255)                                                          tGGGAATGGGTGATCAACGGGAAGAAAAACGGCAGGCCAGCGTCTTTTTTCAGACGCGGTAAGACAAAATGTCGAACACGCTGAACCATACATCCTCCCGGCAAATTGTCCGGCATTATACTCATCGTCAGAAGCGgcgc                                                                           <  1:405243/140‑1 (MQ=255)                                                                  ggcgATCAACGGGAAGAAAAACGGCAGGCCAGCGTCTTTTTTCAGACGCGGTAAGACAAAATGTCGAACACGCTGAACCATACATCCTCCCGGCAAATTGTCCGGCATTATACTCATCGTCAGAAGCGGCGCAAGCAttt                                                                   >  2:356821/4‑140 (MQ=255)                                                                   gTGATCAACGGGAAGAAAAACGGCAGGCCAGCGTCTTTTTTCAGACGCGGTAAGACAAAATGTCGAACACGCTGAACCATACATCCTCCCGGCAAATTGTCCGGCATTATACTCATCGTCAGAAGCGGCGCAAGCATTTg                                                                  >  2:538875/1‑140 (MQ=255)                                                                       tCAACGGGAAGAAAAACGGCAGGCCAGCGTCTTTTTTCAGACGCGGTAAGACAAAATGTCGAACACGCTGAACCATACATCCTCCCGGCAAATTGTCCGGCATTATACTCATCGTCAGAAGCGGCGCAAGCATTTGAtac                                                              <  2:374121/140‑3 (MQ=255)                                                                       tCAACGGGAAGAAAAACGGCAGGCCAGCGTCTTTTTTCAGACGCGGTAAGACAAAATGTCGAACACGCTGAACCATACATCCTCCCGGCAAATTGTCCGGCATTATACTCATCGTCAGAAGCGGCGCAAGCATTTGATGc                                                              <  2:414267/140‑1 (MQ=255)                                                                       tCAACGGGAAGAAAAACGGCAGGCCAGCGTCTTTTTTCAGACGCGGTAAGACAAAATGTCGAACACGCTGAACCATACATCCTCCCGGCAAATTGTCCGGCATTATACTCATCGTCAGAAGCGGCGCAAGCATTTGATGc                                                              <  2:374033/140‑1 (MQ=255)                                                                        cAACGGGAAGAAAAACGGCAGGCCAGCGTCTTTTTTCAGACGCGGTAAGACAAAATGTCGAACACGCTGAACCATACATCCTCCCGGCAAATTGTCCGGCATTATACTCATCGTCAGAAGCGGCGCAAGCATTTGATGCa                                                             >  1:413793/1‑140 (MQ=255)                                                                                         gCAGGCCAGCGTCTTCTTTCAGCCGGGGTAAGACAAAATGTCGAACACGCTGAACCATACATCCTCCCGGCAAATTGTCCGGCATTATACTCATCGTCAGAAGCGGCGCAAGCATTTGATGCAATATTTTGTCAGCGCaa                                            <  1:320390/140‑1 (MQ=255)                                                                                            ggCCAGCGTCTTTTTTCAGACGCGGTAAGACAAAATGTCGAACACGCTGAACCATACATCCTCCCGGCAAATTGTCCGGCATTATACTCATGGTCAGAAGGGGCTCAAGCATTTGATACAATCTTTTGTCAGCGCAAACg                                         >  2:320390/1‑140 (MQ=255)                                                                                            ggCCAGCGTCTTTTTTCAGACGCGGTAAGACAAAATGTCGAACACGCTGAACCATACATCCTCCCGGCAAATTGTCCGGCATTATACTCATCGTCAGAAGCGGCGCAAGCATTTGATGCAATATTTTGTCAGCGCAAACg                                         >  1:439177/1‑140 (MQ=255)                                                                                             gCCAGCGTCTTTTTTCAGACGCGGTAAGACAAAATGTCGAACACGCTGAACCATACATCCTCCCGGCAAATTGTCCGGCATTATACTCATCGTCAGAAGCGGCGCAAGCATTTGATGCAATATTTTGTCAGCGCAAACgg                                        >  2:525977/1‑140 (MQ=255)                                                                                                       tttttCAGACGCGGTAAGACAAAATGTCGAACACGCTGAACCATACATCCTCCCGGCAAATTGTCCGGCATTATACTCATCGTCAGAAGCGGCGCAAGCATTTGATGCAATATTTTGTCAGCGCAAACGGTTTATTTGaa                              >  2:580479/1‑140 (MQ=255)                                                                                                       tttttCAGACGCGGTAAGACAAAATGTCGAACACGCTGAACCATACATCCTCCCGGCAAATTGTCCGGCATTATACTCATCGTCAGAAGCGGCGCAAGCATTTGATGCAATATTTTGTCAGCGCAAACGGCTTATTTGaa                              >  1:650288/1‑140 (MQ=255)                                                                                                         tttCAGACGCGGTAAGACAAAATGTCGAACACGCTGAACCATACATCCTCCCGGCAAATTGTCCGGCATTATACTCATCGTCAGAAGCGGCGCAAGCATTTGATGCAATATTTTGTCAGCGCAAACGGTTTATTTGAAtt                            <  2:106618/140‑1 (MQ=255)                                                                                                           tCAGACGCGGTAAGACAAAATGTCGAACACGCTGAACCATACATCCTCCCGGCAAATTGTCCGGCATTATACTCATCGTCAGAAGCGGCGCAAGCATTTGATGCAATATTTTGTCAGCGCAAACGGTTTATTTGAATTaa                          <  2:48251/140‑1 (MQ=255)                                                                                                               aCGCGGTAAGACAAAATGTCGAACACGCTGAACCATACATGCTCCCGGCAAATTGTCCGGCATTATACTCATCGTCAGAAGCGGCGCAAGCATTTGATGCAATATTTTGTCAGCGCAAACGGTTTATTTGAATTAAAAGt                      >  2:421032/1‑140 (MQ=255)                                                                                                               aCGCGGTAAGACAAAATGTCGAACACGCTGAACCATACATCCTCCCGGCAAATTGTCCGGCATTATACTCATCGTCAGAAGCGGCGCAAGCATTTGATGCAATATTTTGTCAGCGCAAACGGTTTATTTGAATTAAAAGt                      <  1:356821/140‑1 (MQ=255)                                                                                                                 gcgGTAAGACAAAATGTTGAACTCGCTGAACCATACATCCTCCCGGCAGATTGTCCGGCATTATACTCATCGTCAGAAGCGGCGCAAGCATTTGATGCAATATTTTGTCAGCGCAAACGGTTTATTTGAATTAAAAGTCa                    <  2:195805/140‑1 (MQ=255)                                                                                                                 gcgGTAAGACAAAATGTCGAACACGCTGAACCATACATCCTCCCGGCAAATTGTCCGGCATTATACTCATCGTCAGAAGCGGCGCAAGCATTTGATGCAATATTTTGTCAGCGCAAACGGTTTATTTGAATTAAAAGTCa                    <  2:650288/140‑1 (MQ=255)                                                                                                                 gcgGTAAGACAAAATGTCGAACACGCTGAACCATACATCCTCCCGGCAAATTGTCCGGCATTATACTCATCGTCAGAAGCGGCGCAAGCATTTGATGCAATATTTTGTCAGCGCAAACGGTTTATTTGAATTAAAAGTCa                    <  2:573500/140‑1 (MQ=255)                                                                                                                    gTAAGACAAAATGTCGAACACGCTGAACCATACATCCTCCCGGCAAATTGTCCGGCATTATACTCATCGTCAGAAGCGGCGCAAGCATTTGATGCAATATTTTGTCAGCGCAAACGGTTTATTTGAATTAAAAGTCAAgg                 <  2:499279/140‑1 (MQ=255)                                                                                                                       aGACAAAATGTCGAACACGCTGAACCATACATCCTCCCGGCAAATTGTCCGGCATTATACTCATCGTCAGAAGCGGCGCAAGCATTTGATGCAATATTTTGTCAGCGCAAACGGTTTATTTGAATTAAAAGTCAAGGtat              <  1:27421/140‑1 (MQ=255)                                                                                                                          cAAAATGTCGAACACGCTGAACCATACATCCTCCCGGCAAATTGTCCGGCATTATACTCATCGTCAGAAGCGGCGCAAGCATTTGATGCAATATTTTGTCAGCGCAAACGGTTTATTTGAATTAAAAGTCAAGGTATATg           <  1:148796/140‑1 (MQ=255)                                                                                                                             aaTGTCGAACACGCTGAACCATACATCCTCCCGGCAAATTGTCCGGCATTATACTCATCGTCAGAAGCGGCGCAAGCATTTGATGCAATATTTTGTCAGCGCAAACGGTTTATTTGAATTAAAAGTCa                    <  1:85612/128‑1 (MQ=255)                                                                                                                             aaTGTCGAACACGCTGAACCATACATCCTCCCGGCAAATTGTCCGGCATTATACTCATCGTCAGAAGCGGCGCAAGCATTTGATGCAATATTTTGTCAGCGCAAACGGTTTATTTGAATTAAAAGTCa                    >  2:85612/1‑128 (MQ=255)                                                                                                                             aaTGTCGAACACGCTGAACCATACATCCTCCCGGCAAATTGTCCGGCATTATACTCATCGTCAGAAGCGGCGCAAGCATTTGATGCAATATTTTGTCAGCGCAAACGGTTTATTTGAATTAAAAGTCAAGGTATATGCAt        >  2:243011/1‑140 (MQ=255)                                                                                                                             aaTGTCGAACACGCTGAACCATACATCCTCCCGGCAAATTGTCCGGCATTATACTCATCGTCAGAAGCGGCGCAAGCATTTGATGCAATATTTTGTCAGCACAACCGGCTTATTTGATTTAAGAGTCAAGGTATATGCAt        >  2:205778/1‑140 (MQ=255)                                                                                                                                  cGAACACGCTGAACCATACATCCTCCCGGCAAATTGTCCGGCATTATACTCATCGTCAGAAGCGGCGCAAGCATTTGATGCAATATTTTGTCAGCGCAAACGGTTTATTTGAATTAAAAGTCAAGGTATATGCATTTTTa   >  1:167425/1‑140 (MQ=255)                                                                                                                                   gAACACGCTGAACCATACATCCTCCCGGCAAATTGTCCGGCATTATACTCATCGTCAGAAGCGGCGCAAGCATTTGATGCAATATTTTGTCAGCGCAAACGGTTTATTTGAATTAAAAGTCAAGGTATATGCATTTTTAt  >  1:357211/1‑140 (MQ=255)                                                                                                                                         |                                                                                                                                      TTAGCCATCTTTTATCTCCTGACGCCTGGGCAAAAAAAAGCCCCTCGATTGAGGGGCTGGGAATGGGTGATCAACGGGAAGAAAAACGGCAGGCCAGCGTCTTTTTTCAGACGCGGTAAGACAAAATGTCGAACACACTGAACCATACATCCTCCCGGCAAATTGTCCGGCATTATACTCATCGTCAGAAGCGGCGCAAGCATTTGATGCAATATTTTGTCAGCGCAAACGGTTTATTTGAATTAAAAGTCAAGGTATATGCATTTTTAT  >  NC_000913/4472388‑4472657

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