Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,877,286	G→A	A236T (GCG→ACG)	iap →	alkaline phosphatase isozyme conversion protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,877,286	0	G	A	100.0%	74.2 / NA	22	A236T (GCG→ACG)	iap	alkaline phosphatase isozyme conversion protein
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (13/9);  total (13/9)

GATAAATTGTATTTCAACAGCGGTGTAAAAACCCCTGAGGCAGTAAGGAAATTAACGCGCGACAGGGCGCTGGCAATTGCGCGCAGTCACGGAATAGCCGCAACGACCAATCCGGGTTTGAATAAAAATTATCCGAAAGGCACTGGGTGTTGTAA  >  NC_000913/2877208‑2877362                                                                               |                                                                             gATAAATTGTATTTCAACAGCGGTGTAAAAACCCCTGAGGCAGTAAGGAAATTAACGCGCGACAGGGCGCTGGCAATTacg                                                                            >  2:496875/1‑81 (MQ=255)                   aGCGGTGTAAAAACCCCTGAGGCAGTAAGGAAATTAACGCGCGACAGGGCGCTGGCAATTACGCGCAGTCACGGAATAGcc                                                          >  2:496867/1‑81 (MQ=255)                   aGCGGTGTAAAAACCCCTGAGGCAGTAAGGAAATTAACGCGCGACAGGGCGCTGGCAATTACGCGCAGTCACGGAATAGcc                                                          >  2:496873/1‑81 (MQ=255)                       gtgtAAAAACCCCTGAGGCAGTAAGGAAATTAACGCGCGACAGGGCGCTGGCAATTACGCGCAGTCACGGAATAGCCGCaa                                                      >  1:496871/1‑81 (MQ=255)                        tgtAAAAACCCCTGAGGCAGTAAGGAAATTAACGCGCGACAGGGCGCTGGCAATTACGCGCAGTCACGGAATAGCCGCAAc                                                     <  2:496848/81‑1 (MQ=255)                           aaaaaCCCCTGAGGCAGTAAGGAAATTAACGCGCGACAGGGCGCTGGCAATTACGCGCAGTCACGGAATAGCCGCAACGAc                                                  >  1:496877/1‑81 (MQ=255)                           aaaaaCCCCTGAGGCAGTAAGGAAATTAACGCGCGACAGGGCGCTGGCAATTACGCGCAGTCACGGAATAGCCGCAACGAc                                                  >  2:496881/1‑81 (MQ=255)                              aaCCCCTGAGGCAGTAAGGAAATTAACGCGCGACAGGGCGCTGGCAATTACGCGCAGTCACGGAATAGCCGCAACGACCaa                                               >  1:496897/1‑81 (MQ=255)                               aCCCCTGAGGCAGTAAGGAAATTAACGCGCGACAGGGCGCTGGCAATTACGCGCAGTCACGGAATAGCCGCAACGACCAAt                                              <  1:496849/81‑1 (MQ=255)                                     gAGGCAGTAAGGAAATTAACGCGCGACAGGGCGCTGGCAATTACGCGCAGTCACGGAATAGCCGCAACGACCAATCCGGGt                                        >  2:496883/1‑81 (MQ=255)                                       ggCAGTAAGGAAATTAACGCGCGACAGGGCGCTGGCAATTACGCGCAGTCACGGAATAGCCGCAACGACCAATCCGGGttt                                      <  2:496850/81‑1 (MQ=255)                                              aGGAAATTAACGCGCGACAGGGCGCTGGCAATTACGCGCAGTCACGGAATAGCCGCAACGACCAATCCGGGTTTGAATaaa                               <  2:496851/81‑1 (MQ=255)                                                 aaaTTAACGCGCGACAGGGCGCTGGCAATTACGCGCAGTCACGGAATAGCCGCAACGACCAATCCGGGTTTGAATAAAAAt                            >  2:496889/1‑81 (MQ=255)                                                  aaTTAACGCGCGACAGGGCGCTGGCAATTACGCGCAGTCACGGAATAGCCGCAACGACCAATCCGGGTTTGAATAAAAAt                            >  2:496887/1‑80 (MQ=255)                                                   aTTAACGCGCGACAGGGCGCTGGCAATTACGCGCAGTCACGGAATAGCCGCAACGACCAATCCGGGTTTGAATAAAAATTa                          >  2:496874/1‑81 (MQ=255)                                                      aaCGCGCGACAGGGCGCTGGCAATTACGCGCAGTCACGGAATAGCCGCAACGACCAATCCGGGTTTGAATAAAAATTATc                        <  1:496852/80‑1 (MQ=255)                                                              aCAGGGCGCTGGCAATTACGCGCAGTCACGGGATAGCCGCAACGACCAATCCGGGTTTGAATAAAAATTATCCGAAAGGCa               <  2:496853/81‑1 (MQ=255)                                                                  ggCGCTGGCAATTACGCGCAGTCACGGAATAGCCGCAACGACCAATCCGGGTTTGAATAAAAATTATCCGAAAGGCACTgg           <  2:496854/81‑1 (MQ=255)                                                                       tGGCAATTACGCGCAGTCACGGAATAGCCGCAACGACCAATCCGGGTTTGAATAAAAATTATCCGAAAGGCACTGGGtgtt      >  2:496878/1‑81 (MQ=255)                                                                         gCAATTACGCGCAGTCACGGAATAGCCGCAACGACCAATCCGGGTTTGAATAAAAATTATCCGAAAGGCACTGGGtgttg     <  2:496856/80‑1 (MQ=255)                                                                         gCAATTACGCGCAGTCACGGAATAGCCGCAACGACCAATCCGGGTTTGAATAAAAATTATCCGAAAGGCACTGGGtgtt      <  2:496855/79‑1 (MQ=255)                                                                          cAATTACGCGCAGTCACGGAATAGCCGCAACGACCAATCCGGGTTTGAATAAAAATTATCCGAAAGGCACTGGGtgttgt    >  2:496885/1‑80 (MQ=255)                                                                           aaTTACGCGCAGTCACGGAATAGCCGCAACGACCAATCCGGGTTTGAATAAAAATTATCCGAAAGGCACTGGGTGTTGTaa  >  1:496890/1‑81 (MQ=255)                                                                               |                                                                             GATAAATTGTATTTCAACAGCGGTGTAAAAACCCCTGAGGCAGTAAGGAAATTAACGCGCGACAGGGCGCTGGCAATTGCGCGCAGTCACGGAATAGCCGCAACGACCAATCCGGGTTTGAATAAAAATTATCCGAAAGGCACTGGGTGTTGTAA  >  NC_000913/2877208‑2877362

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