Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,026,626	T→C	R287R (CGT→CGC)	guaD →	guanine deaminase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,026,626	0	T	C	100.0%	73.8 / NA	22	R287R (CGT→CGC)	guaD	guanine deaminase
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (14/8);  total (14/8)

AATCGCTTTATCCTGACCATGATGGTTATCTGGATGTTTACCATCAGTACGGCCTGACCGGTAAAAACTGTGTCTTTGCTCACTGCGTCCATCTCGAAGAAAAAGAGTGGGATCGTCTCAGCGAAACCAAATCCAGCATTGCTTTCTGTCCGACCTCCAACCTTTACCTCGGCAGCGGCTTATTCAACTTGAAAAAAGCATGGCAGAAGAAAGTTAAAGTGGGCATGGGAACGGATATCGGT  >  NC_000913/3026511‑3026752                                                                                                                    |                                                                                                                               aaTCGCTTTATCCTGACCATGATGGTTATCTGGATGTTTACCATCAGTACGGCCTGACCGGTAAAAACTGTGTCTTTGCTCACTGCGTCCATCTCGAAGAAAAAGAGTGGGATCGCCTCAGCGAAACCAAATc                                                                                                               <  1:56136/133‑1 (MQ=255)               gACCATGATGGTTATCTGGATGTTTACCATCAGTACGGCCTGACCGGTAAAAACTGTGTCTTTGCTCACTGCGTCCATCTCGAAGAAAAAGAGTGGGATCGCCTCAGCGAAACCAAATCCAGCATTGCTTTCt                                                                                                 >  1:337523/1‑133 (MQ=255)                  cATGATGGTTCTCTGAATGTTTACTATCTGTACGGCCTGACCGGTAAAAACTGTGTCTTTGCTCACTGCGTCCATCTCGAAGAAAAAGAGTGGGATCGCCTCAGCGAAACCAAATCCAGCATTGCTTTCTGTc                                                                                              <  2:153904/133‑1 (MQ=255)                  cATGATGGTTATCTGGATGTTTACCATCAGTACGGCCTGACCGGTAAAAACTGTGTCTTTGCTCACTGCGTCCATCTCGAAGAAAAAGAGTGGGATCACCTCAGCGAAACCAAATCCAGCATTGCTTTCTGTc                                                                                              <  1:462294/133‑1 (MQ=255)                   atgatgGTTAACTGGCTGTTTACCCTCATTACCGCCTGACCTGTAAAAACTGTGTCTTTGCTCACTGCGTCCATCTCGAAGAAAAAGAGTGGGATCGCCTCAGCGAAACCAAATCCAGCATTGCTTTCTGTcc                                                                                             <  2:140175/133‑1 (MQ=255)                                       tACCATCAGTACGGCCTGACCGGTAAAAACTGTGTCTTTGCTCACTGCGTCCATCTCGAAGAAAAAGAGTGGGATCGCCTCAGCGAAACCAAATCCAGCATTGCTTTCTGTCCGACCTCCAACCTTTACCTCg                                                                         >  1:401412/1‑133 (MQ=255)                                        aCCATCAGTACGGCCTGACCGGTAAAAACTGTGTCTTTGCTCACTGCGTCCATCTCGAAGAAAAAGAGTGGGATCGCCTCAGCGAAACCAAATCCAGCATTGATTTCTGTCCGACCTACAACCTGTACCTCgg                                                                        >  1:231993/1‑133 (MQ=255)                                          cATCAGTACGGCCTGACCGGTAAAAACTGTGTCTTTGCTCACTGCGTCCATCTCGAAGAAAAAGAGTGGGATCGCCTCAGCGAAACCAAATCCAGCATTGCTTTCTGTCCGACCTCCAACCTTTACCTCGGCa                                                                      <  2:6256/133‑1 (MQ=255)                                             cAGTACGGCCTGACCGGTAAAAACTGTGTCTTTGCTCACTGCGTCCATCTCGAAGAAAAAGAGTGGGATCGCCTCAGCGAAACCAAATCCAGCATTGCTTTCTGTCCGACCTCCAACCTTTACCTCGGCAGCg                                                                   >  1:497944/1‑133 (MQ=255)                                              aGTACGGCCTGACCGGTAAAAACTGTGTCTTTGCTCACTGCGTCCATCTCGAAGAAAAAGAGTGGGATCGCCTCAGCGAAACCAAATCCAGCATTGCTTTCTGTCCGACCTCCAACCTTTACCTCGGCAGCgg                                                                  <  2:497944/133‑1 (MQ=255)                                               gTACGGCCTGACCGGTAAAAACTGTGTCTTTGCTCACTGCGTCCATCTCGAAGAAAAAGAGTGGGATCGCCTCAGCGAAACCAAATCCAGCATTGCTTTCTGTCCGACCTCCAACCTTTACCTCGGCAGCGGc                                                                 >  1:384261/1‑133 (MQ=255)                                                    gCCTGACCGGTAAAAACTGTGTCTTTGCTCACTGCGTCCATCTCGAAGAAAAAGAGTGGGATCGCCTCAGCGAAACCAAATCCAGCATTGCTTTCTGTCCGACCTCCAACTTTTCCCCCGGCAGCGGCTTAtt                                                            >  1:252037/1‑133 (MQ=255)                                                        gACCGGTAAAAACTGTGTCTTTGCTCACTGCGTCCATCTCGAAGAAAAAGAGTGGGATCGCCTCAGCGAAACCAAATCCAGCATTGCTTTCGGTCCGACCGCCAGCCTGTACCTCGGGAGCGTCTTATTCAAc                                                        >  2:132307/1‑133 (MQ=255)                                                         aCCGGTAAAAACTGTGTCTTTGCTCACTGCGTCCATCTCGAAGAAAAAGAGTGGGATCGCCTCAGCGAAACCAAATCCAGCATTGCTTTCTGTCCGACCTCCAACCTTTACCTCGGCAGCGGCTTATTCAACt                                                       >  1:21692/1‑133 (MQ=255)                                                                               cTCACTGCGTCCATCTCGAAGAAAAAGAGTGGGATCGCCTCAGCGAAACCAAATCCAGCATTGCTTTCTGTCCGACCTCCAACCTTTACCTCGGCAGCGGCTGATTCAACGTGAGAAGAGCATGGCAgaagaa                                 >  1:149348/1‑133 (MQ=255)                                                                                 cACTGCGTCCATCTCGAAGAAAAAGAGTGGGATCGCCTCAGCGAAACCAAATCCAGCATTGCTTTCTGTCCGACCTCCAACCTTTACCTCGGCAGCGGCTTATTCAACTTGAAAAAAGCATGGCAGAAGAAAg                               >  1:329000/1‑133 (MQ=255)                                                                                           aTCTCGAAGAAAAAGAGTGGGATCGCCTCAGCGAAACCAAATCCAGCATTGCTTTCTGTCCGACCTCCAACCTTTACCTCGGCAGCGGCTTATTCAACTTGAAAAAAGCATGGCAGAAGAAAGTTAAAGTggg                     <  2:337523/133‑1 (MQ=255)                                                                                              tcGAAGAAAAAGAGTGGGATCGCCTCAGCGAAACCAAATCCAGCATTGCTTTCTGTCCGACCTCCAACCTTTACCTCGGCAGCGGCTTATTCAACTTGAAAAAAGCATGGCAGAAGAAAGTTAAAGTGGGCAt                  >  2:301243/1‑133 (MQ=255)                                                                                                       aaGAGTGGGATCGCCTCAGCGAAACCAAATCCAGCATTGCTTTCTGTCCGACCTCCAACCTTTACCTCGGCAGCGGCTTATTCAACTTGAAAAAAGCATGGCAGAAGAAAGTTAAAGTGGGCATGGGAACGGa         >  2:48200/1‑133 (MQ=255)                                                                                                        agagTGGGATCGCCTCAGCGAAACCAAATCCAGCATTGCTTTCTGTCCGACCTCCAACCTTTACCTCGGCAGCGGCTTATTCAACTTGAAAAAAGCATGGCAGAAGAAAGTTAAAGTGGGCATGGGAACGGat        >  2:100485/1‑133 (MQ=255)                                                                                                            tGGGATCGCCTCAGCGAAACCAAATCCAGCATTGCTTTCTGTCCGACCTCCAACCTTTACCTCGGCAGCGGCTTATTCAACTTGAAAAAAGCATGGCAGAAGAAAGTTAAGGTGAGCATGGGAATGGATATCg    >  1:152866/1‑133 (MQ=255)                                                                                                              ggATCGCCTCAGCGAAACCAAATCCAGCATTGCTTTCTGTCCGACCTCCAACCTTTACCTCGGCAGCGGCTTATTCAACTTGAAAAAAGCATTGCAGAAGAAAGTTAAAGTGGGCATGGGAACGGATATCGGt  <  2:384261/133‑1 (MQ=255)                                                                                                                    |                                                                                                                               AATCGCTTTATCCTGACCATGATGGTTATCTGGATGTTTACCATCAGTACGGCCTGACCGGTAAAAACTGTGTCTTTGCTCACTGCGTCCATCTCGAAGAAAAAGAGTGGGATCGTCTCAGCGAAACCAAATCCAGCATTGCTTTCTGTCCGACCTCCAACCTTTACCTCGGCAGCGGCTTATTCAACTTGAAAAAAGCATGGCAGAAGAAAGTTAAAGTGGGCATGGGAACGGATATCGGT  >  NC_000913/3026511‑3026752

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