Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,516,617	T→A	I347N (ATC→AAC)	patD →	gamma‑aminobutyraldehyde dehydrogenase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,516,617	0	T	A	100.0%	151.4 / NA	43	I347N (ATC→AAC)	patD	gamma‑aminobutyraldehyde dehydrogenase
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base A (18/25);  total (18/25)

GCTGCGGTGGCAACGTTAAAATCTGGTGCGCCAGATGACGAGTCTACGGAGCTTGGACCTTTAAGCTCGCTGGCGCATCTCGAACGCGTCGGCAAGGCAGTAGAAGAGGCGAAAGCGACAGGGCACATCAAAGTGATCACTGGCGGTGAAAAGCGCAAGGGTAATGGCTATTACTATGCGCCGACGCTGCTGGCTGGCGCATTACAGGACGATGCCATCGTGCAAAAAGAGGTATTTGGTCCAGTAGTGAGTGTTACGCCCTTCG  >  NC_000913/1516481‑1516745                                                                                                                                         |                                                                                                                                 gctgcGGTGGCAACGTTAAAATCTGGTGCGCCAGATGACGAGTCTACGGAGCTTGGACCTTTAAGCTCGCTGGCGCATCTCGAACGCGTCGGCAAGGCAGTAGAAGAGGCGAAAGCGACAGGGCACATCAAAGTGAacac                                                                                                                               <  1:157672/140‑1 (MQ=255)        tGGCAACGTTAAAATCTGGTGCGCCAGATGACGAGTCTACGGAGCTTGGACCTTTAAGCTCGCTGGCGCATCTCGAACGCGTCGGCAAGGCAGTAGAAGAGGCGAAAGCGACAGGGCACATCAAAGTGAACACTGGCGGt                                                                                                                        >  1:246384/1‑140 (MQ=255)           cAACGTTAAAATCTGGTGCGCCAGATGACGAGTCTACGGAGCTTGGACCTTTAAGCTCGCTGGCGCATCTCGAACGCGTCGGCAAGGCAGTAGAAGAGGCGAAAGCGACAGGGCACATCAAAGTGAACACTGGCGGTGaa                                                                                                                     <  2:302417/140‑1 (MQ=255)           cAACGTTAAAATCTGGTGCGCCAGATGACGAGTCTACGGAGCTTGGACCTTTAAGCTCGCTGGCGCATCTCGAACGCGTCGGCAAGGCAGTAGAAGAGGCGAAAGCGACAGGGCACATCAAAGTGAACACTGGCGGTGaa                                                                                                                     <  2:246384/140‑1 (MQ=255)                 tAAAATCTGGTGCGCCAGATGACGAGTCTACGGAGCTTGGACCTTTAAGCTCGCTGGCGCATCTCGAACGCGTCGGCAAGGCAGTAGAAGAGGCGAAAGCGACAGGGCACATCAAAGTGAACACTGGCGGTGAAAAgcgc                                                                                                               <  1:127266/140‑1 (MQ=255)                 tAAAATCTGGTGCGCCAGATGACGAGTCTACGGAGCTTGGACCTTTAAGCTCGCTGGCGCATCTCGAACGCGTCGGCAAGGCAGTAGAAGAGGCGAAAGCGACAGGGCACATCAAAGTGAACACTGGCGGTGAAAAgcgc                                                                                                               <  1:184514/140‑1 (MQ=255)                  aaaaTCTGGTGCGCCAGATGACGAGTCTACGGAGCTTGGACCTTTAAGCTCGCTGGCGCATCTCGAACGCGTCGGCAAGGCAGTAGAAGAGGCGAAAGCGACAGGGCACATCAAAGTGAACACTGGCGGTGAAAAGCGCa                                                                                                              <  1:584479/140‑1 (MQ=255)                        tGGTGCGCCAGATGACGAGTCTACGGAGCTTGGACCTTTAAGCTCGCTGGCGCATCTCGAACGCGTCGGCAAGGCAGTAGAAGAGGCGAAAGCGACAGGGCACATCAAAGTGAACACTGGCGGTGAAAAGCGCAAGGGTa                                                                                                        >  1:843798/1‑140 (MQ=255)                        tGGTGCGCCAGATGACGAGTCTACGGAGCTTGGACCTTTAAGCTCGCTGGCGCATCTCGAACGCGTCGGCAAGGCAGTAGAAGAGGCGAAAGCGACAGGGCACATCAAAGTGAACACTGGCGGTGAAAAGCGCAAGGGTa                                                                                                        <  1:442444/140‑1 (MQ=255)                         ggTGCGCCAGATGACGAGTCTACGGAGCTTGGACCTTTAAGCTCGCTGGCGCATCTCGAACGCGTCGGCAAGGCAGTAGAAGAGGCGAAAGCGACAGGGCACATCAAAGTGAACACTGGCGGTGAAAAgcgc                                                                                                               <  2:103293/132‑1 (MQ=255)                         ggTGCGCCAGATGACGAGTCTACGGAGCTTGGACCTTTAAGCTCGCTGGCGCATCTCGAACGCGTCGGCAAGGCAGTAGAAGAGGCGAAAGCGACAGGGCACATCAAAGTGAACACTGGCGGTGAAAAgcgc                                                                                                               >  1:103293/1‑132 (MQ=255)                         ggTGCGCCAGATGACGAGTCTACGGAGCTTGGACCTTTAAGCTCGCTGGCGCATCTCGAACGCGTCGGCAAGGCAGTAGAAGAGGCGAAAGCGACAGGGCACATCAAAGTGAACACTGGCGGTGAAAAGCGCAAAGGGaa                                                                                                       >  2:325400/1‑140 (MQ=255)                          gTGCGCCAGATGACGAGTCTACGGAGCTTGGACCTTTAAGCTCGCTGGCGCATCTCGAACGCGTCGGCAAGGCAGTAGAAGAGGCGAAAGCGACAGGGCACATCAAAGTGAACACTGGCGGTGAAAAGCGCAAGGGTAAt                                                                                                      <  2:866017/140‑1 (MQ=255)                          gTGCGCCAGATGACGAGTCTACGGAGCTTGGACCTTTAAGCTCGCTGGCGCATCTCGAACGCGTCGGCAAGGCAGTAGAAGAGGCGAAAGCGACAGGGCACATCAAAGTGAACACTGGCGGTGAAAAGCGCAAGGGTAAt                                                                                                      <  1:781772/140‑1 (MQ=255)                            gcgcCAGATGACGAGTCTACGGAGCTTGGACCTTTAAGCTCGCTGGCGCATCTCGAACGCGTCGGCAAGGCAGTAGAAGAGGCGAAAGCGACAGGGCACATCAAAGTGAACACTGGCGGTGAAAAGCGCAAGGGTAATgg                                                                                                    <  2:843798/140‑1 (MQ=255)                             cgcCAGATGACGAGTCTACGGAGCTTGGACCTTTAAGCTCGCTGGCGCATCTCGAACGCGTCGGCAAGGCAGTAGAAGAGGCGAAAGCGACAGGGCACATCAAAGTGAACACTGGCGGTGAAAAGCGCAAGGGTAATGGc                                                                                                   <  2:665111/140‑1 (MQ=255)                                       cGAGTCTACGGAGCTTGGACCTTTAAGCTCGCTGGCGCATCTCGAACGCGTCGGCAAGGCAGTAGAAGAGGCGAAAGCGACAGGGCACATCAAAGTGAACACTGGCGGTGAAAAGCGCAAGGGTAATGGCTATTACTATg                                                                                         >  2:197899/1‑140 (MQ=255)                                             tACGGAGCTTGGACCTTTAAGCTCGCTGGCGCATCTCGAACGCGTCGGCAAGGCAGTAGAAGAGGCGAAAGCGACAGGGCACATCAAAGTGAACACTGGCGGTGAAAAGCGCAAGGGTAATGGCTATTACTATGCGCCGa                                                                                   >  1:138576/1‑140 (MQ=255)                                                    cTTGGACCTTTAAGCTCGCTGGCGCATCTCGAACGCGTCGGCAAGGCAGTAGAAGAGGCGAAAGCGACAGGGCACATCAAAGTGAACACTGGCGGTGAAAAGCGCAAGGGTAATGGCTATTACTATGCGCCGACgctgct                                                                            >  1:208195/1‑140 (MQ=255)                                                    cTTGGACCTTTAAGCTCGCTGGCGCATCTCGAACGCGTCGGCAAGGCAGTAGAAGAGGCGAAAGCGACAGGGCACATCAAAGTGAACACTGGCGGTGAAAAGCGCAAGGGTAATGGCTATTACTATGCGCCGACgctgct                                                                            >  1:271723/1‑140 (MQ=255)                                                            tttAAGCTCGCTGGCGCATCTCGAACGCGTCGGCAAGGCAGTAGAAGAGGCGAAAGCGACAGGGCACATCAAAGTGAACACTGGCGGTGAAAAGCGCAAGGGTAATGGCTATTACTATGCGCCGACGCTGCTGGCTGgcg                                                                    <  2:271723/140‑1 (MQ=255)                                                            tttAAGCTCGCTGGCGCATCTCGAACGCGTCGGCAAGGCAGTAGAAGAGGCGAAAGCGACAGGGCACATCAAAGTGAACACTGGCGGTGAAAAGCGCAAGGGTAATGGCTATTACTATGCGCCGACGCTGCTGGCTGgcg                                                                    <  2:208195/140‑1 (MQ=255)                                                            tttAAGCTCGCTGGCGCATCTCGAACGCGTCGGCAAGGCAGTAGAAGAGGCGAAAGCGACAGGGCACATCAAAGTGAACACTGGCGGTGAAAAGCGCAAGGGTAATGGCTATTACTATGCGCAGACGCTGCTGGCTGgcg                                                                    <  1:197899/140‑1 (MQ=255)                                                              tAAGCTCGCTGGCGCATCTCGAACGCGTCGGCAAGGCAGTAGAAGAGGCGAAAGCGACAGGGCACATCAAAGTGAACACTGGCGGTGAAAAGCGCAAGGGTAATGGCTATTACTATGCGCCGACGCTGCTGGCTGgcg                                                                    >  2:768654/1‑138 (MQ=255)                                                              tAAGCTCGCTGGCGCATCTCGAACGCGTCGGCAAGGCAGTAGAAGAGGCGAAAGCGACAGGGCACATCAAAGTGAACACTGGCGGTGAAAAGCGCAAGGGTAATGGCTATTACTATGCGCCGACGCTGCTGGCTGgcg                                                                    <  1:768654/138‑1 (MQ=255)                                                              tAAGCTCGCTGGCGCATCTCGAACGCGTCGGCAAGGCAGTAGAAGAGGCGAAAGCGACAGGGCACATCAAAGTGAACACTGGCGGTGAAAAGCGCAAGGGTAATGGCTATTACTATGCGCCGACGCTGCTGGCTGGCGCa                                                                  >  1:18289/1‑140 (MQ=255)                                                              tAAGCTCGCTGGCGCATCTCGAACGCGTCGGCAAGGCAGTAGAAGAGGCGAAAGCGACAGGGCACATCAAAGTGAACACTGGCGGTGAAAAGCGCAAGGGTAATGGCTATTACTATGCGCCGACGCTGCTGGCTGGCGCa                                                                  >  2:61976/1‑140 (MQ=255)                                                              tAAGCTCGCTGGCGCATCTCGAACGCGTCGGCAAGGCAGTAGAAGAGGCGAAAGCGACAGGGCACATCAAAGTGAACACTGGCGGTGAAAAGCGCAAGGGGAATGGCTATTACTATGCGCCGACGCCGCTGGCTGgggcc                                                                  >  2:712716/1‑139 (MQ=255)                                                                 gCTCGCTGGCGCATCTCGAACGCGTCGGCAAGGCAGTAGAAGAGGCGAAAGCGACAGGGCACATCAAAGTGAACACTGGCGGTGAAAAGCGCAAGGGTAATGGCTATTACTATGCGCCGACGCTGCTGGCTGGCGCATTa                                                               <  1:867924/140‑1 (MQ=255)                                                                       tGGCGCATCTCGAACGCGTCGGCAAGGCAGTAGAAGAGGCGAAAGCGACAGGGCACATCAAAGTGAACACTGGCGGTGAAAAGCGCAAGGGTAATGGCTATTACTATGCGCCGACGCTGCTGGCTGGCGCATTACAGGAc                                                         >  2:816026/1‑140 (MQ=255)                                                                           gcATCTCGAACGCGTCGGCAAGGCAGTAGAAGAGGCGAAAGCGACAGGGCACATCAAAGTGAACACTGGCGGTGAAAAGCGCAAGGGTAATGGCTATTACTATGCGCCGACGCTGCTGGCTGGCGCATTACAGGACGATg                                                     <  1:816794/140‑1 (MQ=255)                                                                               ctcGAACGCGTCGGCAAGGCAGTAGAAGAGGCGAAAGCGACAGGGCACATCAAAGTGAACACTGGCGGTGAAAAGCGCAAGGGTAATGGCTATTACTATGCGCCGACGCTGCTGGCTGGCGCATTACAGGACGATGCCAt                                                 <  2:18289/140‑1 (MQ=255)                                                                                     cgcgTCGGCAAGGCAGTAGAAGAGGCGAAAGCGACAGGGCACATCAAAGTGAACACTGGCGGTGAAAAGCGCAAGGGTAATGGCTATTACTATGCGCCGACGCTGCTGGCTGGCGCATTACAGGACGATGCCATCGTGCa                                           <  2:138576/140‑1 (MQ=255)                                                                                       cgTCGGCAAGGCAGTAGAAGAGGCGAAAGCGACAGGGCACATCAAAGTGAACACTGGCGGTGAAAAGCGCAAGGGTAATGGCTATTACTATGCGCCGACGCTGCTGGCTGGCGCATTACAGGACGATGCCATCGTGCaaa                                         <  1:712716/140‑1 (MQ=255)                                                                                         tCGGCAAGGCAGTAGAAGAGGCGAAAGCGACAGGGCACATCAAAGTGAACACTGGCGGTGAAAAGCGCAAGGGTAATGGCTATTACTATGCGCCGACGCTGCTGGCTGGCGCATTACAGGACGATGCCATCGTGCaaaaa                                       >  2:641270/1‑140 (MQ=255)                                                                                          cGGCAAGGCAGTAGAAGAGGCGAAAGCGACAGGGCACATCAAAGTGAACACTGGCGGTGAAAAGCGCAAGGGTAATGGCTATTACTATGCGCCGACGCTGCTGGCTGGCGCATTACAGGACGATGCCATCGTGCa                                           <  1:684503/135‑1 (MQ=255)                                                                                          cGGCAAGGCAGTAGAAGAGGCGAAAGCGACAGGGCACATCAAAGTGAACACTGGCGGTGAAAAGCGCAAGGGTAATGGCTATTACTATGCGCCGACGCTGCTGGCTGGCGCATTACAGGACGATGCCATCGTGCa                                           >  2:684503/1‑135 (MQ=255)                                                                                                ggCAGTAGAAGAGGCGAAAGCGACAGGGCACATCAAAGTGAACACTGGCGGTGAAAAGCGCAAGGGTAATGGCTATTACTATGCGCCGACGCTGCTGGCTGGCGCATTACAGGACCATGCCATCGTGCAAAAAGAGGTAt                                >  2:656251/1‑140 (MQ=255)                                                                                                     tAGAAGAGGCGAAAGCGACAGGGCACATCAAAGTGAACACTGGCGGTGAAAAGCGCAAGGGTAATGGCTATTACTATGCGCCGACGCTGCTGGCTGGCGCATTACAGGACGATGCCATCGTGCAAAAAGAGGTATTTGGt                           <  1:641270/140‑1 (MQ=255)                                                                                                        aagaGGCGAAAGCGACAGGGCACATCAAAGTGAACACTGGCGGTGAAAAGCGCAAGGGTAATGGCTATTACTATGCGCCGACGCTGCTGGCTGGCGCATTACAGGACGATGCCATCGTGCAAAAAGAGGTATTTGGTCCa                        <  2:415177/140‑1 (MQ=255)                                                                                                            ggCGAAAGCGACAGGGCACATCAAAGTGAACACTGGCGGTGAAAAGCGCAAGGGTAATGGCTATTACTATGCGCCGACGCTGCTGGCTGGCGCATTACAGGACGATGCCATCGTGCAAAAAGAGGGATTTGGTCCagtag                    >  2:307578/1‑140 (MQ=255)                                                                                                             gCGAAAGCGACAGGGCACATCAAAGTGAACACTGGCGGTGAAAAGCGCAAGGGTAATGGCTATTACTATGCGCCGACGCTGCTGGCTGGCGCATTACAGGACGATGCCATCGTGCAAAAAGAGGTATTTGGTCCagtagt                   <  1:61976/140‑1 (MQ=255)                                                                                                              cGAAAGCGACAGGGCACATCAAAGTGAACACTGGCGGTGAAAAGCGCAAGGGTAATGGCTATTACTATGCGCCGACGCTGCTGGCTGGCGCATTACAGGACGATGCCATCGTGCAAAAAGAGGTATTTGGTCCAGTagtg                  <  1:816026/140‑1 (MQ=255)                                                                                                                              caTCAAAGTGAACACTGGCGGTGAAAAGCGCAAGGGGAATGGCTATTACTATGCGCCGACGCTGCTGGCTGGCGCATTACAGGACGATGCCATCGTGCAAAAAGAGGTATTTGGTCCAGTAGTGAGTGTTACGCCCTTCg  >  2:831123/1‑140 (MQ=255)                                                                                                                                         |                                                                                                                                 GCTGCGGTGGCAACGTTAAAATCTGGTGCGCCAGATGACGAGTCTACGGAGCTTGGACCTTTAAGCTCGCTGGCGCATCTCGAACGCGTCGGCAAGGCAGTAGAAGAGGCGAAAGCGACAGGGCACATCAAAGTGATCACTGGCGGTGAAAAGCGCAAGGGTAATGGCTATTACTATGCGCCGACGCTGCTGGCTGGCGCATTACAGGACGATGCCATCGTGCAAAAAGAGGTATTTGGTCCAGTAGTGAGTGTTACGCCCTTCG  >  NC_000913/1516481‑1516745

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