Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_007793	2,641,715	G→A	S21S (AGC→AGT)	SAUSA300_RS13565 ←	GTP pyrophosphokinase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_007793	2,641,715	0	G	A	100.0%	247.9 / NA	87	S21S (AGC→AGT)	SAUSA300_RS13565	GTP pyrophosphokinase
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (2/85);  total (2/85)

TGACGAATAAATATGATCTAACTCTACGAAACCTAATAACGTATCAAATGCTTCATCACCATTTTCAAATTTACTTAAGCTATTTTTAAAATCATGTCGCAAATCCTCTAAATATAGTGATGGTTTTCGATCTACATACATTTTTATACTAACCTCCGATATATAATCATCTTTA  >  NC_007793/2641637‑2641811                                                                               |                                                                                                 tGACGAATAAATATGATCTAACTCTACGAAACCTAATAACGTATCAAATGCTTCATCCCCATTTTCAAATTTACTTAAACTATTTTTAAAATCATGTCGc                                                                             <  1:547821/100‑1 (MQ=255) tGACGAATAAATATGATCTAACTCTACGAAACCTAATAACGTATCAAATGCTTCATCACCATTTTCAAATTTACTTAAACTATTTTTAAAATCATGTCgt                                                                             <  1:2986560/100‑2 (MQ=255) tGACGAATAAATATGATCTAACTCTACGAAACCTAATAACGTATCAAATGCTTCATCACCATTTTCAAATTTACTTAAACTATTTTTAAAATCATGTCGc                                                                             <  1:988582/100‑1 (MQ=255) tGACGAATAAATATGATCTAACTCTACGAAACCTAATAACGTATCAAATGCTTCATCACCATTTTCAAATTTACTTAAACTATTTTTAAAATCATGTCGc                                                                             <  1:844280/100‑1 (MQ=255) tGACGAATAAATATGATCTAACTCTACGAAACCTAATAACGTATCAAATGCTTCATCACCATTTTCAAATTTACTTAAACTATTTTTAAAATCATGTCGc                                                                             <  1:5091475/100‑1 (MQ=255)  gccGAATAAATATGATCTAACTCTACGAAACCTAATAACGTATCAAATGCTTCATCACCATTTTCAAATTTACTTAAACTATTTTTAAAATCATGTCGCa                                                                            <  1:3262887/98‑1 (MQ=255)  gccGAATAAATATGATCTAACTCTACGAAACCTAATAACGTATCAAATGCTTCATCACCATTTTCAAATTTACTTAAACTATTTTTAAAATCATGTCGCa                                                                            <  1:1658706/98‑1 (MQ=255)  gatgAATAAATATGATCTAACTCTACGAAACCTAATAACGTATCAAATGCTTCATCACCATTTTCAAATTTACTTAAACTATTTTTAAAATCATgcgagc                                                                            >  1:812743/4‑95 (MQ=255)  gACGAATAAATATGATCTAACTCTACGAAACCTAATAACGTATCAAATGCTTCGTCACCATTTTCAAATTTACTTAAACTATTTTTAAAATCATGTCGCa                                                                            <  1:626191/100‑1 (MQ=255)  gACGAATAAATATGATCTAACTCTACGAAACCTAATAACGTATCAAATGCTTCATCACCATTTTCAAATTTACTTAAACTATTTTTAAAATCATGTCGCa                                                                            <  1:4044354/100‑1 (MQ=255)  gACGAATAAATATGATCTAACTCTACGAAACCTAATAACGTATCAAATGCTTCATCACCATTTTCAAATTTACTTAAACTATTTTTAAAATCATGTCGCa                                                                            <  1:5478123/100‑1 (MQ=255)  gACGAATAAATATGATCTAACTCTACGAAACCTAATAACGTATCAAATGCTTCATCACCATTTTCAAATTTACTTAAACTATTTTTAAAATCATGTCGCa                                                                            <  1:3130185/100‑1 (MQ=255)  gACGAATAAATATGATCTAACTCTACGAAACCTAATAACGTATCAAATGCTTCATCACCATTTTCAAATTTACTTAAACTATTTTTAAAATCATGTCGCa                                                                            <  1:3114623/100‑1 (MQ=255)  gACGAATAAATATGATCTAACTCTACGAAACCTAATAACGTATCAAATGCTTCATCACCATTTTCAAATTTACTTAAACTATTTTTAAAATCATGTCGCa                                                                            <  1:1459522/100‑1 (MQ=255)  gACGAATAAATATGATCTAACTCTACGAAACCTAATAACGTATCAAATGCTTCATCACCATATTCAAATTTACTTAAACTATTTTTAAAATCATGTCGCa                                                                            <  1:5457626/100‑1 (MQ=255)   aCGAATAAATATGATCTAACTCTACGAAACCTAATAACGTATCAAATGCTTCATCACCATTTTCAAATTTACTTAAACTATTTTTAAAATCATGTCGCag                                                                           <  1:3822215/100‑2 (MQ=255)   aCGAATAAATATGATCTAACTCTACGAAACCTAATAACGTATCAAATGCTTCATCACCATTTTCAAATTTACTTAAACTATTTTTAAAATCATGTCGCag                                                                           <  1:5101777/100‑2 (MQ=255)   aCGAATAAATATGATCTAACTCTACGAAACCTAATAACGTATCAAATGCTTCATCACCATTTTCAAATTTACTTAAACTATTTTTAAAATCATGTCGCaa                                                                           <  1:2272943/100‑1 (MQ=255)   aCGAATAAATATGATCTAACTCTACGAAACCTAATAACGTATCAAATGCTTCATCACCATTTTCAAATTTACTTAAACTATTTTTAAAATCATGTCGCaa                                                                           <  1:2610060/100‑1 (MQ=255)   aCGAATAAATATGATCTAACTCTACGAAACCTAATAACGTATCAAATGCTTCATCACCATTTTCAAATTTACTTAAACTATTTTTAAAATCATGTCGCaa                                                                           <  1:4996157/100‑1 (MQ=255)   aCGAATAAATATGATCTAACTCTACGAAACCTAATAACGTATCAAATGCTTCATCACCATTTTCAAATTTACTTAAACTATTTTTAAAATCATGTCGCaa                                                                           <  1:1951859/100‑1 (MQ=255)   aCGAATAAATATGATCTAACTCTACGAAACCTAATAACGTATCAAATGCTTCATCACCATTTTCAAATTTACTTAAACTATTTTTAAAATCATGTCGCaa                                                                           <  1:1351775/100‑1 (MQ=255)     gAATAAATATGATCTAACTCTACGAAACCTAATAACGTATCAAATGCTTCATCACCATTTTCAAATTTACTTAAACTATTTTTAAAATCATGTCGCaaga                                                                         <  1:4525273/100‑3 (MQ=255)        taaataTGATCTAACTCTACGAAACCTAATAACGTATCAAATGCTTCATCACCATTTTCAAATTTACTTAAACTATTTTTAAAATCATGTCGCAAATCct                                                                      <  1:1888036/100‑1 (MQ=255)         aaataTGATCTAACTCTACGAAACCTAATATCGTATCAAATGCTTCATCACCATTTTCAAATTTACTTAAACTATTTTTAAAATCATGTCGCAAATCctc                                                                     <  1:601975/100‑1 (MQ=255)          cttccgATCTGACTCTACGAAACCTAATAACGTATCAAATGCTTCATCACCATTTTCAAATTTACTTAAACTATTTTTAAAATCATGTCGCAAATCctct                                                                    <  1:5356083/95‑1 (MQ=255)          aataTGATCTAACTCTACGAAACCTAATAACGTATCAAATGCTTCATCACCATTTTCAAATTTACTTAAACTATTTTTAAAATCATGTCGCAAATCctct                                                                    <  1:5572484/100‑1 (MQ=255)          aataTGATCTAACTCTACGAAACCTAATAACGTATCAAATGCTTCATCACCATTTTCAAATTTACTTAAACTATTTTTAAAATCATGTCGCAAATCctct                                                                    <  1:3738972/100‑1 (MQ=255)          aataTGATCTAACTCTACGAAACCTAATAACGTATCAAATGCTTCATCACCATTTTCAAATTTACTTAAACTATTTTTAAAATCATGTCGCAAATCctca                                                                    <  1:5648240/100‑2 (MQ=255)           ttccgATCTAACTCTACGAAACCTAATAACGTATCAAATGCTTCATCACCATTTTCAAATTTACTTAAACTATTTTTAAAATCATGTCGCAAATCCTCTa                                                                   <  1:2936352/96‑1 (MQ=255)           ctccgATCTAACTCTACGAAACCTAATAACGTATCAAATGCTTCATCACCATTTTCAAATTTACTTAAACTATTTTTAAAATCATGTCGCAAATCCTCTa                                                                   <  1:5754815/96‑1 (MQ=255)           atatGATCTAACTCTACGAAACCTAATAACGTATCAAATGCTTCATCACCATTTTCAAATTTACTTAAACTATTTTTAAAATCATGTCGCAAATCctgct                                                                   <  1:11572/100‑4 (MQ=255)           atatGATCTAACTCTACGAAACCTAATAACGTATCAAATGCTTCATCACCATTTTCAAATTTACTTAAACTATTTTTAAAATCATGTCGCAAATCctctg                                                                   <  1:2698551/100‑2 (MQ=255)           atatGATCTAACTCTACGAAACCTAATAACGTATCAAATGCTTCATCACCATTTTCAAATTTACTTAAACTATTTTTAAAATCATGTCGCAAATCctctg                                                                   <  1:2890453/100‑2 (MQ=255)           atatGATCTAACTCTACGAAACCTAATAACGTATCAAATGCTTCATCACCATTTTCAAATTTACTTAAACTATTTTTAAAATCATGTCGCAAATCCTCTa                                                                   <  1:963341/100‑1 (MQ=255)           atatGATCTAACTCTACGAAACCTAATAACGTATCAAATGCTTCATCACCATTTTCAAATTTACTTAAACTATTTTTAAAAGCATGTCGCAAATCCTCTa                                                                   <  1:2533231/100‑1 (MQ=255)            tatGATCTAACTCTACGAAACCTAATAACGTATCAAATGCTTCATCACCATTTTCAGATTTACTTAAACTATTTTTAAAATCATGTCGCAAATCCTCTaa                                                                  <  1:2052065/100‑1 (MQ=255)             ttccgatctctctACGAAACCTAATAACGTATCAAATGCTTCATCACCATTTTCAAATTTACTTAAACTATTTTTAAAATCATGTCGCAAATCCTCTaag                                                                 <  1:1557234/91‑2 (MQ=255)                aTCTGACTCTACGAAACCTAATAACGTATCAAATGCTTCATCACCATTTTCAAATTTACTTAAACTATTTTTAAAATCATGTCGCAGATCCTCTAAatat                                                              <  1:3050519/100‑1 (MQ=255)                aTCTAACTCTACGAAACCTAATAACGTATCAAATGCTTCATCACCATTTTCAAATTTACTTAAACTATTTTTAAAATCATGTCGCAAATCCTCTAAatat                                                              <  1:4277395/100‑1 (MQ=255)                aTCTAACTCTACGAAACCTAATAACGTATCAAATGCTTCATCACCATTTTCAAATTTACTTAAACTATTTTTAAAATCATGTCGCAAATCCTCTAAatat                                                              <  1:4252789/100‑1 (MQ=255)                aTCTAACTCTACGAAACCTAATAACGTATCAAATGCTTCATCACCATTTTCAAATTTACTTAAACTATTTTTAAAATCATGTCGCAAATCCTCTAAatat                                                              <  1:4463247/100‑1 (MQ=255)                aTCTAACTCTACGAAACCTAATAACGTATCAAATGCTTCATCACCATTTTCAAATTTACTTAAACTATTTTTAAAATCATGTCGCAAATCCTCTAAatat                                                              <  1:4657804/100‑1 (MQ=255)                aTCTAACTCTACGAAACCTAATAACGTATCAAATGCTTCATCACCATTTTCAAATTTACTTAAACTATTTTTAAAATCATGTCGCAAATCCTCTAAatat                                                              <  1:1892102/100‑1 (MQ=255)                aTCTAACTCTACGAAACCTAATAACGTATCAAATGCTTCATCACCATTTTCAAATTTACTTAAACTATTTTTAAAATCATGTCGCAAATCCTCTAAatat                                                              <  1:1769584/100‑1 (MQ=255)                aTCTAACTCTACGAAACCTAATAACGTATCAAATGCTTCATCACCATTTTCAAATTTACTTAAACTATTTTTAAAATCATGTCGCAAATCCTCTAAatat                                                              <  1:2382514/100‑1 (MQ=255)                aTCTAACTCTACGAAACCTAATAACGTATCAAATGCTTCATCACCATTTTCAAATTTACTTAAACTATTTTTAAAATCATGCCGCAAATCCTCTAAatat                                                              <  1:4736220/100‑1 (MQ=255)                aTCTAACTCTACGAAACCTAATAACGTATCAAATGCGTCATCTCCATTTTCAAATTTACTTAAACTATTTTTAAAATCATGTCGCAAATCCTCTAAatat                                                              <  1:4551952/100‑1 (MQ=255)                 tCTAACTCTACGAAACCTAATAACGTATCAAATGCTTCATCACCATTTTCAAATTTACTTAAACTATTTTTAAAATCATGTCGCAAATCCTCTAAatacg                                                             <  1:5130277/100‑3 (MQ=255)                   tAACTCTACGAAACCTAATAACGTATCAAATGCTTCATCACCATTTTCAAATTTACTTAAACTATTTTTAAAATCATGTCGCAAATCCTCTAAATATAGt                                                           <  1:4730520/100‑1 (MQ=255)                   tAACTCTACGAAACCTAATAACGTATCAAATGCTTCATCACCATTTTCAAATTTACTTAAACTATTTTTAAAATCATGTCGCAAATCCTCTAAATATAGt                                                           <  1:2908703/100‑1 (MQ=255)                    aaCTCTACGAAAGCTAATAACGTATCAAATGCTTCATCACCATTTTCAAATTTACTTAAACTATTTTTAAAATCATGTCGCAAATCCTCTAAATATAGTg                                                          <  1:4878611/100‑1 (MQ=255)                    aaCTCTACGAAACCTAATAACGTATCAAATGCTTCATCACCATTTTCAAATTTACTTAAACTATTTTTAAAATCATGTCGCAAATCCTCTAAATATAGtc                                                          <  1:2856939/100‑2 (MQ=255)                    aaCTCTACGAAACCTAATAACGTATCAAATGCTTCATCACCATTTTCAAATTTACTTAAACTATTTTTAAAATCATGTCGCAAATCCTCTAAATATAGTg                                                          <  1:215153/100‑1 (MQ=255)                    aaCTCTACGAAACCTAATAACGTATCAAATGCTTCATCACCATTTTCAAATTTACTTAAACTATTTTTAAAATCATGTCGCAAATCCTCTAAATATAGTg                                                          <  1:2596418/100‑1 (MQ=255)                    aaCTCTACGAAACCTAATAACGTATCAAATGCTTCATCACCATTTTCAAATTTACTTAAACTATTTTTAAAATCATGTCGCAAATCCTCTAAATATAGTg                                                          <  1:1306609/100‑1 (MQ=255)                    aaCTCTACGAAACCTAATAACGTATCAAATGCTTCATCACCATTTTCAAATTTACTTAAACTATTTTTAAAATCATGTCGCAAATCCTCTAAATATAGTg                                                          <  1:4755111/100‑1 (MQ=255)                    aaCTCTACGAAACCTAATAACGTATCAAATGCTTCATCACCATTTTCAAATTTACTTAAACTATTTTTAAAATCATGTCGCAAATCCTCTAAATATAGTg                                                          <  1:3895243/100‑1 (MQ=255)                    aaCTCTACGAAACCTAATAACGTATCAAATGCTTCATCACCATTTTCAAATTTACTTAAACTATTTTTAAAATCATGTCGCAAATCCTCTAAATATAGTg                                                          <  1:4215715/100‑1 (MQ=255)                       tctACGAAACCTAATAACGTATCAAATGCTTCATCCCCATTTTCAAATTTACTTAAACTATTTTTAAAATCATGTCGCAAATCCTCTAAATATAGTGagc                                                       <  1:1964920/100‑3 (MQ=255)                       tctACGAAACCTAATAACGTATCAAATGCTTCATCACCATTTTCAAATTTACTTAAACTATTTTTAAAATCATGTCGCAAATCCTCTAAATATAGTgcag                                                       <  1:939722/100‑4 (MQ=255)                        ctACGAAACCTAATAACGTATCAAATGCTTCATCACCATTTTCAAATTTACTTAAACTATTTTTAAAATCATGTCGCAAATCCTCTAAATATAGTGaggt                                                      <  1:1689697/100‑4 (MQ=255)                          aCGAAACCTAATAACGTATCAAATGCTTCATCACCATTTTCAAATTTACTTAAACTATTTTTAAAATCATGTCGCAAATCCTCTAAATATAGTGATGGtt                                                    <  1:4625638/100‑1 (MQ=255)                           cGAAACCTAATAACGTATCAAATGCTTCATCACCATTTTCAAATTTACTTAAACTATTTTTAAAATCATGTCGCAAATCCTCTAAATATAGTGATGGttg                                                   <  1:3802700/100‑2 (MQ=255)                            gAAACCTAATAACGTATCAAATGCTTCATCACCATTTTCAAATTTACTTAAACTATTTTTAAAATCATGTCGCAAATCCTCTAAATATAGTGATGGttta                                                  <  1:4763654/100‑2 (MQ=255)                            gAAACCTAATAACGTATCAAATGCTTCATCACCATTTTCAAATTTACTTAAACTATTTTTAAAATCATGTCGCAAATCCTCTAAATATAGTGATGGttta                                                  <  1:121266/100‑2 (MQ=255)                            gAAACCTAATAACGTATCAAATGCTTCATCACCATTTTCAAATTTACTTAAACTATTTTTAAAATCATGTCGCAAATCCTCTAAATATAGTGATGGtaca                                                  <  1:2248515/100‑4 (MQ=255)                                        cGTATCAAATGCTTCATCACCATTTTCAAATTTACTTAAACTATTTTTAAAATCATGTCGCAAATCCTCTAAATATAGTGATGGTTTTCGATCtacatac                                      <  1:2154830/100‑1 (MQ=255)                                        cGTATCAAATGCTTCATCACCATTTTCAAATTTACTTAAACTATTTTTAAAATCATGTCGCAAATCCTCTAAATATAGTGATGGTTTTCGATCtacatac                                      <  1:584052/100‑1 (MQ=255)                                        cGTATCAAATGCTTCATCACCATTTTCAAATTCACTTAAACTATTTTTAAAATCATGTCGCAAATCCTCTAAATATAGTGATGGTTTTCGATCtacatac                                      <  1:4585738/100‑1 (MQ=255)                                         gTATCAAATGCTTCATCACCATTTTCAGATTTACTTAAACTATTTTTAAAATCATGTCGCAAATCCTCTAAATATAGTGATGGTTTTCGATCtacataca                                     <  1:3155729/100‑1 (MQ=255)                                         gTATCAAATGCTTCATCACCATTTTCAAATTTACTTAAACTATTTTTAAAATCATGTCGCAAATCCTCTAAATATAGTGATGGTTTTCGATCtacataca                                     <  1:3562778/100‑1 (MQ=255)                                          tATCAAATGCTTCATCACCATTTTCAAATTTACTTAAACTATTTTTAAAATCATGTCGCAAATCCTCTAAATATAGTGATGGTTTTCGATCTacatacat                                    <  1:3958256/100‑1 (MQ=255)                                          tATCAAATGCTTCATCACCATTTTCAAATTTACTTAAACTATTTTTAAAATCATGTCGCAAATCCTCTAAATATAGTGATGGTTTTCGATCTacatacat                                    <  1:2033882/100‑1 (MQ=255)                                           aTCAAATGCTTCATCACCATTTTCAAATTTACTTAAACTATTTTTAAAATCATGTCGCAAATCCTCTAAATATAGTGATGGTTTTCGATCTacatacatg                                   <  1:5750773/100‑2 (MQ=255)                                            tCAAATGCTTCATCACCATTTTCAAATTTACTTAAACTATTTTTAAAATCATGTCGCAAATCCTCTAAATATAGTGATGGTTTTCGATCTACATACAttt                                  <  1:2328823/100‑1 (MQ=255)                                               aaTGCTTCATCACCATTTTCAAATTTACTTAAACTATTTTTAAAATCATGTCGCAAATCCTCTAAATATAGTGATGGTTTTCGATCTACATACAtttttt                               <  1:1149803/100‑2 (MQ=255)                                                   cTTCATCACCATTTTCAAATTTACTTAAATTATTTTTAAAATCATGTCGCAAATCCTCTAAATATAGTGATGGTTTTCGATCTACATACATTTTTATACt                           <  1:177374/100‑1 (MQ=255)                                                   cTTCATCACCATTTTCAAATTTACTTAAACTATTTTTAAAATCATGTCGCAAATCCTCTAAATATAGTGATGGTTTTCGATCTACATACATTTTTATACt                           <  1:3575206/100‑1 (MQ=255)                                                   cTTCATCACCATTTTCAAATTTACTTAAACTATTTTTAAAATCATGTCGCAAATCCTCTAAATATAGTGATGGTTTTCGATCTACATACATTTTTATACt                           <  1:1318480/100‑1 (MQ=255)                                                    ttCATCACCATTTTCAAATTTACTTAAACTATTTTTAAAATCATGTCGCAAATCCTCTAAATATAGTGATGGTTTTCGATCTACATACATTTTTATACTa                          <  1:1840882/100‑1 (MQ=255)                                                    ttCATCACCATTTTCAAATTTACTTAAACTATTTTTAAAATCATGTCGCAAATCCTCTAAATATAGTGATGGTTTTCGATCTACATACATTTTTATACTa                          <  1:4943697/100‑1 (MQ=255)                                                      catcaCCATTTTCAAATTTACTTAAACTATTTTTAAAATCATGTCGCAAATCCTCTAAATATAGTGATGGTTTTCGATCTACATACATTTTTATACTAAc                        <  1:3422338/100‑1 (MQ=255)                                                                 tCAAATTTACTTAAACTATTTTTAAAATCATGTCGCAAATCCTCTAAATATAGTGATGGTTTTCGATCTACATACATTTTTATACTAACCTCCGatatat             <  1:3516472/100‑1 (MQ=255)                                                                 tCAAATTTACTTAAACTATTTTTAAAATCATGTCGCAAATCCTCTAAATATAGTGATGGTTTTCGATCTACATACATTTTTATACTAACCTCCGatatat             <  1:5255495/100‑1 (MQ=255)                                                                   aaaTTTACTTAAACTATTTTTAAAATCATGTCGCAAATCCTCTAAATATAGTGATGGTTTTCGATCTACATACATTTTTATACTAACCTCCGATATATaa           <  1:5627929/100‑1 (MQ=255)                                                                            tAAACTATTTTTAAAATCATGTCGCAAATCCTCTAAATATAGTGATGGTTTTCGATCTACATACATTTTTATACTAACCTCCGATATATAATCATCTTTa  >  1:4137551/1‑100 (MQ=255)                                                                               |                                                                                                 TGACGAATAAATATGATCTAACTCTACGAAACCTAATAACGTATCAAATGCTTCATCACCATTTTCAAATTTACTTAAGCTATTTTTAAAATCATGTCGCAAATCCTCTAAATATAGTGATGGTTTTCGATCTACATACATTTTTATACTAACCTCCGATATATAATCATCTTTA  >  NC_007793/2641637‑2641811

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