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

GCATTAGGGATCATCCCCATAATCGCGGTGGTAGAAAGTACATAGCAAACGGCGGCAATCAATACCCCACCAATGGTGGCGATAGGGACATTGCGTTTCGGGTTTTTCACCACACCTGCGGCAACGGAGGCACTTTCCACACCGATGAACGACCACAGCGTAACGTTAAGGGTACTTTGAATTGCACCGAAGGTGCCCAGGCCGCTGACGTTCCATGCCGCCATATAGGTTTCACCACGGAACCAGAACCAGCCAAATACGG  >  NC_000913/4336274‑4336535                                                                                                                                    |                                                                                                                                   gCATTAGGGATCATCCCCATAATCGCGGTGGTAGAAAGTACATAGCAAACGGCGGCAATCAATACCCCACCAATGGTGGCGATAGGGACATTGCGTTTCGGGTTTTTCACCACACCTGCGGCAACGGAGGCGCTTTCca                                                                                                                             >  2:514376/1‑139 (MQ=255)     tAGGGATCATCCCCATAATCGCGGTGGTAGAAAGTACATAGCAAACGGCGGCAATCAATACCCCACCAATGGTGGCGATAGGGACATTGCGTTTCGGGTTTTTCACCACACCTGCGGCAACGGAGGCGCTTTCCACAcc                                                                                                                         >  1:1422616/1‑139 (MQ=255)            catcCCCATAATCGCGGTGGTAGAAAGTACATAGCAAACGGCGGCAATCAATACCCCACCAATGGTGGCGATAGGGACATTGCGTTTCGGGTTTTTCACCACACCTGCGGCAACGGAGGCGCTTTCCACACCGATGAAc                                                                                                                  >  1:109128/1‑139 (MQ=255)            catcCCCATAATCGCGGTGGTAGAAAGTACATAGCAAACGGCGGCAATCAATACCCCACCAATGGTGGCGATAGGGACATTGCGTTTCGGGTTTTTCACCACACCTGCGGCAACGGAGGCGCTTTCCACACCGATGAAc                                                                                                                  >  2:476678/1‑139 (MQ=255)                 ccATAATCGCGGTGGTAGAAAGTACATAGCAAACGGCGGCAATCAATACCCCACCAATGGTGGCGATAGGGACATTGCGTTTCGGGTTTTTCACCACACCTGCGGCAACGGAGGCGCTTTCCACACCGATGAACGACca                                                                                                             <  1:476678/139‑1 (MQ=255)                  cATAATCGCGGTGGTAGAAAGTACATAGCAAACGGCGGCAATCAATACCCCACCAATGGTGGCGATAGGGACATTGCGTTTCGGGTTTTTCACCACACCTGCGGCAACGGAGGCGCTTTCCACACCGATGAACGACcac                                                                                                            >  2:723129/1‑139 (MQ=255)                    tAATCGCGGTGGTAGAAAGTACATAGCAAACGGCGGCAATCAATACCCCACCAATGGTGGCGATAGGGACATTGCGTTTCGGGTTTTTCACCACACCTGCGGCAACGGAGGCGCTTTCCACACCGATGAACGACCACAg                                                                                                          >  1:655185/1‑139 (MQ=255)                      aTCGCGGTGGTAGAAAGTACATAGCAAACGGCGGCAATCAATACCCCACCAATGGTGGCGATAGGGACATTGCGTTTCGGGTTTTTCACCACACCTGCGGCAACGGAGGCGCTTTCCACACCGATGAACGACCACAGCg                                                                                                        >  1:1093093/1‑139 (MQ=255)                       tCGCGGTGGTAGAAAGTACATAGCAAACGGCGGCAATCAATACCCCACCAATGGTGGCGATAGGGACATTGCGTTTCGGGTTTTTCACCACACCTGCGGCAACGGAGGCGCTTTCCACACCGATGAACGACCACAGCGt                                                                                                       >  2:116256/1‑139 (MQ=255)                                 aGAAAGTACATAGCAAACGGCGGCAATCAATACCCCACCAATGGTGGCGATAGGGACATTGCGTTTCGGGTTTTTCACCACACCTGCGGCAACGGAGGCGCTTTCCACACCGATGAACGACCACAGCGTAACGTTAAgg                                                                                             >  2:16807/1‑139 (MQ=255)                                       tACATAGCAAACGGCGGCAATCAATACCCCACCAATGGTGGCGATAGGGACATTGCGTTTCGGGTTTTTCACCACACCTGCGGCAACGGAGGCGCTTTCCACACCGATGAACGAcc                                                                                                              >  1:1399545/1‑116 (MQ=255)                                       tACATAGCAAACGGCGGCAATCAATACCCCACCAATGGTGGCGATAGGGACATTGCGTTTCGGGTTTTTCACCACACCTGCGGCAACGGAGGCGCTTTCCACACCGATGAACGAcc                                                                                                              <  2:1399545/116‑1 (MQ=255)                                       tACATAGCAAACGGCGGCAATCAATACCCCACCAATGGTGGCGATAGGGACATTGCGTTTCGGGTTTTTCACCACACCTGCGGCAACGGAGGCGCTTTCCACACCGATGAACGACCACAGCGTAACGTTAAgg                                                                                             >  1:533937/1‑133 (MQ=255)                                       tACATAGCAAACGGCGGCAATCAATACCCCACCAATGGTGGCGATAGGGACATTGCGTTTCGGGTTTTTCACCACACCTGCGGCAACGGAGGCGCTTTCCACACCGATGAACGACCACAGCGTAACGTTAAgg                                                                                             <  2:533937/133‑1 (MQ=255)                                         cATAGCAAACGGCGGCAATCAATA‑CCCACCAATGGTGGCGATAGGGACATTGCGTTTCGGGTTTTTCACCACACCTGCGGCAACGGAGGCGCTTTCCACACCGATGaa                                                                                                                   >  2:768138/1‑108 (MQ=255)                                         cATAGCAAACGGCGGCAATCAATA‑CCCACCAATGGTGGCGATAGGGACATTGCGTTTCGGGTTTTTCACCACACCTGCGGCAACGGAGGCGCTTTCCACACCGATGaa                                                                                                                   <  1:768138/108‑1 (MQ=255)                                             gCAAACGGCGGCAATCAATACCCCACCAATGGTGGCGATAGGGACATTGCGTTTCGGGTTTTTCACCACACCTGCGGCAACGGAGGCGCTTTCCACACCGATGAACGACcaca                                                                                                           >  2:266092/1‑113 (MQ=255)                                             gCAAACGGCGGCAATCAATACCCCACCAATGGTGGCGATAGGGACATTGCGTTTCGGGTTTTTCACCACACCTGCGGCAACGGAGGCGCTTTCCACACCGATGAACGACcaca                                                                                                           <  1:266092/113‑1 (MQ=255)                                              cAAACGGCGGCAATCAATACCCCACCAATGGTGGCGATAGGGACATTGCGTTTCGGGTTTTTCACCACACCTGCGGCAACGGAGGCGCTTTCCACACCGATGAACGACCACAGCGTAACGTTAAGGGTACTTTGAATTg                                                                                >  1:1388453/1‑139 (MQ=255)                                                 aCGGCGGCAATCAATACCCCACCAATGGTGGCGATAGGGACATTGCGTTTCGGGTTTTTCACCACACCTGCGGCAACGGAGGCGCTTTCCACACCGATGAACGACCACAGCGTAACGTTAAGGGTACTTTGAATTGCAc                                                                             >  2:1456120/1‑139 (MQ=255)                                                   ggcggcAATCAATACCCCACCAATGGTGGCGATAGGGACATTGCGTTTCGGGTTTTTCACCACACCTGCGGCAACGGAGGCGCTTTCCACACCGATGaa                                                                                                                   >  2:345616/1‑99 (MQ=255)                                                   ggcggcAATCAATACCCCACCAATGGTGGCGATAGGGACATTGCGTTTCGGGTTTTTCACCACACCTGCGGCAACGGAGGCGCTTTCCACACCGATGaa                                                                                                                   <  1:345616/99‑1 (MQ=255)                                                    gcggcAATCAATACCCCACCAATGGTGGCGATAGGGACATTGCGTTTCGGGTTTTTCACCACACCTGCGGCAACGGAGGCGCTTTCCACACCGATGAACGACCACAGCGTAACGTTAAGGGTACTTTGAATTGCACCGa                                                                          <  1:512020/139‑1 (MQ=255)                                                    gcggcAATCAATACCCCACCAATGGTGGCGATAGGGACATTGCGTTTCGGGTTTTTCACCACACCTGCGGCAACGGAGGCGCTTTCCACACCGATGAACGACCACAGCGTAACGTTAAGGGTACTTTGAATTGCACCGa                                                                          >  2:687458/1‑139 (MQ=255)                                                     cggcAATCAATACCCCACCAATGGTGGCGATAGGGACATTGCGTTTCGGGTTTTTCACCACACCTGCGGCAACGGAGGCGCTTTCCACACCGATGAACGACCACAGCGt                                                                                                       <  2:716357/109‑1 (MQ=255)                                                     cggcAATCAATACCCCACCAATGGTGGCGATAGGGACATTGCGTTTCGGGTTTTTCACCACACCTGCGGCAACGGAGGCGCTTTCCACACCGATGAACGACCACAGCGt                                                                                                       >  1:716357/1‑109 (MQ=255)                                                     cggcAATCAATACCCCACCAATGGTGGCGATAGGGACATTGCGTTTCGGGTTTTTCACCACACCTGCGGCAACGGAGGCGCTTTCCACACCGATGAACGACCACAGCGTAACGTTAAGGGTACTTTGAATTGCACCGaa                                                                         <  2:1388453/139‑1 (MQ=255)                                                     cggcAATCAATACCCCACCAATGGTGGCGATAGGGACATTGCGTTTCGGGTTTTTCACCACACCTGCGGCAACGGAGGCGCTTTCCACACCGATGAACGACCACAGCGTAACGTTAAGGGTACTTTGAATTGCACCGaa                                                                         <  1:16807/139‑1 (MQ=255)                                                        caatcaatACCCCACCAATGGTGGCGATAGGGACATTGCGTTTCGGGTTTTTCACCACACCTGTGGCAACGGAGGCGCTTTCCACACCGATGAACGACCACAGCGTAACGTTAAGGGTACTTTGAATTGAACCGAAGGt                                                                      >  1:1009700/1‑139 (MQ=255)                                                        caatcaatACCCCACCAATGGTGGCGATAGGGACATTGCGTTTCGGGTTTTTCACCACACCTGCGGCAACGGAGGCGCTTTCCACACCGATGAACGACCACAGCGTAACGTTAAGGGTACTTTGAATTGCACCGAAGGt                                                                      >  1:109418/1‑139 (MQ=255)                                                                    caccaATGGTGGCGATAGGGACATTGCGTTTCGGGTTTTTCACCACACCTGCGGCAACGGAGGCGCTTTCCACACCGATGAACGACCACAGCGTAACGTTAAGGGTACTTTGAATTGCACCGAAGGTGCCCAGGCCGCt                                                          >  2:925189/1‑139 (MQ=255)                                                                     accaATGGTGGCGATAGGGACATTGCGTTTCGGGTTTTTCACCACACCTGCGGCAACGGAGGCGCTTTCCACACCGATGAACGACCACAGCGTAACGTTAAGGGTACTTTGAATTGCACCGAAGGTGCCCAGGCCGCTg                                                         >  2:1117743/1‑139 (MQ=255)                                                                         aTGGTGGCGATAGGGACATTGCGTTTCGGGTTTTTCACCACACCTGCGGCAACGGAGGCGCTTTCCACACCGATGAACGACCACAGCGTAAc                                                                                                    >  1:273179/1‑92 (MQ=255)                                                                         aTGGTGGCGATAGGGACATTGCGTTTCGGGTTTTTCACCACACCTGCGGCAACGGAGGCGCTTTCCACACCGATGAACGACCACAGCGTAAc                                                                                                    <  2:273179/92‑1 (MQ=255)                                                                          tggtggCGATAGGGACATTGCGTTTCGGGTTTTTCACCACACCTGCGGCAACGGAGGCGCTTTCCACACCGATGAACGACCACAGCGTAACGTTAAGGGTACTTTGAATTGCACCGAAGGTGCCCAGGCCGCTGACGtt                                                    <  1:1117743/139‑1 (MQ=255)                                                                                  aTAGGGACATTGCGTTTCGGGTTTTTCACCACACCTGCGGCAACGGAGGCGCTTTCCACACCGATGAACGACCACAGCGTAACGTTAAGGGTACTTTGAATTGCACCGAAGGTGCCCAGGCCGCTGACGTTCCATgccg                                            <  1:687458/139‑1 (MQ=255)                                                                                    aGGGACATTGCGTTTCGGGTTTTTCACCACACCTGCGGCAACGGAGGCGCTTTCCACACCGATGAACGACCACAGCGTAACGTTAAGGGTACTTTGAATTGCACCGAAGGTGCCCAGGCCGCTGACGTTCCATgccgcc                                          >  2:1285986/1‑139 (MQ=255)                                                                                    aGGGACATTGCGTTTCGGGTTTTTCACCACACCTGCGGCAACGGAGGCGCTTTCCACACCGATGAACGACCACAGCGTAACGTTAAGGGTACTTTGAATTGCACCGAAGGTGCCCAGGCCGCTGACGTTCCATgccgcc                                          >  1:1266014/1‑139 (MQ=255)                                                                                         cATTGCGTTTCGGGTTTTTCACCACACCTGCGGCAACGGAGGCGCTTTCCACACCGATGAACGACCACAGCGTAACGTTaa                                                                                               >  2:974125/1‑81 (MQ=255)                                                                                         cATTGCGTTTCGGGTTTTTCACCACACCTGCGGCAACGGAGGCGCTTTCCACACCGATGAACGACCACAGCGTAACGTTaa                                                                                               <  1:974125/81‑1 (MQ=255)                                                                                           ttGCGTTTCGGGTTTTTCACCACACCTGCGGCAACGGAGGCGCTTTCCACACCGATGAACGACCACAGCGTAACGTTAAGGGTACTTTGAATTGCACCGAAGGTGCCCAGGCCGCTGACGTTCCATGCCGCCATATAgg                                   >  2:51866/1‑139 (MQ=255)                                                                                              cGTTTCGGGTTTTTCACCACACCTGCGGCAACGGAGGCGCTTTCCACACCGATGAACGACCACAGCGTAACGTTAAGGGTACTTTGAATTGTACCGAAGGTGCCCAGGCCGCTGACGTTCCATgccg                                            >  1:1003840/1‑127 (MQ=255)                                                                                              cGTTTCGGGTTTTTCACCACACCTGCGGCAACGGAGGCGCTTTCCACACCGATGAACGACCACAGCGTAACGTTAAGGGTACTTTGAATTGCACCGAAGGTGCCCAGGCCGCTGACGTTCCATgccg                                            <  2:1003840/127‑1 (MQ=255)                                                                                                tttCGGGTTTTTCACCACACCTGCGGCAACGGAGGCGCTTTCCACACCGATGAACGACCACAGCGt                                                                                                       <  2:262240/66‑1 (MQ=255)                                                                                                tttCGGGTTTTTCACCACACCTGCGGCAACGGAGGCGCTTTCCACACCGATGAACGACCACAGCGt                                                                                                       >  1:262240/1‑66 (MQ=255)                                                                                                                    cTGCGGCAACGGAGGCGCTTTCCACACCGATGAACGACCACAGCGTAAc                                                                                                    <  2:232148/49‑1 (MQ=255)                                                                                                                    cTGCGGCAACGGAGGCGCTTTCCACACCGATGAACGACCACAGCGTAAc                                                                                                    >  1:232148/1‑49 (MQ=255)                                                                                                                    cTGCGGCAACGGAGGCGCTTTCCACACCGATGAACGACCACAGCGTAACGTTAAGGGTACTTTGAATTGCACCGAAGGTGCCCAGGc                                                              >  2:615911/1‑87 (MQ=255)                                                                                                                    cTGCGGCAACGGAGGCGCTTTCCACACCGATGAACGACCACAGCGTAACGTTAAGGGTACTTTGAATTGCACCGAAGGTGCCCAGGc                                                              <  1:615911/87‑1 (MQ=255)                                                                                                                         gCAACGGAGGCGCTTTCCACACCGATGAACGACcaca                                                                                                           <  2:714198/37‑1 (MQ=255)                                                                                                                         gCAACGGAGGCGCTTTCCACACCGATGAACGACcaca                                                                                                           >  1:714198/1‑37 (MQ=255)                                                                                                                            aCGGAGGCGCTTTCCACACCGATGAACGACCACAGCGTAACGTTAAGGGTACTTTGAATTGCACCGAAGGTGCCCAGGCCGCTGACGTTCCATGCCGCCATATAGGTTTCACCACGGAACCAGAACCAGCCAAATACgg  >  2:682888/1‑139 (MQ=255)                                                                                                                                    |                                                                                                                                   GCATTAGGGATCATCCCCATAATCGCGGTGGTAGAAAGTACATAGCAAACGGCGGCAATCAATACCCCACCAATGGTGGCGATAGGGACATTGCGTTTCGGGTTTTTCACCACACCTGCGGCAACGGAGGCACTTTCCACACCGATGAACGACCACAGCGTAACGTTAAGGGTACTTTGAATTGCACCGAAGGTGCCCAGGCCGCTGACGTTCCATGCCGCCATATAGGTTTCACCACGGAACCAGAACCAGCCAAATACGG  >  NC_000913/4336274‑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: ‑