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

GCGGTGGTAGAAAGTACATAGCAAACGGCGGCAATCAATACCCCACCAATGGTGGCGATAGGGACATTGCGTTTCGGGTTTTTCACCACACCTGCGGCAACGGAGGCACTTTCCACACCGATGAACGACCACAGCGTAACGTTAAGGGTACTTTGAATTGCACCGAAGGTGCCCAGGCCGCTGACGTTCCATGCCGCCATATAGGTTTCACCACGGAACCAGAACCAGCCAAATACGGCA  >  NC_000913/4336298‑4336537                                                                                                            |                                                                                                                                     gcgGTGGTAGAAAGTACATAGCAAACGGCGACAATCAATACCCCACCAATGGTGGCGATAGGGACATTGCGTTTCGGGTTTTTCACCACACCTGCGGCAACGGAGGCGCTTTCCACACCGATGAACGACCACAGCGTaa                                                                                                       >  2:483299/1‑139 (MQ=255)  cgGTGGTAGAAAGTACATAGCAAACGGCGGCAATCAATACCCCACCAATGGTGGCGATAGGGACATTGCGTTTCGGGTTTTTCACCACACCTGCGGCAACGGAGGCGCTTTCCACACCGATGAACGACca                                                                                                               >  2:906445/1‑130 (MQ=255)  cgGTGGTAGAAAGTACATAGCAAACGGCGGCAATCAATACCCCACCAATGGTGGCGATAGGGACATTGCGTTTCGGGTTTTTCACCACACCTGCGGCAACGGAGGCGCTTTCCACACCGATGAACGACca                                                                                                               <  1:906445/130‑1 (MQ=255)  cgGTGGTAGAAAGTACATAGCAAACGGCGGCAATCAATACCCCACCAATGGTGGCGATAGGGACATTGCGTTTCGGGTTTTTCACCACACCTGCGGCAACGGAGGCGCTTTCCACACCGATGAACGACCACAGCGt                                                                                                         >  2:56937/1‑136 (MQ=255)  cgGTGGTAGAAAGTACATAGCAAACGGCGGCAATCAATACCCCACCAATGGTGGCGATAGGGACATTGCGTTTCGGGTTTTTCACCACACCTGCGGCAACGGAGGCGCTTTCCACACCGATGAACGACCACAGCGt                                                                                                         <  1:56937/136‑1 (MQ=255)  cgGTGGTAGAAAGTACATAGCAAACGGCGGCAATCAATACCCCACCAATGGTGGCGATAGGGACATTGCGTTTCGGGTTTTTCACCACACCTGCGGCAACGGAGGCGCTTTCCACACCGATGAACGACCACAGCGTAAc                                                                                                      >  2:359502/1‑139 (MQ=255)  cgGTGGTAGAAAGTACATAGCAAACGGCGGCAATCAATACCCCACCAATGGTGGCGATAGGGACATTGCGTTTCGGGTTTTTCACCACACCTGCGGCAACGGAGGCGCTTTCCACACCGATGAACGACCACAGCGTAAc                                                                                                      <  1:359502/139‑1 (MQ=255)                 cATAGCAAACGGCGGCAATCAATACCCCACCAATGGTGGCGATAGGGACATTGCGTTTCGGGTTTTTCACCACACCTGCGGCAACGGAGGCGCTTTCCACACCGATGAACGACCACAGCGTAACGTTAAGGGTACTTTg                                                                                       >  2:913751/1‑139 (MQ=255)                    aGCAAACGGCGGCAATCAATACCCCACCAATGGTGGCGATAGGGACATTGCGTTTCGGGTTTTTCACCACACCTGCGGCAACGGAGGCGCTTTCCACACCGATGAAc                                                                                                                    >  1:879241/1‑107 (MQ=255)                    aGCAAACGGCGGCAATCAATACCCCACCAATGGTGGCGATAGGGACATTGCGTTTCGGGTTTTTCACCACACCTGCGGCAACGGAGGCGCTTTCCACACCGATGAAc                                                                                                                    <  2:879241/107‑1 (MQ=255)                    aGCAAACGGCGGCAATCAATACCCCACCAATGGTGGCGATAGGGACATTGCGTTTCGGGTTTTTCACCACACCTGCGGCAACGGAGGCGCTTTCCACACCGATGAACGACCACAGCGTAACGTTAAGGGTACTTTGAAt                                                                                    >  1:70738/1‑139 (MQ=255)                     gCAAACGGCGGCAATCAATACCCCACCAATGGTGGCGATAGGGACATTGCGTTTCGGGTTTTTCACCACACCTGCGGCAACGGAGGCGCTTTCCACACCGATGAACGACCACAGCGTAACGTTAAGGGTACTTTGAAtt                                                                                   >  2:1030692/1‑139 (MQ=255)                      cAAACGGCGGCAATCAATACCCCACCAATGGTGGCGATAGGGACATTGCGTTTCGGGTTTTTCACCACACCTGCGGCAACGGAGGCGCTTTCCACACCGATGAACGACCACAGCGTAACGTTAAGGGTACTTTGAATTg                                                                                  >  1:225973/1‑139 (MQ=255)                        aaCGGCGACAATCAATACCCCACCAATGGTGGCGATAGGGACATTGCGTTTCGGGTTTTTCACCACACCTGCGGCAACGGAGGCGCTTTCCACACCGATGAACGACCACAGCGTAACGTTAAGGGTACTTTGAATTGCa                                                                                <  1:483299/139‑1 (MQ=255)                         aCGGCGGCAATCAATACCCCACCAATGGTGGCGATAGGGACATTGCGTTTCGGGTTTTTCACCACACCTGCGGCAACGGAGGCGCTTTCCACACCGATGAACGACCACAGCGTAACGTTAAGGGTACTTTGAATTGCAc                                                                               >  2:125987/1‑139 (MQ=255)                                   tcaatACCCCACCAATGGTGGCGATAGGGACATTGCGTTTCGGGTTTTTCACCACACCTGCGGCAACGGAGGCGCTTTCCACACCGATGAACGACCACAGCGTAACGTTAAGGGTACTTTGAATTGCACCGAAGGTGcc                                                                     <  1:742296/139‑1 (MQ=255)                                          cccACCAATGGTGGCGATAGGGACATTGCGTTTCGGGTTTTTCACCACACCTGCGGCAACGGAGGCGCTTTCCACACCGATGAACGACCACAGCGTAACGTTAAGGGTACTTTGAATTGCACCGAAGGTGCCCAGGCCg                                                              <  1:606220/139‑1 (MQ=255)                                            caccaTTGGTGGGGATAGGGACATTGCGTTTCGGGTTTTTCACCACACCTGCGGCAACGGAGGCGCTTTCCACACCGATGAACGACCACAGCGTAACGTTAAGGGTACTTTGAATTGCAc                                                                               >  2:783592/1‑120 (MQ=255)                                            caccaTTGGTGGGGATAGGGACATTGCGTTTCGGGTTTTTCACCACACCTGCGGCAACGGAGGCGCTTTCCACACCGATGAACGACCACAGCGTAACGTTAAGGGTACTTTGAATTGCAc                                                                               <  1:783592/120‑1 (MQ=255)                                             accaATGGTGGCGATAGGGACATTGCGTTTCGGGTTTTTCACCACACCTGCGGCAACGGAGGCGCTTTCCACACCGATGAACGACCACAGCGTa                                                                                                        >  2:333386/1‑94 (MQ=255)                                             accaATGGTGGCGATAGGGACATTGCGTTTCGGGTTTTTCACCACACCTGCGGCAACGGAGGCGCTTTCCACACCGATGAACGACCACAGCGTa                                                                                                        <  1:333386/94‑1 (MQ=255)                                             accaATGGTGGCGATAGGGACATTGCGTTTCGGGTTTTTCACCACACCTGCGGCAACGGAGGCGCTTTCCACACCGATGAACGACCACAGCGTAACGTTAAGGGTACTTTGAATTGCACCGAAGGTGCCCAGGCCGCTg                                                           >  2:838127/1‑139 (MQ=255)                                                    gtggCGATAGGGACATTGCGTTTCGGGTTTCTCACCACACCTGCGGCAACGGAGGCGCTTTCCACACCGATGAACGACCACAGCGTAACGTTAAGGGTACTTTGAATTGCACCGAAGGTGCCCAGGc                                                                <  2:234867/127‑1 (MQ=255)                                                    gtggCGATAGGGACATTGCGTTTCGGGTTTCTCACCACACCTGCGGCAACGGAGGCGCTTTCCACACCGATGAACGACCACAGCGTAACGTTAAGGGTACTTTGAATTGCACCGAAGGTGCCCAGGc                                                                >  1:234867/1‑127 (MQ=255)                                                        cGATAGGGACATTGCGTTTCGGGTTTTTCACCACACCTGCGGCAACGGAGGCGCTTTCCACACCGATGaa                                                                                                                     >  1:517176/1‑70 (MQ=255)                                                        cGATAGGGACATTGCGTTTCGGGTTTTTCACCACACCTGCGGCAACGGAGGCGCTTTCCACACCGATGaa                                                                                                                     <  2:517176/70‑1 (MQ=255)                                                         gATAGGGACATTGCGTTTCGGGTTTTTCACCACACCTGCGGCAACGGAGGCGCTTTCCACACCGATGAACGACCACAGCGTAACGTTAAGGGTACTTTGAATTGCACCGAAGGTGCCCAGGCCGCTGACGTTCCATgcc                                               <  1:680948/139‑1 (MQ=255)                                                               gACATTGCGTTTCGGGTTTTTCACCACACCTGCGGCAACGGAGGCGCTTTCCACACCGa                                                                                                                         <  1:191719/59‑1 (MQ=255)                                                               gACATTGCGTTTCGGGTTTTTCACCACACCTGCGGCAACGGAGGCGCTTTCCACACCGa                                                                                                                         >  2:191719/1‑59 (MQ=255)                                                               gACATTGCGTTTCGGGTTTTTCACCACACCTGCGGCAACGGAGGCGCTTTCCACACCGATGAACGACCACAGCGTAACGTTAAGGGTACTTTGAATTGCACCGAAGGTGCCCAGGCCGCTGACGTTCCATGCCGCCata                                         >  1:194856/1‑139 (MQ=255)                                                                 aaTTGCGTTTCGGGTTTTTCACCACACCTGCGGCAACGGAGGCGCTTACCACACCGATCAACAACCACAGGGTAACGTTAAGGGTACTTTGCATTGCACCGAAGGTGTCCAGGCCGCTGACGTTCCATGCCGCCAtata                                       <  2:901746/138‑1 (MQ=255)                                                                    tGCGTTTCGGGTTTTTCACCACACCTGCGGCAACGGAGGCGCTTTCCACACCGATGAACGACCACAGCGTAACGTTAAGGGTACTTTGAATTGCACCGAAGGTGCCCAGGCCGCTGACGTTCCATgccg                                              >  1:1045399/1‑129 (MQ=255)                                                                    tGCGTTTCGGGTTTTTCACCACACCTGCGGCAACGGAGGCGCTTTCCACACCGATGAACGACCACAGCGTAACGTTAAGGGTACTTTGAATTGCACCGAAGGTGCCCAGGCCGCTGACGTTCCATgccg                                              <  2:1045399/129‑1 (MQ=255)                                                                       gTTTCGGGTTTTTCACCACACCTGCGGCAACGGAGGCGCTTTCCACACCGATGAACGACCACAGCGTAACGTTAAGGGTACTTTGAATTGCACCGAAGGTGCCCAGGCCGCTGACGTTCCATGCCGCCATATAGGTTTc                                 <  1:1030692/139‑1 (MQ=255)                                                                        tttCGGGTTTTTCACCACACCTGCGGCAACGGAGGCGCTTTCCACACCGATGAACGACCACAGCGTAACGTTAAGGGTACTTTGAATTGCACCGAAGGTGCCCAGGCCGCTGACGTTCCATGCCGCCATATAGGTTTca                                >  1:994017/1‑139 (MQ=255)                                                                                ttttCACCACACCTGCGGCAACGGAGGCGCTTTCCACACCGATGAACGACCACAGCGTAACGTTAAgg                                                                                               <  2:320474/68‑1 (MQ=255)                                                                                ttttCACCACACCTGCGGCAACGGAGGCGCTTTCCACACCGATGAACGACCACAGCGTAACGTTAAgg                                                                                               >  1:320474/1‑68 (MQ=255)                                                                                 tttCACCACACCTGCGGCAACGGAGGCGCTTTCCACACCGATGAACGACCACAGCGTAACGTTAAGGGTACTTTGAATTGCACCGAAGGTGCCCAGGCCGCTGACGTTCCATGCCGCCATATAGGTTTCACCACGgaac                       <  1:125987/139‑1 (MQ=255)                                                                                            cTGCGGCAACGGAGGCGCTTTCCACACCGATGAACGACCACAGCGTAACGTTAAGGGTACTTTGAATTGCACCGAAGGTGCCCAGGCCGCTGACGTTCCATGCCGCCATATAGGTTTCACCACGGAACCAGAACCAGcc            >  1:413954/1‑139 (MQ=255)                                                                                                    aCGGAGGCGCTTTCCACACCGATGAACGACCACAGCGTAACGTTAAGGGTACTTTGAATTGCACCGAAGGTGCCCAGGCCGCTGACGTTCCATGCCGCCATATAGGTTTCACCACGGAACCAGAACCAGCCAAATACgg    >  2:212164/1‑139 (MQ=255)                                                                                                      ggAGGCGCTTTCCACACCGATGAACGACCACAGCGTAACGTTAAGGGTACTTTGAATTGCACCGAAGGTGCCCAGGCCGCTGACGTTCCATGCCGCCATATAGGTTTCACCACGGAACCAGAACCAGCCAAATACGGCa  >  1:515809/1‑139 (MQ=255)                                                                                                      ggAGGCGCTTTCCACACCGATGAACGACCACAGCGTAACGTTAAGGGTACTTTGAATTGCACCGAAGGTGCCCAGGCCGCTGACGTTCCATGCCGCCATATAGGTTTCACCACGGAACCAGAACCAGCCAAATACGGCa  >  1:52200/1‑139 (MQ=255)                                                                                                            |                                                                                                                                     GCGGTGGTAGAAAGTACATAGCAAACGGCGGCAATCAATACCCCACCAATGGTGGCGATAGGGACATTGCGTTTCGGGTTTTTCACCACACCTGCGGCAACGGAGGCACTTTCCACACCGATGAACGACCACAGCGTAACGTTAAGGGTACTTTGAATTGCACCGAAGGTGCCCAGGCCGCTGACGTTCCATGCCGCCATATAGGTTTCACCACGGAACCAGAACCAGCCAAATACGGCA  >  NC_000913/4336298‑4336537

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