Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913_3_pae_tpiA	2,795,396	C→A	Q382K (CAG→AAG)	gabP →	gamma‑aminobutyrate transporter

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913_3_pae_tpiA	2,795,396	0	C	A	100.0%	160.4 / NA	48	Q382K (CAG→AAG)	gabP	gamma‑aminobutyrate transporter
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base A (26/22);  total (26/22)

CGTGGCGGTGTTACTCTCCACCGGCGCGGCATTCTTAACGGTGGTGGTGAACTATTACGCGCCTGCGAAGGTATTTAAATTTCTGATCGACAGCTCCGGCGCTATCGCCCTGCTGGTTTATTTAGTCATCGCCGTTTCACAGTTGCGGATGCGCAAAATTCTGCGAGCAGAAGGAAGCGAAATTCGCTTGCGAATGTGGCTTTATCCGTGGCT‑CACCTGGCTAGTCATCGGCTTTATTACCTTTGTGTTGGTAGTGATGCTATTCCGTCCGGCGCAACA  >  NC_000913_3_pae_tpiA/2795257‑2795535                                                                                                                                            |                                                                                                                                             cGTGGCGGTGTTACTCTCCACCGGCGCGGCATTCTTAACGGTGGTGGTGAACTATTACGCGCCTGCGAAGGTATTTAAATTTCTGATCGACAGCTCCGGCGCTATCGCCCTGCTGGTTTATTTAGTCATCGCCGTTTCAAAGTTGCGGa                                                                                                                                     >  1:70695/1‑149 (MQ=255)  gTGGCGGTGTTACTCTCCACCGGCGCGGCATTCTTAACGGTGGTGGTGAACTATTACGCGCCTGCGAAGGTATTTAAATTTCTGATCGACAGCTCCGGCGCTATCGCCCTGCTGGTTTATTTAGTCATCGCCGTTTCAAAGTTGCGGAt                                                                                                                                    <  1:50504/149‑1 (MQ=255)    ggCGGTTTTACTTTCGTCCGTCGCGGCATTCTTAACGGTGGTGGTGAACTATTACGCGCCTGCGAAGGTATTTAAATTTCTGATCGACAGCTCCGGCGCTATCGCCCTGCTGGTTTATTTAGTCATCGCCGTTTCAAAGTTGCGGATgc                                                                                                                                  <  1:71266/149‑1 (MQ=255)     gCGGTGTTACTCTCCACCGGCGAGGCATTCATAACGGGGAGCGTGAACTATTAAGCGCCTGCGAAGGTATTTAAATTTCTGATGGACAGCTCCGGCCCTATCGCCGTGCTGGTTTATTTAGTACTGGCGGTGTCAAAGTTGCGACTgcg                                                                                                                                 >  2:144463/1‑149 (MQ=255)       ggTGTTACTCTCCACCGGCGCGGCATTCTTAACGGTGGTGGTGAACTATTACGCGCCTGCGAAGGTATTTAAATTTCTGATCGACAGCTCCGGCGCTATCGCCCTGCTGGTTTATTTAGTCATCGCCGTTTCAAAGTTGCGGATGCGCa                                                                                                                               <  1:188078/149‑1 (MQ=255)               tctcCACCGGCGCGGCATTCTTAACGGTGGTGGTGAACTATTACGCGCCTGCGAAGGTATTTAAATTTCTGATCGACAGCTCCGGCGCTATCGCCCTGCTGGTTTATTTAGTCATCGCCGTTTCAAAGTTGCGGATGCGCAAAATTCTg                                                                                                                       >  1:747497/1‑149 (MQ=255)               tctcCACCGGCGCGGCATTCTTAACGGTGGTGGTGAACTATTACGCGCCTGCGAAGGTATTTAAATTTCTGATCGACAGCTCCGGCGCTATCGCCCTGCTGGTTTATTTAGTCATCGCCGTTTCAAAGTTGCGGATGCGCAAAATTCTg                                                                                                                       >  1:636888/1‑149 (MQ=255)                    aCCGGCGCGGCATTCTTAACGGTGGTGGTGAACTATTACGCGCCTGCGAAGGTATTTAAATTTCTGATCGACAGCTCCGGCGCTATCGCCCTGCTGGTTTATTTAGTCATCGCCGTTTCAAAGTTGCGGATGCGCAAAATTCTGCGAGc                                                                                                                  >  1:609324/1‑149 (MQ=255)                     ccGGCGCGGCATTCTTAACGGTGGTGGTGAACTATTACGCGCCTGCGAAGGTATTTAAATTTCTGATCGACAGCTCCGGCGCTATCGCCCTGCTGGTTTATTTAGTCATCGCCGTTTCAAAGTTGCGGATGCGCAAAATTCTGCGAGCa                                                                                                                 <  1:468082/149‑1 (MQ=255)                             gCATTCTTAACGGTGGTGGTGAACTATTACGCGCCTGCGAAGGTATTTAAATTTCTGATCGACAGCTCCGGCGCTATCGCCCTGCTGGTTTATTTAGTCATCGCCGTTTCAAAGTTGCGGATGCGCAAAATTCTGCGAGCAgaaggaag                                                                                                         >  1:821913/1‑149 (MQ=255)                                   ttAACGGTGGTGGTGAACTATTACGCGCCTGCGAAGGTATTTAAATTTCTGATCGACAGCTCCGGCGCTATCGCCCTGCTGGTTTATTTAGTCATCGCCGTTTCAAAGTTGCGGATGCGCAAAATTCTGCGAGCAGAAGGAAGCGAAAt                                                                                                   >  2:61434/1‑149 (MQ=255)                                   ttAACGGTGGCGGTGAACTATTACGCGCCTGCGAAGGTATTTAAATTTCTGATCGACAGCTCCGGCGCTATCGCCCTGCTGGTTTATTTAGTCATCGCCGTTTCAAAGTTGCGGATGCGAAAAATTCTGCGAGCACAAGGAAGCGAAAt                                                                                                   >  2:318374/1‑149 (MQ=255)                                     aaCGGTGGTGGTGAACTATTACGCGCCTGCGAAGGTATTTAAATTTCTGATCGACAGCTCCGGCGCTATCGCCCTGCTGGTTTATTTAGTCCTAGGAGTCGCAAAGTTGCGGATGCGCAAAATTCTGCGAGCAGAAGGAAGCGAAATTc                                                                                                 <  1:803901/149‑1 (MQ=255)                                                tgAACTATTACGCGCCTGCGAAGGTATTTAAATTTCTGATCGACAGCTCCGGCGCTATCGCCCTGCTGGTTTATTTAGTCATCGCCGTTTCAAAGTTGCGGATGCGCAAAATTCTGCGAGCAGAAGGAAGCGAAATTCGCTTGCGAAtg                                                                                      >  2:493817/1‑149 (MQ=255)                                                  aaCTATTACGCGCCTGCGAAGGTATTTAAATTTCTGATCGACAGCTCCGGCGCTATCGCCCTGCTGGTTTATTTAGTCATCGCCGTTTCAAAGTTGCGGATGCGCAAAATTCTGCGAGCAGAAGGAAGCGAAATTCGCTTGCGAAtgtg                                                                                    >  1:267032/1‑149 (MQ=255)                                                   aCTATTACGGGCCTGCGAAGGTATTTAAATTTCTGATCGACAGCTCCGGCGCTATCGCCCTGCTGGTTTATTTAGTCATCGCCGTTTCAAAGTTGCGGATGCGCAAAATTCTGCGAGCAGAAGGAAGCGAAATTCGCTTGCGAATTTgg                                                                                   <  2:528214/149‑1 (MQ=255)                                                      aTTACGCGCCTGCGAAGGTATTTAAATTTCTGATCGACAGCTCCGGCGCTATCGCCCTGCTGGTTTATTTAGTCATCGCCGTTTCAAAGTTGCGGATGCGCAAAATTCTGCGAGCAGAAGGAAGCGAAATTCGCTTGCGAATGTGGCtt                                                                                >  1:400459/1‑149 (MQ=255)                                                                   gAAGGTATTTAAATTTCTGATCGACAGCTCCGGCGCTATCGCCCTGCTGGTTTATTTAGTCATCGCCGTTTCAAAGTTGCGGATGCGCAAAATTCTGCGAGCAGAAGGAAGCGAAATTCGCTTGCGAATGTGGCTTTATCCGTGGCtaa                                                                   >  1:403431/1‑147 (MQ=255)                                                                      ggTATTTAAATTTCTGATCGACAGCTCCGGCGCTATCGCCCTGCTGGTTTATTTAGTCATCGCCGTTTCAAAGTTGCGGATGCGCAAAATTCTGCGAGCAGAAGGAAGCGAAATTCGCTTGCGAATGTGGCTTTATCCGTGGCT‑CACCt                                                               <  2:728688/149‑1 (MQ=255)                                                                      ggTATTTAAATTTCTGATCGACAGCTCCGGCGCTATCGCCCTGCTGGTTTATTTAGTCATCGCCGTTTCAAAGTTGCGGATGCGCAAAATTCTGCGAGCAGAAGGAAGCGAAATTCGCTTGCGAATGTGGCTTTATCCGTGGCT‑CACCt                                                               <  1:637062/149‑1 (MQ=255)                                                                         aTTTAAAGTTCTGATCGACAGCTCCGGCGCTATCGCCCTGCTGGTTTATTTAGTCATCGCCGTTTCAAAGTTGCGGATGCGCAAAATTCTGCGAGCAGAAGGAAGCGAAATTCGCTTGCGAATGTGGCTTTATCCGTGGCT‑CACCTGGc                                                            <  1:720856/149‑1 (MQ=255)                                                                              aaTTTCTGATCGACAGCTCCGGCGCTATCGCCCTGCTGGTTTATTTAGTCATCGCCGTTTCAAAGTTGCGGATGCGCAAAATTCTGCGAGCAGAAGGAAGCGAAATTCGCTTGCGAATGTGGCTTTATCCGTGGCT‑CACCTGGCTAGTc                                                       >  1:808261/1‑149 (MQ=255)                                                                                  tgcgATCGACAGCTCCGGCGCTCTCCCCCTGCTGGTTTATTTAGTCATCGCCGTTTCAAAGTTGCGGATGCGCAAAATTCTGCGAGCAGAAGGAAGCGAAATTCGCTTGCGAATGTGGCTTTATCCGTGGCT‑CACCGGGCTAGTCATCg                                                   <  1:324313/146‑1 (MQ=255)                                                                                  tCTGATCGACAGCTCCGGCGCTATCGCCCTGCTGGTTTATTTAGTCATCGCCGTTTCAAAGTTGCGGATGCGCAAAATTCTGCGAGCAGAAGGAAGCGAAATTCGCTTGCGAATGTGGCTTTATCCGTGGCT‑CACCTGGCTAGTCATCg                                                   <  2:614949/149‑1 (MQ=255)                                                                                    tGATCGACAGCTCCGGCGCTATCGCCCTGCTGGTTTATTTAGTCATCGCCGTTTCAAAGTTGCGGATGCGCAAAATTCTGCGAGCAGAAGGAAGCGAAATTCGCTTGCGAATGTGGCTTTATCCGTGGCT‑CACCTGGCTAGTCCTCGGc                                                 >  1:33959/1‑149 (MQ=255)                                                                                    tGATCGACAGCTCCGGCGCTATCGCCCTGCTGGTTTATTTAGTCATCGCCGTTTCAAAGTTGCGGATGCGCAAAATTCTGCGAGCAGAAGGAAGCGAAATTCGCTTGCGAATGTGGCTTTATCCGTGGCT‑CACCTGGCTAGTCATCGGc                                                 >  1:560301/1‑149 (MQ=255)                                                                                      aTCGACAGCTCCGGCGCTATCGCCCTGCTGGTTTATTTAGTCATCGCCGTTTCAAAGTTGCGGATGCGCAAAATTCTGCGAGCAGAAGGAAGCGAAATTCGCTTGCGAATGTGGCTTTATCCGTGGCT‑CACCTGGCTAGTCATCGGAtt                                               >  1:349823/1‑149 (MQ=255)                                                                                       tCGACAGCTCCGGCGCTATCGCCCTGCTGGTTTATTTAGTCATCGCCGTTTCAAAGTTGCGGATGCGCAAAATTCTGCGAGCAGAAGGAAGCGAAATTCGCTTGCGAATGTGGCTTTATCCGTGGCT‑CACCTGGCTAGTCATCGGCttt                                              >  1:485469/1‑149 (MQ=255)                                                                                       tCGACAGCTCCGGCGCTATCGCCCTGCTGGTTTATTTAGTCATCGCCGTTTCAAAGTTACGGATGCGCAAAATTATGCGAGCAGAAGCAAGCGAAATTCGCTTGCGAATGTGGCTTTGTCCGTGGCC‑CACCTGGCTAGTCTTCGGCttt                                              >  2:380074/1‑149 (MQ=255)                                                                                          aCAGCTCCGGCGCTATCGCCCTGCTGGTTTATTTAGTCATCGCCGTTTCAAAGTTGCGGATGCGCAAAATTCTGCGAGCAGAAGGAAGCGAAATTCGCTTGCGAATGTGGCTTTATCCGTGGCT‑CACCTGGCTAGTCATCGGCTttatt                                           >  2:90450/1‑149 (MQ=255)                                                                                          aCAGCTCCGGCGCGATCGCCCTGCTGGCTTATTTAGTCATCGCCGTTTCAAAGTTGCGGATGCGCAAAATTCTGCGAGCAGAAGGACGCGAAATGCGCTTGCGAATGTGGCTCTATCCGTGGCT‑CACCTGGCTAGTCGTCGGCTttatt                                           >  2:472898/1‑149 (MQ=255)                                                                                           cAGCTCCGGCGCTATCGCCCTGCTGGTTTATTTAGTCATCGCCGTTTCAAAGTTGCGGATGCGCAAAATTCTGCGAGCAGAAGGAAGCGAAATTCGCTTGCGAATGTGGCTTTATCCGTGGCT‑CACCTGGCTAGTCATCGGCTttatta                                          >  1:397542/1‑149 (MQ=255)                                                                                             ggcccGGCGCTATCGCCCTGCTGGTTTATTTAGTCATCGCCGTTTCAAAGTTGCGGATGCGCAAAATTCTGCGAGCAGAAGGAAGCGAAATTCGCTTGCGAATGTGGCTTTA‑CAGTTAATGCACCTGGCTAGTCATCGGCTTTATTAcc                                        <  2:822550/146‑1 (MQ=255)                                                                                               tCCGGCGCTATCGCCCTGCTGGTTTATTTAGTCATCGCCGTTTCAAAGTTGCGGATGCGCAAAATTCTGCGAGCAGAAGGAAGc                                                                                                        >  1:22361/1‑84 (MQ=255)                                                                                               tCCGGCGCTATCGCCCTGCTGGTTTATTTAGTCATCGCCGTTTCAAAGTTGCGGATGCGCAAAATTCTGCGAGCAGAAGGAAGc                                                                                                        <  2:22361/84‑1 (MQ=255)                                                                                                    cgcTATCGCCCTGCTGGTTTATTTATTCATCGCCGTTTCAAAGTTGCGGATGCGCAAAATTCTGCGAGCAGAAGGAAGCGAAATTCGCTTGCGAATGTGGCTTTATCCGTGGCT‑CACCTGGCTAGTCATCGGCTTTATTACCTTtgttt                                 <  2:529934/149‑3 (MQ=255)                                                                                                      cTATCGCCCTGCTGGTTTATTTAGTCATCGCCGTTTCAAAGTTGCGGATGCGCAAAATTCTGCGAGCAGAAGGAAGCGAAATTCGCTTGCGAATGTGGCTTTATCCGTGGCT‑CACCTGGCTAGTCATCGGCTTTATTACCTTTGTGTTg                               <  2:560301/149‑1 (MQ=255)                                                                                                         tCGCCCTGCTGGTTTATTTAGTCATCGCCGTTTCAAAGTTGCGGATGCGCAAAATTCTGCGAGCAGAAGGAAGCGAAATTCGCTTGCGAATGTGGCTATATCCGTGGCT‑CACCTGGCTAGTCATCGGCTTTATTACCTTTGTGTTGGTa                            <  1:242292/149‑1 (MQ=255)                                                                                                             ccTGCTGGTTTATTTAGTCATCGCCGTTTCAAAGTTGCGGATGCGCAAAATTCTGCGAGCAGAAGGAAGCGAAATTCGCTTGCGAATGTGGCTTTATCCGTGGCT‑CACCTGGCTAGTCATCGGCTTTATTACCTTTGTGTTGGTAGTGa                        <  2:267032/149‑1 (MQ=255)                                                                                                               tgctgGTTTAGTTGGTCATCGCCGTTTCAAAGTTGCGGATGCGCAAAATTCTGCGAGCAGAAGGAAGCGAAATTCGCTTGCGAATGTGGCTTTATCCGTGGCT‑CACCTGGCTAGTCATCGGCTTTATTACCTTTGTGTTGGTAGTGATg                      <  1:29969/149‑1 (MQ=255)                                                                                                               tgcagttttgtttaGTGATCGCGGTTTCAAAGTTGCGGATGCGCAAAATTCTGCGAGCAGAAGGAAGCGAAATTCGCTTGCGAATGTGGCTTTATCCGTGGCT‑CACCTGGCTAGTCATCGGCTTTATTACCTTTGTGTTGGTAGTGATg                      <  2:650342/143‑1 (MQ=255)                                                                                                                 ctgGTTTATTTAGTCATCGCCGTTTCAAAGTTGCGGATGCGCAAAATTCTGCGAGCAGAAGGAAGCGAAATTCGCTTGCGAATGTGGCTTTATCCGTGGCT‑CACCTGGCTAGTCATCGGCTTTATTACCTTTGTGTTGGTAGTGATGCt                    >  1:729606/1‑149 (MQ=255)                                                                                                                  tgGTTTATTTAGTCATCGCCGTTTCAAAGTTGCGGATGCGCAAAATTCTGCGAGCAGAAGGAAGCGAAATTCGCTTGCGAATGTGGCTTTATCCGTGGCT‑CACCTGGCTAGTCATCGGCTTTATTACCTTTGTGTTGGTAGTGATGCTa                   >  1:530361/1‑149 (MQ=255)                                                                                                                     tttatttaGTCATCGCCGTTTCAAAGTTGCGGATGCGCAAAATTCTGCGAGCAGAAGGAAGCGAAATTCGCTTGCGAATGTGGCTTTATCCGTGGCT‑CACCTGGCTAGTCATCGGCTTTATTACCTTTGTGTTGGTAGTGATGCTATtc                >  1:34587/1‑149 (MQ=255)                                                                                                                       tatttaGTCATCGCCGTTTCAAAGTTGCGGATGCGCAAAATTCTGCGAGCAGAAGGAAGCGAAATTCGCTTGCGAATGTGGCTTTATCCGTGGCT‑CACCTGGCTAGTCATCGGCTTTATTACCTTTGTGTTGGTAGTGATTCTTTtccg              >  2:159531/1‑149 (MQ=255)                                                                                                                                 tggctgTTTCAAATTTGGCGATGCGCAAAATTCTTACAGCAGAAGGAAGCGAAATGCGCTTGCGAATGTGGCTTTATCCGTGCCT‑CACCTGGCTAGGCATCGGCGTTATTGCCTTTGTGTTGGTAGTGATGCTATTCCGTCCGGCGCaa    <  1:343005/144‑1 (MQ=255)                                                                                                                                 tCGCCGTTTCAAAGTTGCGGATGCGCAAAATTCTGCGAGCAGAAGGAAGCGAAATTCGCTTGCGAATGTGGCTTTATCCGTGGCT‑CACCTGGCTAGTCATCGGCTTTATTACCTTTGTGTTGGTAGTGATGCTATTCCGTCCGGCGCaa    <  2:808261/149‑1 (MQ=255)                                                                                                                                   gCCGTTTCAAAGTTGCGGATGCGCAAACTTCTGCGCGCAGGAAGAAGCGAAATTCGCTTGCGAATGTGGCTTTATCCGTGGCT‑CACCTGGCTAGTCATCGGCTTTATTACCTTTGTGTTGGTAGTGATGCTATTCCGTCCGGCGCAACa  <  1:232983/149‑1 (MQ=255)                                                                                                                                            |                                                                                                                                             CGTGGCGGTGTTACTCTCCACCGGCGCGGCATTCTTAACGGTGGTGGTGAACTATTACGCGCCTGCGAAGGTATTTAAATTTCTGATCGACAGCTCCGGCGCTATCGCCCTGCTGGTTTATTTAGTCATCGCCGTTTCACAGTTGCGGATGCGCAAAATTCTGCGAGCAGAAGGAAGCGAAATTCGCTTGCGAATGTGGCTTTATCCGTGGCT‑CACCTGGCTAGTCATCGGCTTTATTACCTTTGTGTTGGTAGTGATGCTATTCCGTCCGGCGCAACA  >  NC_000913_3_pae_tpiA/2795257‑2795535

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 13 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑