Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,516,460	A→G	T295A (ACG→GCG)	patD →	gamma‑aminobutyraldehyde dehydrogenase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,516,460	0	A	G	100.0%	141.2 / NA	40	T295A (ACG→GCG)	patD	gamma‑aminobutyraldehyde dehydrogenase
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (21/19);  total (21/19)

AAAGCGCCAGTGATTGTTTTTGATGATGCGGATATTGAAGCAGTGGTCGAAGGTGTACGTACATTTGGCTATTACAATGCTGGACAGGATTGTACTGCGGCTTGTCGGATCTACGCGCAAAAAGGCATTTACGATACGCTGGTGGAAAAACTGGGTGCTGCGGTGGCAACGTTAAAATCTGGTGCGCCAGATGACGAGTCTACGGAGCTTGGACCTTTAAGCTCGCTGGCGCATCTCGAACGCGTCGGCAAGGCAGTAGA  >  NC_000913/1516325‑1516584                                                                                                                                        |                                                                                                                             aaaGCGCCAGTGATTGTTTTTGATGATGCGGATATTGAAGCAGTGGTCGAAGGTGTACGTACATTTGGCTATTACAATGCTGGACAGGATTGTACTGCGGCTTGTCGGATCTACGCGCAAAAAGGCATTTACGATGCGCt                                                                                                                          <  1:748007/140‑1 (MQ=255)                gTTTTTGATGATGCGGATATTGAAGCAGTGGTCGAAGGTGTACGTACATTTGGCTATTACAATGCTGGACAGGATTGTACTGCGGCTTGTCGGATCTACGCGCAAAAAGGCATTTACGATGCGCTGGTGGAAAAACTggg                                                                                                           <  1:150896/140‑1 (MQ=255)                  ttttGATGATGCGGATATTGAAGCAGTGGTCGAAGGTGTACGTACATTTGGCTATTACAATGCTGGACAGGATTGTACTGCGGCTTGTCGGATCTACGCGCAAAAAGGCATTTACGATGCGCTGGTGGAAAAACTGGGtg                                                                                                         <  2:515042/140‑1 (MQ=255)                    ttGATGATGCGGATATTGAAGCAGTGGTCGAAGGTGTACGTACATTTGGCTATTACAATGCTGGACAGGATTGTACTGCGGCTTGTCGGATCTACGCGCAAAAAGGCATTTACGATGCGCTGGTGGAAAAACTGGGtgct                                                                                                       <  1:24281/140‑1 (MQ=255)                     tgatgaTGCGGATATTGAAGCAGTGGTCGAAGGTGTACGTACATTTGGCTATTACAATGCTGGACAGGATTGTACTGCGGCTTGTCGGATCTACGCGCAAAAAGGCATTTACGATGCGCTGGTGGAAAAACTGGGtgctg                                                                                                      >  2:314469/1‑140 (MQ=255)                           tgCGGATATTGAAGCAGTGGTCGAAGGTGTACGTACATTTGGCTATTACAATGCTGGACAGGATTGTACTGCGGCTTGTCGGATCTACGCGCAAAAAGGCATTTACGATGCGCTGGTGGAAAAACTGGGTGCTGAGGTgg                                                                                                >  2:733877/1‑140 (MQ=255)                              ggATATTGAAGCAGTGGTCGAAGGTGTACGTACATTTGGCTATTACAATGCTGGACAGGATTGTACTGCGGCTTGTCGGATCTACGCGCAAAAAGGCATTTACGATGCGCTGGTGGAAAAACTGGGTGCTGCGGTGGCaa                                                                                             >  1:168681/1‑140 (MQ=255)                                  atTGAAGCAGTGGTCGAAGGTGTACGTACATTTGGCTATTACAATGCTGGACAGGATTGTACTGCGGCTTGTCGGATCTACGCGCAAAAAGGCATTTACGATGCGCTGGTGGAAAAACTGGGTGCTGCGGTGGCAACGtt                                                                                         <  2:25240/140‑1 (MQ=255)                                       aGCAGTGGTCGAAGGTGTACGTACATTTGGCTATTACAATGCTGGACAGGATTGTACTGCGGCTTGTCGGATCTACGCGCAAAAAGGCATTTACGATGCGCTGGTGGAAAAACTGGGTGCTGCGGTGGCAACGTTAAAAt                                                                                    >  2:530686/1‑140 (MQ=255)                                              gTCGAAGGTGTACGTACATTTGGCTATTACAATGCTGGACAGGATTGTACTGCGGCTTGTCGGATCTACGCGCAAAAAGGCATTTACGATGCGCTGGTGGAAAAACTGGGTGCTGCGGTGGCAACGTTAAAATCTGGTgc                                                                             >  2:397567/1‑140 (MQ=255)                                              gTCGAAGGTGTACGTACATTTGGCTATTACAATGCTGGACAGGATTGTACTGCGGCTTGTCGGATCTACGCGCAAAAAGGCATTTACGATGCGCTGGTGGAAAAACTGGGTGCTGCGGTGGCAACGTTAAAATCTGGTgc                                                                             <  1:314469/140‑1 (MQ=255)                                                cGAAGGTGTACGTACATTTGGCTATTACAATGCTGGACAGGATTGTACTGCGGCTTGTCGGATCTACGCGCAAAAAGGCATTTACGATGCGCTGGTGGAAAAACTGGGTGCTGCGGTGGCAACGTTAAAATCTTGTgcgc                                                                           <  2:476295/140‑1 (MQ=255)                                                  aaGGTGTACGTACATTTGGCTATTACAATGCTGGACAGGATTGTACTGCGGCTTGTCGGATCTACGCGCAAAAAGGCATTTACGATGCGCTGGTGGAAAAACTGGGTGCTGCGGTGGCAACGTTAAAATCTGGTGCGCCa                                                                         <  2:209269/140‑1 (MQ=255)                                                     gtgtACGTACATTTGGCTATTACAATGCTGGACAGGATTGTACTGCGGCTTGTCGGATCTACGCGCAAAAAGGCATTTACGATGCGCTGGTGGAAAAACTGGGTGCTGCGGTGGCAACGTTAAAATCTGGTGCGCCAGAt                                                                      >  1:505195/1‑140 (MQ=255)                                                       gtacgtacATTTGGCTATTACAATGCTGGACAGGATTGTACTGCGGCTTGTCGGATCTACGCGCAAAAAGGCATTTACGATGCGCTGGTGGAAAAACTGGGTGCTGCGGTGGCAACGTTAAAATCTGGTGCGCCAGATGa                                                                    <  1:733877/140‑1 (MQ=255)                                                       gtacgtacATTTGGCTATTACAATGCTGGACAGGATTGTACTGCGGCTTGTCGGATCTACGCGCAAAAAGGCATTTACGATGCGCTGGTGGAAAAACTGGGTGCTGCGGTGGCAACGTTAAAATCTGGTGCGCCAGATGa                                                                    <  1:397567/140‑1 (MQ=255)                                                        tacgtacATTTGGCTATTACAATGCTGGACAGGATTGTACTGCGGCTTGTCGGATCTACGCGCAAAAAGGCATTTACGATGCGCTGGTGGAAAAACTGGGTGCTGCGGTGGCAACGTTAAAATCTGGTGCGCCAGATGAc                                                                   >  2:641946/1‑140 (MQ=255)                                                                 ttGGCTATTACAATGCTGGACAGGATTGTACTGCGGCTTGTCGGATCTGCGCGCAAAAAGGCATTTACGATGCGCTGGTGGAAAAACTGGGTGCTGCGGTGGCAACGTTAAAATCTGGTGCGCCAGATGACGAGTCTACg                                                          <  1:475524/140‑1 (MQ=255)                                                                 ttGGCTATTACAATGCTGGACAGGATTGTACTGCGGCTTGTCGGATCTACGCGCAAAAAGGCATTTACGATGCGCTGGTGGAAAAACTGGGTGCTGCGGTGGCAACGTTAAAATCTGGTGCGCCAGATGACGAGTCTACg                                                          <  2:685895/140‑1 (MQ=255)                                                                     cTATTACAATGCTGGACAGGATTGTACTGCGGCTTGTCGGATCTACGCGCAAAAAGGCATTTACGATGCGCTGGTGGAAAAACTGGGTGCTGCGGTGGCAACGTTAAAATCTGGTGCGCCAGATGACGAGTCTACGGAGc                                                      >  1:448000/1‑140 (MQ=255)                                                                              tGCTGGACAGGATTGTACTGCGGCTTGTCGGATCTACGCGCAAAAAGGCATTTACGATGCGCTGGTGGAAAAACTGGGTGCTGCGGTGGCAACGTTAAAATCTGGTGCGCCAGATGACGAGTCTACGGAGCTTGGACCtt                                             >  2:77480/1‑140 (MQ=255)                                                                               gCTGGACAGGATTGTACTGCGGCTTGTCGGATCTACGCGCAAAAAGGCATTTACGATGCGCTGGTGGAAAAACTGGGTGCTGCGGTGGCAACGTTAAAATCTGGTGCGCCAGATGACGAGTCTACGGAGCTTGGACCttt                                            >  2:155547/1‑140 (MQ=255)                                                                               gCTGGACAGGATTGTACTGCGGCTTGTCGGATCTACGCGCAAAAAGGCATTTACGATGCGCTGGTGGAAAAACTGGGTGCTGCGGTGGCAACGTTAAAATCTGGTGCGCCAGATGACGAGTCTACGGAGCTTGGACCttt                                            <  2:264054/140‑1 (MQ=255)                                                                                  ggACAGGATTGTACTGCGGCTTGTCGGATCTACGCGCAAAAAGGCATTTACGATGCGCTGGTGGAAAAACTGGGTGCTGCGGTGGCAACGTTAAAATCTGGTGCGCCAGATGACGAGTCTACGGAGCTTGGACCTTTAAg                                         >  1:445988/1‑140 (MQ=255)                                                                                   gACAGGATTGTACTGCGGCTTGTCGGATCTACGCGCAAAAAGGCATTTACGATGCGCTGGTGGAAAAACTGGGTGCTGCGGTGGCAACGTTAAAATCTGGTGCGCCAGATGACGAGTCTACGGAGCTTGGACCTTTAAGc                                        >  2:308154/1‑140 (MQ=255)                                                                                        gATTGTACTGCGGCTTGTCGGATCTACGCGCAAAAAGGCATTTACGATGCGCTGGTGGAAAAACTGGGTGCTGCGGTGGCAACGTTAAAATCTGGTGCGCCAGATGACGAGTCTACGGAGCTTGGACCTTTAAGCTCGCt                                   >  1:776150/1‑140 (MQ=255)                                                                                        gATTGTACTGCGGCTTGTCGGATCTACGCGCAAAAAGGCATTTACGATGCGCTGGTGGAAAAACTGGGTGCTGCGGTGGCAACGTTAAAATCTGGTGCGCCAGATGACGAGTCTACGGAGCTTGGACCTTTAAGCTCGCt                                   >  1:49672/1‑140 (MQ=255)                                                                                          ttGTACTGCGGCTTGTCGGATCTACGCGCAAAAAGGCATTTACGATGCGCTGGTGGAAAAACTGGGTGCTGCGGTGGCAACGTTAAAATCTGGTGCGCCAGATGACGAGTCTACGGAGCTTGGACCTTTAAGCTCGCTgg                                 <  2:448000/140‑1 (MQ=255)                                                                                            gTACTGCGGCTTGTCGGATCTACGCGCAAAAAGGCATTTACGATGCGCTGGTGGAAAAACTGGGTGCTGCGGTGGCAACGTTAAAATCTGGTGCGCCAGATGACGAGTCTACGGAGCTTGGACCTTTAAGCTCGCTGgcg                               <  1:155547/140‑1 (MQ=255)                                                                                            gTACTGCGGCTTGTCGGATCTACGCGCAAAAAGGCATTTACGATGCGCTGGTGGAAAAACTGGGTGCTGCGGTGGCAACGTTAAAATCTGGTGCGCCAGATGACGAGTCTACGGAGCTTGGACCTTTAAGCTCGCTGgcg                               <  1:242034/140‑1 (MQ=255)                                                                                             tACTGCGGCTTGTCGGATCTACGCGCAAAAAGGCATTTACGATGCGCTGGTGGAAAAACTGGGTGCTGCGGTGGCAACGTTAAAATCTGGTGCGCCAGATGACGAGTCTACGGAGCTTGGACCTTTAAGCTCGCTGgcgc                              >  1:122346/1‑140 (MQ=255)                                                                                             tACTGCGGCTTGTCGGATCTACGCGCAAAAAGGCATTTACGATGCGCTGGTGGAAAAACTGGGTGCTGCGGTGGCAACGTTAAAATCTGGTGCGCCAGATGACGAGTCTACGGAGCTTGGACCTTTAAGCTCGCTGgcgc                              >  1:156809/1‑140 (MQ=255)                                                                                                  cGGCTTGTCGGATCTACGCGCAAAAAGGCATTTACGATGCGCTGGTGGAAAAACTGGGTGCTGCGGTGGCAACGTTAAAATCTGGTGCGCCAGATGACGAGTCTACGGAGCTTGGACCTTTAAGCTCGCTGgcg                               >  1:619483/1‑134 (MQ=255)                                                                                                  cGGCTTGTCGGATCTACGCGCAAAAAGGCATTTACGATGCGCTGGTGGAAAAACTGGGTGCTGCGGTGGCAACGTTAAAATCTGGTGCGCCAGATGACGAGTCTACGGAGCTTGGACCTTTAAGCTCGCTGgcg                               <  2:619483/134‑1 (MQ=255)                                                                                                  cGGCTTGTCGGATCTACGCGCAAAAAGGCATTTACGATGCGCTGGTGGAAAAACTGGGTGCTGCGGTGGCAACGTTAAAATCTGGTGCGCCAGATGACGAGTCTACGGAGCTTGGACCTTTAAGCTCGCTGGCGCAtctc                         >  1:41492/1‑140 (MQ=255)                                                                                                  cGGCTTGTCGGATCTACGCGCAAAAAGGCATTTACGATGCGCTGGTGGAAAAACTGGGTGCTGCGGTGGCAACGTTAAAATCTGGTGCGCCAGATGACGAGTCTACGGAGCTTGGACCTTTAAGCTCGCTGGCGCAtctc                         >  1:264555/1‑140 (MQ=255)                                                                                                    gCTTGTCGGATCTACGCGCAAAAAGGCATTTACGATGCGCTGGTGGAAAAACTGGGTGCTGCGGTGGCAACGTTAAAATCTGGTGCGCCAGATGACGAGTCTACGGAGCTTGGACCTTTAAGCTCGCTGGCGCATCTCGa                       >  2:69304/1‑140 (MQ=255)                                                                                                                   gcgcAAAAAGGCATTTACGATGCGCTGGTGGAAAAACTGGGTGCTGCGGTGGCAACGTTAAAATCTGGTGCGCCAGATGACGAGTCTACGGAGCTTGGACCTTTAAGCTCGCTGGCGCATCTCGAACGCGTCGGCAAGGc        <  1:77480/140‑1 (MQ=255)                                                                                                                       aaaaaGGCATTTACGATGCACTGGTGGAAAAACTGGGTGCTGCGGTGGCAACGTTAAAATCTGGTGCGCCAGATGACGAGTCTACGGAGCTTGGACCTTTAAGCTCGCTGGCGCATCTCGAACGCGTCGGCAAGGCAGTa    >  1:162351/1‑140 (MQ=255)                                                                                                                         aaaGGCATTTACGATGCGCTGGTGGAAAAACTGGGTGCTGCGGTGGCAACGTTAAAATCTGGTGCGCCAGATGACGAGTCTACGGAGCTTGGACCTTTAAGCTCGCTGGCGCATCTCGAACGCGTCGGCAAGGCAGTaga  <  2:776150/140‑1 (MQ=255)                                                                                                                                        |                                                                                                                             AAAGCGCCAGTGATTGTTTTTGATGATGCGGATATTGAAGCAGTGGTCGAAGGTGTACGTACATTTGGCTATTACAATGCTGGACAGGATTGTACTGCGGCTTGTCGGATCTACGCGCAAAAAGGCATTTACGATACGCTGGTGGAAAAACTGGGTGCTGCGGTGGCAACGTTAAAATCTGGTGCGCCAGATGACGAGTCTACGGAGCTTGGACCTTTAAGCTCGCTGGCGCATCTCGAACGCGTCGGCAAGGCAGTAGA  >  NC_000913/1516325‑1516584

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