Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_007793	1,781,115	T→G	E338A (GAA→GCA)	SAUSA300_RS08880 ←	gamma‑aminobutyrate permease

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_007793	1,781,115	0	T	G	100.0%	111.1 / NA	39	E338A (GAA→GCA)	SAUSA300_RS08880	gamma‑aminobutyrate permease
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base G (3/36);  total (3/36)

GTAAGGCTTTGTAATACGAAAATTATCACTACAATGATTGCTGTAACGAGTAACGACATATATGGCACACCGTTTTTATTTGTTTTACCAAATGCTTCAAACGCTAATTTATCTTTACTCATTGAATATAGCATTCGAGTTGAAGCATACATACCTGAGTTACCTGCTGATAA  >  NC_007793/1781019‑1781191                                                                                                 |                                                                             gTAAGGCTTTGTAATACGAAAATTATCACTACAATGATTGCTGTAACGAGTAACGACATATATGGCACACCGTTTTTATTTGTTTTACCAAATGCTGCaa                                                                           <  1:1603373/100‑1 (MQ=255) gTAAGGCTTTGTAATACGAAAATTATCACTACAATGATTGCTGTAACGAGTAACGACATATATGGCACACCGTTTTTATTTGTTTTACCAAATGCTGCaa                                                                           <  1:5649497/100‑1 (MQ=255) gTAAGGCTTTGTAATACGAAAATTATCACTACAATGATTGCTGTAACGAGTAACGACATATATGGCACACCGTTTTTATTTGTTTTACCAAATGCTGCaa                                                                           <  1:1695072/100‑1 (MQ=255) gTAAGGCTTTGTAATACGAAAATTATCACTACAATGATTGCTGTAACGAGTAACGACATATATGGCACACCGTTTTTATTTGTTTTACCAAATGCTGCaa                                                                           <  1:2727520/100‑1 (MQ=255)                         aTCACTACAATGATTGCTGTAACGAGTAACGACATATATGGCACACCGTTTTTATTTGTTTTACCAAATGCTGCAAACGCTAATTTATCTTTACTCATTg                                                   <  1:3919323/100‑1 (MQ=255)                              tACAATGATTGCTGTAACGAGTAACGACATATATGGCACACCGTTTTTATTTGTTTTACCAAATGCTGCAAACGCTAATTTATCTTTACTCATTGAatat                                              <  1:472302/100‑1 (MQ=255)                              tACAATGATTGCTGTAACGAGTAACGACATATATGGCACACCGTTTTTATTTGTTTTACCAAATGCTGCAAACGCTAATTTATCTTTACTCATTGAatat                                              <  1:3071131/100‑1 (MQ=255)                                 aaTGATTGCTGTAACGAGTAACGACATATATGGCACACCGTTTTTATTTGTTTCACCAAATGCTGCAAACGCTAATTTATCTTTACTCATTGAATATAGc                                           <  1:1305671/100‑1 (MQ=255)                                  aTGATTGCTGTAACGAGTAACGACATATATGGCACACCGTTTTTATTTGTTTTACCAAATGCTGCAAACGCTAATTTATCTTTACTCATTGAATATAGCa                                          <  1:5827956/100‑1 (MQ=255)                                   tGATTGCTGTAACGAGTAACGACATATATGGCACACCGTTTTTATTTGTTTTACCAAATGCTGCAAACGCTAATTTATCTTTACTCATTGAATATAGCAt                                         <  1:5696227/100‑1 (MQ=255)                                      ttGCTGTAACGAGTAACGACATATATGGCACACCGTTTTTATTTGTTTTACCAAATGCTGCAAACGCTAATTTATCTTTACTCATTGAATATAGCATTct                                      <  1:4519071/100‑2 (MQ=255)                                      ttGCTGTAACGAGTAACGACATATATGGCACACCGTTTTTATTTGTTTTACCAAATGCTGCAAACGCTAATTTATCTTTACTCATTGAATATAGCATTCg                                      <  1:730540/100‑1 (MQ=255)                                      ttGCTGTAACGAGTAACGACATATATGGCACACCGTTTTTATTTGTTTTACCAAATGCTGCAAACGCTAATTTATCTTTACTCATTGAATATAGCATTCg                                      <  1:4535256/100‑1 (MQ=255)                                           gtaacgagtaacgaCATATATGGCACACCGTTTTTATTTGTTTTACCAAATGCTGCAAACGCTAATTTATCTTTACTCATTGAATATAGCATTCGAGTTg                                 >  1:216791/1‑100 (MQ=255)                                                          atatatGGCACACCGTTTTTATTTGTTTTACCAAATGCTGCAAACGCTAATTTATCTTTACTCATTGAATATAGCATTCGAGTTGAAGCATACATACCTg                  <  1:1419288/100‑1 (MQ=255)                                                          atatatGGCACACCGTTTTTATTTGTTTTACCAAATGCTGCAAACGCTAATTTATCTTTACTCATTGAATATAGCATTCGAGTTGAAGCATACATACCTg                  <  1:2696023/100‑1 (MQ=255)                                                          atatatGGCACACCGTTTTTATTTGTTTTACCAAATGCTGCAAACGCTAATTTATCTTTACTCATTGAATATAGCATTCGAGTTGAAGCATACATACCTg                  <  1:5995877/100‑1 (MQ=255)                                                          atatatGGCACACCGTTTTTATTTGTTTTACCAAATGCTGCAAACGCTAATTTATCTTTACTCATTGAATATAGCATTCGAGTTGAAGCATACATACCTg                  <  1:1637061/100‑1 (MQ=255)                                                          atatatGGCACACCGTTTTTATTTGTTTTACCAAATGCTGCAAACGCTAATTTATCTTTACTCATTGAATATAGCATTCGAGTTGAAGCATACATACCTg                  <  1:5572546/100‑1 (MQ=255)                                                          atatatGGCACACCGTTTTTATTTGTTTTACCAAATGCTGCAAACGCTAATTTATCTTTACTCATTGAATATAGCATTCGAGTTGAAGCATACATACCTg                  <  1:5511884/100‑1 (MQ=255)                                                          atatatGGCACACCGTTTTTATTTGTTTTACCAAATGCTGCAAACGCTAATTTATCTTTACTCATTGAATATAGCATTCGAGTTGAAGCATACATACCTg                  <  1:2228023/100‑1 (MQ=255)                                                          atatatGGCACACCGTTTTTATTTGTTTTACCAAATGCTGCAAACGCTAATTTATCTTTACTCATTGAATATAGCATTCGAGTTGAAGCATACATACCTg                  <  1:2418746/100‑1 (MQ=255)                                                           tatatGGCACACCGTTTTTATTTGTTTTACCAAATGCTGCAAACGCTAATTTATCTTTACTCATTGAATATAGCATTCGAGTTGAAGCATACATACCTgc                 <  1:1847928/100‑2 (MQ=255)                                                           tatatGGCACACCGTTTTTATTTGTTTTACCAAATGCTGCAAACGCTAATTTATCTTTACTCATTGAATATAGCATTCGAGTTGAAGCATACATACCTgc                 <  1:2866618/100‑2 (MQ=255)                                                                   acacCGTTTTTATTTGTTTTACCAAATGCTGCAAACGCTAATTTATCTTTACTCATTGAATATAGCATTCGAGTTGAAGCATACATACCTGAGTTACctg         <  1:3939670/100‑1 (MQ=255)                                                                   acacCGTTTTTATTTGTTTTACCAAATGCTGCAAACGCTAATTTATCTTTACTCATTGAATATAGCATTCGAGTTGAAGCATACATACCTGAGTTACctg         >  1:2712626/1‑100 (MQ=255)                                                                   acacCGTTTTTATTTGTTTTACCAAATGCTGCAAACGCTAATTTATCTTTACTCATTGAATATAGCATTCGAGTTGAAGCATACATACCTGAGTTACctg         >  1:5180704/1‑100 (MQ=255)                                                                    cacCGTTTTTATTTGTTTTACCAAATGCTGCAAACGCTAATTTATCTTTACTCATTGAATATAGCATTCGAGTTGAAGCATACATACCTGAGTTACctgt        <  1:4388947/100‑2 (MQ=255)                                                                    cacCGTTTTTATTTGTTTTACCAAATGCTGCAAACGCTAATTTATCTTTACTCATTGAATATAGCATTCGAGTTGAAGCATACATACCTGAGTTACctgc        <  1:3892120/100‑1 (MQ=255)                                                                    cacCGTTTTTATTTGTTTTACCAAATGCTGCAAACGCTAATTTATCTTTACTCATTGAATATAGCATTCGAGTTGAAGCATACATACCTGAGTTACctgc        <  1:4068610/100‑1 (MQ=255)                                                                    cacCGTTTTTATTTGTTTTACCAAATGCTGCAAACGCTAATTTATCTTTACTCATTGAATATAGCATTCGAGTTGAAGCATACATACCTGAGTTACctgc        <  1:142464/100‑1 (MQ=255)                                                                     gccGTTTTTATTTGTTTTACCAAATGCTGCAAACGCTAATTTATCTTTACTCATTGAATTTAGCATTCGAGTTGAAGCGTACATACCTGAGTTACctgcg       <  1:2263053/99‑2 (MQ=255)                                                                     acCGTTTTTATTTGTTTTACCAAATGCTGCAAACGCTAATTTATCTTTACTCATTGAATATAGCATTCGAGTTGAAGCATACATACCTGAGTTACctgct       <  1:3867437/100‑1 (MQ=255)                                                                     acCGTTTTTATTTGTTTTACCAAATGCTGCAAACGCTAATTTATCTTTACTCATTGAACATAGCATTCGAGTTGAAGCATACATACCTGAGTTACctgct       <  1:2821480/100‑1 (MQ=255)                                                                      ccGTTTTTATTTGTTTTACCAAATGCTGCAAACGCTAATTTATCTTTACTCATTGAATATAGCATTCGAGTTGAAGCATACATACCTGAGTTACctgctg      <  1:2917802/100‑1 (MQ=255)                                                                        gTTTTTATTTGTTTTACCAAATGCTGCAAACGCTAATTTATCTTTACTCATTGAATATAGCATTCGAGTTGAAGCATACATACCTGAGTTACCTGCTGAt    <  1:1651294/100‑1 (MQ=255)                                                                        gTTTTTATTTGTTTTACCAAATGCTGCAAACGCTAATTTATCTTTACTCATTGAATATAGCATTCGAGTTGAAGCATACATACCTGAGTTACCTGCTGAt    <  1:2006178/100‑1 (MQ=255)                                                                        gTTTTTATTTGTTTTACCAAATGCTGCAAACGCTAATTTATCTTTACTCATTGAATATAGCATTCGAGTTGAAGCATACATACCTGAGTTACCTGCAGaa    <  1:5967418/100‑2 (MQ=255)                                                                          ttttATTTGTTTTACCAAATGCTGCAAACGCTAATTTATCTTTACTCATTGAATATAGCATTCGAGTTGAAGCATACATACCTGAGTTACCCGCTGATaa  <  1:712929/100‑1 (MQ=255)                                                                                                 |                                                                             GTAAGGCTTTGTAATACGAAAATTATCACTACAATGATTGCTGTAACGAGTAACGACATATATGGCACACCGTTTTTATTTGTTTTACCAAATGCTTCAAACGCTAATTTATCTTTACTCATTGAATATAGCATTCGAGTTGAAGCATACATACCTGAGTTACCTGCTGATAA  >  NC_007793/1781019‑1781191

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