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%	163.4 / NA	46	S209S (AGT→AGC)	adiC	arginine:agmatine antiporter
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (26/20);  total (26/20)

GGGATCATCCCCATAATCGCGGTGGTAGAAAGTACATAGCAAACGGCGGCAATCAATACCCCACCAATGGTGGCGATAGGGACATTGCGTTTCGGGTTTTTCACCACACCTGCGGCAACGGAGGCACTTTCCACACCGATGAACGACCACAGCGTAACGTTAAGGGTACTTTGAATTGCACCGAAGGTGCCCAGGCCGCTGACGTTCCATGCCGCCATATAGGTTTCACCACGGAACCAGAACCAGCCAAA  >  NC_000913/4336280‑4336530                                                                                                                              |                                                                                                                              gggATCATCCCCATAATCGCGGTGGTAGAAAGTACATAGCAAACGGCGGCAATCAATACCCCACCAATGGTGGCGATAGGGACATTGCGTTTCGGGTTTTTCACCACACCTGCGGCAACGGAGGCGCTTTCCACACCGa                                                                                                                  >  1:847083/1‑139 (MQ=255)     tcatcCCCATAATCGCGGTGGTAGAAAGTACATAGCAAACGGCGGCAATCAATACCCCACCAATGGTGGCGATAGGGACATTGCGTTTCGGGTTTTTCACCACACCTGCGGCAACGGAGGCGCTTTCCACACCGATGaa                                                                                                              <  1:803784/139‑1 (MQ=255)            cATAATCGCGGTGGTAGAAAGTACATAGCAAACGGCGGCAATCAATACCCCACCAATGGTGGCGATAGGGACATTGCGTTTCGGGTTTTTCACCACACCTGCGGCAACGGAGGCGCTTTCCACACCGATGAACGACcac                                                                                                       >  1:1382701/1‑139 (MQ=255)              tAATCGCGGTGGTAGAAAGTACATAGCAAACGGCGGCAATCAATACCCCACCAATGGTGGCGATAGGGACATTGCGTTTCGGGTTTTTCACCACACCTGCGGCAACGGAGGCGCTTTCCACACCGATGAACGACCACAg                                                                                                     >  1:1273567/1‑139 (MQ=255)                aTCGCGGTGGTAGAAAGTACATAGCAAACGGCGGCAATCAATACCCCACCAATGGTGGCGATAGGGACATTGCGTTTCGGGTTTTTCACCACACCTGCGGCAACGGAGGCGCTTTCCACACCGATGAACGACCACAGCg                                                                                                   <  2:847083/139‑1 (MQ=255)                    cgGTGGTAGAAAGTACATAGCAAACGGCGGCAATCAATACCCCACCAATGGTGGCGATAGGGACATTGCGTTTCGGGTTTTTCACCACACCTGCGGCAACGGAGGCGCTTTCCACACCGATGAAc                                                                                                             <  1:1145787/125‑1 (MQ=255)                    cgGTGGTAGAAAGTACATAGCAAACGGCGGCAATCAATACCCCACCAATGGTGGCGATAGGGACATTGCGTTTCGGGTTTTTCACCACACCTGCGGCAACGGAGGCGCTTTCCACACCGATGAAc                                                                                                             >  2:1145787/1‑125 (MQ=255)                    cgGTGGTAGAAAGTACATAGCAAACGGCGGCAATCAATACCCCACCAATGGTGGCGATAGGGACATTGCGTTTCGGGTTTTTCACCACACCTGCGGCAACGGAGGCGCTTTCCACACCGATGAACGACCACAGCGTAAc                                                                                               >  2:655886/1‑139 (MQ=255)                      gtggtAGAAAGTACATAGCAAACGGCGGCAATCAATACCCCACCAATGGTGGCGATAGGGACATTGCGTTTCGGGTTTTTCACCACACCTGCGGCAACGGAGGCGCTTTCCACACCGa                                                                                                                  <  2:539807/118‑1 (MQ=255)                      gtggtAGAAAGTACATAGCAAACGGCGGCAATCAATACCCCACCAATGGTGGCGATAGGGACATTGCGTTTCGGGTTTTTCACCACACCTGCGGCAACGGAGGCGCTTTCCACACCGa                                                                                                                  >  1:539807/1‑118 (MQ=255)                                 tACATAGCAAACGGCGGCAATCAATACCCCACCAATGGTGGCGATAGGGACATTGCGTTTCGGGTTTTTCACCACACCTGCGGCAACGGAGGCGCTTTCCACACCGATGAACGACCACAGCGTAACGTTAAGGGTACtt                                                                                  >  1:846548/1‑139 (MQ=255)                                   cATAGCAAACGGCGGCAATCAATACCCCACCAATGGTGGCGATAGGGACATTGCGTTTCGGGTTTTTCACCACACCTGCGGCAACGGAGGCGCTTTCCACAcc                                                                                                                    <  2:797256/103‑1 (MQ=255)                                   cATAGCAAACGGCGGCAATCAATACCCCACCAATGGTGGCGATAGGGACATTGCGTTTCGGGTTTTTCACCACACCTGCGGCAACGGAGGCGCTTTCCACAcc                                                                                                                    >  1:797256/1‑103 (MQ=255)                                                    atcaatACCCCACCAATGGTGGCGATAGGGACATTGCGTTTCGGGTTTTTCACCACACCTGCGGCAACGGAGGCGCTTTCCACACCGATGAACGAcc                                                                                                         >  2:787236/1‑97 (MQ=255)                                                    atcaatACCCCACCAATGGTGGCGATAGGGACATTGCGTTTCGGGTTTTTCACCACACCTGCGGCAACGGAGGCGCTTTCCACACCGATGAACGAcc                                                                                                         <  1:787236/97‑1 (MQ=255)                                                              caccaATGGTGGCGATAGGGACATTGCGTTTCGGGTTTTTCACCACACCTGCGGCAACGGAGGCGCTTTCCACACCGATGAACGACCACAGCGTAACGTTAAGGGTACTTTGAATTGCACCGAAGGTGCCCAGGCCGCt                                                     >  2:939807/1‑139 (MQ=255)                                                                   aTGGTGGCGATAGGGACATTGCGTTTCGGGTTTTTCACCACACCTGCGGCAACGGAGGCGCTTTCCacac                                                                                                                     <  1:27312/70‑1 (MQ=255)                                                                   aTGGTGGCGATAGGGACATTGCGTTTCGGGTTTTTCACCACACCTGCGGCAACGGAGGCGCTTTCCacac                                                                                                                     >  2:27312/1‑70 (MQ=255)                                                                   aTGGTGGCGATAGGGACATTGCGTTTCGGGTTTTTCACCACACCTGCGGCAACGGAGGCGCTTTCCACACCGATGAACGACCACAGCGTAACGTTAAGGGTACTTTGAATTGCACCGAAGGTGCCCAGGCCGCTGACGt                                                >  2:383600/1‑139 (MQ=255)                                                                      gtggCGATAGGGACATTGCGTTTCGGGTTTTTCACCACACCTGCGGCAACGGAGGCGCTTTCCACACCGATGAACGACCACAGCGt                                                                                                  >  1:1226111/1‑86 (MQ=255)                                                                      gtggCGATAGGGACATTGCGTTTCGGGTTTTTCACCACACCTGCGGCAACGGAGGCGCTTTCCACACCGATGAACGACCACAGCGt                                                                                                  <  2:1226111/86‑1 (MQ=255)                                                                      gtggCGATAGGGACATTGCGTTTCGGGTTTTTCACCACACCTGCGGCAACGGAGGCGCTTTCCACACCGATGAACGACCACAGCGTAACGTTAAGGGTACTTTGAATTGCACCGAAGGTGCCCAGGCCGCTGACGTTcc                                             >  2:1046533/1‑139 (MQ=255)                                                                         gCGATAGGGACATTGCGTTTCGGGTTTTTCACCACACCTGCGGCAACGGAGGCGCTTTCCACACCGATGAACGACCACAGCGTAACGTTAAGGGTACTTTGAATTGTACCGAAGGTGCACAGGCCGCTGACGTTACATg                                          >  2:383541/1‑139 (MQ=255)                                                                          cGATAGGGACATTGCGTTTCGGGTTTTTCACCACACCTGCGGCAACGGAGGCGCTTTCCACACCGATGAACGACCACAGCGTAACGTTAAGGGTACTTTGAATTGCACCGAAGGTGCCCAGGCCGCTGACGTTCCATgc                                         <  1:383600/139‑1 (MQ=255)                                                                          cGATAGGGACATTGCGTTTCGGGTTTTTCACCACACCTGCGGCAACGGAGGCGCTTTCCACACCGATGAACGACCACAGCGTAACGTTAAGGGTACTTTGAATTGCACCGAAGGTGCCCAGGCCGCTGACGTTCCATgc                                         <  1:1532030/139‑1 (MQ=255)                                                                            aTAGGGACATTGCGTTTCGGGTTTTTCACCACACCTGCGGCAACGGAGGCGCTTTCCACACCGATGaa                                                                                                              <  2:728456/68‑1 (MQ=255)                                                                            aTAGGGACATTGCGTTTCGGGTTTTTCACCACACCTGCGGCAACGGAGGCGCTTTCCACACCGATGaa                                                                                                              >  1:728456/1‑68 (MQ=255)                                                                              aGGGACATTGCGTTTCGGGTTTTTCACCACACCTGCGGCAACGGAGGCGCTTTCCACACCGATGaa                                                                                                              <  1:43216/66‑1 (MQ=255)                                                                              aGGGACATTGCGTTTCGGGTTTTTCACCACACCTGCGGCAACGGAGGCGCTTTCCACACCGATGaa                                                                                                              >  2:43216/1‑66 (MQ=255)                                                                                     ttGCGTTTCGGGTTTTTCACCACACCTGCGGCAACGGAGGCGCTTTCCACACCGATGAACGACCACAGCGTAACGTTAAGGGTACTTTGAATTGCACCGAAGGTGcc                                                              <  1:882283/107‑1 (MQ=255)                                                                                     ttGCGTTTCGGGTTTTTCACCACACCTGCGGCAACGGAGGCGCTTTCCACACCGATGAACGACCACAGCGTAACGTTAAGGGTACTTTGAATTGCACCGAAGGTGcc                                                              >  2:882283/1‑107 (MQ=255)                                                                                      tGCGTTTCGGGTTTTTCACCACACCTGCGGCAACGGAGGCGCGTTCCACACCGATGAACGACCACAGCGTAACGTTAAGGGTACTTTGAATTGCACCGAAGGTGCCCAGGCCGCTGACGTTCCATGCCGCCATATAGGt                             >  1:1001851/1‑139 (MQ=255)                                                                                           ttCGGGTTTTTCACCACACCTGCGGCAACGGAGGCGCTTTCCACACCGATGAACGACCACAGCGTAATGTTAAGGGTACTTTGAATTGCACCGAAGGTGCCCAGGc                                                         >  2:897998/1‑106 (MQ=255)                                                                                           ttCGGGTTTTTCACCACACCTGCGGCAACGGAGGCGCTTTCCACACCGATGAACGACCACAGCGTAACGTTAAGGGTACTTTGAATTGCACCGAAGGTGCCCAGGc                                                         <  1:897998/106‑1 (MQ=255)                                                                                            tCGGGTTTTTCACCACACCTGCGGCAACGGAGGCGCTTTCCACACCGATGAACGACCACAGCGTAACGTTAAGGGTACt                                                                                   <  2:264930/79‑1 (MQ=255)                                                                                            tCGGGTTTTTCACCACACCTGCGGCAACGGAGGCGCTTTCCACACCGATGAACGACCACAGCGTAACGTTAAGGGTACt                                                                                   >  1:264930/1‑79 (MQ=255)                                                                                              gggTTTTTCACCACACCTGCGGCAACGGAGGCGCTTTCCACACCGATTAACGACCACAGCGTAACGTTAAGGGTACTTTGAATTGCACCGAAGGTGCCCAGGCCGCTGACGTTCCATGCCGCCATATAGGTTTcaccac                     <  1:383541/139‑1 (MQ=255)                                                                                                  ttttCACCACACCTGCGGCAACGGAGGCGCTTTCCACACCGATGAACGACCACAGCGTAACGTTAAGGGTACTTTGAATTGCACCGAAGGTGCCCAGGCCGCTGACGTTCCATGCCGCCATATAGGTTTCACCACGgaa                 >  2:375771/1‑139 (MQ=255)                                                                                                     tCACCACACCTGCGGCAACGGAGGCGCTTTCCACACCGATGAACGAc                                                                                                          >  2:688862/1‑47 (MQ=255)                                                                                                     tCACCACACCTGCGGCAACGGAGGCGCTTTCCACACCGATGAACGAc                                                                                                          <  1:688862/47‑1 (MQ=255)                                                                                                       accacACCTGCGGCAACGGAGGCGCTTTCCACACCGATGAACGACCACAGCGTAACGTTAAGGGTACTTTGAATTGCACCGAAGGTGCCCAGGCCGCTGACGTTCCATGCCGCCATATAGGTTTCACCACGgaaccaga            <  1:375771/139‑1 (MQ=255)                                                                                                           cacCTGCGGCAACGGAGGCGCTTTCCACACCGATGAACGACCACAGCGTAACGTTAAGGGTACTTTGAATTGCACCGAAGGTGCCCAGGCCGCTGACGTTCCATGCCGCCATATAGGTTTCACCACGgaaccagaacca        >  2:138949/1‑139 (MQ=255)                                                                                                              cTGCGGCAACGGAGGCGCTTTCCACACCGATGAACGACCACAGCGTAACGTTAAGGGTACTTTGAATTGCACCGAAGGTGCCCAGGCCGCTGACGTTCCATGCCGCCATATAGGTTTCACCACGGAACCAGAACCAGcc     >  2:1495164/1‑139 (MQ=255)                                                                                                                gCGGCAACGGAGGCGCTTTCCACACCGATGAACGACCACAGCGt                                                                                                  <  2:616563/44‑1 (MQ=255)                                                                                                                gCGGCAACGGAGGCGCTTTCCACACCGATGAACGACCACAGCGt                                                                                                  >  1:616563/1‑44 (MQ=255)                                                                                                                 cGGCAACGGAGGCGCTTTCCACACCGATGAACGACCACAGCGTAACGTTAAGGGTACTTTGAATTGCACCGAAGGTGCCCAGGCCGCTGACGTTCCATGCCGCCATATAGGTTTCACCACGGAACCAGAACCAGCCaaa  <  1:655886/139‑1 (MQ=255)                                                                                                                              |                                                                                                                              GGGATCATCCCCATAATCGCGGTGGTAGAAAGTACATAGCAAACGGCGGCAATCAATACCCCACCAATGGTGGCGATAGGGACATTGCGTTTCGGGTTTTTCACCACACCTGCGGCAACGGAGGCACTTTCCACACCGATGAACGACCACAGCGTAACGTTAAGGGTACTTTGAATTGCACCGAAGGTGCCCAGGCCGCTGACGTTCCATGCCGCCATATAGGTTTCACCACGGAACCAGAACCAGCCAAA  >  NC_000913/4336280‑4336530

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: ‑