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	98.0%	177.2 / NA	50	S209S (AGT→AGC)	adiC	arginine:agmatine antiporter
Reads supporting (aligned to +/- strand):  ref base A (0/0);  major base G (23/26);  minor base T (0/1);  total (23/27)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GCATTAGGGATCATCCCCATAATCGCGGTGGTAGAAAGTACATAGCAAACGGCGGCAATCAATACCCCACCAATGGTGGCGATAGGGACATTGCGTTTCGGGTTTTTCACCACACCTGCGGCAACGGAGGCACTTTCCACACCGATGAACGACCACAGCGTAACGTTAAGGGTACTTTGAATTGCACCGAAGGTGCCCAGGCCGCTGACGTTCCATGCCGCCATATAGGTTTCACCACGGAACCAGAACCAGCCAAATACGGCAAT  >  NC_000913/4336274‑4336539                                                                                                                                    |                                                                                                                                       gCATTAGGGATCATCCCCATAATCGCGGTGGTAGAAAGTACATAGCAAACGGCGGCAATCAATACTCCACCAATGGTGGCGATAGGGACATTGCGTTTCGGCTTTTTCACCACACCTGCGGCAACTGATGCGCTTTCca                                                                                                                                 <  2:1134127/139‑1 (MQ=255)  cATTAGGGATCATCCCCATAATCGCGGTGGTAGAAAGTACATAGCAAACGGCGGCAATCAATACCCCACCAATGGTGGCGATAGGGACATTGCGTTTCGGGTTTTTCACCACACCTGCGGCAACGGAGGCGCTTTCcac                                                                                                                                >  1:1259753/1‑139 (MQ=255)           tcatcCCCATAATCGCGGTGGTAGAAAGTACATAGCAAACGGCGGCAATCAATACCCCACCAATGGTGGCGATAGGGACATTGCGTTTCGGGTTTTTCACCACACCTGCGGCAACGGAGGCGCTTTCCACACCGATGaa                                                                                                                       <  2:242867/139‑1 (MQ=255)            catcCCCATAATCGCGGTGGTAGAAAGTACATAGCAAACGGCGGCAATCAATTCCCCACCAATGGTGGCGATAGGGACATTGCGTTTCGGGTTTTTCACCACACCTGCGGCAACGGAGGCGCTTTCCACACCGATGAAc                                                                                                                      >  1:1143384/1‑139 (MQ=255)              tcCCCATAATCGCGGTGGTAGAAAGTACATAGCAAACGGCGGCAATCAATACCCCACCAATGGTGGCGATAGGGACATTGCGTTTCGGGTTTTTCACCACACCTGCGGCAACGGAGGCGCTTTCCACACCGATGAACGa                                                                                                                    >  2:752311/1‑139 (MQ=255)                 ccATAATCGCGGTGGTAGAAAGTACATAGCAAACGGCGGCAATCAATACCCCACCAATGGTGGCGATAGGGACATTGCGTTTCGGGTTTTTCACCACACCTGCGGCAACGGAGGCGCTTTCCACACCGATGAACGACca                                                                                                                 <  2:1259753/139‑1 (MQ=255)                 ccATAATCGCGGTGGTAGAAAGTACATAGCAAACGGCGGCAATCAATACCCCACCAATGGTGGCGATAGGGACATTGCGTTTCGGGTTTTTCACCACACCTGCGGCAACGGAGGCGCTTTCCACACCGATGAACGACca                                                                                                                 <  2:1143384/139‑1 (MQ=255)                       tCGCGGTGGTAGAAAGTACATAGCAAACGGCGGCAATCAATACCCCACCAATGGTGGCGATAGGGACATTGCGTTTCGGGTTTTTCACCACACCTGCGGCAACGGAGGCGCTTTCCACACCGATGAACGACCACAGCGt                                                                                                           <  1:752311/139‑1 (MQ=255)                          cgGTGGTAGAAAGTACATAGCAAACGGCGGCAATCAATACCCCACCAATGGTGGCGATAGGGACATTGCGTTTCGGGTTTTTCACCACACCTGCGGCAACGGAGGCGCTTTCCACACCGATGAACGACCACAGCGTAAc                                                                                                        <  1:1349199/139‑1 (MQ=255)                                         cATAGCAAACGGCGGCAATCAATACCCCACCAATGGTGGCGATAGGGACATTGCGTTTCGGGTTTTTCACCACACCTGCGGCAACGGAGGCGCTTTCCACACCGATGAACGACCACAGCGTAACGTTaa                                                                                                   <  2:944505/129‑1 (MQ=255)                                         cATAGCAAACGGCGGCAATCAATACCCCACCAATGGTGGCGATAGGGACATTGCGTTTCGGGTTTTTCACCACACCTGCGGCAACGGAGGCGCTTTCCACACCGATGAACGACCACAGCGTAACGTTaa                                                                                                   >  1:944505/1‑129 (MQ=255)                                            aGCAAACGGCGGCAATCAATACCCCACCAATGGTGGCGATAGGGACATTGCGTTTCGGGTTTTTCACCACACCTGCGGCAACGGAGGCGCTTTCCACACCGATGAACGACCACAGCGTAACGTTAAgg                                                                                                 >  1:517795/1‑128 (MQ=255)                                            aGCAAACGGCGGCAATCAATACCCCACCAATGGTGGCGATAGGGACATTGCGTTTCGGGTTTTTCACCACACCTGCGGCAACGGAGGCGCTTTCCACACCGATGAACGACCACAGCGTAACGTTAAgg                                                                                                 <  2:517795/128‑1 (MQ=255)                                              cAAACGGCGGCAATCAATACCCCACCAATGGTGGCGATAGGGACATTGCGTTTCGGGTTTTTCACCACACCTGCGGCAACGGAGGCGCTTTCCACACCGATGAACGACCACAGCGTAACGTTAAGGGTACTTTGAATTg                                                                                    >  1:1570338/1‑139 (MQ=255)                                                      ggcAATCAATACCCCACCAATGGTGGCGATAGGGACATTGCGTTTCGGGTTTTTCACCACACCTGCGGCAACGGAGGCGCTTTCCACACCGATGAACGACCACAGCGTAACGTTAAGGGTACTTTGAATTGCACCGAAg                                                                            >  1:66949/1‑139 (MQ=255)                                                           tcaatACCCCACCAATGGTGGCGATAGGGACATTGCGTTTCGGGTTTTTCACCACACCTGCGGCAACGGAGGCGCTTTCCACACCGATGAACGACCACAGCGTAACGTTAAGGGTACTTTGAATTGCACCGAAGGTGcc                                                                       <  2:66949/139‑1 (MQ=255)                                                             ctatcgccaccaATGGTGGCGATAGGGACATTGCGTTTCGGGTTTTTCACCACACCTGCGGCAACGGAGGCGCTTTCCACACCGATGAACGACCACAGCGTAACGTTAAGGGTACTTTGAATTGCACCGAAGGTGCCCa                                                                     >  1:539724/7‑139 (MQ=255)                                                               tACCCCACCAATGGTGGCGATAGGGACATTGCGTTTCGGGTTTTTCACCACACCTGCGGCAACGGAGGCGCTTTCCACACCGATGAACGACCACAGCGTAACGTTAAGGGTACTTTGAATTGCACCGAAg                                                                            <  1:396645/130‑1 (MQ=255)                                                               tACCCCACCAATGGTGGCGATAGGGACATTGCGTTTCGGGTTTTTCACCACACCTGCGGCAACGGAGGCGCTTTCCACACCGATGAACGACCACAGCGTAACGTTAAGGGTACTTTGAATTGCACCGAAg                                                                            >  2:396645/1‑130 (MQ=255)                                                                 ccccACCAATGGTGGCGATAGGGACATTGCGTTTCGGGTTTTTCACCACACCTGCGGCAACGGAGGCGCTTTCca                                                                                                                                 <  2:455524/75‑1 (MQ=255)                                                                 ccccACCAATGGTGGCGATAGGGACATTGCGTTTCGGGTTTTTCACCACACCTGCGGCAACGGAGGCGCTTTCca                                                                                                                                 >  1:455524/1‑75 (MQ=255)                                                                    caccaATGGTGGCGATAGGGACATTGCGTTTCGGGTTTTTCACCACACCTGCGGCAACGGAGGCGCTTTCCACACCGATGAACGACCACAGCGTAACGTTAAGGGTACTTTGAATTGCACCGAAGGTGcc                                                                       <  1:43090/130‑1 (MQ=255)                                                                    caccaATGGTGGCGATAGGGACATTGCGTTTCGGGTTTTTCACCACACCTGCGGCAACGGAGGCGCTTTCCACACCGATGAACGACCACAGCGTAACGTTAAGGGTACTTTGAATTGCACCGAAGGTGcc                                                                       >  2:43090/1‑130 (MQ=255)                                                                     accaATGGTGGCGATAGGGACATTGCGTTTCGGGTTTTTCACCACACCTGCGGCAACGGAGGCGCTTTCCACACCGa                                                                                                                           <  1:735376/77‑1 (MQ=255)                                                                     accaATGGTGGCGATAGGGACATTGCGTTTCGGGTTTTTCACCACACCTGCGGCAACGGAGGCGCTTTCCACACCGa                                                                                                                           >  2:735376/1‑77 (MQ=255)                                                                     accaATGGTGGCGATAGGGACATTGCGTTTCGGGTTTTTCACCACACCTGCGGCAACGGAGGCGCTTTCCACACCAATGAACGACCACAGCGTAACGTTAAgg                                                                                                 >  1:883751/1‑103 (MQ=255)                                                                     accaATGGTGGCGATAGGGACATTGCGTTTCGGGTTTTTCACCACACCTGCGGCAACGGAGGCGCTTTCCACACCAATGAACGACCACAGCGTAACGTTAAgg                                                                                                 <  2:883751/103‑1 (MQ=255)                                                                      ccaATGGTGGCGATAGGGACATTGCGTTTCGGGTTTTTCACCACACCTGCGGCAACGGAGGCGCTTTCCACACCGATGAACGACCACAGCGTAACGTTAAGGGTACTTTGAATTGCACCGAAGGTGCCCAGGCCGCTGa                                                            <  2:836238/139‑1 (MQ=255)                                                                                  aTAGGGTCATTGCGTTTCGGGTTTTTCACTACACCTGCGGCAACGGAGTCTCTTTCCACACCGATGAACGACCGCAGCATAGCGGTAATGGTACTTTGAATTGCACCGAATGTGCCCAGGCCGCTGACGTTCCATgccg                                                <  2:28457/139‑1 (MQ=255)                                                                                    aGGGACATTGCGTTTCGGGTTTTTCACCACACCTGCGGCAACGGAGGCGCTTTCCACACCGATGAACGACCACAGCGTAACGTTAAGGGTACTTTGAATTGCACCGAAGGTGCCCAGGCCGCTGACGTTCCATgccgcc                                              <  2:539724/139‑1 (MQ=255)                                                                                            tGCGTTTCGGGTTTTTCACCACACCTGCGGCAACGGAGGCGCTTTCCACACCGATGaa                                                                                                                       >  2:935872/1‑58 (MQ=255)                                                                                            tGCGTTTCGGGTTTTTCACCACACCTGCGGCAACGGAGGCGCTTTCCACACCGATGaa                                                                                                                       <  1:935872/58‑1 (MQ=255)                                                                                            tGCGTTTCGGGTTTTTCACCACACCTGCGGCAACGGAGGCGCTTTCCACACCGATGAACGACCACAGCGt                                                                                                           >  1:799340/1‑70 (MQ=255)                                                                                            tGCGTTTCGGGTTTTTCACCACACCTGCGGCAACGGAGGCGCTTTCCACACCGATGAACGACCACAGCGt                                                                                                           <  2:799340/70‑1 (MQ=255)                                                                                            tGCGTTTCGGGTTTTTCACCACACCTGCGGCAACGGAGGCGCTTTCCACACCGATGAACGACCACAGCGTAACGTTAAGGGTACTTTGAATTGCACCGAAGGTGCCCAGGCCGCTGACGTTCCATGCCGCCATATAGGt                                      >  1:896797/1‑139 (MQ=255)                                                                                                tttCGGGTTTTTCACCACACCTGCGGCAGCGGAGGCGCTTTCCACACCGATGAACGACCACAGCGTAACGTTAAGGGTACTTTGAATTGCACCGAAGGTGCCCAGGCCGCTGACGtt                                                        >  1:1048800/1‑117 (MQ=255)                                                                                                tttCGGGTTTTTCACCACACCTGCGGCAACGGAGGCGCTTTCCACACCGATGAACGACCACAGCGTAACGTTAAGGGTACTTTGAATTGCACCGAAGGTGCCCAGGCCGCTGACGtt                                                        <  2:1048800/117‑1 (MQ=255)                                                                                                 ttCGGGTTTTTCACCACACCTGCGGCAACGGAGGCGCTTTCCACACCGATGAACGACCACAGCGTAACGTTAAGGGTACTTTGAATTGCACCGAAGGTGCCCAGGCCGCTGACGTTCCATGCCGCCATATAGGTTTcac                                 >  1:795549/1‑139 (MQ=255)                                                                                                  tCGGGTTTTTCACCACACCTGCGGCAACGGAGGCGCTTTCCACACCGATGAACGACCACAGCGTAACGTTAAgg                                                                                                 >  1:1353783/1‑74 (MQ=255)                                                                                                  tCGGGTTTTTCACCACACCTGCGGCAACGGAGGCGCTTTCCACACCGATGAACGACCACAGCGTAACGTTAAgg                                                                                                 <  2:1353783/74‑1 (MQ=255)                                                                                                   cGGGTTTTTCACCACACCTGCGGCAACGGAGGCGCTTTCCACACCGATGAACGACCACAGCGTAACGTTAAGGGTACTTTGAATTGCACCGAAGGTGCCCAGGCCGCTGACGTTCCATGCCGCCATATAGGTTTcacca                               <  2:795549/139‑1 (MQ=255)                                                                                                    gggTTTTTCACCACACCTGCGGCAACGGAGGCGCTTTCCACACCGATGAACGACCACAGCGTAACGTTAAGGGTACTTTGAATTGCACCGAAGGTGCCCAGGCCGCTGACGTTCCATGCCGCCATATAGGTTTcaccac                              <  2:896797/139‑1 (MQ=255)                                                                                                                 cacCTGCGGCAACGGAGGCGCTTTCCACACCGATGAACGACCACAGCGTAACGTTAAGGGTACTTTGAATTGCACCGAAGGTGCCCAGGCCGCTGACGTTCCATGCCGCCATATAGGTTTCACCACg                             <  2:62212/127‑1 (MQ=255)                                                                                                                 cacCTGCGGCAACGGAGGCGCTTTCCACACCGATGAACGACCACAGCGTAACGTTAAGGGTACTTTGAATTGCACCGAAGGTGCCCAGGCCGCTGACGTTCCATGCCGCCATATAGGTTTCACCACg                             >  1:62212/1‑127 (MQ=255)                                                                                                                  acCTGCGGCAACGGAGGCGCTTTCCACACCGATGAACGACCACAGCGTAACGTTAAGGGTACTTTGAATTGCACCGAAGGTGCCCAGGCCGCTGACGTTCCATGCCGCCAtata                                         >  2:837457/1‑114 (MQ=255)                                                                                                                  acCTGCGGAAACGGAGGCGCTTTCCACACCGATGAACGACCACAGCGTAACGTTAAGGGTACTTTGAATTGCACCGAAGGTGCCCAGGCCGCTGACGTTCCATGCCGCCAtata                                         <  1:837457/114‑1 (MQ=255)                                                                                                                       cGGCAACGGAGGCGCTTTCCACACCGATGAACGACCACAGCGTAACGTTAAGGGTACTTTGAATTGCACCGAAGGTGCCCAg                                                                    <  2:1629029/82‑1 (MQ=255)                                                                                                                       cGGCAACGGAGGCGCTTTCCACACCGATGAACGACCACAGCGTAACGTTAAGGGTACTTTGAATTGCACCGAAGGTGCCCAg                                                                    >  1:1629029/1‑82 (MQ=255)                                                                                                                              ggAGGCGCTTTCCACACCGATGAACGACCACAGCGTAACGTTAAGGGTACTTTGAATTGCACCGAAGGTGCCCAGGCCGCTGACGTTCCATGCCGCCATATAGGTTTCACCACGGAACCAGAACCAGCCAAATACGGCa    <  2:1570338/139‑1 (MQ=255)                                                                                                                                aGGCGCTTTCCACACCGATGAACGACCACAGCGTAACGTTAAGGGTACTTTGAATTGCACCGAAGGTGCCCAGGCCGCTGACGTTCCATGCCGCCATATAGGTTTCACCACGGAACCAGAACCAGCCAAATACGGCAAt  >  1:1662181/1‑139 (MQ=255)                                                                                                                                    |                                                                                                                                       GCATTAGGGATCATCCCCATAATCGCGGTGGTAGAAAGTACATAGCAAACGGCGGCAATCAATACCCCACCAATGGTGGCGATAGGGACATTGCGTTTCGGGTTTTTCACCACACCTGCGGCAACGGAGGCACTTTCCACACCGATGAACGACCACAGCGTAACGTTAAGGGTACTTTGAATTGCACCGAAGGTGCCCAGGCCGCTGACGTTCCATGCCGCCATATAGGTTTCACCACGGAACCAGAACCAGCCAAATACGGCAAT  >  NC_000913/4336274‑4336539

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