Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,336,405	A→G	S209S (AGT→AGC)	adiC ←	arginine:agmatine antiporter

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,336,405	0	A	G	100.0%	169.7 / NA	47	S209S (AGT→AGC)	adiC	arginine:agmatine antiporter
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (25/22);  total (25/22)

ATTAGGGATCATCCCCATAATCGCGGTGGTAGAAAGTACATAGCAAACGGCGGCAATCAATACCCCACCAATGGTGGCGATAGGGACATTGCGTTTCGGGTTTTTCACCACACCTGCGGCAACGGAGGCACTTTCCACACCGATGAACGACCACAGCGTAACGTTAAGGGTACTTTGAATTGCACCGAAGGTGCCCAGGCCGCTGACGTTCCATGCCGCCATATAGGTTTCACCACGGAACCAGAACCAGCCAAATACGG  >  NC_000913/4336276‑4336535                                                                                                                                  |                                                                                                                                   aTTAGGGATCATCCCCATAATCGCGGTGGTAGAAAGTACATAGCAAACGGCGGCAATCAATACCCCACCAATGGTGGCGATAGGGACATTGCGTTTCGGGTTTTTCACCACACCTGCGGCAACGGAGGCGCTTTCcaca                                                                                                                           <  1:500182/139‑1 (MQ=255)     gggATCATCCCCATAATCGCGGTGGTAGCAAGTACATAGCAAACGGCGGCAATCAATACCCCACCAATGGTGGCGATAGGGACATTGCGTTTCGGGTTTTTCACCACACCTGCGGCAACGGAGGCGCTTTCCACACCGa                                                                                                                       >  2:754458/1‑139 (MQ=255)     gggATCATCCCCATAATCGCGGTGGTAGAAAGTACATAGCAAACGGCGGCAATCAATACCCCACCAATGGTGGCGATAGGGACATTGCGTTTCGGGTTTTTCACCACACCTGCGGCAACGGAGGCGCTTTCCACACCGa                                                                                                                       >  1:1192536/1‑139 (MQ=255)     gggATCATCCCCATAATCGCGGTGGTAGAAAGTACATAGCAAACGGCGGCAATCAATACCCCACCAATGGTGGCGATAGGGACATTGCGTTTCGGGTTTTTCACCACACCTGCGGCAACGGAGGCGCTTTCCACACCGa                                                                                                                       >  2:354400/1‑139 (MQ=255)                cATAATCGCGGTGGTAGAAAGTACATAGCAAACGGCGGCAATCAATACCCCACCAATGGTGGCGATAGGGACATTGCGTTTCGGGTTTTTCACCACACCTGCGGCAACGGAGGCGCTTTCCACACCGATGAACGACcac                                                                                                            >  2:819484/1‑139 (MQ=255)                     tCGCGGTGGTAGAAAGTACATAGCAAACGGCGGCAATCAATACCCCACCAATGGTGGCGATAGGGACATTGCGTTTCGGGTTTTTCACCACACCTGCGGCAACGGAGGCGCTTTCCACACCGATGAACGACCACAGCGt                                                                                                       <  1:819484/139‑1 (MQ=255)                          gtggtAGAAAGTACATAGCAAACGGCGGCAATCAATACCCCACCAATGGTGGCGATAGGGACATTGCGTTTCGGGTTTTTCACCACACCTGCGGCAACGGAGGCGCTTTCCACACCGATGAACGACCACAGCGTAACGt                                                                                                  >  1:431689/1‑139 (MQ=255)                             gtAGAAAGTACATAGCAAACGGCGGCAATCAATACCCCACCAATGGTGGCGATAGGGACATTGCGTTTCGGGTTTTTCACCACACCTGCGGCAACGGAGGCGCTTTCCACACCGATGAACGACCACAGCGTAACGTTaa                                                                                               <  2:431689/139‑1 (MQ=255)                                   aGTACATAGCAAACGGCGGCAATCAATACCCCACCAATGGTGGCGATAGGGACATTGCGTTTCGGGTTTTTCACCACACCTGCGGCAACGGAGGCGCTTTCCACACCGATGAACGACCACAGCGTAACGTTAAGGGTAc                                                                                         >  1:876562/1‑139 (MQ=255)                                     tACATAGCAAACGGCGGCAATCAATACCCCACCAATGGTGGCGATAGGGACATTGCGTTTCGGGTTTTTCACCACACCTGCGGCAACGGAGGCGCTTTCCACACCGATg                                                                                                                     <  2:620248/109‑1 (MQ=255)                                     tACATAGCAAACGGCGGCAATCAATACCCCACCAATGGTGGCGATAGGGACATTGCGTTTCGGGTTTTTCACCACACCTGCGGCAACGGAGGCGCTTTCCACACCGATg                                                                                                                     >  1:620248/1‑109 (MQ=255)                                      aCATAGCAAACGGCGGCAATCTATACCCCACCAATGGTGGCGATAGGGACATTTCGTTTCGGGTTTTTCACCACACCTGCGGCAACGGAGGCGCTTTCCACACCGATGAACGACCACAGCGTAACGTTAAGGGTACttt                                                                                      <  1:560961/139‑1 (MQ=255)                                       cATAGCAAACGGCGGCAATCAATACCCCACCAATGGTGGCGATAGGGACATTGCGTTTCGGGTTTTTCACCACACCTGCGGCAACGGAGGCGCTTTCCACACCGATGAACGACCACAGCGTAACGTTAAGGGTACTTTg                                                                                     >  2:533121/1‑139 (MQ=255)                                            cAAACGGCGGCAATCAATACCCCACCAATGGTGGCGATAGGGACATTGCGTTTCGGGTTTTTCACCACACCTGCGGCAACGGAGGCGCTTTCCACACCGATGAACGACCACAGCGTAACGTTAAGGGTACTTTGAATTg                                                                                >  2:429122/1‑139 (MQ=255)                                               aCGGCGGCAATCAATACCCCACCAATGGTGGCGATAGGGACATTGCGTTTCGGGTTTTTCACCACACCTGCGGCAACGGAGGCGCTTTCCACACCGATGAACGACCACAGCGTAACGTTAAGGGTACTTTGAATTGCAc                                                                             >  1:1169211/1‑139 (MQ=255)                                               aCGGCGGCAATCAATACCCCACCAATGGTGGCGATAGGGACATTGCGTTTCGGGTTTTTCACCACACCTGCGGCAACGGAGGCGCTTTCCACACCGATGAACGACCACAGCGTAACGTTAAGGGTACTTTGAATTGCAc                                                                             >  2:873455/1‑139 (MQ=255)                                                      caatcaatACCCCACCAATGGTGGCGATAGGGACATTGCGTTTCGGGTTTTTCACCACACCTGCGGCAACGGAGGCGCTTTCCACACCGATGAACGACCACAGCGTAACGTTAAg                                                                                              <  2:897361/115‑1 (MQ=255)                                                      caatcaatACCCCACCAATGGTGGCGATAGGGACATTGCGTTTCGGGTTTTTCACCACACCTGCGGCAACGGAGGCGCTTTCCACACCGATGAACGACCACAGCGTAACGTTAAg                                                                                              >  1:897361/1‑115 (MQ=255)                                                                    ccaATGGTGGCGATAGGGACATTGCGTTTCGGGTTTTTCACCACACCTGCGGCAACGGAGGCGCTTTCCACACCGATGAACGACCACAGCGt                                                                                                       >  2:373696/1‑92 (MQ=255)                                                                    ccaATGGTGGCGATAGGGACATTGCGTTTCGGGTTTTTCACCACACCTGCGGCAACGGAGGCGCTTTCCACACCGATGAACGACCACAGCGt                                                                                                       <  1:373696/92‑1 (MQ=255)                                                                    ccaATGGTGGCGATAGGGACATTGCGTTTCGGGTTTTTCACCACACCTGCGGCAACGGAGGCGCTTTCCACACCGATGAACGACCACAGCGTAACGTTAAGGGTACTTTGAATTGCACCGAAGGTGCCCAGGCCGCTGa                                                        >  2:616796/1‑139 (MQ=255)                                                                        tggtggCGATAGGGACATTGCGTTTCGGGTTTTTCACCACACCTGCGGCAACGGAGGCGCTTTCCACACCGATGAACGACCACAGCGTAACGTTAAGGGTACTTTGAATTGCACCGAAGGTGCCCAGGCCGCTGACGtt                                                    <  2:876562/139‑1 (MQ=255)                                                                         ggtggCGATAGGGACATTGCGTTTCGGGTTTTTCACCACACCTGCGGCAACGGAGGCGCTTTCCACACCGATGAACGACCACAGCGTAACGTTAAg                                                                                              >  2:599467/1‑96 (MQ=255)                                                                         ggtggCGATAGGGACATTGCGTTTCGGGTTTTTCACCACACCTGCGGCAACGGAGGCGCTTTCCACACCGATGAACGACCACAGCGTAACGTTAAg                                                                                              <  1:599467/96‑1 (MQ=255)                                                                          gtggCGATAGGGACATTGCGTTTCGGGTTTTTCACCACACCTGCGGCAACGGAGGCGCTTTCCACACCGATGAACGACCACAGCGTAACGTTAAGGGTACTTTGAATTGCACCGAAGGTGCCCAGGCCGCTGACGTTcc                                                  >  2:640411/1‑139 (MQ=255)                                                                           tggCGATAGGGACATTGCGTTTCGGGTTTTTCACCACACCTGCGGCAACGGAGGCGCTTTCCACACCGATGAACGACCACAGCGTAACGTTAAGGGTACTTTGAATTGCACCGAAGGTGCCCAGGCCGCTGACGTTCCa                                                 <  2:163563/139‑1 (MQ=255)                                                                           tggCGATAGGGACATTGCGTTTCGGGTTTTTCACCACACCTGCGGCAACGGAGGCGCTTTCCACACCGATGAACGACCACAGCGTAACGTTAAGGGTACTTTGAATTGCACCGAAGGTGCCCAGGCCGCTGACGTTCCa                                                 <  2:326047/139‑1 (MQ=255)                                                                             gCGATAGGGACATTGCGTTTCGGGTTTTTCACCACACCTGCGGCAACGGAGGCGCTTTCCACACCGATGAACGACCAAAGCGTAACGTTAAGGGTACTTTGAATTGCACCGAAGGTGCCCAGGCCGCTGACGTTCCATg                                               >  1:38560/1‑139 (MQ=255)                                                                                aTAGGGACATTGCGTTTCGGGTTTTTCACCACACCTGCGGCAACGGAGGCGCTTTCCACACCGATGAACGACCACAGCGTAACGTTAAGGGTACTTTGAATTGCACCGAAGGTGCCCAGGCCGCTGACGTTCCATgccg                                            <  1:533121/139‑1 (MQ=255)                                                                                  aGGGACATTGCGTTTCGGGTTTTTCACCACACCTGCGGCAACGGAGGCGCTTTCCACACCGATGAACGACCACAGCGTAACGTTAAGGGTACTTTGAATTGCACCGAAg                                                                        >  1:669794/1‑109 (MQ=255)                                                                                  aGGGACATTGCGTTTCGGGTTTTTCACCACACCTGCGGCAACGGAGGCGCTTTCCACACCGATGAACGACCACAGCGTAACGTTAAGGGTACTTTGAATTGCACCGAAg                                                                        <  2:669794/109‑1 (MQ=255)                                                                                    ggACATTGCGTTTCGGGTTTTTCACCACACCTGCGGCAACGGAGGCGCTTTCCACACCGATGAACGACCACAGCGTAACGTTAAGGGTACTTTGAATTGCACCGAAGGTGCCCAGGCCGCTGACGTTCCATGCCGCCat                                        >  2:659915/1‑139 (MQ=255)                                                                                     gACATTGCGTTTCGGGTTTTTCACCACACCTGCGGCAACGGAGGCGCTTTCCACACCGATGAACGACCACAGCGTAACGTTAAGGGTACTTTGAATTGCACCGAAGGTGCCCAGGCCGCTGACGTTCCATGCCGCCata                                       <  1:984943/139‑1 (MQ=255)                                                                                     gACATTGCGTTTCGGGTTTTTCACCACACCTGCGGCAACGGAGGCGCTTTCCACACCGATGAACGACCACAGCGTAACGTTAAGGGTACTTTGAATTGCACCGAAGGTGCCCAGGCCGCTGACGTTCCATGCCGCCata                                       <  1:873455/139‑1 (MQ=255)                                                                                     gACATTGCGTTTCGGGTTTTTCACCACACCTGCGGCAACGGAGGCGCTTTCCACACCGATGAACGACCACAGCGTAACGTTAAGGGTACTTTGAATTGCACCGAAGGTGCCCAGGCCGCTGACGTTCCATGCCGCCata                                       <  2:819693/139‑1 (MQ=255)                                                                                          tGCGTTTCGGGTTTTTCACCACACCTGCGGCAACGGAGGCGCTTTCCACACCGATGAACGACCACAGCGTAACGTTAAGGGTACTTTGAATTGCACCGAAGGTGCCCAGGCCGCTGACGTTCCATGCCGCCATATAGGt                                  >  1:23267/1‑139 (MQ=255)                                                                                               ttCGGGTTTTTCACCACACCTGCGGCAACGGAGGCGCTTTCCACACCGATGAACGACCACAGCGTAACGTTAAGGGTACTTTGAATTGCACCGAAGGTGCCCa                                                                 >  2:915005/1‑103 (MQ=255)                                                                                               ttCGGGTTTTTCACCACACCTGCGGCAACGGAGGCGCTTTCCACACCGATGAACGACCACAGCGTAACGTTAAGGGTACTTTGAATTGCACCGAAGGTGCCCa                                                                 <  1:915005/103‑1 (MQ=255)                                                                                                   ggTTTTTCACCACACCTGCGGCAACGGAGGCGCTTTCCACACCGATGAACGACCACAGCGTAACGTTAAGGGTACTTTGAATTGCACCGAAGGTGCCCAGGCCGCTGACGTTCCATGCCGCCATATAGGTTTCACCACg                         <  1:754458/139‑1 (MQ=255)                                                                                                            ccacACCTGCGGCAACGGAGGCGCTTTCCACACCGATGAACGACCACAGCGTAACGTTAAGGGTACTTTGAATTGCACCGAAGGTGCCCAGGCCGCTGACGTTCCATGCCGCCATATAGGTTTCACCACGgaaccagaa                <  1:616796/139‑1 (MQ=255)                                                                                                             cacaCCTGCGGCAACGGAGGCGCTTTCCACACCGATGAACGACCACAGCGTAACGTTAAGGGTACTTTGAATTGCACCGAAGGTGCCCAGGCCGCTGACGTTCCATGCCGCCATATAGGTTTCACCACGgaaccagaac               <  1:1181351/139‑1 (MQ=255)                                                                                                                  cTGCGGCAACGGAGGCGCTTTCCACACCGATGAACGACCACAGTGTAACGTTAAGGGTACTTTGAATTGCACCGAAGGTGCCCAGGCCGCTGACGTTCCATGCCGCCATATAGGTTTCACCACGGAACCAGAACCAGcc          >  2:887369/1‑139 (MQ=255)                                                                                                                     cGGCAACGGAGGCGCTTTCCACACCGATGAACGACCACAGCGTAACGtt                                                                                                 >  2:599222/1‑49 (MQ=255)                                                                                                                     cGGCAACGGAGGCGCTTTCCACACCGATGAACGACCACAGCGTAACGtt                                                                                                 <  1:599222/49‑1 (MQ=255)                                                                                                                          aCGGAGGCGCTTTCCACACCGATGAACGACCACAGCGTAACGTTAAGGGTACTTTGAATTGCACCGAAGGTGCCCAGGCCGCTGACGTTCCATGCCGCCATATAGGTTTCACCACGGAACCAGAACCAGCCAAATACgg  >  2:77281/1‑139 (MQ=255)                                                                                                                              aGGCGCTTTCCACACCGATGAACGACCACAGCGTAAc                                                                                                    <  2:648585/37‑1 (MQ=255)                                                                                                                              aGGCGCTTTCCACACCGATGAACGACCACAGCGTAAc                                                                                                    >  1:648585/1‑37 (MQ=255)                                                                                                                                  |                                                                                                                                   ATTAGGGATCATCCCCATAATCGCGGTGGTAGAAAGTACATAGCAAACGGCGGCAATCAATACCCCACCAATGGTGGCGATAGGGACATTGCGTTTCGGGTTTTTCACCACACCTGCGGCAACGGAGGCACTTTCCACACCGATGAACGACCACAGCGTAACGTTAAGGGTACTTTGAATTGCACCGAAGGTGCCCAGGCCGCTGACGTTCCATGCCGCCATATAGGTTTCACCACGGAACCAGAACCAGCCAAATACGG  >  NC_000913/4336276‑4336535

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