Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,515,985	C→T	R136R (CGC→CGT)	patD →	gamma‑aminobutyraldehyde dehydrogenase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,515,985	0	C	T	100.0%	119.1 / NA	34	R136R (CGC→CGT)	patD	gamma‑aminobutyraldehyde dehydrogenase
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (19/15);  total (19/15)

TTCGCGGGTGCGGCGCGCTGTCTGAATGGTCTGGCGGCAGGTGAATATCTTGAAGGTCATACTTCGATGATCCGTCGCGATCCGTTGGGGGTCGTGGCTTCTATCGCACCGTGGAATTATCCGCTGATGATGGCCGCGTGGAAACTTGCTCC  >  NC_000913/1515908‑1516059                                                                              |                                                                           ttcgcGGGTGCGGCGCGCTGTCTGAATGGTCTGGCGGCAGGTGAATATCTTGAAGGTCATACTTCGATGATCCGTCGTGAt                                                                         >  2:266865/1‑81 (MQ=255)    gcgGGTGCGGCGCGCTGTCTGAATGGTCTGGCGGCAGGTGAATATCTTGAAGGTCATACTTCGATGATCCGTCGTGATCCg                                                                      >  2:266848/1‑81 (MQ=255)     cgGGTGCGGCGCGCTGTCTGAATGGTCTGGCGGCAGGTGAATATCTTGAAGGTCATACTTCGATGATCCGTCGTGATCCGt                                                                     <  1:266805/81‑1 (MQ=255)     cgGGTGCGGCGCGCTGTCTGAATGGTCTGGCGGCAGGTGAATATCTTGAAGGTCATACTTCGATGATCCGTCGTGATCCGt                                                                     >  2:266857/1‑81 (MQ=255)      gggTGCGGCGCGCTGTCTGAATGGTCTGGCGGCAGGTGAATATCTTGAAGGTCATACTTCGATGATCCGTCGTGATCCGtt                                                                    <  2:266806/81‑1 (MQ=255)        gTGCGGCGCGCTGTCTGAATGGTCTGGCGGCAGGTGAATATCTTGAAGGTCATACTTCGATGATCCGTCGTGATCCGTTgg                                                                  >  2:266855/1‑81 (MQ=255)            ggcgCGCTGTCTGAATGGTCTGGCGGCAGGTGAATATCTTGAAGGTCATACTTCGATGATCCGTCGTGATCCGTTGGGGGt                                                              >  2:266847/1‑81 (MQ=255)            ggcgCGCTGTCTGAATGGTCTGGCGGCAGGTGAATATCTTGAAGGTCATACTTCGATGATCCGTCGTGATCCGTTGGGGGt                                                              >  1:266854/1‑81 (MQ=255)            ggcgCGCTGTCTGAATGGTCTGGCGGCAGGTGAATATCTTGAAGGTCATACTTCGATGATCCGTCGTGATCCGTTGGGGGt                                                              >  1:266835/1‑81 (MQ=255)               gcgcTGTCTGAATGGTCTGGCGGCAGGTGAATATCTTGAAGGTCATACTTCGATGATCCGTCGTGATCCGTTGGGGGTCGt                                                           >  2:266858/1‑81 (MQ=255)                cgcTGTCTGAATGGTCTGGCGGCAGGTGAATATCTTGAAGGTCATACTTCGATGATCCGTCGTGATCCGTTGGGGGTCGTg                                                          <  2:266807/81‑1 (MQ=255)                 gcTGTCTGAATGGTCTGGCGGCAGGTGAATATCTTGAAGGTCATACTTCGATGATCCGTCGTGATCCGTTGGGGGTCGTgg                                                         >  1:266843/1‑81 (MQ=255)                     tCTGAATGGTCTGGCGGCAGGTGAATATCATGAAGGTCATACTTCGATGATCCGTCGTGATCCGTTGGGGGTCGTGGCt                                                       <  2:266808/79‑1 (MQ=255)                         aaTGGTCTGGCGGCAGGTGAATATCTTGAAGGTCATACTTCGATGATCCGTCGTGATCCGTTGGGGGTCGTGGCTTCTATc                                                 <  2:266809/81‑1 (MQ=255)                           tGGTCTGGCGGCAGGTGAATATCTTGAAGGTCATACTTCGATGATCCGTCGTGATCCGTTGGGGGTCGTGGCTTCTATCGc                                               <  2:266810/81‑1 (MQ=255)                                 ggcggcAGGTGAATATCTTGAAGGTCATACTTCGATGATCCGTCGTGATCCGTTGGGGGTCGTGGCTTCTATCGCACCGTg                                         >  1:266869/1‑81 (MQ=255)                                    ggcAGGTGAATATCTTGAAGGTCATACTTCGATGATCCGTCGTGATCCGTTGGGGGTCGTGGCTTCTATCGCACCGTGGaa                                      >  1:266864/1‑81 (MQ=255)                                            aaTATCTTGAAGGTCATACTTCGATGATCCGTCGTGATCCGTTGGGGGTCGTGGCTTCTATCGCACCGTGGAATTATCCGc                              <  1:266811/81‑1 (MQ=255)                                            aaTATCTTGAAGGTCATACTTCGATGATCCGTCGTGATCCGTTGGGGGTCGTGGCTTCTATCGCACCGTGGAATTATCCGc                              <  2:266812/81‑1 (MQ=255)                                             atatCTTGAAGGTCATACTTCGATGATCCGTCGTGATCCGTTGGGGGTCGTGGCTTCTATCGCACCGTGGAATTATCCGCt                             <  2:266813/81‑1 (MQ=255)                                             atatCTTGAAGGTCATACTTCGATGATCCGTCGTGATCCGTTGGGGGTCGTGGCTTCTATCGCACCGTGGAATTATCCGCt                             <  1:266814/81‑1 (MQ=255)                                                 cTTGAAGGTCATACTTCGATGATCCGTCGTGATCCGTTGGGGGTCGTGGCTTCTATCGCACCGTGGAATTATCCGCtgat                          >  1:266874/1‑80 (MQ=255)                                                 cTTGAAGGTCATACTTCGATGATCCGTCGTGATCCGTTGGGGGTCGTGGCTTCTATCGCACCGTGGAATTATCCGCtgat                          >  2:266849/1‑80 (MQ=255)                                                   tGAAGGTCATACTTCGATGATCCGTCGTGATCCGTTGGGGGTCGTGGCTTCTATCGCACCGTGGAATTATCCGCTgattat                       >  2:266868/1‑81 (MQ=255)                                                   tGAAGGTCATACTTCGATGATCCGTCGTGATCCGTTGGGGGTCGTGGCTCCTATCGCACCGTGGAATTATCCGCTgatgat                       >  2:266851/1‑81 (MQ=255)                                                    gtaGGTCATACTTCGATGATCCGTCGTGATCCGTTGGGGGTCGTGGCTTCTATCGCACCGTGGAATTATCCGCtgat                          <  2:266815/75‑1 (MQ=255)                                                    gAAGGTCATACTTCGATGATCCGTCGTGATCCGTTGGGGGTCGTGGCTTCTATCGCACCGTGGAATTATCCGCtgat                          <  1:266816/77‑1 (MQ=255)                                                      aGGTCATACTTCGATGATCCGTCGTGATCCGTTGGGGGTCGTGGCTTCTATCGCACCGTGGAATTATCCGCTGATGATGGc                    <  2:266817/81‑1 (MQ=255)                                                         tCATACTTCGATGATCCGTCGTGATCCATTGGGGGTCGTGGCTTCTATCGCACCGTGGAATTATCCGCTGATGATGGCcgc                 <  1:266818/81‑1 (MQ=255)                                                          cATACTTCGATGATCCGTCGTGATCCGTTGGGGGTCGTGGCTTCTATCGCACCGTGGAATTATCCGCTGATGATGGCcgcg                >  2:266852/1‑81 (MQ=255)                                                                 cGATGATCCGTCGTGATCCGTTGGGGGTCGTGGCTTCTATCGCACCGTGGAATTATCCGCTGATGATGGCCGCGTGGAAAc         >  1:266863/1‑81 (MQ=255)                                                                     gatCCGTCGTGATCCGTTGGGGGTCGTGGCTTCTATCGCACCGTGGAATTATCCGCTGATGATGGCCGCGTGGAAACTGGc     >  2:266853/1‑81 (MQ=255)                                                                      atCCGTCGTGATCCGTTGGGGGTCGTGGCTTCTATCGCACCGTGGAATTATCCGCTGATGATGGCCGCGTGGAAACTTGCt    <  2:266819/81‑1 (MQ=255)                                                                        ccGTCGTGATCCGTTGGGGGTCGTGGCTTCTATCGCACCGTGGAATTATCCGCTGATGATGGCCGCGTGGAAACTTGCTcc  >  2:266850/1‑81 (MQ=255)                                                                              |                                                                           TTCGCGGGTGCGGCGCGCTGTCTGAATGGTCTGGCGGCAGGTGAATATCTTGAAGGTCATACTTCGATGATCCGTCGCGATCCGTTGGGGGTCGTGGCTTCTATCGCACCGTGGAATTATCCGCTGATGATGGCCGCGTGGAAACTTGCTCC  >  NC_000913/1515908‑1516059

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