Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,146,079	C→A	V173F (GTT→TTT)	yceF ←	m(7)GTP pyrophosphatase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,146,079	0	C	A	100.0%	73.3 / NA	27	V173F (GTT→TTT)	yceF	m(7)GTP pyrophosphatase
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base A (16/11);  total (16/11)

GAAGCTCCGATGGATGAAGGTTTGAAGCGTTGTTTGCAAAAGCTGCGTAACGCGCGCTAATACAAGCATATAAACCTGATAATGGGCGGTTGGCCTGATAAGGCGTTTACGCCGCATCAGGCCGCCAGCACCGATTGCCGGATGCGACGTAACCACATCCGGCGCAAAAAACTATCCCATCAGCGGGTTTTTCCCTTCCCGACGTAACATCTGACACAGCGCGATTAGCGGTAAACCAACCAGCGTGTTAGGATCACGCCCCTCTAAACGCTCAAACAGCGTAATGCCAAATCCTTCACTCTTAAAGCTACCCGCGCAGTGCAGGGGATGCTCTTTACGCACGTAATTATCAATCTCCGCCTCGCTCAGGTGACGGAAATGGACGTCAAAAGGCTCCACTTCTGTTTGCAGATGCCCATTCGCCGAATTAAACAACGCCAGTCCAGTATAGAAGGTGACGATATTGCCGCTGG  >  NC_000913/1145840‑1146312                                                                                                                                                                                                                                                |                                                                                                                                                                                                                                          gcaGCTCCGGGGGATGAAGGTTTGAAGCGTTGTTTGCAAAAGCTGCGTATCGCGCGCTAATACAAGCATATAAACCTGATAATGGGCGGTTGGCTTGATAAGGCTTTTACGCCGCATCAGGCCGCCAGCACCGATTGCCGGATGCGACGTAACTACATCCGGCGCAAAAAACTATCCCATCAGCGGGTTTTTCCCTTCCCGACGTAACATCTGACACAGCGCGATTAGCGGTAAACCAAACAGCGTGTTAg                                                                                                                                                                                                                                <  2:71039/249‑1 (MQ=255)                 aaGGTTTGAAGCGTTGTTTGCAAAAGCTGCGTAACGCGCGCTAATACAAGCATATAAACCTGATAATGGGCGGTTGGCCTGATAAGGCGTTTACGCCGCATCAGGCCGCCAGCACCGATTGCCGGATGCGACGTAACCACATCCGGCGCAAAAAACTATCCCATCAGCGGGTTTTTCCCTTCCCGACGTAACATCTGACACAGCGCGATTAGCGGTAAACCAAACAGCGTGTTAGGATCACGCCCCTCTaa                                                                                                                                                                                                                <  1:259569/251‑1 (MQ=255)                 aaGGATTGAAGCGTTGTTTGCAAAAGCTGCGTAACGCGCGCTAATACAAGCATATAAACCTGATAATGGGCGGTTGGCCTGATAAGGCGTTTACGCCGCATCAGGCCGCCAGCACCGATTGCCGGATGCGACGTAACCACATCCGGCGCAAAAAACTATCCCATCAGCGGGTTTTTCCCTTCCCGACGTAACATCTGACACAGCGCGATTAGCGGTAAACCAAACAGCGTGTTAGGATCACGCCCCTCTa                                                                                                                                                                                                                 >  1:298681/1‑250 (MQ=255)                    gTTTGAAGCGTTGTTTGCAAAAGCTGCGTAACGCGCGCTAATACAAGCATATAAACCTGCTAATGGGCGGTTGGCCTGATAAGGCGTTTACGCCGCATCAGGCCGCCAGCACCGATTGCCGGATGCGACGTAACCACATCCGGCGCAAAAAACTATCCCATCAGCGGGGTTTTCCCTTCCCGACGTAACATCTGACACAGCGCGCTTAGCGGTAAACCAAACAGCGTGTTAGGATAACGCCCCTCTAAACg                                                                                                                                                                                                             >  1:214503/1‑251 (MQ=255)                        gAAGCGTTGTTTGCAAAAGCTGCGTAACGCGCGCTAATACAAGCATATAACCCTGATAATGGGCGGTTGGCCTGATAAGGCGTTTACGCCGCATCAGGCCGCCAGCACCGATTGCCGGATGCGACGTAACCACATCCGGCGCAAAAAACTATCCCATCAGCGGGTTTTTCCCTTCCCGACGTAACATCTGACACAGCGCGATTAGCGGTAAACCAAACAGCGTGTTAGGATCACGCCCCTCTAAACGCTCa                                                                                                                                                                                                         <  1:45596/251‑1 (MQ=255)                               tgttTGCAAAAGCTGCGTAACGCGCGCTAATACAAGCATATAAACCTGATAATGGGCGGTTGGCCTGATAAGGCGTTTACGCCGCATCAGGCCGCCAGCACCGATTGCCGGATGCGACGTAACCACATCCGGCGCAAAAAACTATCCCATCAGCGGGTTTTTCCCTTCCCGACGTAACATCTGACACAGCGCGATTAGCGGTAAACCAAACAGCGTGTTAGGATCACGCCCCTCTAAACGCTCAAACAGCg                                                                                                                                                                                                  >  1:163170/1‑251 (MQ=255)                                    gCAAAAGCTGCGTAACGCGCGCTAATAAAAGCATATAAACCTGATAATGGGCGGTTGGCCTGATAAGGCGTTTACGCCGCATCAGGCCGCCAGCACCGATTGCCGGATGCGACGTAACCACATCCGGCGCAAAAAACTATCCCATCAGCGGGTTTTTCCCTTCCCGACGTAACATCTGACACAGCGCGATTAGAGGTACACCAAACAGCGTGTTAGGATCACGCCCCTCTAAACGCTCAAACAGCGTactt                                                                                                                                                                                             >  1:239826/1‑248 (MQ=255)                                    ccAAAAGCTGCGTAACGCGCGCTAATACAAGCATATAAACCTGATAATGGGCGGTTGGCCTGATAAGGCGTTTACGCCGCCTCAGGCCGCCAGCACCGATTGCCGGATGCGACGTAACCACATCCGGCGCAAAAAACTCTCCCATCAGCGGGTTTTTCCCTTCCCGACGTAACATCTCACACAGCGCGATTAGCCGTAAACCAAACAGCGTGTTAGGATCACGCCCCTCTAAACGCTCAAACAGCGTAAtt                                                                                                                                                                                             >  1:239835/2‑250 (MQ=255)                                                 aaCGCGCGCTAATACAAGCATATAAACCTGATAATGGGCGGTTGGCCTGATAAGGCGTTTACGCCGCATCAGGCCGCCAGCACCGATTGCCGGATGCGACGTAACCACATCCGGCGCAAAAAACTATCCCATCAGCGGGTTTTTCCCTTCCCGACGTAACCACTGACAAAGCGCGATTAGCGGTAAACCAAACAGCGTGTTAGGATCACGCCCCTCTAAACGCTCACAAAGCGTAATGCCAAATCCTTCAc                                                                                                                                                                                >  1:333620/1‑251 (MQ=255)                                                   cgcgcgCTAATACAAGCATATAAACCTGATAATGGGCGGTTGGCCTGATAAGGCGTTTACGCCGCATCAGGCCGCCAGCACCGATTGCCGGATGCGACGTAACCACATCCGGCGCAAAAAACTATCCCATCAGCGGGTTTTTCCCTTCCCGACGTAACATCTGACACAGCGCGATTAGCGGTAAACCAAACAGCGTGTTAGGATCACGCCCCTCTAAACGCTCAAACAGCGTAATGCCAAATCATTCActc                                                                                                                                                                              >  1:209318/1‑251 (MQ=255)                                                                            cTGATAATGGGCGGTTGGCCTGATAAGGCGTTTACGCCGCATCAGGCCGCCAGCACCGATTGCCGGATGCGACGTAACCACATCCGGCGCAAAAAACTATCCCATCAGCGGGGTTTTCCCTTCCCGACGTAACATCTGACACAACGCGATTAGCGGTAAACCAAACAGCGTGTTAGGATCACGCCCCTCTAAACGCTCAAACAGCGTAATGCCAAATCCTTCACTCTTAAAGCTACCCGCGCAGTGCAggg                                                                                                                                                     >  1:93182/1‑251 (MQ=255)                                                                               atcaTGGGCGATTGGCCTGATACGGCGTCTACGCCGCATCAGGCCGCCAGCACCGATTGCCGGATGCGACGTAACCACATCCGGCGCACAAAACTATCCCATCAGCGGGTTTTTCCCTTCCCGAGGTAACATATGACCCAGCGCGATTAGCGGTACACCAAACAGCATGTTAAGATCAAGCCCCTCTAAACGCTCAAACAGCGTAATTGCAACTCCTTAACTCTTCAAGCTACCcgccccgttacggggat                                                                                                                                                  >  2:60662/4‑237 (MQ=255)                                                                               aTAATGGGCGGTTGGCCTGATAAGGCGTTTACGCCGCATCAGGCCGCCAGCACCGATTGCCGGATGCGACGTAACCACATCCGGCGCAAAAAACTATCCCATCAGCGGGTTTTTCCCTTCCCGACGTAACATCTGACACAGCGCGATTAGCGGTAAACCAAACAGCGTGTTAGGATCACGCCCCTCTAAACGCTCAAACAGCGTAATGCCAAATCCTTCACTCTTAAAGCTACCCGCGCAGTGCAGGGGAt                                                                                                                                                  >  1:77440/1‑251 (MQ=255)                                                                                        ggTTGGCCTGATAAGGCGTTTACGCCGCATCAGGCCGCCAGCACCGATTGCCGGATGCGACGTAACCACATCCGGCGCAAAAAACTATCCCATCAGCGGGTTTTTCCCTTCCCGACGTAACATCTGACACAGCGCGATTAGCGGTAAACCAAACAGCGTGTTAGGATCACGCCCCTCTAAACGCTCAAACAGCGTAATGCCAAATCCTTCACTCTTAAAGCTACCCGCGCAGTGCAGGGGATGCTCTTTAc                                                                                                                                         <  1:77893/251‑1 (MQ=255)                                                                                                    aaGGCTTTTACGCCGCACCTGGCCGCCAGCACCGATTGCCGGATGCGTCGTAACCACACCCGGCGCAAAAAACGATCCCATCAGCGGTTTTTTCCCTTCCCGACGTAACATCTGACACAGCGCTATTAGCGGTAAACCAAACAGCGTGTTAGGATCACGCCCCTCTAAACGCTCAAACAGCGTAATGCCAAATCCTTCACTCTTAAAGCTACCCGCGCAGTGCAGGGGATGCTCGTTACGCACGTAATTa                                                                                                                              <  2:219459/250‑1 (MQ=255)                                                                                                                 cacaccAGGCCGCCAGCGCGGCTTGCCGGATGCGACGCATCCCCAACCGGAGCAAAAAACTAATCCATCAGCGGGTTTTTCCCTTCCCGACGTAACATCGGACACAGCGCGATTAGCGGTAAACCAAACAGCGTGTGAGGATCACGCCCCGCTAAACGCTCAAACAGCGTAATGTCAAATCCTTCACTCTTAAAGCTACCCGCGCAGTGCAGGGGATGCTCTTTCCGCACGTAATTATCAATCTACGCCTc                                                                                                                <  2:118767/246‑1 (MQ=255)                                                                                                                                   ccGATTGCCGGATGCGACGTAACCACATCCGGCGCAAAAAACTATCCCATCAGCGGGTTTTTCCCTTCCCGACGTAACATCTGACACAGCGCGATTAGCGGTAAACCAAACAGCGTGTTAGGATCACGCCCCTCTAAACGCTCAAACAGCGTAATGCCAAATCCTTCACTCTTAAAGCTACCCGCGCAGTGCAGGGGATGCTCTTTACGCACGTAATTATCAATCTCCCCCTCGCTCAGGTGACGGAAATg                                                                                              >  2:298063/1‑251 (MQ=255)                                                                                                                                                                     caaaaaccttacccAACAGCGGTTTGTTCCCTTCCCGACGTATCATCTGCCACAGCGCGATTAGCGTTAAACCAAACGGCGTGTTAGGATCACGCACCTATAAACTCTCAAACCGCGTAATGCCAAATCCTTCGCTCTTAAAGCTACCCGCGCAGTGCAGGGGATGCTCTTTACGCACGTCATTATCATTCTCCGCCGCGCTCAGGTGTCGGAAATGTACGTCAAAAGTCTCCACTTCTGTTTGCAGATGc                                                            <  2:238253/240‑1 (MQ=255)                                                                                                                                                                       aaaaaCTATCCCATCAGCGGGTTTTTCCCTTCCCGACGTAACATCTGACACAGCGCGATTAGCGGTAAACCAAACAGCGTGTTAGGATCACGCCCCTCTAAACGCTCAAACAGCGTAATGCCAAATCCTTCACTCTTAAAGCTACCCGCGCAGTGCAGGGGATGCTCTTTACGCACGTAATTATCAATCTCCGCCTCGCTCAGGTGACGGAAATGGACGTCAAAAGGCTCCACTTCTGTTTGCAGATGccc                                                          <  1:219977/251‑1 (MQ=255)                                                                                                                                                                                     cAGCGGGATTTTCCCTTCCCGACGTAACATCTGACACAGCGCGATTAGCGGTAAACCAAACAGCGTGTTAGGATCACGCCCCTCTAAACGCTCAAACAGCGTAATGCCAAATCCTTCACTCTTAAAGCTACCCGCGCAGTGCAGGGGATGCTCTTTACGCACGTAATTATCAATCTCCGCCTCGCTCAGGTGACGGAAATGGACGTCAAAAGGCTCAACTTCTGTTTGCAGATGCCCATTCGCCGAATTa                                             >  1:219622/1‑250 (MQ=255)                                                                                                                                                                                          ggTTTTTCCCTTCCCGACGTAACATCTGACACAGCGCGATTAGCGGTAAACCAAACAGCGTGTTAGGATCACGCCCCTCTAAACGCTCAAACAGCGTAATGCCAAATCCTTCACTCTTAAAGCTACCCGCGCAGTGCAGGGGATGCTCTTTACGCACGTAATTATCAATCTCCGCCTCGCTCAGGTGACGGAAATGGACGTCAAAAGGCTCCACTTCTGTTTGCAGATGCCCATTCGCCGAATTAaacaac                                       >  1:197466/1‑251 (MQ=255)                                                                                                                                                                                                      cccGACGTAACATCTGACACAGCGCGATTAGCGGTAAACCAAACAGCGTGTTAGGATCACGCCCCTCTAAACGCTCAAACAGCGTAATGCCAAATCCTTCACTCTTAAAGCt                                                                                                                                                                      <  1:7307/112‑1 (MQ=255)                                                                                                                                                                                                      cccGACGTAACATCTGACACAGCGCGATTAGCGGTAAACCAAACAGCGTGTTAGGATCACGCCCCTCTAAACGCTCAAACAGCGTAATGCCAAATCCTTCACTCTTAAAGCt                                                                                                                                                                      >  2:7307/1‑112 (MQ=255)                                                                                                                                                                                                                cATCTGACACAGCGCGATTAGCGGTAAACCAAACAGCGTGTTAGGATCACGCCCCTCTAAACGCTCAAACAGCGTAATGCCAAATCCTTCACTCTTAAAGCTACCCGCGCAGTGCAGGGGATGCTCTTTACGCACGTAATTATCAATCTCCGCCTCGCTCAGGTGACGGAAATGGACGTCAAAAGGCTCCACTTCTGTTTGCAGATGCCCATTCGCCGAATTAAACAACGCCAGTCCAGTATAGAAGGTg                  >  2:115519/1‑250 (MQ=255)                                                                                                                                                                                                                     gACACAGCGCGATGAACGGTAAACCAAACAGCGTGCTAGGAACCCGCCCCTCTAAACGCTCAAACAGCGTAATGCCAAATCCTTCACTCTTAAAGCTTCCCGCGCAGTGCAGGGGATGCTCTTTACGCACGTAATTATCAATCTCCGCCTCGCTCAGGTGACGGAAATGGACGTCAAAAGGCTCCACTTCTGTTTGCAGATGCCCATTCGCCGAATTAAACAACGCCAGTCCAGTATAGAAGGTGACGat             <  2:77440/250‑1 (MQ=255)                                                                                                                                                                                                                              cgATTAGCGGTAAACCAAACAGCGTGTTAGGATCACGCCCCTCTAAACGCTCAAACAGCGTAATGCCAAATCCTTCACTCTTAAAGCTACCCGCGCAGTGCAGGGGATGCTCTTTACGCACGTAATTATCAATCTCCGCCTCGCTCAGGTGACGGAAATGGACGTCAAAAGGCTCCACTTCTGTTTGCAGATGCCCATTCGCCGAATTAAACAACGCCAGTCCAGTATAGAAGGTGACGATATTGCCGCTg   >  1:104040/1‑251 (MQ=255)                                                                                                                                                                                                                               gATTAGCGGTAAACCAAACAGCGTGTTAGGATCACGCCCCTCTAAACGCTCAAACAGCGTAATGCCAAATCCTTCACTCTTAAAGCTACCCGCGCAGTGCAGGGGATGCTCTTTACGCACGTAATTATCAATCTCCGCCTCGCTCAGGTGACGGAAATGGACGTCAAAAGGCTCCACTTCTGTTTGCAGATGCCCATTCGCCGAATTAAACAACGCCAGTCCAGTATAGAAGGTGACGATATTGCCGCTgg  <  1:115519/251‑1 (MQ=255)                                                                                                                                                                                                                                                |                                                                                                                                                                                                                                          GAAGCTCCGATGGATGAAGGTTTGAAGCGTTGTTTGCAAAAGCTGCGTAACGCGCGCTAATACAAGCATATAAACCTGATAATGGGCGGTTGGCCTGATAAGGCGTTTACGCCGCATCAGGCCGCCAGCACCGATTGCCGGATGCGACGTAACCACATCCGGCGCAAAAAACTATCCCATCAGCGGGTTTTTCCCTTCCCGACGTAACATCTGACACAGCGCGATTAGCGGTAAACCAACCAGCGTGTTAGGATCACGCCCCTCTAAACGCTCAAACAGCGTAATGCCAAATCCTTCACTCTTAAAGCTACCCGCGCAGTGCAGGGGATGCTCTTTACGCACGTAATTATCAATCTCCGCCTCGCTCAGGTGACGGAAATGGACGTCAAAAGGCTCCACTTCTGTTTGCAGATGCCCATTCGCCGAATTAAACAACGCCAGTCCAGTATAGAAGGTGACGATATTGCCGCTGG  >  NC_000913/1145840‑1146312

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