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%	151.7 / NA	56	S21S (AGC→AGT)	SAUSA300_RS13565	GTP pyrophosphokinase
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (1/55);  total (1/55)

GCTGACGAATAAATATGATCTAACTCTACGAAACCTAATAACGTATCAAATGCTTCATCACCATTTTCAAATTTACTTAAGCTATTTTTAAAATCATGTCGCAAATCCTCTAAATATAGTGATGGTTTTCGATCTACATACATTTTTATACTAACCTCCGATATATAATCATC  >  NC_007793/2641635‑2641807                                                                                 |                                                                                             gCTGACGAATAAATATGATCTAACTCTACGAAACCTAATAACGTATCAAATGCTTCATCACCATTTTCAAATTTACTTAAACTATTTTTAAAATCATGTc                                                                           <  1:5576744/100‑1 (MQ=255)   tGACGAATAAATATGATCTAACTCTACGAAACCTAATAACGTATCAAATGCTTCATCGCCATTTTCAAATTTACTTAAACTATTTTTAAAATCATGTCgt                                                                         <  1:1713712/100‑2 (MQ=255)   tGACGAATAAATATGATCTAACTCTACGAAACCTAATAACGTATCAAATGCTTCATCACCATTTTCAAATTTACTTAAACTATTTTTAAAATCATGTCgg                                                                         <  1:47963/100‑2 (MQ=255)   tGACGAATAAATATGATCTAACTCTACGAAACCTAATAACGTATCAAATGCTTCATCACCATTTTCAAATTTACTTAAACTATTTTTAAAATCATGTCGc                                                                         <  1:1025082/100‑1 (MQ=255)   tGACGAATAAATATGATCTAACTCTACGAAACCTAATAACGTATCAAATGCTTCATCACCATTTTCAAATTTACTTAAACTATTTTTAAAATCATGTCGc                                                                         <  1:1628569/100‑1 (MQ=255)   tGACGAATAAATATGATCTAACTCTACGAAACCTAATAACGTATCAAATGCTTCATCACCATTTTCAAATTTACTTAAACTATTTTTAAAATCATGTCGc                                                                         <  1:5273623/100‑1 (MQ=255)    gacgacTAAATATGATCTAACTCTACGAAACCTAATAACGTATCAAATGCTTCATCACCATTTTCAAATTTACTTAAACTATTTTTAAAATCATGTCGCa                                                                        >  1:3452006/1‑100 (MQ=255)     gcGAATAAATATGATCTAACTCTACGAAACCTAATAACGTATCAAATGCTTCATCACCATTTTCAAATTTACTTAAACTATTTTTAAAATCATGTCGCaa                                                                       <  1:3424902/99‑1 (MQ=255)     aCGAATAAATATGATCTAACTCTACGAAACCTAATAGCGTATCAAATGCTTCATCACCATTTTCAAATTTACTTAAACTATTTTTAAATTCATGTCGCaa                                                                       <  1:6367645/100‑1 (MQ=255)     aCGAATAAATATGATCTAACTCTACGAAACCTAATAACGTATCAAATGCTTCATCACCATTTTCAAATTTACTTAAACTATTTTTAAAATCATGTCGCaa                                                                       <  1:2702031/100‑1 (MQ=255)          taaataTGATCTAACTCTACGAAACCTAATAACGTATCAAATGCTTCATCACCATTTTCAAATTTACTTAAACTATTTTTAAAATCATGTCGCAAATCct                                                                  <  1:2963796/100‑1 (MQ=255)            aataTGATCTAACTCTACGAAACCTAATAACGTATCAAATGCTTCATCACCATTGTCAAATTTACTTAAACTATTTTTAAAATCATGTCGCAAATCctca                                                                <  1:5556897/100‑2 (MQ=255)                  aTCTAACTCTACGAAACCTAATAACGTATCAAATGCTTCATCACCATTTTCAAATTTACTTAAACTATTTTTAAAATCATGTCGCAAATCCTCTAAatat                                                          <  1:819309/100‑1 (MQ=255)                  aTCTAACTCTACGAAACCTAATAACGTATCAAATGCTTCATCACCATTTTCAAATTTACTTAAACTATTTTTAAAATCATGTCGCAAATCCTCTAAatat                                                          <  1:3959419/100‑1 (MQ=255)                  aTCTAACTCTACGAAACCTAATAACGTATCAAATGCTTCATCACCATTTTCAAATTTACTTAAACTATTTTTAAAATCATGTCGCAAATCCTCTAAatat                                                          <  1:640702/100‑1 (MQ=255)                  aTCTAACTCTACGAAACCTAATAACGTATCAAATGCTTCATCACCATTTTCAAATTTACTTAAACTATTTTTAAAATCATGTCGCAAATCCTCTAAatat                                                          <  1:1850011/100‑1 (MQ=255)                  aTCTAACTCTACGAAACCTAATAACGTATCAAATGCTTCATCACCATTTTCAAATTTACTTAAACTATTGTTAAAATCATGTCGCAAATCCTCTAAatat                                                          <  1:5893140/100‑1 (MQ=255)                      aaCTCTACGAAACCTAATAACGTATCAAATGCTTCATCACCATTTTCAAATTTACTTAAACTATTTTTAAAATCATGTCGCAAATCCTCTAAATATAggc                                                      <  1:4614986/100‑3 (MQ=255)                      aaCTCTACGAAACCTAATAACGTATCAAATGCTTCATCACCATTTTCAAATTTACTTAAACTATTTTTAAAATCATGTCGCAAATCCTCTAAATATAGTg                                                      <  1:5586449/100‑1 (MQ=255)                      aaCTCTACGAAACCTAATAACGTATCAAATGCTTCATCACCATTTTCAAATTTACTTAAACTATTTTTAAAATCATGTCGCAAATCCTCTAAATATAGTg                                                      <  1:2892281/100‑1 (MQ=255)                      aaCTCTACGAAACCTAATAACGTATCAAATGCTTCATCACCATTTTCAAATTTACTTAAACTATTTTTAAAATCATGTCGCAAATCCTCTAAATATAGTg                                                      <  1:4636172/100‑1 (MQ=255)                          ctACGAAACCTAATAACGTATCAAATGCTTCATCACCATTTTCAAATTTACTTAAACTATTTTTAAAATCATGTCGCAAATCCTCTAAATATAGTGATgg                                                  <  1:2280977/100‑1 (MQ=255)                          ctACGAAACCTAATAACGTATCAAATGCTTCATCACCATTTTCAAATTTACTTAAACTATTTTTAAAATCATGTCGCAAATCCTCTAAATATAGTGATgg                                                  <  1:5841432/100‑1 (MQ=255)                          ctACGAAACCTAATAACGTATCAAATGCTTCATCACCATTTTCAAATTTACTTAAACTATTTTTAAAATCATGTCGCAAATCCTCTAAATATAGTGATgg                                                  <  1:3787302/100‑1 (MQ=255)                           tACGAAACCTAATAACGTATCAAATGCTTCATCACCATTTTCAAATTTACTTAAACTATTTTTAAAATCATGTCGCAAATCCTCTAAATATAGTGagaga                                                 <  1:4734322/100‑5 (MQ=255)                           tACGAAACCTAATAACGTATCAAATGCTTCATCACCATTTTCAAATTTACTTAAACTATTTTTAAAATCATGTCGCAAATCCTCTAAATATAGTGATGGt                                                 <  1:4541688/100‑1 (MQ=255)                            aCGAAACCTAATAACGTATCAAATGCTTCATCACCATTTTCAAATTTACTTAAACTATTTTTAAAATCATGTCGCAAAGCCTCTAAATATAGTGATGGta                                                <  1:454891/100‑2 (MQ=255)                            aCGAAACCTAATAACGTATCAAATGCTTCATCACCATTTTCAAATTTACTTAAACTATTTTTAAAATCATGTCGCAAAACCTCTAAATATAGTGATgggg                                                <  1:5998597/100‑3 (MQ=255)                            aCGAAACCTAATAACGTATCAAATGCTCCATCACCATTTTCAAATTTACTTAAACTATTTTTAAAATCATGTCGCAAATCCTCTAAATATAGTGATGGtt                                                <  1:4125855/100‑1 (MQ=255)                              gAAACCTAATAACGTATCAAATGCTTCATCACCATTTTCAAATTTACTTAAACTATTTTTAAAATCATGTCGCAAATCCTCTAAATATAGTGATGGtttt                                              <  1:397933/100‑1 (MQ=255)                              gAAACCTAATAACGTATCAAATGCTTCATCACCATTTTCAAATTTACTTAAACTATTTTTAAAATCATGTCGCAAATCCTCTAAATATAGTGATGGttat                                              <  1:1742314/100‑3 (MQ=255)                               aaaCCTAATAACGTATCAAATGCTTCATCACCATTTTCAAATTTACTTAAACTATTTTTAAAATCATGTCGCAAATCCTCTAAATATAGTGATGGTTTTc                                             <  1:5744925/100‑1 (MQ=255)                                        aaCGTATCAAATGCTTCATCACCATTTTCAAATTTACTTAAACTATTTTTAAAATCATGTCGCAAATCCTCTAAATATAGTGATGGTTTTCGATCtacat                                    <  1:2148622/100‑1 (MQ=255)                                          cGTATCAAATGCTTCATCACCATTTTCAAATTTACTTAAACTATTTTTAAAATCATGTCGCAAATCCTCTAAATATAGTGATGGTTTTCGATCtacatac                                  <  1:3907169/100‑1 (MQ=255)                                          cGTATCAAATGCTTCATCACCATTTTCAAATTTACTTAAACTATTTTTAAAATCATGTCGCAAATCCTCTAAATATAGTGATGGTTTTCGATCtacatac                                  <  1:1711171/100‑1 (MQ=255)                                           gTATCAAATGCTTCATCACCATTTTCAAATTTACTTAAACTATTTTTAAAATCATGTCGCAAATCCTCTAAATATAGTGATGGTTTTCGATCtacataca                                 <  1:2494596/100‑1 (MQ=255)                                            tATCAAATGCTTCATCACCATTTTCAAATTTACTTAAACTATTTTTAAAATCATGTCGCAAATCCTCTAAATATAGTGATGGTTTTCGATCTacatacat                                <  1:2042441/100‑1 (MQ=255)                                             aTCAAATGCTTCATCACCATTTTCAAATTTACTTAAACTATTTTTAAAATCATGTCGCAAATCCTCTAAATATAGTGATGGTTTTCGATCtacatacagc                               <  1:4381467/100‑3 (MQ=255)                                              tCAAATGCTTCATCACCATTTTCAAATTTACTTAAACTATTTTTAAAATCATGTCGCAAATCCTCTAAATATAGTGATGGTTTTCGATCTACATACAttt                              <  1:6332641/100‑1 (MQ=255)                                                     cTTCATCACCATTTTCAAATTTACTTAAACTATTTTTAAAATCATGTCGCAAATCCTCTAAATATAGTGATGGTTTTCGATCTACATACATTTTTATACt                       <  1:4375659/100‑1 (MQ=255)                                                     cTTCATCACCATTTTCAAATTTACTTAAACTATTTTTAAAATCATGTCGCAAATCCTCTAAATATAGTGATGGTTTTCGATCTACATACATTTTTATACt                       <  1:6346281/100‑1 (MQ=255)                                                     cTTCATCACCATTTTCAAATTTACTTAAACTATTTTTAAAATCATGTCGCAAATCCTCTAAATATAGTGATGGTTTTCGATCTACATACATTTTTATACt                       <  1:1766024/100‑1 (MQ=255)                                                     cTTCATCACCATTTTCAAATTTACTTAAACTATTTTTAAAATCATGTCGCAAATCCTCTAAATATAGTGATGGTTTTCGATCTACATACATTTTTATACt                       <  1:2021739/100‑1 (MQ=255)                                                     cTTCATCACCATTTTCAAATTTACTTAAACTATTTTTAAAATCATGTCGCAAATCCTCTAAATATAGTGATGGTTTTCGATCTACATACATTTTTATACt                       <  1:5839494/100‑1 (MQ=255)                                                      ttCATCACCATTTTCAAATTTACTTAAACTATTTTTAAAATCATGTCGCAAATCCTCTAAATATAGTGATGGTTTTCGATCTACATACATTTTTATACTa                      <  1:2107450/100‑1 (MQ=255)                                                       tcatcaCCATTTTCAAATTTACTTAAACTATTTTTAAAATCATGTCGCAAATCCTCTAAATATAGTGATGGTTTTCGATCTACATACATTTTTATACTac                     <  1:2617526/100‑2 (MQ=255)                                                       tcatcaCCATTTTCAAATTTACTTAAACTATTTTTAAAATCATGTCGCAAATCCTCTAAATATAGTGATGGTTTTCGATCTACATACATTTTTATACTac                     <  1:767210/100‑2 (MQ=255)                                                        catcaCCATTTTCAAATTTACTTAAACTATTTTTAAAATCATGTCGCAAATCCTCTAAATATAGTGATGGTTTTCGATCTACATACATTTTTATACTagg                    <  1:1310903/100‑3 (MQ=255)                                                          tcaCCATTTTCAAATTTACTTAAACTATTTTTAAAATCATGTCGCAAATCCTCTAAATATAGTGATGGTTTTCGATCTACATACATTTTTATACTAAcca                  <  1:1741253/100‑2 (MQ=255)                                                             ccATTTTCAAATTTACTTAAACTATTTTTAAAATCATGTCGCAAATCCTCTAAATCTAGTGATGGTTTTCGATCTACATACATTTTTATACTAACCTCCg               <  1:6054488/100‑1 (MQ=255)                                                              cATTTTCAAATTTACTTAAACTATTTTTAAAATCATGTCGCAAATCCTCTAAATATAGTGATGGTTTTCGATCTACATACATTTTTATACTAACCTCCgt              <  1:3242082/100‑2 (MQ=255)                                                                   tCAAATTTACTTAAACTATTTTTAAAATCATGTCGCAAATCCTCTAAATATAGTGATGGTTTTCGATCTACATACATTTTTATACTAACCTCCGatatat         <  1:3058756/100‑1 (MQ=255)                                                                   tCAAATTTACTTAAACTATTTTTAAAATCATGTCGCAAATCCTCTAAATATAGTGATGGTTTTCGATCTACATACATTTTTATACTAACCTCCGatatat         <  1:1480702/100‑1 (MQ=255)                                                                    cAAATTTACTTAAACTATTTTTAAAATCATGTCGCAAATCCTCTAAATATAGTGATGGTTTTCGATCTACATACATTTTTATACTAACCTCCGatatacc        <  1:3554630/100‑3 (MQ=255)                                                                    cAAATTTACTTAAACTATTTTTAAAATCATGTCGCAAATCCTCTAAATATAGTGATGGTTTTCGATCTACATACATTTTTATACTAACCTCCGAtatata        <  1:2665872/100‑1 (MQ=255)                                                                          tACTTAAACTATTTTTAAAATCATGTCGCAAATCCTCTAAATATAGTGATGGTTTTCGATCTACATACATTTTTATACTAACCTCCGATATATAatcagg  <  1:3428604/100‑3 (MQ=255)                                                                                 |                                                                                             GCTGACGAATAAATATGATCTAACTCTACGAAACCTAATAACGTATCAAATGCTTCATCACCATTTTCAAATTTACTTAAGCTATTTTTAAAATCATGTCGCAAATCCTCTAAATATAGTGATGGTTTTCGATCTACATACATTTTTATACTAACCTCCGATATATAATCATC  >  NC_007793/2641635‑2641807

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