Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,877,592	G→A	V338M (GTG→ATG)	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,592	0	G	A	100.0%	118.2 / NA	33	V338M (GTG→ATG)	iap	alkaline phosphatase isozyme conversion protein
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (22/11);  total (22/11)

AACATATTGATAAGGCTCTTCCTGGAAGAATAGAACGTCGCTGCCGTGACGTTATGCGGATAATGCTACCTCTGGTGAAGGAGTTGGCGAAGGCGTCTTGATGGGTTTGAAAATGGGAGCTGGGAGTTCTACCGCAGAGGCGGGGGAA  >  NC_000913/2877518‑2877665                                                                           |                                                                          aaCATATTGATAAGGCTCTTCCTGGAAGAATAGAACGTCGCTGCCGTGACGTTATGCGGATAATGCTACCTCTGATGAAgg                                                                     <  2:522078/81‑1 (MQ=255)  aCATATTGATAAGGCTCTTCCTGGAAGAATAGAACGTCGCTGCCGTGACGTTATGCGGATAATGCTACCTCTGATGAAGGa                                                                    <  2:522079/81‑1 (MQ=255)   cATATTGATAAGGCTCTTCCTGGAAGAATAGAACGTCGCTGCCGTGACGTTATGCGGATAATGCTACCTCTGATGAAGGAg                                                                   >  1:522104/1‑81 (MQ=255)   cATATTGATAAGGCTCTTCCTGGAAGAATAGAACGTCGCTGCCGTGACGTTATGCGGATAATGCTACCTCTGATGAAGGAg                                                                   >  2:522133/1‑81 (MQ=255)          aTAAGGCTCTTCCTGGAAGAATAGAACGTCGCTGCCGTGACGTTATGCGGATAATGCTACCTCTGATGAAGGAGTTGGCg                                                             <  2:522080/80‑1 (MQ=255)             aGGCTCTTCCTGGAAGAATAGAACGTCGCTGCCGTGACGTTATGCGGATAATGCTACCTCTGATGAAGGAGTTGGCGAAgg                                                         >  2:522110/1‑81 (MQ=255)                  ctTCCTGGAAGAATAGAACGTCGCTGCCGTGACGTTATGCGGATAATGCTACCTCTGATGAAGGAGTTGGCGAAGGCGTCt                                                    >  1:522122/1‑81 (MQ=255)                     ccTGGAAGAATAGAACGTCGCTGCCGTGACGTTATGCGGATAATGCTACCTCTGATGAAGGAGTTGGCGAAGGCGTCtt                                                   <  1:522081/79‑1 (MQ=255)                         gaagaaTAGAACGTCGCTGCCGTGACGTTATGCGGATAATGCTACCTCTGATGAAGGAGTTGGCGAAGGCGTCTTGATgg                                              <  1:522082/80‑1 (MQ=255)                            gaaTAGAACGTCGCTGCCGTGACGTTATGCGGATAATGCTACCTCTGATGAAGGAGTTGGCGAAGGCGTCTTGATGGGttt                                          >  1:522117/1‑81 (MQ=255)                               tAGAACGTCGCTGCCGTGACGTTATGCGGATAATGCTACCTCTGATGAAGGAGTTGGCGAAGGCGTCTTGATGGGTTTGa                                        >  1:522107/1‑80 (MQ=255)                                 gAACGTCGCTGCCGTGACGTTATGCGGATAATGCTACCTCTGATGAAGGAGTTGGCGAAGGCGTCTTGATGGGTTTGaaaa                                     >  1:522116/1‑81 (MQ=255)                                      tCGCTGCCGTGACGTTATGCGGATAATGCTACCTCTGATGAAGGAGTTGGCGAAGGCGTCTTGATGGGTTTGAAAATGGGa                                <  2:522083/81‑1 (MQ=255)                                      tCGCTGCCGTGACGTTATGCGGATAATGCTACCTCTGATGAAGGAGTTGGCGAAGGCGTCTTGATGGGTTTGAAAATGGGa                                <  1:522084/81‑1 (MQ=255)                                          tGCCGTGACGTTATGCGGATAATGCTACCTCTGATGAAGGAGTTGGCGAAGGCGTCTTGATGGGTTTGAAAATGGGAGCTg                            >  1:522126/1‑81 (MQ=255)                                          tGCCGTGACGTTATGCGGATAATGCTACCTCTGATGAAGGAGTTGGCGAAGGCGTCTTGATGGGTTTGAAAATGGGAGCTg                            >  2:522121/1‑81 (MQ=255)                                           gCCGTGACGTTATGCGGATAATGCTACCTCTGATGAAGGAGTTGGCGAAGGCGTCTTGATGGGTTTGAAAATGGGAGCTgg                           >  1:522125/1‑81 (MQ=255)                                                gACGTTATGCGGATAATGCTACCTCTGATGAAGGAGTTGGCGAAGGCGTCTTGATGGGTTTGAAAATGGGAGCTGGGAGtt                      >  1:522129/1‑81 (MQ=255)                                                   gTTATGCGGATAATGCTACCTCTGATGAAGGAGTTGGCGAAGGCGTCTTGATGGGTTTGAAAATGGGAGCTGGGAGTTCTa                   >  2:522130/1‑81 (MQ=255)                                                      aTGCGGATAATGCTACCTCTGATGAAGGAGTTGGCGAAGGCGTCTTGATGGGTTTGAAAATGGGAGCTGGGAGTTCTACCg                >  1:522143/1‑81 (MQ=255)                                                       tGCGGATAATGCTACCTCTGATGAAGGAGTTGGCGAAGGCGTCTTGATGGGTTTGAAAATGGGAGCTGGGAGTTCTACCg                <  1:522085/80‑1 (MQ=255)                                                        gCGGATAATGCTACCTCTGATGAAGGAGTTGGCGAAGGCGTCTTGATGGGTTTGAAAATGGGAGCTGGGAGTTCTACCg                >  2:522127/1‑79 (MQ=255)                                                          ggATAATGCTACCTCTGATGAAGGAGTTGGCGAAGGCGTCTTGATGGGTTTGAAAATGGGAGCTGGGAGTTCTACCGCaga            <  1:522086/81‑1 (MQ=255)                                                          ggATAATGCTACCTCTGATGAAGGAGTTGGCGAAGGCGTCTTGATGGGTTTGAAAATGGGAGCTGGGAGTTCTACCGCa              >  2:522134/1‑79 (MQ=255)                                                           gATAATGCTACCTCTGATGAAGGAGTTGGCGAAGGCGTCTTGATGGGTTTGAAAATGGGAGCTGGGAGTTCTACCGCagag           >  1:522153/1‑81 (MQ=255)                                                           gATAATGCTACCTCTGATGAAGGAGTTGGCGAAGGCGTCTTGATGGGTTTGAAAATGGGAGCTGGGAGTTCTACCGCagag           <  1:522087/81‑1 (MQ=255)                                                            aTAATGCTACCTCTGATGAAGGAGTTGGCGAAGGCGTCTTGATGGGTTCGAAAATGGGAGCTGGGAGTTCTACCGCAGAgg          >  1:522131/1‑81 (MQ=255)                                                             tAATGCTACCTCTGATGAAGGAGTTGGCGAAGGCGTCTTGATGGGTTTGAAAATGGGAGCTGGGAGTTCTACCGCAGAGGc         >  2:522118/1‑81 (MQ=255)                                                              aaTGCTACCTCTGATGAAGGAGTTGGCGAAGGCGTCTTGATGGGTTTGAAAATGGGAGCTGGGAGTTCTACCGCAGAGGCg        >  1:522135/1‑81 (MQ=255)                                                               aTGCTACCTCTGATGAAGGAGTTGGCGAAGGCGTCTTGATGGGTTTGAAAATGGGAGCTGGGAGTTCTACCGCAGAGGCgg       >  1:522132/1‑81 (MQ=255)                                                                  cTACCTCTGATGAAGGAGTTGGCGAAGGCGTCTTGATGGGTTTGAAAATGGGAGCTGGGAGTTCTACCGCAGAGGCgggg     >  1:522137/1‑80 (MQ=255)                                                                    aCCTCTGATGAAGGAGTTGGCGAAGGCGTCTTGATGGGTTTGAAAATGGGAGCTGGGAGTTCTACCGCAGAGGCGGGGGaa  >  2:522123/1‑81 (MQ=255)                                                                     ccTCTGATGAAGGAGTTGGCGAAGGCGTCTTGATGGGTTTGAAAATGGGAGCTGGGAGTTCTACCGCAGAGGCGGGGGa   <  2:522088/79‑1 (MQ=255)                                                                           |                                                                          AACATATTGATAAGGCTCTTCCTGGAAGAATAGAACGTCGCTGCCGTGACGTTATGCGGATAATGCTACCTCTGGTGAAGGAGTTGGCGAAGGCGTCTTGATGGGTTTGAAAATGGGAGCTGGGAGTTCTACCGCAGAGGCGGGGGAA  >  NC_000913/2877518‑2877665

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