Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	433,896	T→C	Y148H (TAT→CAT)	ribD →	fused diaminohydroxyphosphoribosylaminopyrimidine deaminase/5‑amino‑6‑(5‑phosphoribosylamino)uracil reductase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	433,896	0	T	C	100.0%	90.1 / NA	25	Y148H (TAT→CAT)	ribD	fused diaminohydroxyphosphoribosylaminopyrimidine deaminase/5‑amino‑6‑(5‑phosphoribosylamino)uracil reductase
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (11/14);  total (11/14)

CCTATGTCACACTCGAACCCTGTAGCCATCATGGTCGTACGCCACCGTGCTGTGACGCACTCATCGCCGCTGGCGTAGCGCGCGTGGTTGCCTCGATGCAAGATCCTAACCCGCAGGTCGCTGGGCGTGGACTTTACCGTCTGCAACAGGCTGGCATTGACGTCAGCCACGGCCTGATGATGAGTGAAGCCGAGCAATTGAATAAAGGCTTTCTCAAGCGGATGCGCACCGGCTTTCCTTATATTCAGTTAAAACTTGGCGCATCGCTTGATGGTCGCACGGCGATGGCGAGCGGCGAAAGCCAGTGGATCACTTCGCCCCAGGCGCGGCGTGATGTACAACTACTGCGCGCGCAAAGTCATGCCATTTTAACCAGCAGCGCCACGGTGCTGGCGGATGATCCTGCCTTAACGGTGCGTTGGTCTGAACTGGATGAACAAACTCAGGCGCTCTATCCGCAACAAAATCTCCGTCAGC  >  NC_000913/433657‑434133                                                                                                                                                                                                                                                |                                                                                                                                                                                                                                              ccTATGTCACACGCGAACCCTGTAGCCATCATGGTCGTACGCCACCGTGCTGTGACGCACTCATCGCCGCTGGCGTAGCGCGCGTGGTTGCCTCGATGCAAGAACCTAACCCGCAGGTCGCTGGGCGTGGACTTTACCGTCTGCAACAGGCTGGCATTGACGTCAGCCACGGCCTGATGATGAGTGAAGCCGAGCAATTGAATAAAGGCTTGCTCAAGCGGATGCGCACCGGCTTTCCTCATAGTCAGTTa                                                                                                                                                                                                                                    <  2:189686/251‑1 (MQ=255)   tATGTCACACTCGAACCCTGTAGCCATCATGGTCGTACGCCACCGTGCTGTGACGCACTCATCGCCGCTGGCGTAGCGCGCGTGGTTGCCTCGATGCAAGATCCTAACCCGCAGGTCGCTGGGCGTGGACTTTACCGTCTGCAACAGGCTGGCATTGACGTCAGCCACGGCCTGATGATGAGTGAAGCCGAGCAATTGAATAAAGGCTTTCTCAAGCGGATGCGCACCGGCTTTCCTCATATTCAGTTaa                                                                                                                                                                                                                                   >  2:110255/1‑250 (MQ=255)                   ccTGTAGCCATCATGGTCGTACGCCACCGTGCTGTGACGCACTCATCGCCGCTGGCGTAGCGCGCGTGGTTGCCTCGATGCAAGATCCTAACCCGCAGGTCGCTGGGCGTGGACTTTACCGTCTGCAACAGGCTGGCATTGACGTCAGCCACGGCCTGATGATGAGTGAAGCCGAGCAATTGAATAAAGGCTTTCTCAAGCGGATGCGCACCGGCTTTCCTCATATTCAGTTAAAACTTGGCGCATCGCtt                                                                                                                                                                                                                  >  1:68753/1‑251 (MQ=255)                                 ggTCGTACGCCACCGTGCTGTGACGCACTCATCGCCGCTGGCGTAGCGCGCGTGGTTGCCTCGATGCAAGATCCTAACCCGCAGGTCGCTGGGCGTGGACTTTACCGTCTGCAACAGGCTGGCATTGACGTCAGCCACGGCCTGATGATGAGTGAAGCCGAGCAATTGAATAAAGGCGTTCTCAAGCGGATGCGCACCGGCTTTCCTCATATTCAGTTAAAACTTGGCGCATCGCTTGATGGTCGCACGGc                                                                                                                                                                                                    <  1:117288/251‑1 (MQ=255)                                                            cTCATCGCCGCTGGCGTAGCGCGCGTGGTTGCCTCGATGCAAGATCCTAACCCGCAGGTCGCTGGGCGTGGACTTTACCGTCTGCAACAGGCTGGCATTGACGTCAGCCACGGCCTGATGATGAGTGAAGCCGAGCAATTGAGTAAAGGCTTTCTCAAGCGGATGCGCACCGGCTTTCCTCATATTCAGTTAAAACTTGTCGCATCGCTTGATGGTCGCACGGCGATGGCGAGCGGCGAAAGCCAGTGGAt                                                                                                                                                                         >  2:99072/1‑251 (MQ=255)                                                                    cgcTGGCGTAGCGCGCGTGGTTGCCTCGATGCAAGATCCTAACCCGCAGGTCGCTGGGCGTGGACTTTACCGTCTGCAACAGGCTGGCATTGACGTCAGCCACGGCCTGATGATGAGTGAAGCCGAGCAATTGAATAAAGGCTTTCTCAAGCGGATGCGCACCGGCTTTCCTCATATTCAGTTAAAACTTGGCGCATCGCTTGATGGTCGCACGGCGATGGCGAGCGGCGAAAGCCAGTGGATCACTTCGc                                                                                                                                                                 >  1:46105/1‑251 (MQ=255)                                                                    cgcTGGCGTAGCGCGCGTGGTTGCCTCGATGCAAGATCCTAACCCGCAGGTCGCTGGGCGTGGACTTTACCGTCTGCAACAGGCTGGCATTGACGTCAGCCACGGCCTGATGATGAGTGAAGCCGAGCAATTGAATAAAGGCTTTCTCAAGCGGATGCGCACCGGCTTTCCTCATATTCAGTTAAAACTTGGCGCATCGCTTGATGGTCGCACGGCGATGGCGAGCGGCGAAAGCCAGGGGATCACTTCGc                                                                                                                                                                 >  2:191595/1‑251 (MQ=255)                                                                     gcTGGCGTAGCGCGCGTGGTTGCCTCGATGCAAGATCCTAACCCGCAGGTCGCTGGGCGTGGACTTTACCGTCTGCAACAGGCTGGCATTGACGTCAGCCACGGCCTGATGATGAGTGAAGCCGAGCAATTGAATAAAGGCTTTCTCAAGCGGATGCGCACCGGCTTTCCTCATATTCAGTTAAAACTTGGCGCATCGCTTGATGGTCGCACGGCGATGGCGAGCGGCGAAAGCCAGTGGATCACTTCGcc                                                                                                                                                                <  1:223538/251‑1 (MQ=255)                                                                     gcTGGCGTAGCGCGCGTGGTTGCCTCGATGCAAGATCCTAACCCGCAGGTCGCTGGGCGTGGACTTTACCGTCTGCAACAGGCTGGCATTGACGTCAGCCACGGCCTGATGATGAGTGAAGCCGAGCAATTGAATAAAGGCTTTCTCAAGCGGATGCGCACCGGCTTTCCTCATATTCAGTTAAAACTTGGCGCATCGCTTGATGGTCGCACGGCGATGGCGAGCGGCGAAAGCCAGTGGATCACTTCGcc                                                                                                                                                                >  2:181610/1‑251 (MQ=255)                                                                         gCGTAGCGCGCGTGGTTGCCTCGATGCAAGATCCTAACCCGCAGGTCGCTGGGCGTGGACTTTACCGTCTGCAACAGGCTGGCATTGACGTCAGCCACGGCCTGATGATGAGTGAAGCCGAGCAATTGAATAAAGGCTTTCTCAAGCGGATGCGCACCGGCTTTCCTCATATTCAGTTAAAACTTGGCGCATCGCTTGATGGTCGCACGGCGATGGCGAGCGGCGAAAGCCAGTGGATCACTTCGCCCCAg                                                                                                                                                            <  2:73542/251‑1 (MQ=255)                                                                                                     aGATCCTAACCCGCAGGTCGCTGGGCGTGGACTTTACCGTCTGCAACAGGCTGGCATTGACGTCAGCCACGGCCTGATGATGAGTGAAGCCGAGCAATTGAATAAAGGCTTTCTCAAGCGGATGCGCACCGGCTTTCCTCATATTCAGTTAAAACTTGGCGCATCGCTTGATGGTCGCACGGCGATGTCGAGCGGCGAACGCCAGTGGATCACTTCGCCCCAGGCGCGGCGGGATGTACAACTACTgcgcg                                                                                                                                >  2:183870/1‑251 (MQ=255)                                                                                                                cGCAGGTCGCTGGGCGTGGACTTTACCGTCTGCAACAGGCTGGCATTGACGTCAGCCACGGCCTGATGATGAGTGAAGCCGAGCAATTGAATAAAGGCTTTCTCAAGCGGATGCGCACCGGCTTTCCTCATATTCAGTTAAAACTTGGCGCATCGCTTGATGGTCGCACGGCGATGGCGAGCGGCGAAAGCCAGTGGATCACTTCGCCCCAGGCGCGGCGTGATGTACAACTACTGCGCGCGCAAAGTCAt                                                                                                                     <  2:205133/251‑1 (MQ=255)                                                                                                                                               gCAACAGGCTGGCATTGACGTCAGCCACGGCCTGATGATGAGTGAAGCCGAGCAATTGAATAAAGGCTTTCTCAAGCGGATGCGCACCGGCTTTCCTCATATTCAGTTAAAACTTGGCGCATCGCTTGATGGTCGCACGGCGATGGCGAGCGGCGAAAGCCAGTGGATCACTTCGCCCCAGGCGCGGCGTGATGTACAACTACTGCGCGCGCAAAGTCATGCCATTTTAACCAGCAGCGCCACGGTGCt                                                                                        <  2:155300/249‑1 (MQ=255)                                                                                                                                                                aCGTCAGCCACGTCCTGATGAGGAGTGAAGCCGAGCAATTGAATAAAGGCTTTCTCAAGCGGATGCGCACCGGCTTTCCTCATATTCAGTTAAAACTTGGCGCATCGCTTGATGGTCTCACGGCGATGGCGAGCGGCGAAAGCCAGTGGATCACTTCGCCCCAGGCGCGGCGTGATGTACAACTACTGCGCGCGCAAAGTCATGCCATTTTAACCAGCAGCGCCACGGTGCTGGCGGATGATCCTGCCTTa                                                                     <  2:176121/251‑1 (MQ=255)                                                                                                                                                                        cACGGCCTGATGATGAGTGAAGCCGAGCAATTGAATAAAGGCTTTCTCAAGCGGATGCGCACCGGCTTTCCTCATATTCAGTTAAAACTTGGCGCATCGCTTGATGGTCGCACGGCGATGGCGAGCGGCGAAAGCCAGTGGATCACTTCGCCCCAGGCGCGGCGTGATGTACAACTACTGCGCGCGCAAAGTCATGCCATTTTAACCAGCAGCGCCACGGTGCTGGCGGATGATCCTGCCTTAACGGTGCg                                                             <  1:235045/251‑1 (MQ=255)                                                                                                                                                                          cGGCCTGATGATGAGTGAAGCCGAGCAATTGAATAAAGGCTTTCTCAAGCGGATGCGCACCGGCTTTCCTCATATTCAGTTAAAACTTGGCGCATCGCTTGATGGTCGCACGGCGATGGCGAGCGGCGAAAGCCAGTGGATCACTTCGCCCCAGGCGCGGCGTGATGTACAACTACTGCGCGCGCAAAGTCATGCCATTTTAACCAGCAGCGCCACGGTGCTGGCGGATGATCCTGCCTTAACGGTGCGtt                                                           >  1:95357/1‑251 (MQ=255)                                                                                                                                                                                          gAAGCCGATCAATTGAATAAAGGCTTTCTCAAGCGGATGCGCACCGGCTTTCCGCAGATTCAGTTAAAACTTGGCGCATCGCTTGATGGTCGCACGGCGATGGCGAGCGGCGAAAGCCAGTGGATCACTTCGCCCCAGGCGCGGCGTGATGTACAACTACTGCGCGCGCAAAGTCATGCCATTTTAACCAGCAGCGCCACGGTGCTGGCGGATGATCCTGCCTTAACGGTGCGTTGGTCTGAACTGGATgc                                           <  2:203966/251‑2 (MQ=255)                                                                                                                                                                                               cGAGCAATTGAATAAAGGCTTTCTCAAGCGGATGCGCACCGGCTTTCCTCATATTCAGTTAAAACTTGGCGCATCGCTTGATGGTCGCACGGCGATGGCGAGCGGCGAAAGCCAGTGGATCACTTCGCCCCAGGCGCGGCGTGATGTACAACTACTGCGCGCGCAAAGTCATGCCATTTTAACCAGCAGCGCCACGGTGCTGGCGGATGATCCTGCCTTAACGGTGCGTTGGTCTGAACTGGATGAACaa                                       >  1:191631/1‑250 (MQ=255)                                                                                                                                                                                                    aaTTGAATAAAGGCTTTCTCAAGCGGATGCGCACCGGCTTTCCTCATATTCAGTTAAAACTTGGCGCATCGCTTGATGGTCGCACGGCGATGGCGAGCGGCGAAAGCCAGTGGATCACTTCGCCCCAGGCGCGGCGTGATGTACAACTACTGCGCGCGCAAAGTCATGCCATTTTAACCAGCAGCGCCACGGTGCTGGCGGATGATCCTGCCTTAACGGTGCGTTGGTCTGAACTGGATGAACAAACTCAg                                 <  1:24280/251‑1 (MQ=255)                                                                                                                                                                                                                 cTTTCTCAAGCGGATGCGCACCGGCTTTCCTCATATTCAGTTAAAACTTGGCGCATCGCTTGATGGTCGCACGGCGATGGCGAGCGGCGAAAGCCAGTGGATCACTTCGCCCCAGGCGCGGCGTGATGTACAACTACTGCGCGCGCAAAGTCATGCCATTTTAACCAGCAGCGCCACGGTGCTGGCGGATGATCCTGCCTTAACGGTGCGTTGGTCTGAACTGGATGAACAAACTCAGGCGCTCTATCCGc                    >  1:163132/1‑251 (MQ=255)                                                                                                                                                                                                                   ttCTCAAGCGGATGCGCACCGGCTTTCCTCATATTCAGTTAAAACTTGGCGCATCGCTTGATGGTCGCACGGCGATGGCGAGCGGCGAAAGCCAGTGGATCACTTCGCCCCAGGCGCGGCGTGATGTACAACTACTGCGCGCGCAAAGTCATGCCATTTTAACCAGCAGCGCCACGGTGCTGGCGGATGATCCTGCCTTAACGGTGCGTTGGTCTGAACTGGATGAACAAACTCAGGCGCTCTATCCGc                    <  1:176921/249‑1 (MQ=255)                                                                                                                                                                                                                       cAAGCGGATGCGCACCGGCTTTCCTCATATTCAGTTAAAACTTGGCGCATCGCTTGATGGTCGCACGGCGATGGCGAGCGGCGAAAGCCAGTGGATCACTTCGCCCCAGGCGCGGCGTGATGTACAACTACTGCGCGCGCAAAGTCATGCCATTTTAACCAGCAGCGCCACGGTGCTGGCGGATGATCCTGCCTTAACGGTGCGTTGGTCTGAACTGGATGAACAAACTCAGGCGCTCTATCCGccacaa               >  2:204126/1‑250 (MQ=255)                                                                                                                                                                                                                                gcgcACCGGCTTTCCGCATATTCAGTTAAAACTTGGCGCATCGCTTGATGGTCGCACGGCGATGGCGAGCGGCGAAAGCCAGTGGATCACTTCGCCCCAGGCGCGGCGTGATGTACAACTACTGCGCGCGCAAAGTCATGCCATTTTAACCAGCAGCGCCACGGTGCTGGCGGATGATCCTGCCTTAACGGTGCGTTGGTCTGAACTGGATGAACAAACTCAGGCGCTCTATCCGCAACAAAATCTCCGTc     <  2:183339/251‑1 (MQ=255)                                                                                                                                                                                                                                 cgcACCGGCTTTCCTCATATTCAGTTAAAACTTGGCGCATCGCTTGATGGTCGCACGGCGATGGCGAGCGGCGAAAGCCAGTGGATCACTTCGCCCCAGGCGCGGCGTGATGTACAACTACTGCGCGCGCAAAGTCATGCCATTTTAACCAGCAGCGCCACGGTGCTGGCGGATGATCCTGCCTTAACGGTGCGTTGGTCTGAACTGGATGAACAAACTCAGGCGCTCTATCCGCAACAAAATCTCCGTCa    <  2:130380/251‑1 (MQ=255)                                                                                                                                                                                                                                   cgccGGCTTTCCTCATATTCAGTTAAACCTTGGCGCATCGCTTGATGGTCGCACGGCGATGGCGAGCGGCGAAAGCCAGTGGATCACTTCGCCCCAGGCGCGGCGTGATGTACAACTACTGCGCGCGCAAAGTCATGCCATTTTAACCAGCAGCGCCACGGTGCTGGCGGATGATCCTGCCTTAACGGTGCGTTGGTCTGAACTGGATGAACAAACTCAGGCGCTCTATCCGCAACAAAATCTCCGTCAGc  <  2:78062/249‑1 (MQ=255)                                                                                                                                                                                                                                                |                                                                                                                                                                                                                                              CCTATGTCACACTCGAACCCTGTAGCCATCATGGTCGTACGCCACCGTGCTGTGACGCACTCATCGCCGCTGGCGTAGCGCGCGTGGTTGCCTCGATGCAAGATCCTAACCCGCAGGTCGCTGGGCGTGGACTTTACCGTCTGCAACAGGCTGGCATTGACGTCAGCCACGGCCTGATGATGAGTGAAGCCGAGCAATTGAATAAAGGCTTTCTCAAGCGGATGCGCACCGGCTTTCCTTATATTCAGTTAAAACTTGGCGCATCGCTTGATGGTCGCACGGCGATGGCGAGCGGCGAAAGCCAGTGGATCACTTCGCCCCAGGCGCGGCGTGATGTACAACTACTGCGCGCGCAAAGTCATGCCATTTTAACCAGCAGCGCCACGGTGCTGGCGGATGATCCTGCCTTAACGGTGCGTTGGTCTGAACTGGATGAACAAACTCAGGCGCTCTATCCGCAACAAAATCTCCGTCAGC  >  NC_000913/433657‑434133

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