Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	992,825	A→G	Y735C (TAC→TGC)	pepN →	aminopeptidase N

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	992,825	0	A	G	100.0%	176.7 / NA	49	Y735C (TAC→TGC)	pepN	aminopeptidase N
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (18/31);  total (18/31)

TGGTGAAACGCATCTGGCTGATGTGCTGGTGAGCAAGCAGTTCCACGAAGCAAACAATATGACTGATGCGCTGGCGGCGCTTTCTGCGGCGGTTGCCGCACAGCTGCCTTGCCGTGACGCGCTGATGCAGGAGTACGACGACAAGTGGCATCAGAACGGTCTGGTGATGGATAAATGGTTTATCCTGCAAGCCACCAGCCCGGCGGCGAATGTGCTGGAGACGGTGCGCGGCCTGTTGCAGCATCGCTCATTTACCATGAGCAACCCGAA  >  NC_000913/992691‑992960                                                                                                                                       |                                                                                                                                        tGGTGAAACGCATCTGGCTGATGTGCTGGTGAGCAAGCAGTTCCACGAAGCAAACAATATGACTGATGCGCTGGCGGCGCTTTCTGCGGCGGTTGCCGCACAGCTGCCTTGCCGTGACGCGCTGATGCAGGAGTgcgac                                                                                                                                     <  1:861110/139‑1 (MQ=255)      aaaCGCATCTGGCTGATGTGCTGGTGAGCAAGCAGTTCCACGAAGCAAACAATATGACTGATGCGCTGGCGGCGCTTTCTGCGGCGGTTGCCGCACAGCTGCCTTGCCGTGACGCGCTGATGCAGGAGTGCGACGACaa                                                                                                                                >  2:1364689/1‑139 (MQ=255)       aaCGCATCTGGCTGATGTGCTGGTGAGCAAGCAGTTCCACGAAGCAAACAATATGACTGATGCGCTGGCGGCGCTTTCTGCGGCGGTTGCCGCACAGCTGCCTTGCCGTGACGCGCTGATGCAGGAGTGCGACGACAAg                                                                                                                               <  2:253751/139‑1 (MQ=255)        aCGCATCTGGCTGATGTGCTGGTGAGCAAGCAGTTCCACGAAGCAAACAATATGACTGATGCGCTGGCGGCGCTTTCTGCGGCGGTTGCCGCACAGCTGCCTTGCCGTGACGCGCTGATGCAGGAGTGCGACGACAAGt                                                                                                                              <  2:1014176/139‑1 (MQ=255)        aCGCATCTGGCTGATGTGCTGGTGAGCAAGCAGTTCCACGAAGCAAACAATATGACTGATGCGCTGGCCGCGCTTTCTGCGGCGGTTGCCGCACAGCTGCCTTGCCGTGACGCGCTGATGCAGGAGTGCGACGACAAGt                                                                                                                              >  1:1321777/1‑139 (MQ=255)            aTCTGGCTGATGTGCTGGTGAGCAAGCAGTTCCACGAAGCAAACAATATGACTGATGCGCTGGCGGCGCTTTCTGCGGCGGTTGCCGCACAGCTGCCTTGCCGTGACGCGCTGATGCAGGAGTGCGACGACAAGTGGCa                                                                                                                          >  1:1048404/1‑139 (MQ=255)              cTGGCTGATGTGCTGGTGAGCAAGCAGTTCCACGAAGCAAACAATATGACTGATGCGCTGGCGGCGCTTTCTGCGGCGGTTACCGCACAGCTGCCTTGCCGTGACGCGCTGATGCAGGAGTGcgacga                                                                                                                                   >  1:1063715/1‑128 (MQ=255)              cTGGCTGATGGGCTGGTGAGCAAGCAGTTCCACGAAGCAAACAATATGACTGATGCGCTGGCGGCGCTTTCTGCGGCGGTTGCCGCACAGCTGCCTTGCCGTGACGCGCTGATGCAGGAGTGcgacga                                                                                                                                   <  2:1063715/128‑1 (MQ=255)                  cTGATGTGCTGGTGAGCAAGCAGTTCCACGAAGCAAACAATATGACTGATGCGCTGGCGGCGCTTTCTGCGGCGGTTGCCGCACAGCTGCCTTGCCGTGACGCGCTGATGCAGGAGTGCGACGACAAGTGGCATCAGaa                                                                                                                    <  2:676358/139‑1 (MQ=255)                  cTGATGTGCTGGTGAGCAAGCAGTTCCACGAAGCAAACAATATGACTGATGCGCTGGCGGCGCTTTCTGCGGCGGTTGCCGCACAGCTGCCTTGCCGTGACGCGCTGATGCAGGAGTGCGACGACAAGTGGCATCAGaa                                                                                                                    <  2:1321777/139‑1 (MQ=255)                  cTGATGTGCTGGTGAGCAAGCAGTTCCACGAAGCAAACAATATGACTGATGCGCTGGCGGCGCTTTCTGCGGCGGTTGCCGCACAGCTGCCTTGCCGTGACGCGCTGATGCAGGAGTGCGACGACAAGTGGCATCAGaa                                                                                                                    <  1:426342/139‑1 (MQ=255)                      tgtgCTGGTGAGCAAGCAGTTCCACGAAGCAAACAATATGACTGATGCGCTGGCGGCGCTTTCTGCGGCGGTTGCCGCACAGCTGCCTTGCCGTGACGCGCTGATGCAGGAGTGCGACGACAAGt                                                                                                                              >  2:946941/1‑125 (MQ=255)                      tgtgCTGGTGAGCAAGCAGTTCCACGAAGCAAACAATATGACTGATGCGCTGGCGGCGCTTTCTGCGGCGGTTGCCGCACAGCTGCCTTGCCGTGACGCGCTGATGCAGGAGTGCGACGACAAGt                                                                                                                              <  1:946941/125‑1 (MQ=255)                       gtgCTGGTGAGCAAGCAGTTCCACGAAGCAAACAATATGACTGATGCGCTGGCGGCGCTTTCTGCGGCGGTTGCCGCACAGCTGCCTTGCCGTGACGCGCTGATGCAGGAGTGCGACGACAAGTGGCATCAGAACGGTc                                                                                                               <  1:106237/139‑1 (MQ=255)                           tGGTGAGCAAGCAGTTCCACGAAGCAAACAATATGACTGATGCGCTGGCGGCGCTTTCTGCGGCGGTTGCCGCACAGCTGCCTTGCCGTGACGCGCTGATGCAGGAGTGCGGCGACAAGTGGCATCAGAACGGTCTGGt                                                                                                           <  1:1618009/139‑1 (MQ=255)                           tGGTGAGCAAGCAGTTCCACGAAGCAAACAATATGACTGATGCGCTGGCGGCGCTTTCTGCGGCGGTTGCCGCACAGCTGCCTTGCCGTGACGCGCTGATGCAGGAGTGCGACGACAAGTGGCAt                                                                                                                         >  1:1256842/1‑125 (MQ=255)                           tGGTGAGCAAGCAGTTCCACGAAGCAAACAATATGACTGATGCGCTGGCGGCGCTTTCTGCGGCGGTTGCCGCACAGCTGCCTTGCCGTGACGCGCTGATGCAGGAGTGCGACGACAAGTGGCAt                                                                                                                         <  2:1256842/125‑1 (MQ=255)                           tGGTGAGCAAGCAGTTCCACGAAGCAAACAATATGACTGATGCGCTGGCGGCGCTTTCTGCGGCGGTTGCCGCACAGCTGCCTTGCCGTGACGCGCTGATGCAGGAGTGCGACGACAAGTGGCATCAGAACGGt                                                                                                                <  1:633194/134‑1 (MQ=255)                           tGGTGAGCAAGCAGTTCCACGAAGCAAACAATATGACTGATGCGCTGGCGGCGCTTTCTGCGGCGGTTGCCGCACAGCTGCCTTGCCGTGACGCGCTGATGCAGGAGTGCGACGACAAGTGGCATCAGAACGGt                                                                                                                >  2:633194/1‑134 (MQ=255)                             gTGAGCAAGCAGTTCCACGAAGCAAACAATATGACTGATGCGCTGGCGGCGCTTTCTGCGGCGGTTGCCGCACAGCTGCCTTGCCGTGACGCGCTGATGCAGGAGTGCGACGACAAGTGGCATCAGAACGGTCTGGTGa                                                                                                         <  2:1019712/139‑1 (MQ=255)                             gTGAGCAAGCAGTTCCACGAAGCAAACAATATGACTGATGCGCTGGCGGCGCTTTCTGCGGCGGTTGCCGCACAGCTGCCTTGCCGTGACGCGCTGATGCAGGAGTGCGACGACAAGTGGCATCAGAACGGTCTGGTGa                                                                                                         <  1:99703/139‑1 (MQ=255)                               gAGCAAGCAGTTCCACGAAGCAAACAATATGACTGATGCGCTGGCGGCGCTTTCTGCGGCGGTTGCCGCACAGCTGCCTTGCCGTGACGCGCTGATGCAGGAGTGCGACGACAAGTGGCATCAGaa                                                                                                                    >  1:154426/1‑126 (MQ=255)                               gAGCAAGCAGTTCCACGAAGCAAACAATATGACTGATGCGCTGGCGGCGCTTTCTGCGGCGGTTGCCGCACAGCTGCCTTGCCGTGACGCGCTGATGCAGGAGTGCGACGACAAGTGGCATCAGaa                                                                                                                    <  2:154426/126‑1 (MQ=255)                                  caagcaGTTCCACGAAGCAAACAATATGACTGATGCGCTGGCGGCGCTTTCTGCGGCGGTTGCCGCACAGCTGCCTTGCCGTGACGCGCTGATGCAGGAGTGCGACGACAAGTGGCATCAGAACGGTCTGGTGATGGAt                                                                                                    >  2:119851/1‑139 (MQ=255)                                    agcaGTTCCACGAAGCAAACAATATGACTGATGCGCTGGCGGCGCTTTCTGCGGCGGTTGCCGCACAGCTGCCTTGCCGTGACGCGCTGATGCAGGAGTGCGACGACAAGTGGCATCAGAACGGTCTGGTGATGGATaa                                                                                                  <  1:539775/139‑1 (MQ=255)                                                   cAAACAATATGACTGATGCGCTGGCGGCGCTTTCTGCGGCGGTTGCCGCACAGCTGCCTTGCCGTGACGCGCTGATGCAGGAGTGCGACGACAAGTGGCATCAGaa                                                                                                                    >  1:959215/1‑106 (MQ=255)                                                   cAAACAATATGACTGATGCGCTGGCGGCGCTTTCTGCGGCGGTTGCCGCACAGCTGCCTTGCCGTGACGCGCTGATGCAGGAGTGCGACGACAAGTGGCATCAGaa                                                                                                                    <  2:959215/106‑1 (MQ=255)                                                           atGACTGATGCGCTGGCGGCGCTTTCTGCGGCGGTTGCCGCACAGCTGCCTTGCCGTGACGCGCTGATGCAGGAGTGCGACGACAAGt                                                                                                                              >  2:678461/1‑88 (MQ=255)                                                           atGACTGATGCGCTGGCGGCGCTTTCTGCGGCGGTTGCCGCACAGCTGCCTTGCCGTGACGCGCTGATGCAGGAGTGCGACGACAAGt                                                                                                                              <  1:678461/88‑1 (MQ=255)                                                                  aTGCGCTGGCGGCGCTTTCTGCGGCGGTTGCCGCACAGCTGCCTTGCCGTGACGCGCTGATGCAGGAGTGCGACGACAAGTGGCATCAGaa                                                                                                                    <  1:1130667/91‑1 (MQ=255)                                                                  aTGCGCTGGCGGCGCTTTCTGCGGCGGTTGCCGCACAGCTGCCTTGCCGTGACGCGCTGATGCAGGAGTGCGACGACAAGTGGCATCAGaa                                                                                                                    >  2:1130667/1‑91 (MQ=255)                                                                  aTGCGCTGGCGGCGCTTTCTGCGGCGGTTGCCGCACAGCTGCCTTGCCGTGACGCGCTGATGCAGGAGTGCGACGACAAGTGGCATCAGAACGGTCTGGTGATGGATaa                                                                                                  >  2:823411/1‑109 (MQ=255)                                                                  aTGCGCTGGCGGCGCTTTCTGCGGCGGTTGCCGCACAGCTGCCTTGCCGTGACGCGCTGATGCAGGAGTGCGACGACAAGTGGCATCAGAACGGTCTGGTGATGGATaa                                                                                                  <  1:823411/109‑1 (MQ=255)                                                                       cTGGCGGCGCTTTCTGCGGCGGTTGCCGCACAGCTGCCTTGCCGTGACGCGCTGATGCAGGAGTGCGACGACAAGTGGCATCAGAACGGTCTGGTGATGGATAAATGGTTTATCCTGCAAGCCACCAGCCCGGCGGCGa                                                               <  1:1364689/139‑1 (MQ=255)                                                                        tGGCGGCGCTTTCTGCGGCGGTTGCCGCACAGCTGCCTTGCCGTGACGCGCTGATGCAGGAGTGCGACGACAAGTGGCATCAGAACGGTCTGGTGATGGATAAATGGTTTATCCTGCAAGCCACCAGCCCGGCGGCGaa                                                              <  1:749311/139‑1 (MQ=255)                                                                                          cggTTGCCGCACAGCTGCCTTGCCGTGACGCGCTGATGCAGGAGTGCGACGACAAGTGGCATCAGAACGGTCTGGTGATGGATAAATGGTTTATCCTGCAAGCCACCAGCCCGGCGGCGAATGTGCTGGAGACGGTgcg                                            <  1:554576/139‑1 (MQ=255)                                                                                            gTTGCCGCACAGCTGCCTTGCCGTGACGCGCTGATGCAGGAGTGCGACGACAAGTGGCATCAg                                                                                                                      <  1:627513/63‑1 (MQ=255)                                                                                            gTTGCCGCACAGCTGCCTTGCCGTGACGCGCTGATGCAGGAGTGCGACGACAAGTGGCATCAg                                                                                                                      >  2:627513/1‑63 (MQ=255)                                                                                              tGCCGCACAGCTGCCTTGCCGTGACGCGCTGATGCAGGAGTGCGACGACAAGTGGCATCAGAACGGTCTGGTGATGGATAAATGGTTTATCCTGCAAGCCACCAGCCCGGCGGCGAATGTGCTGGAGACGGTGCGCGGc                                        <  1:119851/139‑1 (MQ=255)                                                                                              tGCCGCACAGCTGCCTTGCCGTGACGCGCTGATGCAGGAGTGCGACGACAAGTGGCATAAGAACGGTCTGGTGATGGATAAATGGTTTATCCTGCAAGCCACCAGCCCGGCGGCGAATGTGCTGGAGACGGTGCGCGGc                                        <  2:1048404/139‑1 (MQ=255)                                                                                               gCCGCACAGCTGCCTTGCCGTGACGCGCTGATGCAGGAGTGCGACGACAAGTGGCATCAGAACGGTCTGGTGATGGATAAATGGTTTATCCTGCAAGCCACCAGCCCGGCGGCGAATGTGCTGGAGACGGTGCGCGGcc                                       <  1:910399/139‑1 (MQ=255)                                                                                                       gCTGCCTTGCCGTGACGCGCTGATGCAGGAGTGCGAAGACAAGTGGCATCAGAACGGTCTGGTGATGGATAAATGGTTTATCCTGCAAGCCACCAGCCCGGCGGCGAATGTGCTGGAGACGGTGCGCGGCCTGTTgcag                               >  2:1143076/1‑139 (MQ=255)                                                                                                        cTGCCTTGCCGTGACGCGCTGATGCAGGAGTGCGACGACAAGTGGCATCAGaa                                                                                                                    <  1:857876/53‑1 (MQ=255)                                                                                                        cTGCCTTGCCGTGACGCGCTGATGCAGGAGTGCGACGACAAGTGGCATCAGaa                                                                                                                    >  2:857876/1‑53 (MQ=255)                                                                                                            cTTGCCGTGACGCGCTGATGCAGGAGTGCGACGACAAGTGGCATCAGAACGGTCTGGTGATGGATAAATGGTTTATCCTGCAAGCCACCAGCCCGGCGGCGAATGTGCTGGAGACGGTGCGCGGCCTGTTGCAGCATCg                          <  1:1149565/139‑1 (MQ=255)                                                                                                                          cTGATGCAGGAGTGCGACGACAAGTGGCATCAGaa                                                                                                                    >  1:124766/1‑35 (MQ=255)                                                                                                                          cTGATGCAGGAGTGCGACGACAAGTGGCATCAGaa                                                                                                                    <  2:124766/35‑1 (MQ=255)                                                                                                                          cTGATGCAGGAGTGCGACGACAAGTGGCATCAGAACGGTCTGGTGATGGATAAATGGTTTATCCTGCAAGCCACCAGCCcggc                                                                    <  2:1644733/83‑1 (MQ=255)                                                                                                                          cTGATGCAGGAGTGCGACGACAAGTGGCATCAGAACGGTCTGGTGATGGATAAATGGTTTATCCTGCAAGCCACCAGCCcggc                                                                    >  1:1644733/1‑83 (MQ=255)                                                                                                                                    aGTGCGACGACAAGTGGCATCAGAACGGTCTGGTGATGGATAAATGGTTTATCCTGCAAGCCACCAGCCCGGCGGCGAATGTGCTGGAGACGGTGCGCGGCCTGTTGCAGCATCGCTCATTTACCATGAGCAACCCGaa  <  1:728400/139‑1 (MQ=255)                                                                                                                                       |                                                                                                                                        TGGTGAAACGCATCTGGCTGATGTGCTGGTGAGCAAGCAGTTCCACGAAGCAAACAATATGACTGATGCGCTGGCGGCGCTTTCTGCGGCGGTTGCCGCACAGCTGCCTTGCCGTGACGCGCTGATGCAGGAGTACGACGACAAGTGGCATCAGAACGGTCTGGTGATGGATAAATGGTTTATCCTGCAAGCCACCAGCCCGGCGGCGAATGTGCTGGAGACGGTGCGCGGCCTGTTGCAGCATCGCTCATTTACCATGAGCAACCCGAA  >  NC_000913/992691‑992960

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