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%	113.7 / NA	32	K66E (AAA→GAA)	rhtC	L‑threonine exporter
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (14/18);  total (14/18)

ATGGGCGTGCTGGGCATTACCTGCGGCGTAATGGTTTGGGCTGGGATTGCGCTGCTTGGCCTGCATTTGATTATCGAAAAAATGGCCTGGCTGCATACGCTGATTATGGTGGGCGGTGGCCTGTATCTCTGCTGGATGGGTTACCAGATGCTA  >  NC_000913/4007874‑4008026                                                                               |                                                                           atgGGCGTGCTGGGCATTACCTGCGGCGTAATGGTTTGGGCTGGGATTGCGCTGCTTGGCCTGCATTTGATTATCgaagaa                                                                          <  2:788015/81‑1 (MQ=255)    ggCGTGCTGGGCATTACCTGCGGCGTAATGGTTTGGGCTGGGATTGCGCTGCTTGGCCTGCATTTGATTATCGAAGAAATg                                                                       <  2:788016/81‑1 (MQ=255)    ggCGTGCTGGGCATTACCTGCGGCGTAATGGATTGGGCTGGGATTGCGCTGCTTGGCCTGCATTTGATTATCGAAGAAATg                                                                       <  2:788017/81‑1 (MQ=255)     gCGTGCTGGGCATTACCTGCGGCGTAATGGTTTGGGCTGGGATTGCGCTGCTTGGCCTGCATTTGATTATCGAAGAAATgg                                                                      <  2:788018/81‑1 (MQ=255)         gCTGGGCATTACCTGCGGCGTAATGGTTTGGGCTGGGATTGCGCTGCTTGGCCTGCATTTGATTATCGAAGAAATGGCCTg                                                                  >  2:788075/1‑81 (MQ=255)         gCTGGGCATTACCTGCGGCGTAATGGTTTGGGCTGGGATTGCGCTGCTTGGCCTGCATTTGATTATCGAAGAAATGGCCTg                                                                  >  2:788062/1‑81 (MQ=255)             ggCATTACCTGCGGCGTAATGGTTTGGGCTGGGATTGCGCTGCTTGGCCTGCATTTGATTATCGAAGAAATGGCCTGGCt                                                               <  1:788019/80‑1 (MQ=255)               cATTACCTGCGGCGTAATGGTTTGGGCTGGGATTGCGCTGCTTGGCCTGCATTTGATTATCGAAGAAATGGCCTGGCTGCa                                                            >  2:788083/1‑81 (MQ=255)                aTTACCTGCGGCGTAATGGTTTGGGCTGGGATTGCGCTGCTTGGCCTGCATTTGATTATCGAAGAAATGGCCTGGCTGCAt                                                           <  2:788020/81‑1 (MQ=255)                   aCCTGCGGCGTAATGGTTTGGGCTGGGATTGCGCTGCTTGGCCTGCATTTGATTATCGAAGAAATGGCCTGGCTGCATAc                                                         <  2:788021/80‑1 (MQ=255)                   aCCTGCGGCGTAATGGTTTGGGCTGGGATTGCGCTGCTTGGCCTGCATTTGATTATCGAAGAAATGGCCTGGCTGCATACg                                                        >  2:788065/1‑81 (MQ=255)                     cTGCGGCGTAATGGTTTGGGCTGGGATTGCGCTGCTTGGCCTGCATTTGATTATCGAAGAAATGGCCTGGCTGCATACGCt                                                      <  2:788022/81‑1 (MQ=255)                        cggcgTAATGGTTTGGGCTGGGATTGCGCTGCTTGGCCTGCATTTGATTATCGAAGAAATGGCCTGGCTGCATACGCTGAt                                                   <  1:788023/81‑1 (MQ=255)                        cggcaTAATGGTTTGGGCTGGGATTGCGCTGCTTGGCCTGCATTTGATTATCGAAGAAATGGCCTGGCTGCATACGCTg                                                     <  1:788024/79‑1 (MQ=255)                         ggcgTAATGGTTTGGGCTGGGATTGCGCTGCTTGGCCTGCATTTGATTATCGAAGAAATGGCCTGGCTGCATACGCTGAtt                                                  <  2:788025/81‑1 (MQ=255)                           cgTAATGGTTTGGGCTGGGATTGCGCTGCTTGGCCTGCATTTGATTATCGAAGAAATGGCCTGGCTGCATACGCTGATTAt                                                <  2:788026/81‑1 (MQ=255)                              aaTGGTTTGGGCTGGGATTGCGCTGCTTGGCCTGCATTTGATTATCGAAGAAATGGCCTGGCTGCATACGCTGATTAtggt                                             >  1:788098/1‑81 (MQ=255)                               aTGGTTTGGGCTGGGATTGCGCTGCTTGGCCTGCATTTGATTATCGAAGAAATGGCCTGGCTGCATACGCTGATTAtggtg                                            >  1:788066/1‑81 (MQ=255)                                  gTTTGGGCTGGGATTGCGCTGCTTGGCCTGCATTTGATTATCGAAGAAATGGCCTGGCTGCATACGCTGATTATGGTGGGc                                         >  2:788067/1‑81 (MQ=255)                                          tGGGATTGCGCTGCTTGGCCTGCATTTGATTATCGAAGAAATGGCCTGGCTGCATACGCTGATTATGGTGGGCGGTGGCCt                                 >  1:788073/1‑81 (MQ=255)                                               ttGCGCTGCTTGGCCTGCATTTGATTATCGAAGAAATGGCCTGGCTGCATACGCTGATTATGGTGGGCGGTGGCCTGTAtc                            >  1:788074/1‑81 (MQ=255)                                                  cgcTGCTTGGCCTGCATTTGATTATCGAAGAAATGGCCTGGCTGCATACGCTGATTATGGTGGGCGGTGGCCTGTAtctc                          <  2:788027/80‑1 (MQ=255)                                                  cgcTGCTTGGCCTGCATTTGATTATCGAAGAAATGGCCTGGCTGCATACGCTGATTATGGTGGGCGGTGGCCTGTATctct                         <  2:788028/81‑1 (MQ=255)                                                       ctTGGCCTGCATTTGATTATCGAAGAAATGGCCTGGCTGCATACGCTGATTATGGTGGGCGGTGGCCTGTATCTctgctg                     >  1:788063/1‑80 (MQ=255)                                                       ctTGGCCTGCATTTGATTATCGAAGAAATGGCCTGGCTGCATACGCTGATTATGGTGGGCGGTGGCCTGTATCTCTGCTgg                    <  2:788029/81‑1 (MQ=255)                                                       ctTGGCCTGCATTTGATTATCGAAGAAATGACCTGGCTGCATACGCTGATTATGGTGGGCGGTGGCCTGTATCTctgctg                     <  2:788030/80‑1 (MQ=255)                                                          ggCCTGCATTTGATTATCGAAGAAATGGCCTGGCTGCATACGCTGATTATGGTGGGCGGTGGCCTGTATCTCTGCTGGATg                 >  1:788079/1‑81 (MQ=255)                                                                cATTTGATTATCGAAGAAATGGCCTGGCTGCATACGCTGATTATGGTGGGCGGTGGCCTGTATCTCTGCTGGATGGGtt             <  2:788031/79‑1 (MQ=255)                                                                   ttGATTATCGAAGAAATGGCCTGGCTGCATACGCTGATTATGGTGGGCGGTGGCCTGTATCTCTGCTGGATGGGTTAc           >  1:788077/1‑78 (MQ=255)                                                                    tGATTATCGAAGAAATGGCCTGGCTGCATACGCTGATTATGGTGGGCGGTGGCCTGTATCTCTGCTGGATGGGTTAc           >  2:788087/1‑77 (MQ=255)                                                                      aTTATCGAAGAAATGGCCTGGCTGCATACGCTGATTATGGTGGGCGGTGGCCTGTATCTCTGCTGGATGGGTTACCAGATg     <  1:788032/81‑1 (MQ=255)                                                                      aTTATCGAAGAAATGGCCTGGCTGCATACGCTGATTATGGTGGGCGGTGGCCTGTATCTCTGCTGGATGGGTTACCAGATg     >  2:788082/1‑81 (MQ=255)                                                                          tCGAAGAAATGGCCTGGCTGCATACGCTGATTATGGTGGGCGGTGGCCTGTATCTCTGCTGGATGGGTTACCAGATGCTa  <  1:788033/80‑1 (MQ=255)                                                                               |                                                                           ATGGGCGTGCTGGGCATTACCTGCGGCGTAATGGTTTGGGCTGGGATTGCGCTGCTTGGCCTGCATTTGATTATCGAAAAAATGGCCTGGCTGCATACGCTGATTATGGTGGGCGGTGGCCTGTATCTCTGCTGGATGGGTTACCAGATGCTA  >  NC_000913/4007874‑4008026

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