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

CCCGTACGTGGCGGTGTTACTCTCCACCGGCGCGGCATTCTTAACGGTGGTGGTGAACTATTACGCGCCTGCGAAGGTATTTAAATTTCTGATCGACAGCTCCGGCGCTATCGCCCTGCTGGTTTATTTAGTCATCGCCGTTTCACAGTTGCGGATGCGCAAAATTCTGCGAGCAGA‑AGGAAGCGAAATTCGCTTGCGAATGTGGCTTTATCCGTGGCTCACCTGGCTAGTCATCGGCTTTATTACCTTTGTGTTGGTAGTGATGCTATTCCGTCCGGCGCAACAGTTAG  >  NC_000913_3_pae_tpiA/2795251‑2795540                                                                                                                                                  |                                                                                                                                                  cccGTACGTGGCGGTGTTACTCTCCACCGGCGCGGCATTCTTAACGGTGGTGGTGAACTATTACGCGCCTGCGAAGGTATTTAAATTTCTGATCGACAGCTCCGGCGCTATCGCCCTGCTGGTTTATTTAGTCATCGCCGTTTCAAAGt                                                                                                                                                >  2:582792/1‑149 (MQ=255)   cGTACGTGGCGGTGTTACTCTCCACCGGCGCGGCATTCTTAACGGTGGTGGTGAACTATTACGCGCCTGCGAAGGTATTTAAATTTCTGATCGACAGCTCCGGCGCTATCGCCCTGCTGGTTTATTTAGTCATCGCCGTTTCAAAGTTg                                                                                                                                              >  2:452663/1‑149 (MQ=255)    gTACGTGGCGGTGTTACTCTCCACAGGCGCGGCATTCTTAACGGTGGTGGTGAACTATTACCCGCCTGCGAAGGTATTTAAATTTCTGATCGACAGCTCTGGCGGTATCGCGCTGTTGGTTTATTTGGTCTTCTGGGTTTGAAAGTTgg                                                                                                                                             >  1:105183/1‑148 (MQ=255)             ggTGTTACTCTCCACCGGCGCGGCATTCTTAACGGTGGTGGTGAACAATTCAGCTCTTGCGAAGGTATTTAAATTTCTGATCGACAGCTCCGGCGCTATCGCCCTGCTGGTTTCTTTGGTCATCGCCGTTTCAAACTTGGGGaacggca                                                                                                                                    >  1:140329/1‑143 (MQ=255)                     tctcCACCGGCGCGGCATTCTTAACGGTGGTGGTGAACTATTACGCGCCTGCGAAGGTATTTAAATTTCTGATCGACAGCTCCGGCGCTATCGCCCTGCTGGTTTATTTAGTCATCGCCGTTTCAAAGTTGCGGATGCGCAAAATTCTg                                                                                                                            <  1:397559/149‑1 (MQ=255)                     tctcCACCGGCGCGGCATTCTTAACGGTGGTGGTGAACTATTACGCGCCTGCGAAGGTATTTAAATTTCTGATCGACAGCTCCGGCGCTATCGCCCTGCTGGTTTATTTAGTCATCGCCGTTTCAAAGTTGCGGATGCGCAAAATTCTg                                                                                                                            <  2:401664/149‑1 (MQ=255)                        ccACCGGCGCGGCATTCTTAACGGGGGTGGTGAACTATTACGCGCCTGCGAAGGTATTTAAATTTCTGATCGACAGCTCCGGCGCTATCGCCCTGCTGGTTTATTTAGTCATCGCCGTTTCAAAGTCGCGGATGCGCAAAATTCTGCGa                                                                                                                         <  2:617693/149‑1 (MQ=255)                            cGGCGCGGCATTCTTAACGGTGGTGGTGAACTATTACGCGCCTGCGAAGGTATTTAAATTTCTGATCGACAGCTCCGGCGCTATCGCCCTGCTGGTTTATTTAGTCATCGCCGTTTCAAAGTTGCGGATGCGCAAAATTCTGCGAGCAg                                                                                                                     >  2:683769/1‑149 (MQ=255)                             ggCGCGGCATTCTTAACGGTGGTGGTGAACTATTACGCGCCTGCGAAGGTATTTAAATTTCTGATCGACAGCTCCGGCGCTATCGCCCTGCTGGTTTATTTAGTCATCGCCGTTTCAAAGTTGCGGATGCGCAAAATTCTGCGAGCAga                                                                                                                    <  1:512621/149‑1 (MQ=255)                                gcgGCATTCTTAACGGTGGTGGTGAACTATTACGCGCCTGCGAAGGTATTTAAATTTCTGATCGACAGCTCCGGCGCTATCGCCCTGCTGGTTTATTTAGTCATCGCCGTTTCAAAGTTGCGAATGCGCAAAATTCTGCGAGCAga‑agg                                                                                                                >  1:186156/1‑149 (MQ=255)                                  ggCATTCTTAACGGTGGTGGTGAACTATTACGCGCCTGCGAAGGTATTTAAATTTCTGATCGACAGCTCCGGCGCTATCGCCCTGCTGGTTTATTTAGTCATCGCCGTTTCAAAGTTGCGGATGCGCAAAATTCTGCGAGCAga‑aggaa                                                                                                              >  1:104689/1‑149 (MQ=255)                                  ggCATGCTTAACGGTGGTGGTGAACTATTACGCGCCTGCGAAGGTATTTAAATTTCTGATCGACAGCTCCGGCGCTATCGCCCTGCTGGTTTATTTAGTCATCGCCGTTTCAAAGTTGCGGATGCGCAAAATTCTGCGAGCAga‑aggaa                                                                                                              <  2:195732/149‑1 (MQ=255)                                   gCATTCTTAACGGTGGTGGTGAACTATTACGCGCCTGCGAAGGTATTTAAATTTCTGATCGACAGCTCCGGCGCTATCGCCCTGCTGGTTTATTTAGTCATCGCCGTTTCAAAGTTGCGGATGCGCAAAATTCTGCGAGCAga‑aggaag                                                                                                             >  1:297058/1‑149 (MQ=255)                                             tggtggtggtGAACTATTACGCGCCTGCGAAGGTATTTAAATTTCTGATCGACAGCTCCGGCGCTATCGCCCTGCTGGTTTATTTAGTCATCGCCGTTTCAAAGTTGCGGATGCGCAAAATTCTGCGAGCAGA‑AGGAAGCGAAATTCGc                                                                                                   <  1:708228/148‑1 (MQ=255)                                                  gtggtgAACTATTACGCGCCTGCGAAGGTATTTAAATTTCTGATCGACAGCTCCGGCGCTATCGCCCTGCTGGTTTATTTAGTCATCGCCGTTTCAAAGTTGCGGATGCGCAAAATTCTGCGGGCGGA‑CGGAAGCGAAATTCGCTTGCg                                                                                              >  1:152743/1‑149 (MQ=255)                                                     gtgAACTATTACGCGCCTGCGAAGGTATTTAAATTTCTGATCGACAGCTCCGGCGCTATCGCCCTGCTGGTTTATTTAGTCATCGCCGTTTCAAAGTTGc                                                                                                                                             >  2:535201/1‑100 (MQ=255)                                                     gtgAACTATTACGCGCCTGCGAAGGTATTTAAATTTCTGATCGACAGCTCCGGCGCTATCGCCCTGCTGGTTTATTTAGTCATCGCCGTTTCAAAGTTGc                                                                                                                                             <  1:535201/100‑1 (MQ=255)                                                           tATTACGCGCCTGCGAAGGTATTTAAATTTCTGATCGACAGCTCCGGCGCTATCGCCCTGCTGGTTTATTTAGTCATCGCCGTTTCAAAGTTGCGGATGCGCAAAATTCTGCGAGCAGA‑AGGAAGCGAAATTCGCTTGCGAATGTGGCt                                                                                     <  1:84054/149‑1 (MQ=255)                                                            aTTACGCGCCTGCGAAGGTATTTAAATTTCTGATCGACAGCTCCGGCGCTATCGCCCTGCTGGTTTATTTAGTCATCGCCGTTTCAAAGTTGCGGATGCGCAAAATTCTGCGAGCAGA‑AGGAAGCGAAATTCGCTTGCGAATGTGGCtt                                                                                    >  1:706440/1‑149 (MQ=255)                                                            aTTACGCGCCTGCGAAGGTATTTAAATTTCTGATCGACAGCTCCGGCGCTATCGCCCTGCTGGTTTATTTAGTCATCGCCGTTTCAAAGTTGCGGATCCGCAAAATTCTGCGAGCAGA‑AGGAAGCGAAATTCGCTTGCGAATGTGGCtt                                                                                    <  1:670847/149‑1 (MQ=255)                                                                   gcCTGCGAAGGTATTTAAATTTCTGATCGACAGCTCCGGCGCTATCGCCCTGCTGGTTTATTTAGTCATCGCCGTTTCAAAGTTGCGGATGCGCAAAATTCTGCGAGCAGA‑AGGAAGCGAAATTCGCTTGCGAATGTGGCTTTATCCGt                                                                             >  1:214216/1‑149 (MQ=255)                                                                    ccTGCGAAGGTATTTAAATTTCTGATCGACAGCTCCGGCGCTATCGCCCTGCTGGTTTATTTAGTCATCGCCGTTTCAAAGTTGCGGATGCGCAAAATTCTGCGAGCAGA‑AGGAAGCGAAATTCGCTTGCGAATGTGGCTTTATCCGTg                                                                            <  1:452663/149‑1 (MQ=255)                                                                     cTGCGAAGGTATTTAAATTTCTGATCGACAGCTCCGGCGCTATCGCCCTGCTGGTTTATTTAGTCATCGCCGTTTCAAAGTTGCGGATGCGCAAAATTCTGCGAGCAGA‑AGGAAGCGAAATTCGCTTGCGAATGTGGCTTTATCCGTgg                                                                           <  2:297058/149‑1 (MQ=255)                                                                                tttAAATTTCTGATCGACAGCTCCGGCGCTATCGCCCTGCTGGTTTATTTAGTCATCGCCGTTTCAAAGTTGCGGATGCGCAAAATTCTGCGAGCa                                                                                                                      >  1:344190/1‑96 (MQ=255)                                                                                tttAAATTTCTGATCGACAGCTCCGGCGCTATCGCCCTGCTGGTTTATTTAGTCATCGCCGTTTCAAAGTTGCGGATGCGCAAAATTCTGCGAGCa                                                                                                                      <  2:344190/96‑1 (MQ=255)                                                                                 ttAAATTTCTGATCGACAGCTCCGGCGCTATCGCCCTGCTGGTTTATTTAGTCATCGCCGTTTCAAAGTTGCGGATGCGCAAAATTCTGCGAGCAga‑aggaag                                                                                                             <  2:581685/103‑1 (MQ=255)                                                                                 ttAAATTTCTGATCGACAGCTCCGGCGCTATCGCCCTGCTGGTTTATTTAGTCATCGCCGTTTCAAAGTTGCGGATGCGCAAAATTCTGCGAGCAga‑aggaag                                                                                                             >  1:581685/1‑103 (MQ=255)                                                                                 ttAAATTTCTGATCGACAGCTCCGGCGCTATCGCCCTGCTGGTTTATTTAGTCATCGCCGTTTCAAAGTTGCGGATGCGCAAAATTCTGCGAGCAGA‑AGGAAGCGAAATTCGCTTGCGAATGTGGCTTTATCCGTGGCTCACCt                                                                    <  1:543324/144‑1 (MQ=255)                                                                                 ttAAATTTCTGATCGACAGCTCCGGCGCTATCGCCCTGCTGGTTTATTTAGTCATCGCCGTTTCAAAGTTGCGGATGCGCAAAATTCTGCGAGCAGA‑AGGAAGCGAAATTCGCTTGCGAATGTGGCTTTATCCGTGGCTCACCt                                                                    >  2:543324/1‑144 (MQ=255)                                                                                 ttAAATTTCTGATCGACAGCTCCGGCGCTATCGCCCTGCTGGTTTATTTAGTCATCGCCGTTTCAAAGTTGCGGATGCGCAAAATTCTGCGAGCAGA‑AGGAAGCGAAATTCGCTTGCGAATGTGGCTTTATCCGTGGCTCACCTGGCTa                                                               >  1:106173/1‑149 (MQ=255)                                                                                  tAAATTTCTGATCGACAGCTCCGGCGCTATCGCCCTGCTGGTTTATTTAGTCATCGCCGTTTCAAAGTTGCGGATGCGCAAAATTCTGCGAGCAGA‑AGGAAGCGAAATTTGCTTGCGAATGTGGCTTTATCCGTGGCTAACCTGGCTAg                                                              >  1:432401/1‑149 (MQ=255)                                                                                        tCTGATCGACAGCTCCGGCGCTATCGCCCTGCTGGTTTATTTAGTCATCGCCGTTTCAAAGTTGCGGATGCGCAAAATTCTGCGAGCAGA‑AGGAAGCGAAATTCGCTTGCGAATGTGGCTTTATCCGTGGCTCACCTGGCTAGTCATCg                                                        >  2:393651/1‑149 (MQ=255)                                                                                            aTCGACAGCTCCGGCGCTATCGCCCTGCTGGTTTATTTAGTCATCGCCGTTTCAAAGTTGCGGATGCGCAAAATTCTGCGAGCAGA‑AGGAAGCGAAATTCGCTTGCGAATGTGGCTTTATCCGTGGCTCACCTGGCTAGTCATCGGCtt                                                    <  1:110532/149‑1 (MQ=255)                                                                                              cGACAGCTCCGGCGCTATCGCCCTGCTGGTTTATTTAGTCATCGCCGTTTCAAAGTTGCGGATGCGCAAAATTCTGCGAGCAGA‑AGGAAGCGAAATTCGCTTGCGAATGTGGCTTTATCCGTGGCTCACCTGGCTAGTCATCGGCTtta                                                  <  2:104689/149‑1 (MQ=255)                                                                                               gACAGCTCCGGCGCTATCGCCCTGCTGGTTTATTTAGTCATCGCCGTTTCAAAGTTGCGGATGCGCAAAATTCTGCGAGCAGA‑AGGAAGCGAAATTCGCTTGCGAATGTGGCTTTATCCGTGGCTCACCTGGCTAGTCATCGGCGttat                                                 >  2:692917/1‑149 (MQ=255)                                                                                                aCAGCTCCGGCGCTATCGCCCTGCTGGTTTATTTAGTCATCGCCGTTTCAAAGTTGCGGATGCGCACAATTCTGCGAGCAGA‑AGGAAGCGAAATTCGATTGCGAATGTAGCTTTATCCGTGGCACACCTGGCTCGCCATCGGCGttatt                                                >  2:20806/1‑149 (MQ=255)                                                                                                   gCTCCGGCGCTATCGCCCTGCTGGTTTATTTAGTCATCGCCGTTTCAAAGTTGCGGATGCGCAAAATTCTGCGAGCAGA‑AGGAAGCGAAATTCGCTTGCGAATGTGGCTTTATCCGTGGCTCACCTGGCTAGTCATCGGCTTTATTAcc                                             <  1:582792/149‑1 (MQ=255)                                                                                                            cTATCGCCCTGCTGGTTTATTTAGTCATCGCCGGTTCAAAGTTGCGGATGCGCAAAATTCTGCGAGCAGA‑AGGAAGCGAAATTCGCTTGCGAATGTGGCTTTATCCGTGGCTCACCTGGCTAGTCATCGGCTTTATTACCTTTGTGTTg                                    <  1:22535/149‑1 (MQ=255)                                                                                                             tATCGCGCTGATGGTTTATTTAGTCATCGCCGTTTCAAAGTGGCGGATGCGCAAAATTCTGCGAGCAGA‑AGGAAGCGAAATTCGCTTGCGAATGTGGCTTTATCCGTGGCTTACCTGGCTAGTCATCGGCTTTATTACCTTTGTGTTgg                                   <  2:279452/149‑1 (MQ=255)                                                                                                                   ccTGCTGGTTTATTTAGTCATCGCCGTTTCAAAGTTGCGGATGCGCAAAATTCTGCGAGCAGA‑AGGAAGCGAAATTCGCTTGCGAATGTGGCTTTATCCGTGGCTCACCTGGc                                                                 <  2:315619/113‑1 (MQ=255)                                                                                                                   ccTGCTGGTTTATTTAGTCATCGCCGTTTCAAAGTTGCGGATGCGCAAAATTCTGCGAGCAGA‑AGGAAGCGAAATTCGCTTGCGAATGTGGCTTTATCCGTGGCTCACCTGGc                                                                 >  1:315619/1‑113 (MQ=255)                                                                                                                          gTTTATTTAGTCATCGCCGTTTCAAAGTTGCGGATGCGCAAAATTCTGCGAGCAGA‑AGGAAGCGAAATTCGCTTGCGAATGTGGCTTTATCCGTGGCTCACCTGGCTAGTCATCGGCTTTATTACCTTTGTGTTGGTAGTGATGCTAtt                      <  2:106173/149‑1 (MQ=255)                                                                                                                          gTTTATTTAGGCATCGCCGTTTCAAAGTTGCGGATGCGCAAAATTATGC‑AGAAGACGGGAAGCGAAATTCGCTTGCGAATGTGGCTTTATCCGTGGCTCACCTGGCTAGTCATCGGCTTTATTACCTTTGTGTTGGTAGTGATGCTAtt                      <  1:20806/149‑1 (MQ=255)                                                                                                                               tttaGTCATCGCCGTTTCAAAGTTGCGAATGCGCAAAATTCTGCGAGCAGA‑AGGAAGCGAAATTCGCTTGCGAATGTGGCTTTATCCGTGGCTCAACTGGCTAGTCATCGGCTTTATTACCTTTGTGTTGGTAGTGATGCTATtccgtc                 >  2:65917/1‑149 (MQ=255)                                                                                                                                ttaGTCATCGCCGTTTCAAAGTTGCGGATGCGCAAAATTCTGCGAGCAGA‑AGGAAGCGAAATTCGCTTGCGAATGTGGCTTTATCCGTGGCTCACCTGGCTAGTCATCGGCTTTATTACCTTTGTGTTGGTAGTGATGCTATtccgtcc                >  1:650244/1‑149 (MQ=255)                                                                                                                                  aGTCATCGCCGTTTCAAAGTTGCGGATGCGCAAAATTCTGCGAGCAGA‑AGGAAGCGAAATTCGCTTGCGAATGTGGCTTTATCCGTGGCTCACCTGGCTAGTCATCGGCTTTATTACCTTTGTGTTGGTAGTGATGCTATTCCGTCCgg              >  1:391544/1‑149 (MQ=255)                                                                                                                                   gTCATCGCCGTTTCAAAGTTGCGGAGGCGCAAAATTCTGCGAGCAGA‑AGGAAGCGAAATTCGCTTGCGAATGTGGCTTTATCCGTGGCTCACCTGGCTAGTCATCGGCTTTATTACCTTTGTGTTGGTAGTGATGCTATTCCATCCGgc             >  1:644733/1‑149 (MQ=255)                                                                                                                                      aTCGCCGTTTCAAAGTTGCGGATGCGCAAAATTCTTCGAGCAGA‑AGGAAGCGAAATTCGCTTGCGAATGTGGCTTTATCCGTGGCTCACCTGGCTAGTCATCGGCTTTATTACCTTTGTGTTGGTAGTGATGCTATTCCGTCCGGAGCa          >  1:425726/1‑149 (MQ=255)                                                                                                                                      aTCGCCGTTTCAAAGTTGCGGATGCGCAAAATTCTGCGAGCAGA‑AGGAAGCGAAATTCGCTTGCGAATGTGGCTTTATCCGTGGCTCACCTGGCTAGTCATCGGCTTTATTACCTTTGTGTTGGTAGTGATGCTATTCCGTCCGGCGCa          >  1:343587/1‑149 (MQ=255)                                                                                                                                             tttCAAAGTTGCGGATGCGCAAAATTCTGCGAGCAGA‑AGGAAGCGAAATTCGCTTGCGAATGTGGCTTTATCCGTGGCTCACCTGGCTAGTCATCGGCt                                                     >  1:119456/1‑99 (MQ=255)                                                                                                                                             tttCAAAGTTGCGGATGCGCAAAATTCTGCGAGCAGA‑AGGAAGCGAAATTCGCTTGCGAATGTGGCTTTATCCGTGGCTCACCTGGCTAGTCATCGGCTTTATTACCTTTGTGTTGGTAGTGATGCTATTCCGTCCGGCCCAACCGTTa   >  1:684973/1‑149 (MQ=255)                                                                                                                                             tttCAAAGTTGCGGATGCGCAAAATTCTGCGAGCAGA‑AGGAAGCGAAATTCGCGTGCGAATGTGGCTTTATCCGTGGCTCACCTGGCTAGTCATCGGCt                                                     <  2:119456/99‑1 (MQ=255)                                                                                                                                              ttCAAAGTTGCGGATGCGCAAAATTCTGCGAGCAGA‑AGGAAGCGAAATTCGCTTGCGAATGTGGCTTTATCCGTGGCTCACCTGGCTAGTCATCGGCTTTATTACCTTTGTGTTGGTAGTGATGCTATTCCGTCCGGCGCAACAGTTAg  <  1:221319/149‑1 (MQ=255)                                                                                                                                                  |                                                                                                                                                  CCCGTACGTGGCGGTGTTACTCTCCACCGGCGCGGCATTCTTAACGGTGGTGGTGAACTATTACGCGCCTGCGAAGGTATTTAAATTTCTGATCGACAGCTCCGGCGCTATCGCCCTGCTGGTTTATTTAGTCATCGCCGTTTCACAGTTGCGGATGCGCAAAATTCTGCGAGCAGA‑AGGAAGCGAAATTCGCTTGCGAATGTGGCTTTATCCGTGGCTCACCTGGCTAGTCATCGGCTTTATTACCTTTGTGTTGGTAGTGATGCTATTCCGTCCGGCGCAACAGTTAG  >  NC_000913_3_pae_tpiA/2795251‑2795540

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