Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,932,077	T→G	intergenic (+18/+38)	purT → / ← eda	phosphoribosylglycinamide formyltransferase 2/KHG/KDPG aldolase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,932,077	0	T	G	100.0%	89.2 / NA	25	intergenic (+18/+38)	purT/eda	phosphoribosylglycinamide formyltransferase 2/KHG/KDPG aldolase
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base G (12/13);  total (12/13)

TTGTTGACGCCATTGAACGCGCGAAGCACGCCGCCGGACAGGTAAAAGTACAGGGTTAAACCCGGCAAAAAAACGCTACAAAAATGCCCGATCCTCGATCGGGCATTTTGACTTTTACAGCTTAGCGCCTTCTACAGCTTCACGCGCCAGCTT  >  NC_000913/1932001‑1932153                                                                             |                                                                             ttgttgACGCCATTGAACGCGCGAAGCACGCCGCCGGACAGGTAAAAGTACAGGGTTAAACCCGGCAAAAAAACGCGACa                                                                           >  1:319495/1‑80 (MQ=255)            aTTGAACGCGCGAAGCACGCCGCCGGACAGGTAAAAGTACAGGGTTAAACCCGGCAAAAAAACGCGACAAAAATGCCCGAt                                                               <  2:319468/81‑1 (MQ=255)                aaCGCGCGAAGCACGCCGCCGGACAGGTAAAAGTACAGGGTTAAACCCGGCAAAAAAACGCGACAAAAATGCCCGATCCTc                                                           <  2:319469/81‑1 (MQ=255)                   gcgcgAAGCACGCCGCCGGACAGGTAAAAGTACAGGGTTAAACCCGGCAAAAAAACGCGACAAAAATGCCCGATCCTCGAt                                                        <  1:319471/81‑1 (MQ=255)                   gcgcgAAGCACGCCGCCGGACAGGTAAAAGTACAGGGTTAAACCCGGCAAAAAAACGCGACAAAAATGCCCGATCCTCGAt                                                        <  1:319470/81‑1 (MQ=255)                     gcgAAGCACGCCGCCGGACAGGTAAAAGTACAGGGTTAAACCCGGCAAAAAAACGCGACAAAAATGCCCGATCCTCGATCg                                                      <  1:319472/81‑1 (MQ=255)                       gAAGCACGCCGCCGGACAGGTAAAAGTACAGGGTTAAACCCGGCAAAAAAACGCGACAAAAATGCCCGATCCTCGATCggg                                                    >  2:319496/1‑81 (MQ=255)                        aaGCACGCCGCCGGACAGGTAAAAGTACAGGGTTAAACCCGGCAAAAAAACGCGACAAAAATGCCCGATCCTCAATCggg                                                    >  2:319508/1‑80 (MQ=255)                           cACGCCGCCGGACAGGTAAAAGTACAGGGTTAAACCCGGCAAAAAAACGCGACAAAAATGCCCGATCCTCGATCGGGCAtt                                                >  2:319500/1‑81 (MQ=255)                            aCGCCGCCGGACAGGTAAAAGTACAGGGTTAGACCCGGCAAAAAAACGCGACAAAAATGCCCGATCCTCGATCGGGCAttt                                               <  2:319473/81‑1 (MQ=255)                              gccgccGGACAGGTAAAAGTACAGGGTTAAACCCGGCAAAAAAACGCGACAAAAATGCCCGATCCTCGATCGGGCATTTTg                                             <  1:319474/81‑1 (MQ=255)                                      aCAGGTAAAAGTACAGGGTTAAACCCGGCAAAAAAACGCGACAAAAATGCCCGATCCTCGATCGGGCATTTTGACTTTTa                                      >  2:319510/1‑80 (MQ=255)                                      aCAGGTAAAAGTACAGGGTTAAACCCGGCAAAAAAACGCGACAAAAATGCCCGATCCTCGATCGGGCATTTTGACTTTTa                                      >  1:319503/1‑80 (MQ=255)                                      aCAGGTAAAAGTACAGGGTTAAACCCGGCAAAAAAACGCGACAAAAATGCCCGATCCTCGATCGGGCATTTTGACTTTTa                                      >  1:319507/1‑80 (MQ=255)                                          gTAAAAGTACAGGGTTAAACCCGGCAAAAAAACGCGACAAAAATGCCCGATCCTCGATCGGGCATTTTGACTTTTACAGCt                                 <  2:319475/81‑1 (MQ=255)                                             aaaGTACAGGGTTAAACCCGGCAAAAAAACGCGACAAAAATGCCCGATCCTCGATCGGGCATTTTGACTTTTACAGCTTAg                              >  1:319505/1‑81 (MQ=255)                                                  aCAGGGTTAAACCCGGCAAAAAAACGCGACAAAAATGCCCGATCCTCGATCGGGCATTTTGACTTTTACAGCTTAGCGCCt                         >  2:319502/1‑81 (MQ=255)                                                         tAAACCCGGCAAAAAAACGCGACAAAAATGCCCGATCCTCGATCGGGCATTTTGACTTTTACAGCTTAGCGCCTTCTACAg                  >  2:319499/1‑81 (MQ=255)                                                           aaCCCGGCAAAAAAACGCGACAAAAATGCCCGATCCTCGATCGGGCATTTTGACTTTTACAGCTTAGCGCCTTCTACAGCt                <  1:319476/81‑1 (MQ=255)                                                             cccGGCAAAAAAACGCGACAAAAATGCCCGATCCTCGATCGGGCATTTTGACTTTTACAGCTTAGCGCCTTCTACAGCCTc              <  2:319477/81‑1 (MQ=255)                                                               cGGCAAAAAAACGCGACAAAAATGCCCGATCCTCGATCGGGCATTTTGACTTTTACAGCTTAGCGCCTTCTACAGCTTCAc            <  1:319478/81‑1 (MQ=255)                                                                    aaaaaaCGCGACAAAAATGCCCGATCCTCGATCGGGCATTTTGACTTTTACAGCTTAGCGCCTTCTACAGCTTCACgcgc        >  2:319497/1‑80 (MQ=255)                                                                    aaaaaaCGCGACAAAAATGCCCGATCCTCGATCGGGCATTTTGACTTTTACAGCTTAGCGCCTTCTACAGCTTCACgcgc        >  2:319504/1‑80 (MQ=255)                                                                     aaaaaCGCGACAAAAATGCCCGATCCTCGATCGGGCATTTTGACTTTTACAGCTTAGCGCCTTCTACAGCTTCACGCGCCa      <  1:319479/81‑1 (MQ=255)                                                                         aCGCGACAAAAATGCCCGATCCTCGATCGGGCATTTTGACTTTTACAGCTTAGCGCCTTCTACAGCTTCACGCGCCAGCtt  <  2:319480/81‑1 (MQ=255)                                                                             |                                                                             TTGTTGACGCCATTGAACGCGCGAAGCACGCCGCCGGACAGGTAAAAGTACAGGGTTAAACCCGGCAAAAAAACGCTACAAAAATGCCCGATCCTCGATCGGGCATTTTGACTTTTACAGCTTAGCGCCTTCTACAGCTTCACGCGCCAGCTT  >  NC_000913/1932001‑1932153

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑