Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,139,431	A→G	L93P (CTG→CCG)	fsaB ←	fructose‑6‑phosphate aldolase 2

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,139,431	0	A	G	100.0%	125.2 / NA	34	L93P (CTG→CCG)	fsaB	fructose‑6‑phosphate aldolase 2
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (26/8);  total (26/8)

TAACCCTTGTGCGGCGCTATATACAGCAGTGCCAAGTGTAGTAATACCCTCTTTTTTCAGTATTTTAATTGCTGCCAGACCTTCGGAAGTCACCGGGA‑TTTTCACCACAATACCCGGAATAGCGTCGCGCAGGCGCTTCGCTTCTTCCACCAT  >  NC_000913/4139355‑4139507                                                                             |                                                                              taaCCCTTGTGCGGCGCTATATACAGCAGTGCCAAGTGTAGTAATACCCTCTTTTTTCAGTATTTTAATTGCTGCCGGAcc                                                                           >  2:819295/1‑81 (MQ=255)      cTTGTGCGGCGCTATATACAGCAGTGCCAAGTGTAGTAATACCCTCTTTTTTCAGTATTTTAATTGCTGCCGGACCTTCgg                                                                      <  1:819264/81‑1 (MQ=255)         gtgCGGCGCTATATACAGCAGTGCCAAGTGTAGTAATACCCTCTTTTTTCAGTATTTTAATTGCTGCCGGACCTTCGGaa                                                                    <  2:819265/80‑1 (MQ=255)          tgCGGCGCTATATACAGCAGTGCCAAGTGTAGTAATACCCTCTTTTTTCAGTATTTTAATTGCTGCCGGACCTTCGGAAGt                                                                  >  1:819318/1‑81 (MQ=255)           gcggcgCTATATACAGCAGTGCCAAGTGTAGTAATACCCTCTTTTTTCAGTATTTTAATTGCTGCCGGACCTTCGGAAGt                                                                  >  2:819302/1‑80 (MQ=255)            cggcgCTATATACAGCAGTGCCAAGTGTAGTAATACCCTCTTTTTTCAGTATTTTAATTGCTGCCGGACCTTCGGAAGTCa                                                                >  2:819293/1‑81 (MQ=255)              gcgcTATATACAGCAGTGCCAAGTGTAGTAATACCCTCTTTTTTCAGTATTTTAATTGCTGCCGGACCTTCGGAAGTCAcc                                                              <  1:819266/81‑1 (MQ=255)              gcgcTATATACAGCAGTGCCAAGTGTAGTAATACCCTCTTTTTTCAGTATTTTAATTGCTGCCGGACCTTCGGAAGTCAcc                                                              <  2:819267/81‑1 (MQ=255)                  tatataCAGCAGTGCCAAGTGTAGTAATACCCTCTTTTTTCAGTATTTTAATTGCTGCCGGACCTTCGGAAGTCACCGGGa                                                          >  1:819317/1‑81 (MQ=255)                     ataCAGCAGTGCCAAGTGTAGTAATACCCTCTTTTTTCAGTATTTTAATTGCTGCCGGACCTTCGGAAGTCACCGGGA‑ttt                                                      >  2:819285/1‑81 (MQ=255)                       aCAGCAGTGCCAAGTGTAGTAATACCCTCTTTTTTCAGTATTTTAATTGCTGCCGGACCTTCGGAAGTCACCGGGA‑tttt                                                     >  1:819301/1‑80 (MQ=255)                       aCAGCAGTGCCAAGTGTAGTAATACCCTCTTTTTTCAGTATTTTAATTGCTGCCGGACCTTCGGAAGTCACCGGGA‑tttt                                                     >  1:819314/1‑80 (MQ=255)                          gcagTGCCAAGTGTAGTAATACCCTCTTTTTTCAGTATTTTAATTGCTGCCGGACCTTCGGAAGTCACCGGGA‑TTTTcac                                                  >  2:819290/1‑80 (MQ=255)                              tGCCAAGTGTAGTAATACCCTCTTTTTTCAGTATTTTAATTGCTGCCGGACCTTCGGAAGTCACCGGGA‑TTTTCACCACaa                                             >  2:819326/1‑81 (MQ=255)                                 cAAGTGTAGTAATACCCTCATTATTCAGTATTTTAATTGCTGCCGAACCTTCGAAAGTCACCGGGA‑TTTTCACCACAATAc                                          <  1:819268/81‑1 (MQ=255)                                  aaGTGTAGTAATACCCTCTTTTTTCAGTATTTTAATTGCTGCCGGACCTTCGGAAGTCACCGGGA‑TTTTCACCACAATAcc                                         >  2:819306/1‑81 (MQ=255)                                        agtaATACCCTCTTTTTTCAGTATTTTAATTGCTGCCGGACCTTCGGAAGTCACCGGGA‑TTTTCACCACAATACCCGGaa                                    >  1:819320/1‑80 (MQ=255)                                         gtaATACCCTCTTTTTTCAGTATTTTAATTGCTGCCGGACCTTCGGAAGTCACCGGGA‑TTTTCACCACAATACCCGGAAt                                   <  2:819269/80‑1 (MQ=255)                                          taATACCCTCTTTTTTCAGTATTTTAATTGCTGCCGGACCTTCGGAAGTCACCGGGA‑TTTTCACCACAATACCCGGAATAg                                 >  2:819307/1‑81 (MQ=255)                                           aaTACCCTCTTTTTTCAGTATTTTAATTGCTGCCGGACCTTCGGAAGTCACCGGGA‑TTTTCACCACAATACCCGGAATAGc                                >  1:819305/1‑81 (MQ=255)                                            aTACCCTCTTTTTTCAGTATTTTAATTGCTGCCGGACCTTCGGAAGTCACCGGGA‑TTTTCACCACAATACCCGGAATAGc                                <  1:819271/80‑1 (MQ=255)                                            aTACCCTCTTTTTTCAGTATTTTAATTGCTGCCGGACCTTCGGAAGTCACCGGGA‑TTTTCACCACAATACCCGGAATAGCg                               <  2:819270/81‑1 (MQ=255)                                             tACCCTCTTTTTTCAGTATTTTAATTGCTGCCGGACCTTCGGAAGTCACCGGGA‑TTTTCACCACAATACCCGGAATAGCGt                              >  2:819297/1‑81 (MQ=255)                                              aCCCTCTTTTTTCAGTATTTTAATTGCTGCCGGACCTTCGGAAGTCACCGGGA‑TTTTCACCACAATACCCGGAATAGCGTc                             >  1:819304/1‑81 (MQ=255)                                                ccTCTTTTTTCAGTATTTTAATTGCTGCCGGACCTTCGGAAGTCACCGGGA‑TTTTCACCACAATACCCGGAATAGCGTcgc                           >  1:819321/1‑81 (MQ=255)                                                 ctctTTTTTCAGTATTTTAATTGCTGCCGGACCTTCGGAAGTCACCGGGA‑TTTTCACCACAATACCCGGAATAGCGTcgcg                          >  2:819315/1‑81 (MQ=255)                                                         tCAGTATTTTAATTGCTGCCGGACCTTCGGAAGTCACCGGGA‑TTTTCACCACAATACCCGGAATAGCGTCGCGCAGGcgct                  >  1:819319/1‑81 (MQ=255)                                                           aGTATTTTAATTGCTGCCGGACCTTCGGAAGTCACCGGGA‑TTTTCACCACAATACCCGGAATAGCGTCGCGCAGGcgcttc                >  1:819353/1‑81 (MQ=255)                                                            gTATTTTAATTGCTGCCGGACCTTCGGAAGTCACCGGGATTTTTCACCACAATACCCGGAATAGCGTCGCGCAGGcgcttc                >  1:819308/1‑81 (MQ=255)                                                                 ttAATTGCTGCCGGACCTTCGGAAGTCACCGGGA‑TTTTCACCACAATACCCGGAATAGCGTCGCGCAGGCGCTTCGcttct          >  2:819312/1‑81 (MQ=255)                                                                  tAATTGCTGCCGGACCTTCGGAAGTCACCGGGA‑TTTTCACCACAATACCCGGAATAGCGTCGCGCAGGCGCTTCGcttcta         >  1:819329/1‑80 (MQ=255)                                                                       gctgcCGGACCTTCGGAAGTCACCGGGA‑TTTTCACCACAATACCCGGAATAGCGTCGCGCAGGCGCTTCGCTTCTTccac     >  1:819303/1‑80 (MQ=255)                                                                       gctgcCGGACCTTCGGAAGTCACCGGGA‑TTTTCACCACAATACCCGGAATAGCGTCGCGCAGGCGCTTCGCTTCTTccac     >  1:819346/1‑80 (MQ=255)                                                                         tgcCGGACCTTCGGAAGTCACCGGGA‑TTTTCACCACAATACCCGGAATAGCGTCGCGCAGGCGCTTCGCTTCTTCCACCAt  >  2:819310/1‑81 (MQ=255)                                                                             |                                                                              TAACCCTTGTGCGGCGCTATATACAGCAGTGCCAAGTGTAGTAATACCCTCTTTTTTCAGTATTTTAATTGCTGCCAGACCTTCGGAAGTCACCGGGA‑TTTTCACCACAATACCCGGAATAGCGTCGCGCAGGCGCTTCGCTTCTTCCACCAT  >  NC_000913/4139355‑4139507

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: ‑