Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,361,590	T→C	H157H (CAT→CAC)	puuD →	gamma‑glutamyl‑gamma‑aminobutyrate hydrolase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,361,590	0	T	C	100.0%	69.1 / NA	21	H157H (CAT→CAC)	puuD	gamma‑glutamyl‑gamma‑aminobutyrate hydrolase
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (12/9);  total (12/9)

TACAAGAACTGGTGGTGGCAACCGGTGGTTCGTTGCATCGCAAGCTGTGCGAACAGCCTGAATTGCTGGAACATCGGGAAGATCCCGAACTGCCGGTGGAACAGCAATATGCACCGTCGCATGAAGTTCAGGTTGAAGAGGGGGGATTACTGTCTGCGTTGTTACCTGAATGTAGCAACTTTTGGGTAAACTCTCTACATGGACAAGGGGCGAAGGTCGTTAGCCCACGGTTGCGTGTTGAAGCTCGCTCG  >  NC_000913/1361469‑1361719                                                                                                                          |                                                                                                                                  tACAAGAACTGGTGGTGGCAACCGGTGGTTCGTTGCATCGCAAGCTGTGCGAACAGCCTGAATTGCTGGAACATCGGGAAGATCCCGAACTGCCGGTGGAACAGCAATATACACTGTCGCACGAAGTTCAGGt                                                                                                                        >  1:186466/1‑133 (MQ=255)       aaCTGGTGGTGGCAACCGGTGGTTCGTTGCATCGCAAGCTGTGCGAACAGCCTGAATTGCTGGAACATCGGGAAGATCCCGAACTGCCGGTGGAACAGCAATATGCACCGTCGCACGAAGTTCAGGTTGAaga                                                                                                                  >  2:475554/1‑133 (MQ=255)         cTGGTGGTGGCAACCGGTGGTTCGTTGCATCGCAAGCTGTGCGAACAGCCTGAATTGCTGGAACATCGGGAAGATCCCGAACTGCCGGTGGAACAGCAATATGCACCGTCGCACGAAGTTCAGGTTGAAGAgg                                                                                                                <  1:40023/133‑1 (MQ=255)         cTGGTGGTGGCAACCGGTGGTTCGTTGCATCGCAAGCTGTGCGAACAGCCTGAATTGCTGGAACATCGGGAAGATCCCGAACTGCCGGTGGAACAGCAATATGCACCGTCGCACGAAGTTCAGGTTGAAGAgg                                                                                                                <  1:475554/133‑1 (MQ=255)                            gtTCGTTGCATCGCAAGCTGTGCGAACAGCCTGAATTGCTGGAACATCGGGAAGATCCCGAACTGCCGGTGGAACAGCAATATGCACCGTCGCACGAAGTTCAGGTTGAAGAGGGGGGATTACTGTCTGCgtt                                                                                             <  2:417663/133‑1 (MQ=255)                                 ttGCATCGCAAGCTGTGCGAACAGCCTGAATTGCTGGAACATCGGGAAGATCCCGAACTGCCGGTGGAACAGCAATATGCACCGTCGCACGAAGTTCAGGTTGAAGAGGGGGGATTACTGTCTGCGTTGTTAc                                                                                        >  1:76307/1‑133 (MQ=255)                                                 gCGAACAGCCTGAATTGCTGGAACATCGGGAAGATCCCGAACTGCCGGTGGAACAGCAATATGCACCGTCGCACGAAGTTCAGGTTGAAGAGGGGGGATTACTGTCTGCGTTGTTACCTGAATGTAGCAACtt                                                                        >  2:284171/1‑133 (MQ=255)                                                          cTGAATTGCTGGAACATCGGGAAGATCCCGAACTGCCGGTGGAACAGCAATATGCACCGTCGCACGAAGTTCAGGTTGAAGAGGGGGGATTACTGTCTGCGTTGTTACCTGAATGTAGCAACTTTTGGGTaaa                                                               >  1:321574/1‑133 (MQ=255)                                                          cTGAATTGCTGGAACATCGGGAAGATCCCGAACTGCCGGTGGAACAGCAATATGCACCGTCGCACGAAGTTCAGGTTGAAGAGGGGGGATTACTGTCTGCGTTGTTACCTGAATGTAGCAACTTTTGGGTaaa                                                               >  1:401505/1‑133 (MQ=255)                                                              aTTGCTGGAACATCGGGAAGATCCCGAACTGCCGGTGGAACAGCAATATGCACCGTCGCACGAAGTTCAGGTTGAAGAGGGGGGATTACTGTCTGCGTTGTTACCTGAATGTAGCAACTTTTGGGTAAActct                                                           <  2:401505/133‑1 (MQ=255)                                                                tGCTGGAACATCGGGAAGATCCCGAACTGCCGGTGGAACAGCAATATGCACCGTCGCACGAAGTTCAGGTTGAAGAGGGGGGATTACTGTCTGCGTTGTTACCTGAATGTAGCAACTTTTGGGTAAActctct                                                         <  2:438694/133‑1 (MQ=255)                                                                  cTGGAACATCGGGAAGATCCCGAACTGCCGGTGGAACAGCAATATGCACCGTCGCACGAAGTTCAGGTTGAAGAGGGGGGATTACTGTCTGCGTTGTTACCTGAATGTAGCAACTTTTGGGt                                                                  <  1:31118/122‑1 (MQ=255)                                                                  cTGGAACATCGGGAAGATCCCGAACTGCCGGTGGAACAGCAATATGCACCGTCGCACGAAGTTCAGGTTGAAGAGGGGGGATTACTGTCTGCGTTGTTACCTGAATGTAGCAACTTTTGGGt                                                                  >  2:31118/1‑122 (MQ=255)                                                                          tCGGGAAGATCCCGAACTGCCGGTGGAACAGCAATATGCACCGTCGCACGAAGTTCAGGTTGAAGAGGGGGGATTACTGTCTGCGTTGTTAACTGAATGTCGCATCTTATGGGTAAAGTCTCGACATGGACaa                                               >  2:155683/1‑133 (MQ=255)                                                                                aGATCCCGAACTGCCTGTGGAACACCAATATGAACCGTCGCTCGAAGTTCAGGTTGAAGAGGGGGGATTACTGTCTGCGTTGTTACCTGAATGTAGCAACTTTTGGGTAAACTCTCTACATGGACAAGGGGCg                                         <  2:185856/133‑1 (MQ=255)                                                                                               ggTGGAACAGCAATATGCACCGTCGCACGAAGTTCAGGTTGAAGAGGGGGGATTACTGTCTGCGTTGTTACCTGAATGTAGCAACTTTTGGGTAAACTCTCTACATGGACAAGGGGCGAAGGTCGTTAGCCCa                          <  2:76307/133‑1 (MQ=255)                                                                                                                 aCCGTCACACGAAGTTCAGGTTGAAGAGGGGGGATTACTGTCTGCGTTGTTACCTGAATGTAGCAACTTTTGGGTAAACTCTCTACATGGACAAGGGGCGAAGGTCGTTAGCCCTCGATTGCGTGTTGAAgct        >  1:140179/1‑133 (MQ=255)                                                                                                                   cGTCGCACGAAGTTCAGGTTGAAGAGGGGGGATTACTGTCTGCGTTGTTACCTGAATGTAGCAACTTTTGGGTAAACTCTCTACATGGACACGGGGCGAAAGTCCTTAGCCCACGGTCGCGTGTTGAAgctcg      >  2:176669/1‑133 (MQ=255)                                                                                                                   cGTCGCACGAAGTTCAGGTTGAAGAGGGGGGATTACTGTCTGCGTTGTTACCTGAATGTAGCAACTTTTGGGTAAACTCTCTACATGGACAAGGGGCGAAGGTCGTTAGCCCACGGTTGCGTGTTGAAgctcg      >  2:277914/1‑133 (MQ=255)                                                                                                                   cGTCGCACGAAGTTCAGGTTGAAGAGGGGGGATTACTGTCTGCGTTGTTACCTGAATGTAGCAACTTTTGGGTAAACTCTCTACATGGACAAGGGGCGAAGGTCGTTAGCCCACGGTTGCGTGTTGAAgctcg      >  2:420939/1‑133 (MQ=255)                                                                                                                       gCACGAAGTTCAGGTTGAAGAGGGGGGATTACTGTCTGCGTTGTTACCTGAATGTAGCAACTTTTGGGTAAACTCTCTACATGGACAAGGGGCGAAGGTCGTTAGCCCACGGTTGCGTGTTGAAGctcgctcg  <  2:321574/133‑1 (MQ=255)                                                                                                                          |                                                                                                                                  TACAAGAACTGGTGGTGGCAACCGGTGGTTCGTTGCATCGCAAGCTGTGCGAACAGCCTGAATTGCTGGAACATCGGGAAGATCCCGAACTGCCGGTGGAACAGCAATATGCACCGTCGCATGAAGTTCAGGTTGAAGAGGGGGGATTACTGTCTGCGTTGTTACCTGAATGTAGCAACTTTTGGGTAAACTCTCTACATGGACAAGGGGCGAAGGTCGTTAGCCCACGGTTGCGTGTTGAAGCTCGCTCG  >  NC_000913/1361469‑1361719

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑