Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	W3110S.gb	2,875,612	G→A	D126N (GAT→AAT)	iap →	aminopeptidase in alkaline phosphatase isozyme conversion

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	W3110S.gb	2,875,612	0	G	A	100.0%	190.9 / NA	63	D126N (GAT→AAT)	iap	aminopeptidase in alkaline phosphatase isozyme conversion
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (25/38);  total (25/38)

ATGAAGGCAAAGCGCCGCAGCAGATCATCATTATGGCGCATCTGGATACTTACGCCCCGCTGAGCGATGCTGACGCCGATGCCAATCTCGGCGGGCTGACGTTACAAGGAATGGATGATAACGCCGCAGGTT  >  W3110S.gb/2875547‑2875678                                                                  |                                                                   aTGAAGGCAAAGCGCCGCAGCAGATCATCATTATGGCGCATCTGGATACTTACGCCCCGCTGAGCAATGCt                                                               <  1:2808255/71‑1 (MQ=255) aTGAAGGCAAAGCGCCGCAGCAGATCATCATTATGGCGCATCTGGATACTTACGCCCCGCTGAGCAATGCt                                                               <  1:3312871/71‑1 (MQ=255) aTGAAGGCAAAGCGCCGCAGCAGATCATCATTATGGCGCATCTGGATACTTACGCCCCGCTGAGCAATGCt                                                               <  1:4432278/71‑1 (MQ=255) aTGAAGGCAAAGCGCCGCAGCAGATCATCATTATGGCGCATCTGGATACTTACGCCCCGCTGAGCAATGCt                                                               <  1:1610594/71‑1 (MQ=255)  tGAAGGCAAAGCGCCGCAGCAGATCATCATTATGGCGCATCTGGATACTTACGCCCCGCTGAGCAATGCt                                                               <  1:150013/70‑1 (MQ=255)    aaGGCAAAGCGCCGCAGCAGATCATCATTATGGCGCATCTGGATACTTACGCCCCGCTGAGCAATGCTGa                                                             >  1:3923440/1‑70 (MQ=255)    aaGGCAAAGCGCCGCAGCAGATCATCATTATGGCGCATCTGGATACTTACGCCCCGCTGAGCAATGCTGa                                                             >  1:1546456/1‑70 (MQ=255)       gCAAAGCGCCGCAGCAGATCATCATTATGGCGCATCTGGATACTTACGCCCCGCTGAGCAATGCTGACGcc                                                         <  1:4632932/71‑1 (MQ=255)            gcgcCGCAGCAGATCATCATTATGGCGCATCTGGATACTTACGCCCCGCTGAGCAATGCTGACGCCGATGc                                                    >  1:615411/1‑71 (MQ=255)             cgccgcAGCAGATCATCATTATGGCGCATCTGGATACTTACGCCCCGCTGAGCAATGCTGACGCCGATGcc                                                   <  1:1310983/71‑1 (MQ=255)                      agATCATCATTATGGCGCATCTGGATACTTACGCCCCGCTGAGCAATGCTGACGCCGATGCCAATCTcggc                                          >  1:1929491/1‑71 (MQ=255)                       gATCATCATTATGGCGCATCTGGATACTTACGCCCCGCTGAGCAATGCt                                                               <  1:3884526/49‑1 (MQ=255)                       gATCATCATTATGGCGCATCTGGATACTTACGCCCCGCTGAGCAATGCt                                                               <  1:2941716/49‑1 (MQ=255)                         tcatcaTTATGGCGCATCTGGATACTTACGCCCCGCTGAGCAATGCTGACGcc                                                         <  1:3325069/53‑1 (MQ=255)                          catcatTATGGCGCATCTGGATACTTACGCCCCGCTGAGCAATGCTGACGCCGATGCCAATCTCGGCggg                                       >  1:930946/1‑70 (MQ=255)                          catcatTATGGCGCATCTGGATACTTACGCCCCGCTGAGCAATGCTGACGCCGATGCCAATCTCGGCGGGc                                      >  1:1697169/1‑71 (MQ=255)                          catcatTATGGCGCATCTGGATACTTACGCCCCGCTGAGCAATGCTGACGCCGATGCCAATCTCGGCGGGc                                      >  1:4264549/1‑71 (MQ=255)                          catcatTATGGCGCATCTGGATACTTACGCCCCGCTGAGCAATGCTGACGCCGATGCCAATCTCGGCGGGc                                      >  1:242067/1‑71 (MQ=255)                          catcatTATGGCGCATCTGGATACTTACGCCCCGCTGAGCAATGCTGACGCCGATGCCAATCTCGGCGGGc                                      >  1:954356/1‑71 (MQ=255)                           atcatTATGGCGCATCTGGATACTTACGCCCCGCTGAGCAATGCTGACGCCGATGCCAATCTCGGCGGGCt                                     <  1:1482605/71‑1 (MQ=255)                           atcatTATGGCGCATCTGGATACTTACGCCCCGCTGAGCAATGCTGACGCCGATGCCAATCTCGGCGGGCt                                     <  1:4653775/71‑1 (MQ=255)                              atTATGGCGCATCTGGATACTTACGCCCCGCTGAGCAATg                                                                 >  1:3802520/1‑40 (MQ=255)                              atTATGGCGCATCTGGATACTTACGCCCCGCTGAGCAATGCTGACGCCGATGCCAATCTCGGCGGGCTGAc                                  >  1:3499766/1‑71 (MQ=255)                              atTATGGCGCATCTGGATACTTACGCCCCGCTGAGCAATGCTGACGCCGATGCCAATCTCGGCGGGCTGAc                                  >  1:359173/1‑71 (MQ=255)                                    gcgcATCTGGATACTTACGCCCCGCTGAGCAATGCTGACGCCGATGCCAATCTCGGCGGGCTGACGTTAc                             <  1:3020171/70‑1 (MQ=255)                                     cgcATCTGGATACTTACGCCCCGCTGAGCAATGCTGACGCCGATGCCAATCTCGGCGGGCTGACGTTaagg                           <  1:4157419/71‑4 (MQ=255)                                       cATCTGGATACTTACGCCCCGCTGAGCATTGCTGACGCCGATGCCAATCTCGGCGGGCTGACGTTACAAgg                         <  1:2930205/71‑1 (MQ=255)                                       cATCTGGATACTTACGCCCCGCTGAGCAATGCTGACGCCGATGCCAATCTCGGCGGGCTGACGTTACAAgg                         <  1:2470903/71‑1 (MQ=255)                                       cATCTGGATACTTACGCCCCGCTGAGCAATGCTGACGCCGATGCCAATCTCGGCGGGCTGACGTTACAAgg                         <  1:3526989/71‑1 (MQ=255)                                       cATCTGGATACTTACGCCCCGCTGAGCAATGCTGACGCCGATGCCAATCTCGGCGGGCTGACGTTACAAgg                         <  1:1101333/71‑1 (MQ=255)                                       cATCTGGATACTTACGCCCCGCTGAGCAATGCTGACGCCGATGCCAATCTCGGCGGGCTGACGTTACAAgg                         <  1:4670025/71‑1 (MQ=255)                                       cATCTGGATACTTACGCCCCGCTGAGCAATGCTGACGCCGATGCCAATCTCGGCGGGCTGACGTTACAAgg                         <  1:4660465/71‑1 (MQ=255)                                       cATCTGGATACTTACGCCCCGCTGAGCAATGCTGACGCCGATGCCAATCTCGGCGGGCTGACGTTACAAgg                         <  1:1723963/71‑1 (MQ=255)                                       cATCTGGATACTTACGCCCCGCTGAGCAATGCTGACGCCGATGCCAATCTCGGCGGGCTGACGTTACAAgg                         <  1:1794266/71‑1 (MQ=255)                                       cATCTGGATACTTACGCCCCGCTGAGCAATGCTGACGCCGATGCCAATCTCGGCGGGCTGACGTTACAAgg                         <  1:2935169/71‑1 (MQ=255)                                       cATCTGGATACTTACGCCCCGCTGAGCAATGCTGACGCCGATGCCAATCTCGGCGGGCTGACGTTACAAgg                         <  1:415520/71‑1 (MQ=255)                                       cATCTGGATACTTACGCCCCGCTGAGCAATGCCGACGCCGATGCCAATCTCGGCGGGCTGACGTTACAAgg                         <  1:3688152/71‑1 (MQ=255)                                        aTCTGGATACTTACGCCCCGCTGAGCAATGCTGACGCCGATGCCAATCTCGGCGGGCTGACGTTACAAgg                         <  1:253116/70‑1 (MQ=255)                                        aTCTGGATACTTACGCCCCGCTGAGCAATGCTGACGCCGATGCCAATCTCGGCGGGCTGACGTTACAAgg                         <  1:3541400/70‑1 (MQ=255)                                           tGGATACTTACGCCCCGCTGAGCAATGCTGACGCCGATGCCAATCTCGGCGGGCTGACGTTACAAGGAATg                     >  1:3836379/1‑71 (MQ=255)                                           tGGATACTTACGCCCCGCTGAGCAATGCTGACGCCGATGCCAATCTCGGCGGGCTGACGTTACAAGGAATg                     >  1:300/1‑71 (MQ=255)                                           tGGATACTTACGCCCCGCTGAGCAATGCTGACGCCGATGCCAATCTCGGCGGGCTGACGTTACAAGGAATg                     >  1:1016100/1‑71 (MQ=255)                                           tGGATACTTACGCCCCGCTGAGCAATGCTGACGCCGATGCCAATCTCGGCGGGCTGACGTTACAAGGAATg                     >  1:120059/1‑71 (MQ=255)                                           tGGATACTTACGCCCCGCTGAGCAATGCTGACGCCGATGCCAATCTCGGCGGGCTGACGTTACAAGGAATg                     >  1:1579783/1‑71 (MQ=255)                                           tGGATACTTACGCCCCGCTGAGCAATGCTGACGCCGATGCCAATCTCGGCGGGCTGACGTTACAAGGAATg                     >  1:2401731/1‑71 (MQ=255)                                           tGGATACTTACGCCCCGCTGAGCAATGCTGACGCCGATGCCAATCTCGGCGGGCTGACGTTACAAGGAATg                     >  1:1734393/1‑71 (MQ=255)                                           tGGATACTTACGCCCCGCTGAGCAATGCTGACGCCGATGCCAATCTCGGCGGGCTGACGTTACAAGGAATg                     >  1:2127811/1‑71 (MQ=255)                                             gATACTTACGCCCCGCTGAGCAATGCTGACGCCGATGCCAATCTCGGCGGGCTGACGTTACAAGGAATgg                    <  1:1442740/70‑1 (MQ=255)                                                aCTTATGCCCCGCTGAGCAATGCTGACGCCGATGCCAATCTCGGCGGGCTGACGTTACAAGGAATGgatga                >  1:3636732/1‑71 (MQ=255)                                                aCTTACGCCCCGCTGAGCAATGCTGACGCCGATGCCAATCTCGGCGGGCTGACGTTACAAGGAATGgatga                >  1:2148041/1‑71 (MQ=255)                                                 cTTACGCCCCGCTGAGCAATGCTGACGCCGATGCCAATCTCGGCGGGCTGACGTTACAAGGAATGgatgat               <  1:3138540/71‑1 (MQ=255)                                                 cTTACGCCCCGCTGAGCAATGCTGACGCCGATGCCAATCTCGGCGGGCTGACGTTACAAGGAATGgatgat               <  1:1169024/71‑1 (MQ=255)                                                 cTTACGCCCCGCTGAGCAATGCTGACGCCGATGCCAATCTCGGCGGGCTGACGTTACAAGGAATGgatgat               <  1:1022903/71‑1 (MQ=255)                                                 cTTACGCCCCGCTGAGCAATGCTGACGCCGATGCCAATCTCGGCGGGCTGACGTTACAAGGAATGgatga                >  1:3571503/1‑70 (MQ=255)                                                 cTTACGCCCCGCTGAGCAATGCTGACGCCGATGCCAATCTCGGCGGGCTGACGTTACAAGGAATGgat                  <  1:1994277/68‑1 (MQ=255)                                                 cTTACGCCCCGCTGAGCAATGCTGACGCCGATGCCAATCTCGGCGGGCTGACGTTACAAGGAATGgat                  <  1:3192676/68‑1 (MQ=255)                                                 cTTACGCCCCGCTGAGCAATGCTGACGCCGATGCCAATCTCGGCGGGCTGACCTTACAAGGAATGgatgat               <  1:1002178/71‑1 (MQ=255)                                                   tACGCCCCGCTGAGCAATGCTGACGCCGATGCCAATCTCGGCGGGCTGACGTTACAAgg                         <  1:3962129/59‑1 (MQ=255)                                                     cGCCCCGCTGAGCAATGCTGACGCCGATGCCAATCTCGGCGGGCTGACGTTACAAGGAATGGATGATAAc            <  1:4032843/70‑1 (MQ=255)                                                     cGCCCCGCTGAGCAATGCTGACGCCGATGCCAATCTCGGCGGGCTGACGTTACAAGGAATGGATGATAAc            <  1:4482659/70‑1 (MQ=255)                                                           gCTGAGCAATGCTGACGCCGATGCCAATCTCGGCGGGCTGACGTTACAAGGAATg                     <  1:1069497/55‑1 (MQ=255)                                                              gAGCAATGCTGACGCCGATGCCAATCTCGGCGGGCTGACGTTACAAGGAATGGATGATAACGCCGCAGGtt  >  1:4150184/1‑71 (MQ=255)                                                              gAGCAATGCTGACGCCGATGCCAATCTCGGCGGGCTGACGTTACAAGGAATGGATGATAACGCCGCAGGtt  >  1:1998950/1‑71 (MQ=255)                                                                  |                                                                   ATGAAGGCAAAGCGCCGCAGCAGATCATCATTATGGCGCATCTGGATACTTACGCCCCGCTGAGCGATGCTGACGCCGATGCCAATCTCGGCGGGCTGACGTTACAAGGAATGGATGATAACGCCGCAGGTT  >  W3110S.gb/2875547‑2875678

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