Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,361,613	(G)6→7	coding (494/765 nt)	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,607	1	.	G	100.0%	125.0 / NA	35	E163G (GAG→GGG)	puuD	gamma‑glutamyl‑gamma‑aminobutyrate hydrolase
Reads supporting (aligned to +/- strand):  ref base . (0/0);  new base G (20/15);  total (20/15)

CAAGAACTGGTGGTGGCAACCGGTGGTTCGTTGCATCGCAAGCTGTGCGAACAGCCTGAATTGCTGGAACATCGGGAAGATCCCGAACTGCCGGTGGAACAGCAATATGCACCGTCGCATGAAGTTCAGGTTGAAGA‑GGGGGGATTACTGTCTGCGTTGTTACCTGAATGTAGCAACTTTTGGGTAAACTCTCTACATGGACAAGGGGCGAAGGTCGTTAGCCCACGGTTGCGTGTTGAAGCTCGCTCGCCGGATGGCTTAGTTGAGGCGG  >  NC_000913/1361471‑1361741                                                                                                                                          |                                                                                                                                       cAAGAACTGGTGGTGGCAACCGGTGGTTCGTTGCATCGCAAGCTGTGCGAACAGCCTGAATTGCTGGAACATCGGGAAGATCCCGAACTGCCGGTGGAACAGCAATATGCACCGTCGCATGAAGTTCAGGTTGAAGA‑ggg                                                                                                                                     <  2:877316/140‑1 (MQ=255)       cTGGTGGTGGCAACCGGTGGTTCGTTGCATCGCAAGCTGTGCGAACAGCCTGAATTGCTGGAACATCGGGAAGATCCCGAACTGCCGGTGGAACAGCAATATGCACCGTCGCATGAAGTTCAGGTTGAAGA‑gggggggat                                                                                                                               <  1:438730/140‑4 (MQ=255)       cTGGTGGTGGCAACCGGTGGTTCGTTGCATCGCAAGCTGTGCGAACAGCCTGAATTGCTGGAACATCGGGAAGATCCCGAACTGCCGGTGGAACAGCAATATGCACCGTCGCATGAAGTTCAGGTTGAAGA‑gggggggat                                                                                                                               >  1:837544/1‑137 (MQ=255)          gtggtggCAACCGGTGGTTCGTTGCATCGCAAGCTGTGCGAACAGCCTGAATTGCTGGAACATCGGGAAGATCCCGAACTGCCGGTGGAACAGCAATATGCACCGTCGCATGAAGTTCAGGTTGAAGAGGGGGGGATTAc                                                                                                                             >  2:812412/1‑140 (MQ=255)               ggCAACCGGTGGTTCGTTGCATCGCAAGCTGTGCGAACAGCCTGAATTGCTGGAACATCGGGAAGATCCCGAACTGCCGGTGGAACAGCAATATGCACCGTCGCATGAAGTTCAGGTTGAAGAGGGGGGGATTACTGTCt                                                                                                                        <  2:107715/140‑1 (MQ=255)                    ccGGTGGTTCGTTGCATCGCAAGCTGTGCGAACAGCCTGAATTGCTGGAACATCGGGAAGATCCCGAACTGCCGGTGGAACAGCAATATGCACCGTCGCATGAAGTTCAGGTTGAAGAGGGGGGGATTACTGTCTGc                                                                                                                      <  2:53065/137‑1 (MQ=255)                    ccGGTGGTTCGTTGCATCGCAAGCTGTGCGAACAGCCTGAATTGCTGGAACATCGGGAAGATCCCGAACTGCCGGTGGAACAGCAATATGCACCGTCGCATGAAGTTCAGGTTGAAGAGGGGGGGATTACTGTCTGc                                                                                                                      >  1:53065/1‑137 (MQ=255)                     cGGTGGTTCGTTGCATCGCAAGCTGTGCGAACAGCCTGAATTGCTGGAACATCGGGAAGATCCCGAACTGCCGGTGGAACAGCAATATGCACCGTCGCATGAAGTTCAGGTTGAAGAGGGGGGGATTACTGTCTGCgttg                                                                                                                  <  1:102073/140‑1 (MQ=255)                            tCGTTGCATCGCAAGCTGTGCGAACAGCCTGAATTGCTGGAACATCGGGAAGATCCCGAACTGCCGGTGGAACAGCAATATGCACCGTCGCATGAAGTTCAGGTTGAAGAGGGGGGGATTACTGTCTGCGTTGTTACCTg                                                                                                           <  2:196550/140‑1 (MQ=255)                               ttGCATCGCAAGCTGTGCGAACAGCCTGAATTGCTGGAACATCGGGAAGATCCCGAACTGCCGGTGGAACAGCAATATGCACCGTCGCATGAAGTTCAGGTTGAAGAGGGGGGGATTACTGTCTGCGTTGTTACCTGAAt                                                                                                        >  2:17892/1‑140 (MQ=255)                                       cAAGCTGTGCGAACAGCCTGAATTGCTGGAACATCGGGAAGATCCCGAACTGCCGGTGGAACAGCAATATGCACCGTCGCATGAAGTTCAGGTTGAAGAGGGGGGGATTACTGTCTGCGTTGTTACCTGAATGTAGCAAc                                                                                                >  1:858359/1‑140 (MQ=255)                                              tgCGAACAGCCTGAATTGCTGGAACATCGGGAAGATCCCGAACTGCCGGTGGAACAGCAATATGCACCGTCGCATGAAGTTCAGGTTGAAGAGGGGGGGATTACTGTCTGCGTTGTTACCTGAATGTAGCAACTTTTggg                                                                                         >  2:76402/1‑140 (MQ=255)                                              tgCGAACAGCCTGAATTGCTGGAACATCGGGAAGATCCCGAACTGCCGGTGGAACAGCAATATGCACCGTCGCATGAAGTTCAGGTTGAAGAGGGGGGGATTACTGTCTGCGTTGTTACCTGAATGTAGCAACTTTTggg                                                                                         >  1:762762/1‑140 (MQ=255)                                               gCGAACAGCCTGAATTGCTGGAACATCGGGAAGATCCCGAACTGCCGGTGGAACAGCAATATGCACCGTCGCATGAAGTTCAGGTTGAAGAGGGGGGGATTACTGTCTGCGTTGTTACCTGAATGTAGCAACTTTTGGGt                                                                                        <  2:1004176/140‑1 (MQ=255)                                                cGAACAGCCTGAATTGCTGGAACATCGGGAAGATCCCGAACTGCCGGTGGAACAGCAATATGCACCGTCGCATGAAGTTCAGGTTGAAGAGGGGGGGATTACTGTCTGCGTTGTTACCTGAATGTAGCAACTTTTGGGTa                                                                                       <  1:1026047/140‑1 (MQ=255)                                                cGAACAGCCTGAATTGCTGGAACATCGGGAAGATCCCGAACTGCCGGTGGAACAGCAATATGCACCGTCGCATGAAGTTCAGGTTGAAGAGGGGGGGATTACTGTCTGCGTTGTTACCTGAATGTAGCAACTTTTGGGTa                                                                                       >  1:162465/1‑140 (MQ=255)                                                 gAACAGCCTGAATTGCTGGAACATCGGGAAGATCCCGAACTGCCGGTGGAACAGCAATATGCACCGTCGCATGAAGTTCAGGTTGAAGAGGGGGGGATTACTGTCTGCGTTGTTACCTGAATGTAGCAACTTTTGGGt                                                                                        >  1:512276/1‑138 (MQ=255)                                                 gAACAGCCTGAATTGCTGGAACATCGGGAAGATCCCGAACTGCCGGTGGAACAGCAATATGCACCGTCGCATGAAGTTCAGGTTGAAGAGGGGGGGATTACTGTCTGCGTTGTTACCTGAATGTAGCAACTTTTGGGt                                                                                        <  2:512276/138‑1 (MQ=255)                                                  aaCAGCCTGAATTGCTGGAACATCGGGAAGATCCCGAACTGCCGGTGGAACAGCAATATGCACCGTCGCATGAAGTTCAGGTTGAAGAGGGGGGGATTACTGTCTGCGTTGTTACCTGAATGTAGCAACTTTTGGGTaaa                                                                                     <  2:762762/140‑1 (MQ=255)                                                  aaCAGCCTGAATTGCTGGAACATCGGGAAGATCCCGAACTGCCGGTGGAACAGCAATATGCACCGTCGCATGAAGTTCAGGTTGAAGAGGGGGGGATTACTGTCTGCGTTGTTACCTGAATGTAGCAACTTTTGGGTaaa                                                                                     <  1:17892/140‑1 (MQ=255)                                                    cAGCCTGAATTGCTGGAACATCGGGAAGATCCCGAACTGCCGGTGGAACAGCAATATGCACCGTCGCATGAAGTTCAGGTTGAAGAGGGGGGGATTACTGTCTGCGTTGTTACCTGAATGTAGCAACTTTTGGGTAAAct                                                                                   <  2:268714/140‑1 (MQ=255)                                                     aGCCTGAATTGCTGGAACATCGGGAAGATCCCGAACTGCCGGTGGAACAGCAATATGCACCGTCGCATGAAGTTCAGGTTGAAGAGGGGGGGATTACTGTCTGCGTTGTTACCTGAATGTAGCAACTTTTGGGTAAActc                                                                                  <  1:812412/140‑1 (MQ=255)                                                      gCCTGAATTGCTGGAACATCGGGAAGATCCCGAACTGCCGGTGGAACAGCAATATGCACCGTCGCATGAAGTTCAGGTTGAAGAGGGGGGGATTACTGTCTGCGTTGTTACCTGAATGTAGCAACTTTTGGGTAAActct                                                                                 <  1:973351/140‑1 (MQ=255)                                                        cTGAATTGCTGGAACATCGGGAAGATCCCGAACTGCCGGTGGAACAGCAATATGCACCGTCGCATGAAGTTCAGGTTGAAGAGGGGGGGATTACTGTCTGCGTTGTTACCTGAATGTAGCAACTTTTGGGTAAActctct                                                                               >  1:970506/1‑140 (MQ=255)                                                                 tGGAACATCGGGAAGATCCCGAACTGCCGGTGGAACAGCAATATGCACCGTCGCATGAAGTTCAGGTTGAAGAGGGGGGGATTACTGTCTGCGTTGTTACCTGAATGTAGCGACTTTTGGGTAAACTCTCTACATGGACa                                                                      >  2:617750/1‑140 (MQ=255)                                                                 tGGAACATCGGGAAGATCCCGAACTGCCGGTGGAACAGCAATATGCACCGTCGCATGAAGTTCAGGTTGAAGAGGGGGGGATTACTGTCTGCGTTGTTACCTGAATGTAGCAACTTTTGGGTAAACTCTCTACATGGACa                                                                      >  1:831694/1‑140 (MQ=255)                                                                                       aCTGCCGGTGGAACAGCAATATGCACCGTCGCATGAAGTTCAGGTTGAAGAGGGGGGGATTACTGTCTGCGTTGTTACCTGAATGTAGCAACTTTTGGGTAAACTCTCTACATGGACAAGGGGCGAAGGTCGTTAGCCCa                                                <  2:323499/140‑1 (MQ=255)                                                                                                            tGCACCGTCGCATGAAGTTCAGGTTGAAGAGGGGGGGATTACTGTCTGCGTTGTTACCTGAATGTAGCAACTTTTGGGTAAACTCTCTACATGGACAAGGGGCGAAGGTCGTTAGCCCACGGTTGCGTGTTGAAgctcgc                           >  1:157413/1‑140 (MQ=255)                                                                                                             gCACCGTCGCATGAAGTTCAGGTTGAAGAGGGGGGGATTACTGTCTGCGTTGTTACCTGAATGTAGCAACTTTTGGGTAAACTCTCTACATGGACAAGGGGCGAAGGTCGTTAGCCCACGGTTGCGTGTTGAAgctcgct                          >  1:454431/1‑140 (MQ=255)                                                                                                                 cGTCGCATGAAGTTCAGGTTGAAGAGGGGGGGATTACTGTCTGCGTTGTTACCTGAATGTAGCAACTTTTGGGTAAACTCTCTACATGGACAAGGGGCGAAGGTCGTTAGCCCACGGTTGCGTGTTGAAGCTCGCTCGcc                      >  2:513780/1‑140 (MQ=255)                                                                                                                 cGTCGCATGAAGTTCAGGTTGAAGAGGGGGGGATTACTGTCTGCGTTGTTACCTGAATGTAGCAACTTTTGGGTAAACTCTCTACATGGACAAGGGGCGAAGGTCGTTAGCCCACGGTTGCGTGTTGAAGCTCGCTCGcc                      >  2:824174/1‑140 (MQ=255)                                                                                                                   tCGCATGAAGTTCAGGTTGAAGAGGGGGGGGTTACTGTCTGCGTTGTTACCTGAATGTAGCAACTTTTGGGTAAACTCTCTACATGGACAAGGGGCGAAGGTCGTTAGCCCACGGTTGCGTGTTGAAGCTCGCTCGCCgg                    >  2:264240/1‑140 (MQ=255)                                                                                                                      cATGAAGTTCAGGTTGAAGAGGGGGGGATTACTGTCTGCGTTGTTACCTGAATGTAGCAACTTTTGGGTAAACTCTCTACATGGACAAGGGGCGAAGGTCGTTAGCCCACGGTTGCGTGTTGAAGCTCGCTCGCCGGATg                 >  1:314136/1‑140 (MQ=255)                                                                                                                      cATGAAGTTCAGGTTGAAGAGGGGGGGATTACTGTCTGCGTTGTTACCTGAATGTAGCAACTTTTGGGTAAACTCTCTACATGGACAAGGGGCGAAGGTCGTTAGCCCACGGTTGCGTGTTGAAGCTCGCTCGCCGGATg                 >  1:326011/1‑140 (MQ=255)                                                                                                                      cATGAAGTTCAGGTTGAAGAGGGGGGGATTACTGTCTGCGTTGTTACCTGAATGTAGCAACTTTTGGGTAAACTCTCTACATGGACAAGGGGCGAAGGTCGTTAGCCCACGGTTGCGTGTTGAAGCTCGCTCGCCGGATg                 >  2:58891/1‑140 (MQ=255)                                                                                                                                aGGTTGAAGAGGGGGGGATTACTGTCTGCGTTGTTACCTGAATGTAGCAACTTTTGGGTAAACTCTCTACATGGACAAGGGGCGAAGGTCGTTAGCCCACGGTTGCGTGTTGAAGCTCGCTCGCCGGATGGCTTAGTTGa       <  2:454431/140‑1 (MQ=255)                                                                                                                                 ggTTGAAGAGGGGGGGATTACTGTCTGCGTTGTTACCTGAATGTAGCAACTTTTGGGTAAACTCTCTACATGGACAAGGGGCGAAGGTCGTTAGCCCACGGTTGCGTGTTGAAGCTCGCTCGCCGGATGGCTTAGTTGAg      <  2:168263/140‑1 (MQ=255)                                                                                                                                     gAAGAGGGGGGGATTACTGTCTGCGTTGTTACCTGAATGTAGCAACTTTTGGGTAAACTCTCTACATGGACAAGGGGCGAAGGTCGTTAGCCCACGGTTGCGTGTTGAAGCTCGCTCGCCGGATGGCTTAGTTGAGGCgg  >  1:146967/1‑140 (MQ=255)                                                                                                                                          |                                                                                                                                       CAAGAACTGGTGGTGGCAACCGGTGGTTCGTTGCATCGCAAGCTGTGCGAACAGCCTGAATTGCTGGAACATCGGGAAGATCCCGAACTGCCGGTGGAACAGCAATATGCACCGTCGCATGAAGTTCAGGTTGAAGA‑GGGGGGATTACTGTCTGCGTTGTTACCTGAATGTAGCAACTTTTGGGTAAACTCTCTACATGGACAAGGGGCGAAGGTCGTTAGCCCACGGTTGCGTGTTGAAGCTCGCTCGCCGGATGGCTTAGTTGAGGCGG  >  NC_000913/1361471‑1361741

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