Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,007,952	A→G	K66E (AAA→GAA)	rhtC →	L‑threonine exporter

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,007,952	0	A	G	100.0%	126.8 / NA	35	K66E (AAA→GAA)	rhtC	L‑threonine exporter
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (18/17);  total (18/17)

GGCGTGCTGGGCATTACCTGCGGCGTAATGGTTTGGGCTGGGATTGCGCTGCTTGGCCTGCATTTGATTATCGAAAAAATGGCCTGGCTGCATACGCTGATTATGGTGGGCGGTGGCCTGTATCTCTGCTGGATGGGTTACCAGATGCTAC  >  NC_000913/4007877‑4008027                                                                            |                                                                            ggCGTGCTGGGCATTACCTGCGGCGTAATGGTTTGGGCTGGGATTGCGCTGCTTGGCCTGCATTTGATTATCGAAGAAAt                                                                         <  1:737345/80‑1 (MQ=255) ggCGTGCTGGGCATTACCTGCGGCGTAATGGTTTGGGCTGGGATTGCGCTGCTTGGCCTGCATTTGATTATCGAAGAAATg                                                                        <  2:737346/81‑1 (MQ=255)  gCGTGCTGGGCATTACCTGCGGCGTAATGGTTTGGGCTGGGATTGCGCTGCTTGGCCTGCATTTGATTATCGAAGAAATgg                                                                       <  1:737348/81‑1 (MQ=255)  gCGTGCTGGGCATTACCTGCGGCGTAATGGTTTGGGCTGGGATTGCGCTGCTTGGCCTGCATTTGATTATCGAAGAAATgg                                                                       <  2:737347/81‑1 (MQ=255)    gTGCTGGGCATTACCTGCGGCGTAATGGTTTGGGCTGGGATTGCGCTGCTTGGCCTGCATTTGATTATCGAAGAAATGGcc                                                                     >  2:737404/1‑81 (MQ=255)          ggCATTACCTGCGGCGTAATGGTTTGGGCTGGGATTGCGCTGCTTGGCCTGCATTTGATTATCGAAGAAATGGCCTGGCt                                                                <  1:737349/80‑1 (MQ=255)             aTTACCTGCGGCGTAATGGTTTGGGCTGGGATTGCGCTGCTTGGCCTGCATTTGATTATCGAAGAAATGGCCTGGCTGCa                                                             >  2:737402/1‑80 (MQ=255)                  cTGCGGCGTAATGGTTTGGGCTGGGATTGCGCTGCTTGGCCTGCATTTGATTATCGAAGAAATGGCCTGGCTGCATACGCt                                                       <  1:737350/81‑1 (MQ=255)                   tGCGGCGTAATGGTTTGGGCTGGGATTGCGCTGCTTGGCCTGCATTTGATTATCGAAGAAATGGCCTGGCTGCATACGCt                                                       >  1:737388/1‑80 (MQ=255)                     cggcgTAATGGTTTGGGCTGGGATTGCGCTGCTTGGCCTGCATTTGATTATCGAAGAAATGGCCTGGCTGCATACGCTGAt                                                    >  1:737415/1‑81 (MQ=255)                      ggcgTAATGGTTTGGGCTGGGATTGCGCTGCTTGGCCTGCATTTGATTATCGAAGAAATGGCCTGGCTGCATACGCTGAtt                                                   <  1:737351/81‑1 (MQ=255)                        cgTAATGGTTTGGGCTGGGATTGCGCTGCTTGGCCTGCATTTGATTATCGAAGAAATGGCCTGGCTGCATACGCCGATTAt                                                 >  1:737390/1‑81 (MQ=255)                         gTAATGGTTTGGGCTGGGATTGCGCTGCTTGGCCTGCATTTGATTATCGAAGAAATGGCCTGGCTGCATACGCTGATTAtg                                                >  2:737424/1‑81 (MQ=255)                          tAATGGTTTGGGCTGGGATTGCGCTGCTTGGCCTGCATTTGATTATCGAAGAAATGGCCTGGCTGCATACGCTGATTAtgg                                               <  1:737352/81‑1 (MQ=255)                               gTTTGGGCTGGGATTGCGCTGCTTGGCCTGCATTTGATTATCGAAGAAATGGCCTGGCTGCATACGCTGATTATGGTggg                                           <  1:737355/80‑1 (MQ=255)                               gTTTGGGCTGGGATTGCGCTGCTTGGCCTGCATTTGATTATCGAAGAAATGGCCTGGCTGCATACGCTGATTATGGTggg                                           <  1:737354/80‑1 (MQ=255)                               gTTTGGGCTGGGATTGCGCTGCTTGGCCTGCATTTGATTATCGAAGAAATGGCCTGGCTGCATACGCTGATTATGGTGGGc                                          >  1:737398/1‑81 (MQ=255)                               gTTTGGGCTGGGATTGCGCTGCTTGGCCTGCATTTGATTATCGAAGAAATGGCCTGGCTGCATACGCTGATTATGGTGGGc                                          <  1:737353/81‑1 (MQ=255)                                   gggCTGGGATTGCGCTGCTTGGCCTGCATTTGATTATCGAAGAAATGGCCTGGCTGCATACGCTGATTATGGTGGGCGGTg                                      >  1:737419/1‑81 (MQ=255)                                   gggCTGGGATTGCGCTGCTTGGCCTGCATTTGATTATCGAAGAAATGGCCTGGCTGCATACGCTGATTATGGTGGGCGGTg                                      >  2:737410/1‑81 (MQ=255)                                    ggCTGGGATTGCGCTGCTTGGCCTGCATTTGATTATCGAAGAAATGGCCTGGCTGCATACGCTGATTATGGTGGGCGGTgg                                     >  2:737396/1‑81 (MQ=255)                                     gCTGGGATTGCGCTGCTTGGCCTGCATTTGATTATCGAAGAAATGGCCTGGCTGCATACGCTGATTATGGTGGGCGGTGGc                                    >  1:737432/1‑81 (MQ=255)                                     gCTGGGATTGCGCTGCTTGGCCTGCATTTGATTATCGAAGAAATGGCCTGGCTGCATACGCTGATTATGGTGGGCGGTGGc                                    >  2:737407/1‑81 (MQ=255)                                           aTTGCGCTGCTTGGCCTGCATTTGATTATCGAAGAAATGGCCTGGCTGCATACGCTGATTATGGTGGGCGGTGGCCTGTAt                              >  1:737397/1‑81 (MQ=255)                                             tGCGCTGCTTGGCCTGCATTTGATTATCGAAGAAATGGCCTGGCTGCATACGCTGATTATGGTGGGCGGTGGCCTGTAtct                            >  2:737408/1‑81 (MQ=255)                                                gctgctTGGCCTGCATTTGATTATCGAAGAAATGGCCTGGCTGCATACGCTGATTATGGTGGGCGGTGGCCTGTATCTctg                         <  2:737356/81‑1 (MQ=255)                                                    ctTGGCCTGCATTTGATTATCGAAGAAATGGCCTGGCTGCATACGCTGATTATGGTGGGCGGTGGCCTGTATCTctgctg                      <  2:737357/80‑1 (MQ=255)                                                    ctTGGCCTGCATTTGATTATCGAAGAAATGGCCTGGCTGCATACGCTGATTATGGTGGGCGGTGGCCTGTATCTctgct                       <  1:737358/79‑1 (MQ=255)                                                     ttGGCCTGCATTTGATTATCGAAGAAATGGCCTGGCTGCATACGCTGATTATGGTGGGCGGTGGCCTGTATCTCTGCTGGa                    >  1:737423/1‑81 (MQ=255)                                                      tGGCCTGCATTTGATTATCGAAGAAATGGCCTGGCTGCATACGCTGATTATGGTGGGCGGTGGCCTGTATCTCTGCTGGAt                   >  1:737412/1‑81 (MQ=255)                                                               tttGATTATCGAAGAAATGGCCTGGCTGCATACGCTGATTATGGTGGGCGGTGGCCTGTATCTCTGCTGGATGGGTTACCa          <  1:737359/81‑1 (MQ=255)                                                                 tGATTATCGAAGAAATGGCCTGGCTGCATACGCTGATTATGGTGGGCGGTGGCCTGTATCTCTGCTGGATGGGTTAc            >  1:737411/1‑77 (MQ=255)                                                                   aTTATCGAAGAAATGGCCTGGCTGCATACGCTGATTATGGTGGGCGGTGGCCTGTATCTCTGCTGGATGGGTTACCAGATg      >  2:737422/1‑81 (MQ=255)                                                                    ttATCGAAGAAATGGCCTGGCTGCATACGCTGATTATGGTGGGCGGTGGCCTGTATCTCTGCTGGATGGGTTACCAGATGc     <  1:737360/81‑1 (MQ=255)                                                                       tCGAAGAAATGGCCTGGCTGCATACGCTGATTATGGTGGGCGGTGGCCTGTATCTCTGCTGGATGGGTTACCAGATGCTAc  <  1:737361/81‑1 (MQ=255)                                                                            |                                                                            GGCGTGCTGGGCATTACCTGCGGCGTAATGGTTTGGGCTGGGATTGCGCTGCTTGGCCTGCATTTGATTATCGAAAAAATGGCCTGGCTGCATACGCTGATTATGGTGGGCGGTGGCCTGTATCTCTGCTGGATGGGTTACCAGATGCTAC  >  NC_000913/4007877‑4008027

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