Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	148,936	G→A	T222T (ACC→ACT)	panB ←	3‑methyl‑2‑oxobutanoate hydroxymethyltransferase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	148,936	0	G	A	100.0%	121.1 / NA	33	T222T (ACC→ACT)	panB	3‑methyl‑2‑oxobutanoate hydroxymethyltransferase
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (16/17);  total (16/17)

CATATACTGCCGCACAGCCGCGCGGATGTCGCCCGTTTCGGCGAGGAAATTTTTAGCGAATTTAGGAATGTGACCGCCGGTAATACCAAAGGCGTCGTGCATCACGAGGATCTGCCCGTCAGTGACGTTGCCTGCGCCAATGCCAATAACCGG  >  NC_000913/148858‑149010                                                                               |                                                                           cATATACTGCCGCACAGCCGCGCGGATGTCGCCCGTTTCGGCGAGGAAATTTTTAGCGAATTTAGGAATGTGACCGCCAGt                                                                          <  1:33406/81‑1 (MQ=255)  atatACTGCCGCACAGCCGCGCGGATGTCGCCCGTTTCGGCGAGGAAATTTTTAGCGAATTTAGGAATGTGACCGCCAGTa                                                                         <  1:33407/81‑1 (MQ=255)     taCTGCCGCACAGCCGCGCGGATGTCGCCCGTTTCGGCGAGGAAATTTTTAGCGAATTTAGGAATGTGACCGCCAGTAAt                                                                       <  1:33408/80‑1 (MQ=255)     taCTGCCGCACAGCCGCGCGGATGTCGCCCGTTTCGGCGAGGAAATTTTTAGCGAATTTAGGAATGTGACCGCCAGTAAt                                                                       <  1:33409/80‑1 (MQ=255)          ccGCACAGCCGCGCGGATGTCGCCCGTTTCGGCGAGGAAATTTTTAGCGAATTTAGGAATGTGACCGCCAGTAATACCaa                                                                  <  2:33410/80‑1 (MQ=255)           cGCACAGCCGCGCGGATGTCGCCCGTTTCGGCGAGGAAATTTTTAGCGAATTTAGGAATGTGACCGCCAGTAATACCAAAg                                                                <  2:33412/81‑1 (MQ=255)           cGCACAGCCGCGCGGATGTCGCCCGTTTCGGCGAGGAAATTTTTAGCGAATTTAGGAATGTGACCGCCAGTAATACCAAAg                                                                <  2:33411/81‑1 (MQ=255)               caGCCGCGCGGATGTCGCCCGTTTCGGCGAGGAAATTTTTAGCGAATTTAGGAATGTGACCGCCAGTAATACCAAAGGcgt                                                            <  1:33413/81‑1 (MQ=255)                aGCCGCGCGGATGTCGCCCGTTTCGGCGAGGAAATTTTTAGCGAATTTAGGAATGTGACCGCCAGTAATACCAAAGGcgtc                                                           >  2:716207/1‑81 (MQ=255)                   cgcgcgGATGTCGCCCGTTTCGGCGAGGAAATTTTTAGCGAATTTAGGAATGTGACCGCCAGTAATACCAAAGGCGTCGTg                                                        >  2:33454/1‑81 (MQ=255)                    gcgcgGATGTCGCCCGTTTCGGCGAGGAAATTTTTAGCGAATTTAGGAATGTGACCGCCAGTAATACCAAAGGCGTCGTg                                                        >  1:64406/1‑80 (MQ=255)                        ggATGTCGCCCGTTTCGGCGAGGAAATTTTTAGCGAATTTAGGAATGTGACCGCCAGTAATACCAAAGGCGTCGTGCATCa                                                   >  1:33460/1‑81 (MQ=255)                         gATGTCGCCCGTTTCGGCGAGGAAATTTTTAGCGAATTTAGGAATGTGACCGCCAGTAATACCAAAGGCGTCGTGCATCAc                                                  <  1:33414/81‑1 (MQ=255)                         gATGTCGCCCGTTTCGGCGAGGAAATTTTTAGCGAATTTAGGAATGTGACCGCCAGTAATACCAAAGGCGTCGTGCATCAc                                                  <  1:33415/81‑1 (MQ=255)                            gTCGCCCGTTTCGGCGAGGAAATTTTTAGCGAATTTAGGAATGTGACCGCCAGTAATACCAAAGGCGTCGTGCATCACGAg                                               >  2:33470/1‑81 (MQ=255)                            gTCGCCCGTTTCGGCGAGGAAATTTTTAGCGAATTTAGGAATGTGACCGCCAGTAATACCAAAGGCGTCGTGCATCACGAg                                               >  1:33455/1‑81 (MQ=255)                               gCCCGTTTCGGCGAGGAAATTTTTAGCGAATTTAGGAATGTGACCGCCAGTAATACCAAAGGCGTCGTGCATCACGAGGAt                                            >  1:33469/1‑81 (MQ=255)                                cccGTTTCGGCGAGGAAATTTTTAGCGAATTTAGGAATGTGACCGCCAGTAATACCAAAGGCGTCGTGCATCACGAGGATc                                           <  2:33416/81‑1 (MQ=255)                                   gTTTCGGCGAGGAAATTTTTAGCGAATTTAGGAATGTGACCGCCAGTAATACCAAAGGCGTCGTGCATCACGAGGATCTGc                                        >  2:33472/1‑81 (MQ=255)                                     ttCGGCGAGGAAATTTTTAGCGAATTTAGGAATGTGACCGCCAGTAATACCAACGGCGTCGTGCATCACGAGGATCTGccc                                      >  2:33463/1‑81 (MQ=255)                                        ggCGAGGAAATTTTTAGCGAATTTAGGAATGTGACCGCCAGTAATACCAAAGGCGTCGTGCATCACGAGGATCTGCCCGTc                                   <  2:33417/81‑1 (MQ=255)                                                aaTTTTTAGCGAATTTAGGAATGTGACCGCCAGTAATACCAAAGGCGTCGTGCATCACGAGGATCTGCCCGTCAGTGACGt                           <  1:33418/81‑1 (MQ=255)                                                  tttttAGCGAATTTAGGAATGTGACCGCCAGTAATACCAAAGGCGTCGTGCATCACGAGGATCTGCCCGTCAGTGACGTTg                         <  2:33419/81‑1 (MQ=255)                                                        gCGAATTTAGGAATGTGACCGCCAGTAATACCAAAGGCGTCGTGCATCACGAGGATCTGCCCGTCAGTGACGTTGCCTgc                    <  2:33420/80‑1 (MQ=255)                                                          gAATTTAGGAATGTGACCGCCAGTAATACCAAAGGCGTCGTGCATCACGAGGATCTGCCCGTCAGTGACGTTGCCTgcgc                  >  1:33466/1‑80 (MQ=255)                                                             tttAGGAATGTGACCGCCAGTAATACCAAAGGCGTCGTGCATCACGAGGATCTGCCCGTCAGTGACGTTGCCTGCgccaat              >  2:33468/1‑81 (MQ=255)                                                              ttAGGAATGTGACCGCCAGTAATACCAAAGGCGTCGTGCATCACGAGGATCTGCCCGTCAGTGACGTTGCCTGCgccaatg             >  1:33459/1‑81 (MQ=255)                                                                aGGAATGTGACCGCCAGTAATACCAAAGGCGTCGTGCATCACGAGGATCTGCCCGTCAGTGACGTTGCCTGCgccaatgc            >  2:33461/1‑80 (MQ=255)                                                                    aTGTGACCGCCAGTAATACCAAAGGCGTCGTGCATCACGAGGATCTGCCCGTCAGTGACGTTGCCTGCgccaatgccaat        >  1:881295/1‑80 (MQ=255)                                                                      gtgACCGCCAGTAATACCAAAGGCGTCGTGCATCACGAGGATCTGCCCGTCAGTGACGTTGCCTGCGCCAATGCCAATAAc     >  2:33462/1‑81 (MQ=255)                                                                       tgACCGCCAGTAATACCAAAGGCGTCGTGCATCACGAGGATCTGCCCGTCAGTGACGTTGCCTGCGCCAATGCCAATAAcc    <  1:33421/81‑1 (MQ=255)                                                                        gACCGCCAGTAATACCAAAGGCGTCGTGCATCACGAGGATCTGCCCGTCAGTGACGTTGCCTGCGCCAATGCCAATAACCg   >  2:33465/1‑81 (MQ=255)                                                                         aCCGCCAGTAATACCAAAGGCGTCGTGCATCACGAGGATCTGCCCGTCAGTGACGTTGCCTGCGCCAATGCCAATAACCg   <  2:33422/80‑1 (MQ=255)                                                                          ccgccaGTAATACCAAAGGCGTCGTGCATCACGAGGATCTGCCCGTCAGTGACGTTGCCTGCGCCAATGCCAATAACCgg  >  1:33464/1‑80 (MQ=255)                                                                               |                                                                           CATATACTGCCGCACAGCCGCGCGGATGTCGCCCGTTTCGGCGAGGAAATTTTTAGCGAATTTAGGAATGTGACCGCCGGTAATACCAAAGGCGTCGTGCATCACGAGGATCTGCCCGTCAGTGACGTTGCCTGCGCCAATGCCAATAACCGG  >  NC_000913/148858‑149010

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