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%	178.0 / NA	52	Q382K (CAG→AAG)	gabP	gamma‑aminobutyrate transporter
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base A (23/29);  total (23/29)

ACTCTCCACCGGCGCGGCATTCTTAACGGTGGTGGTGAACTATTACGCGCCTGCGAAGGTATTTAAATTTCTGATCGACAGCTCCGGCGCTATCGCCCTGCTGGTTTATTTAGTCATCGCCGTTTCACAGTTGCGGATGCGCAAAATTCTGCGAGCAGAAGGAAGCGAAATTCGCTTGCGAATGTGGCTTTATCCGTGGCTCACCTGGCTAGTCATCGGCTTTATTACCTTTGTGTTGGTAGTGATGCTATTCCGTCCGGCGCAACAGTTA  >  NC_000913_3_pae_tpiA/2795269‑2795539                                                                                                                                |                                                                                                                                                aCTCTCCACCGGCGCGGCATTCTTAACGGTGGTGGTAAACTATTACGCGCCTGCGAAGGTATTTAAATTTCTGATCGACAGCTCCGGCGCTATCGCCCTGCTGGTTTATTTAGTACTCGCCGTTTCAACGTTGCGGATGCGCAAAATTc                                                                                                                            >  1:107960/1‑149 (MQ=255)   tctcCACCGGCGCGGCATTCTTAACGGTGGTGGTGAACTATTACGCGCCTGCGAAGGTATTTAAATTTCTGATCGACAGCTCCGGCGCTATCGCCCTGCTGGTTTATTTAGTCATCGCCGTTTCAAAGTTGCGGATGCGCAAAATTCTg                                                                                                                          <  1:169829/149‑1 (MQ=255)          cGGCGCGGCATTTTTAACGGTGGTGGTGAACTACTACGCGCCGGCGAAGTTATTTAAATTTCTGATCGACAGCTCCGGGGCTATCGCCCTGCTGGTTTATTTAGTCATCGCCGTTTCAAAGTCGCGGATGCGCTTAATTCTGCGAGCAg                                                                                                                   <  2:989970/149‑1 (MQ=255)             cgcgGGATTCTTAACGGTGGTGGTGAACTATTACGCGCCCGCGAAGGTATTTAAATTTCTGATCGACAGCTCCGGCGCTACCGCGCTGCTGGTTTATTTAGTCATCGCCGTTTCAAAGTTGCGGATGCGCAAAATTCTGCGAGCAgaag                                                                                                                <  1:954296/149‑1 (MQ=255)                     tCTTAACGGTGGTGGTGAACTATTCCGCGCCTGCGAAGGTATTTAAATTTCTGATCGACAGCTCCGGCGCTATCGCCCTGCTGGTTTATTTAGTCATCGCCGTTTCAAAGTTGCGGATGCGCAAAATTCTGCGAGCAGAAGGAAGCGaa                                                                                                        <  1:617145/149‑1 (MQ=255)                      tttAACGGTGGTGGTGAACCATTACGCGCCTGCGAAGGTATTTAAATTTCTGATCGACAGCTCCGGCGCTATCGCCCTGCTGGTTTATTTAGTCATCGCCGTTTCAAAGTTGCGGATGCGCAAAATTCTGCGAGCAGAAGGAAGCGaaa                                                                                                       <  2:725052/148‑1 (MQ=255)                      cTTAACGGTGGTGGTGAACTATTACGCGCCTGCGAAGGTATTTAAATTTCTGATCGACAGCTCCGGCGCTATCGCCCTGCTGGTTTATTTAGTCATCGCCGTTTCAAAGTTGCGGATGCGCAAAATTCTGCGAGCAGAAGGAAGCGaaa                                                                                                       >  1:880345/1‑149 (MQ=255)                      cTTAACGGTGGTGGTGAACTATTACGCGCCTGCGAAGGTATTTAAATTTCTGATCGACAGCTCCGGCGCTATCGCCCTGCTGGTTTATTTAGTCATCGCCGTTTCAAAGTTGCGGATGCGCAAAATTCTGCGAGCAGAAGGAAGCGaaa                                                                                                       >  1:425001/1‑149 (MQ=255)                         aaCGGTGGTGGTGAACTATTACGCGCCTGCGAAGGTATTTAAATTTCTGATCGACAGCTCCGGCGCTATCGCCCTGCTGGTTTATTTAGTCATCGCCGTTTCAAAGTTGCGGATGCGCAAAATTCTGCGAGCAGAAGGAAGCGAAATTc                                                                                                    <  2:88768/149‑1 (MQ=255)                           cGGTGGTGGTGAACTATTACGCGCCTGCGAAGGTATTTAAATTTCTGATCGACAGCTCCGGCGCTATCGCCCTGCTGGTTTATTTAGTCATCGCCGTTTCAAAGTTGCGGATGCGCAAAATTCTGCGAGCAGAAGGAAGCGAAATTCGc                                                                                                  <  2:496583/149‑1 (MQ=255)                            ggtggtggtGAACTATTACGCGCCTGCGAAGGTATTTAAATTTCTGATCGACAGCTCCGGCGCTATCGCCCTGCTGGTTTATTTAGTCCTCGCCGTTTCAAAGCTGCGGATGCGCAAATTTTTGCGAGCAGAAGGAAGCGAAATTCGCt                                                                                                 <  2:674091/149‑1 (MQ=255)                            ggtgggtgtGGACTATTACGCGCCGGCTAAGGTATTTAAATTTCTGATCGACAGCTCCGGCGCTATCGCCCTGCTGGTTTATTTAGTCATCGCCGTTTCAAAGTTGCGGATGCGCAAAATTCTGCGAGCAGAAGGAAGCGAAATTCGCt                                                                                                 <  2:404956/149‑1 (MQ=255)                             gtgttggggAACTTTTACGGCGCTCCGAAGTCATTTAAAGTTCTGATCGACAGCTCCGGCTCCATGTCCCTGCTGTTTTACTTAGTCAGCTCCGTTTCAAAGTTGCGCATGCGCAAAATTCTGCGAGCAGAAGGAAGCGAAATTCGCtc                                                                                                <  2:668003/149‑2 (MQ=255)                                gtggtgAACTATTACGCGCCTGCGAAGGTATTTAAATTTCTGATCGACAGCTAAGGCGAACTCCCCCAGCTGGTTTATTTAGACAATCCCGTCTTAAAGCCCCGCATGCGCAAACATCAGCGAACAGAAGGAAGCGAAAATCGCGTGCg                                                                                             >  2:99867/1‑149 (MQ=255)                                 tggtgAACTATTACGCGCCTGCGAAGGTATTTAAATTTCTGATCGACAGCTCCGGCGCTATCGCCCTGCTGGTTTATTTAGTCATCGCCGTTTCAAAGTTGCGGATGCGCAAAATTCTGCGAGCAGAAGGAAGCGAAATTCGCTTGCGa                                                                                            >  1:121026/1‑149 (MQ=255)                                     gAACTATTTCGCGCCTGCAAAGGTATTTAAATTTCTGATCGACAGCTCCGGCGCTATCGCCCTGCTGGTTTATTTAGTCATCGCCGTTTCAAAGTTGCGGATGCGCAAAATTCTGCGAGCAGAAGGAAGCGAAATTCGCTTGCGAAtgt                                                                                        <  2:419676/149‑1 (MQ=255)                                      aaCTATTACGCGCCTGCGAAGGTATTTAAATTTCTGATCGACAGCTCCGGCGCTATCGCCCTGCTGGTTTATTTAGTCATCGCCGTTTCAAAGTTGCGGATGCGCAAAATTCTGCGAGCAGAAGGAAGCGAAATTCGCTTGCGAAtgtg                                                                                       >  2:43121/1‑149 (MQ=255)                                      aaCTATTACGCGCCTGCGAAGGTATTTAAATTTCTGATCGACAGCTCCGGCGCTATCGCCCTGCTGGTTTATTTAGTCATCGCCGTTTCAAAGTTGCGGATGCGCAAAATTCTGCGAGCAGAAGGAAGCGAAATTCGCAAGCGAAtttc                                                                                       <  2:880345/149‑4 (MQ=255)                                       aCTAGTACGCGCCTGCGAAGGTATTTAAATTTCTGATCGACAGCTCCGGCGCTATCGCCCTGCTGGTTTATTTAGTCATCGCCGTTTCAAAGTTGCGGATGCGCAAAATTCTGCGAGCAGAAGGAAGCGAAATTCGCTTGCGAATGTgg                                                                                      <  1:296238/149‑1 (MQ=255)                                            tACGCGCCTGCGAAGGTATTTAAATTTCTGATCGACAGCTCCGGCGCTATCGCCCTGCTGGTTTATTTAGTCATCGCCGTTTCAAAGTTGCGGATGCGCAAAATTCTGCGAGCAGAAGGAAGCGGAATTCGCTTGCGAATGTGGCTTTa                                                                                 >  1:365446/1‑149 (MQ=255)                                               gcgcCTGCGAAGGTATTTAAATTTCTGATCGACAGCTCCGGCGCTATCGCCCTGCTGGTTTATTTAGTCATCGCCGTTTCAAAGTTGCGGATGCGCAAAATTCTGCGAGCAGAAGGAAGCGAAATTCg                                                                                                   >  2:958275/1‑128 (MQ=255)                                               gcgcCTGCGAAGGTATTTAAATTTCTGATCGACAGCTCCGGCGCTATCGCCCTGCTGGTTTATTTAGTCATCGCCGTTTCAAAGTTGCGGATGCGCAAAATTCTGCGAGCAGAAGGAAGCGAAATTCg                                                                                                   <  1:958275/128‑1 (MQ=255)                                               gcgcCTGCGAAGGTATTTAAATTTCTGATCGACAGCTCCGGCGCTATCGCCCTGCTGGTTTATTTAGTCATCGCCGTTTCAAAGTTGCGGATGCGCAAAATTCTGCGAGCAGAAGGAAGCGAAATTCGCTTGCGAATGTGGCTTTATcc                                                                              <  1:915556/149‑1 (MQ=255)                                                   cTGCGAAGGTATTTAAATTTCTGATCGACAGCTCCGGCGCTATCGCCCTGCTGGTTTATTTAGTCATCGCCGTTTCAAAGTTGCGGATGCGCAAAATTCTGCGAGCAGAAGGAAGCGAAATTCGCTTGCGAATGTGGCTTTATCCGTgg                                                                          <  1:810366/149‑1 (MQ=255)                                                         aGGTATTTAAATTTCTGATCGACAGCTCCGGCGCTATCGCCCTGCTGGTTTATTTAGTCATCGCCGTTTCAAAGTTGCGGATGCGCAAAATTCTGCGAGCAGAAGGAAGCGAAATTCGCTTGCGAATGTGGCTTTATCCGTGGCTCAcc                                                                    <  1:1025159/149‑1 (MQ=255)                                                           gTATTTAAATTTCTGATCGACAGCTCCGGCGCTATCGCCCTGCTGGTTTATTTAGTCATCGCCGTTTCAAAGTTGCGGATGCGCAAAATTCTGCGAGCAGAAGGAAGCGAAATTCGCTTGCGAATGTGGCTTTATCCGTGGCTCACCTg                                                                  >  1:1133538/1‑149 (MQ=255)                                                             aTTTAAATTTCTGATCGACAGCTCCGGCGCTATCGCCCTGCTGGTTTATTTAGTCATCGCCGTTTCAAAGTTGCGGATGCGCAAAATTCTGCGAGCAGAAGGAAGCGAAATTCGCTTGCGAATGTGGCTTTATCCGTGGCTCACCTGGc                                                                <  2:121026/149‑1 (MQ=255)                                                               ttAAATTTCTGATCGACAGCTCCGGCGCTATCGCCCTGCTGGTTTATTTAGTCATCGCCGTTTCAAAGTTGCGGATGCGCAAAATTCTGCGGGCGGAAGGAGGCGAAATTCGCTTGCGAATGTGGCTTTATCCGTGGCTCACCTGGCTa                                                              >  2:860198/1‑149 (MQ=255)                                                               ctAAGTTTCTGATCGACAGCTCCGGCGCTATCGCCCTGCTGGTTTATTTAGTCATCGCCGTTTCAAAGTTGCGGATGCGCAAAATTCTGCGAGCAGAAGGAAGCGAAATTCGCTTGCGAATGTGGCTTTATCCGTGGCTCACCTGGCTa                                                              >  2:1107556/2‑149 (MQ=255)                                                                tacaTTTCTGATCGACAGCTCCGGCGCTATCGCCCTGCTGGTTTATTTAGTCATCGCCGTTTCAAAGTTGCGGATGCGCAAAATTCTGCGAGCAGAAGGAAGCGAAATTCGCTTGCGAATGTGGCTTTATCCGTGGCTCACCTGGCTAg                                                             <  1:698243/146‑1 (MQ=255)                                                                    tttCTGATCGACAGCTCCGGCGCTATCGCCCTGCTGGTTTATTTAGTCCTCGCCCTTTCAAAGTTGCGGCTGCGCAACATTCTGCGAGCAGAAGGCAGCGAATTTCGCTTGTGATTGTGGCTGTATCCGTGGCTGCCCTGGATAGTCAt                                                         >  2:386246/1‑149 (MQ=255)                                                                    tttCTGATCGACAGCTCCGGCGCTATCGCCCTGCTGGTTTATTTAGTCATCGCCGTTTCAAAGTTGCGGATGCGCAAAATTCTGCGAGCAGAAGGAAGCGAAATTCGCTTGCGAATGTGGCTTTATCCGTGGCTCACCTGGCTAGTCAt                                                         <  1:812689/149‑1 (MQ=255)                                                                     ttCTGATCGACAGCTCCGGCGCTATCGCCCTGCTGGTTTATTTAGTCATCGCCGTTTCAAAGTTGCGGATGCGCAAAATTCTGCGAGCAGAAGGAAGCGAAATTCACTTGCGAATGTGGCTTTATCCGTGGCTCACCTGGCTACTCATc                                                        >  2:968860/1‑149 (MQ=255)                                                                      tCTGTTCGACGGCTGCGGCGCTATCGCCCTTCTGGTTTATTTAGTCTTCGCCGTTTCAAAGTTGCGGATGCGCAAAATTCTGCGAGCAGAAGGAAGCGAAATTCGCTTGCGAATGTGGCTTTATCCGTGGCTCACCTGGCTAGTCATCg                                                       <  2:337855/149‑1 (MQ=255)                                                                          aTCGACAGCTCCGGCGCTATCGCCCTGCTGGTTTATTTAGTCATCGCCGTTTCAAAGTTGCGGATGCGCAAAATTCTGCGAGCAGAAGGAAGCGAAATTCGCTTGCGAATGTGGCTTTATCCGTGGCTCACCTGGCTAGTCATCGGCtt                                                   >  2:837601/1‑149 (MQ=255)                                                                                aGCTCCGGCGCTATCGCCCTGCTGGTTTATTTAGTCATCGCCGTTTCAAAGTTGCGGATGCGCAAAATTCTGCGAGCAGAAGGAAGCGAAATTCGCTTGCGAATGTGGCTTTATCCGTGGCTCACCTGGCTAGTCATCGGCTTTATTAc                                             >  1:187654/1‑149 (MQ=255)                                                                                    ccGGGGCTATCGCCCTGCTGGTTTGTTTAGTCATCGCCGTTTCAAAGTTGCGGATGCGCAAAATTCTGCGAGCAGAAGGAAGCGCAATTCGCTTGCGAATGTGGCTTTATCCGTGGCTCACCTGGCTAGTCATCGGCTTTATTACCttt                                         <  2:133930/149‑1 (MQ=255)                                                                                           tATCGCCCTGCTGGTTTATTTAGTCATCGCCGTTTCAAAGTTGCGGATGCGCAAAATTCTGCGAGCAGAAGGAAGCGAAATTCGCTTGCGAATGTGGCTTTATCCGTgg                                                                          <  1:689264/109‑1 (MQ=255)                                                                                           tATCGCCCTGCTGGTTTATTTAGTCATCGCCGTTTCAAAGTTGCGGATGCGCAAAATTCTGCGAGCAGAAGGAAGCGAAATTCGCTTGCGAATGTGGCTTTATCCGTgg                                                                          >  2:689264/1‑109 (MQ=255)                                                                                           tATCGCCCTGCTGGTTTATTTAGTCATCGCCGTTTCAAAGTTGCGGATGCGCAAAATTCTGCGAGCAGAAGGAAGCGAAATTCGCTTGCGAATGTGGCTTTATCCGTGGCTCACCTGGCTAGTCATCGGCTTTATTACCTTTGTGTTgg                                  >  2:380036/1‑149 (MQ=255)                                                                                             ttggccTGCTGGTTTATTTAGTCATGGGCGTTTCAAAGTTGCGGATGCGCAAAATTCTGCGAGCAGAAGGAAGCGAAATTCGCTTGCGAATGTGGCTTTATCCGTGGCTCACCTGGCTAGTCATCGGCTTTATTACCTTTGTGTTGGTa                                <  2:149246/145‑1 (MQ=255)                                                                                                cccTGCTGGTTTATTTAGTCATCGCCGTTTCAAAGTTGCGGATGCGCAAAATTCTGCGAGCAGAAGGAAGCGAAATTCGCTTGCGAATGTGGCTTTATCCGTGGCTCACCTGGCTAGTCATCGGCTTTATTACCTTTGTGTTGGTAGTg                             >  2:114128/1‑149 (MQ=255)                                                                                                     ctgGTTTATTTAGTCATCGCCGTCTCAAAGTTGCGGATGCGCAAAATTCTGCGAGCAGAAGGAAGCGAAATTCGCTTGCGAATGTGGCTTTATCCGTGGCTCACCTGGCTAGTCATCGGCTTTATTACCTTTGTGTTGGTAGTGATGCt                        >  2:131242/1‑149 (MQ=255)                                                                                                        gTTTATTTAGTCATCGCTGTTTCAAAGTTGCGGATGCGCAAAATTCTGCGAGCAGAAGGAAGCGAAATTCGCTTGCGAATGTGGTTTTATCCGTGGCTCACCTGGCTAGTCATCGGCTTTATTAGCTTTGTGTTGGTAGTGATGCTAtt                     <  1:331997/149‑1 (MQ=255)                                                                                                        gTTTATTTAGTCATCGCCGTTTCAAAGTTGCGGATGCGTAAAATTTTGGGAGTAGAAGGAAGCGAAATTCGCTTGCGAATGTGGCTTTATCCGTGGCTCACCTGGCTAGGCATCGGCTTTATTACCTTTGTGTTGGTAGTGATGCTAtt                     <  1:363333/149‑1 (MQ=255)                                                                                                         tttatttaGTCATCGCCGTTTCAAAGTTGCGGATGCGCAAAATTCTGCGAGCAGAAGGAAGCGAAATTCGCTTGCGAATGTGGCTTTATCCGTGGCTCACCTGGCTAGTCATCGGCTTTATTACCTTTGTGTTGGTAGTGATGCTAt                      <  2:579152/147‑1 (MQ=255)                                                                                                         tttatttaGTCATCGCCGTTTCAAAGTTGCGGATGCGCAAAATTCTGCGAGCAGAAGGAAGCGAAATTCGCTTGCGAATGTGGCTTTATCCGTGGCTCACCTGGCTAGTCATCGGCTTTATTACCTTTGTGTTGGTAGTGATGCTAt                      >  1:579152/1‑147 (MQ=255)                                                                                                          ttatttaGTCATCGGGGTTTCAAAGTTGCGGATGCGCAAAATTCTGCGAGCAGAAGGAAGCGAAATTCGCTTGCGAATGTGGCTTTATCCGTGGCTCACCTGGCTAGTCATCGGCTTTATTACCTTTGTGTTGGTAGTGATGCTATtcc                   <  1:367161/149‑1 (MQ=255)                                                                                                             tttaGTCATCGCCGTTTCAAAGTTGCGGATGCGCAAAATTCTGCGAGCAGAAGGAAGCGAAATTCGCTTGCGAATGTGGCTTTATCCGTGGCTCACCTGGCTAGTCATCGGCTTTATTACCTTTGTGTTGGTAGTGATGCTATtccgtc                <  1:131242/149‑1 (MQ=255)                                                                                                             tttaGTCATCGCCGTTTCAAAGTTGCGGATGCGCAAAATTCTGCGAGCAGAAGGAAGCGAAATTCGCTTGCGAATGTGGCTTTATCCGTGGCTCACCTGGCTAGTCATCGGCTTTATTACCTTTGTGTTGGTAGTCATGCTATtccgtc                >  2:606306/1‑149 (MQ=255)                                                                                                                     tCGCCGTTTCAAAGTTGCGGATGCGCAAAATTCTGCGAGCAGAAGGAAGCGAAATTCGCTTGCGAATGTGGCTTTATCCGTGGCTCACCTGGCTAGTCATCGGCTTTATTACCTTTGTGTTGGTGGTGATGCTATTCCGTCCGgcgcga        >  1:948465/1‑147 (MQ=255)                                                                                                                           tttCAAAGTTGCGGATGCGCAAAATTCTGCGAGCAGAAGGAAGCGAAATTCGCTTGCGAATGTGGCTTTATCCGTGGCTCACCTGGCTAGTCATCGGCTTTATTACCTTTGTGTTGGTAGTGATGCTATTCCGTCCGGCGCAACAGTTa  >  1:682716/1‑149 (MQ=255)                                                                                                                                |                                                                                                                                                ACTCTCCACCGGCGCGGCATTCTTAACGGTGGTGGTGAACTATTACGCGCCTGCGAAGGTATTTAAATTTCTGATCGACAGCTCCGGCGCTATCGCCCTGCTGGTTTATTTAGTCATCGCCGTTTCACAGTTGCGGATGCGCAAAATTCTGCGAGCAGAAGGAAGCGAAATTCGCTTGCGAATGTGGCTTTATCCGTGGCTCACCTGGCTAGTCATCGGCTTTATTACCTTTGTGTTGGTAGTGATGCTATTCCGTCCGGCGCAACAGTTA  >  NC_000913_3_pae_tpiA/2795269‑2795539

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