Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913_3_vch_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_vch_tpiA	2,795,396	0	C	A	100.0%	165.8 / NA	50	Q382K (CAG→AAG)	gabP	gamma‑aminobutyrate transporter
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base A (24/26);  total (24/26)

TACGTGGCGGTGTTACTCTCCACCGGCGCGGCATTCTTAACGGTGGTGGTGAACTATTACGCGCCTGCGAAGGTATTTAAATTTCTGATCGACAGCTCCGGCGCTATCGCCCTGCTGGTTTATTTAGTCATCGCCGTTTCACAGTTGCGGATGCGCAAAATTCTGCGAGCAGAAGGAAGCGAAATTCGCTTGCGAATGTGGCTTTATCCGTGGCTCACCTGGCTAGTCATCGGCTTTATTACCTTTGTGTTGGTAGTGATGCTATTCCGTCCGGCGCAACAG  >  NC_000913_3_vch_tpiA/2795255‑2795536                                                                                                                                              |                                                                                                                                             tACGTGGCGGTGGTACTCTCCACCGGCGCGGCATTCTTAACGGTGGTGGTGAACTATTACGCGCCTGCGAAGGTATTTAAATTTCTGATCGACAGCTCCGGCGCTTTCGCCCTGCTGGTTTATTTAGTCATCGCCGTTTCAAAGTTGCg                                                                                                                                       <  1:310025/149‑1 (MQ=255)     tGGCGGTGTTACTCTCCACCGGCGCGGCATTCTTAACGGTGGTGGTGAACTATTACGCGCCTGCGAAGGTATTTAAATTTCTGATCGACAGCTCCGGCGCTATCGCCCTGCTGCTTTATTTAGTCATCGCCGTTTCAAAGTTGCGGATg                                                                                                                                   <  2:399807/149‑1 (MQ=255)            gtTACTCTCCACCGGCGCGGCATTCTTAACGGTGGTGGTGAACTATTACGCGCCTGCGAAGGTATTTAAATTTCTGATCGACAGCTCCGGCGCTATCGCCCTGCTGGTTTATTTAGTCATCGCCGTTTCAAAGTTGCGGATGCGCaaaa                                                                                                                            >  2:528227/1‑149 (MQ=255)             ttACTCTCCACCGGCGCGGAATTCTTAACGGTGGTGGTGAACTATTACGCGCCTGCGAAGGTATTTAAATTTCTGATCGACAGCTCCGGCGCTATCGCCCTGCTGGTTTATTTAGTCATCGCCGTTTCAAAGTTGCGGATGCGCaaact                                                                                                                           >  2:662679/1‑147 (MQ=255)                    tcACCGGCGCGGCATTCTTAACGGTCGTGGTGAACTATTACGCGCCTGCGAAGGTATTTAAATTTCTGATCGACAGCTCCGGCGCTCTCGCGGTGCTGGATTCGATATTCAGCGCCGTTTCAAAGTTGCGGATGCGCAAAATTCTGCGa                                                                                                                    <  1:773368/148‑1 (MQ=255)                      aCCGGCGCGGCATTCTTAACGGTGGTGGTGAACTATTACGCGCCTGCGAAGGTATTTAAATTTCTGCTCGACAGCTCCGGCGCTATCGCCCTGCTGGTTTATTTAGTTTTCGCCGTTTCAAAGTTGCGGATGCGGAAAATTCTGCGAgt                                                                                                                  >  2:773368/1‑148 (MQ=255)                            gcgGCATTCTTAACGGTGGTGGTGAACAATTACGCTCCTTCGAAGGTATTCAAATTTCTGATCGACAGCTCCGGCGCTATCGCCCTGCTGGTTTATTTAGTCATCGCCGTTTCAAAGTTGCGGATGCGCAAAATTCTGCGAGCAgaagg                                                                                                            <  1:302069/149‑1 (MQ=255)                              ggCATTCTTAACGGTGGTGGTGAACTATTACGCGCCTGCGAAGGTATTTAAATTTCTGATCGACAGCTCCGGCGCTATCGCCCTGCTGGTTTATTTAGTCATCGCCGTTTCAAAGTTGCGGATGCGCAAAATTCTGCGAGCAgaaggaa                                                                                                          >  1:651795/1‑149 (MQ=255)                               gtatacTTATCGGTGGTGGTGAACTATTACGTGCCTGCGAAAGTATTTTAATTTCTGAGCGACAGCGCCGTCGCTATCGCCCTGCTGGTTTATTTAGTCATCGCCGATTCAAAGTTGCGGATGCGCAAAATTCTGCGAGCAgaaggaag                                                                                                         <  2:485766/144‑1 (MQ=255)                                 aTTCTTAACGGTGGTGGTGAACTATTACGCGCCTGCGAAGGTATTTAAATTTCTGATCGACAGCTCCGGCGCTATCGCCCTGCTGGTTTATTTAGTCATCGCCGTTTCAAAGTTGCGGATGCGCAAAATTCTGCGAGCAGAAGGAAGCg                                                                                                       >  2:702932/1‑149 (MQ=255)                                 aTTCTTAACGGTGGTGGTGAACTATTACGCGCCTGAGAAGGTATTTAAATTTCTGATCGACAGCTCCGGCGCTATCGCCCTGCTGGTTTATTTAGTCATCGACGTTTCAAAGTTGCGGATGCGCAAAATTCTGCGCGCAGAAGGAAGCg                                                                                                       >  2:686534/1‑149 (MQ=255)                                  ttCTTAACGGTGGTGGTGAACTATTACGCGCCTGCGAAGGTATTTAAATTTCTGATCGACAGCTCCGGCGCTATCGCCCTGCTGGTTTATTTAGTCATCGCCGTTTCAAAGTTGCGCATGCGCAAAATTCTGCGAGCAGAAGGAAGcca                                                                                                      >  2:637218/1‑147 (MQ=255)                                     gtAACGGTGGTGGTTAACTAGTACGCGTCTGCGAAGGTATTTAAATTTCTGATCGACAGCTCCGGTGCTATCGCCTTGCTGGTTTATTTAGTCATGGCCGTTTCAAAGTTGCGGATGCGCAAAAGTCTGCGAGCAGAAGGAAGCGAAAt                                                                                                   <  2:61099/148‑1 (MQ=255)                                      gaaCGGTGGTGGTGAACTATTACGCTGCTGGGAAGGTCTTTAAATTTCTGATCGACAGCTGCGGCGCTATCGGCCTGGTGGTTTATTTAGTCATCGCCGTTTCAAAGTTGCGGATGCGCAAAATTCTGCGAGGAGAAGGAAGCGAAAtt                                                                                                  <  2:730529/148‑1 (MQ=255)                                        aCGGTGGTGGTGAACTATTACGCGCCTGCGAAGGTATTTAAATTTCTGATCGACAGCTCCGGCGCTATCGCCCTGCTGGTTTATTTAGTCATCGCCGTTTCAAAGTTGCGGATGCGCAAAATTCTGCGAGCAGAAGGAAGCGAAATTCg                                                                                                <  2:241582/149‑1 (MQ=255)                                            tggtggtgAACTATTACGCGCCTGCGAAGGTATTTAAATTTCTGATCGACAGCTCCGGCGCTATCGCCCTGCTGGTTTATTTAGTCATCGCCGTTTCAAAGTTGCGGATGCGCAAAATTCTGCGAGCAGAAGGAAGCGAAATTCGCTTg                                                                                            <  2:604805/149‑1 (MQ=255)                                            gcgtggtgCATTATTACGCGCCTGCGAAGGTATTGAAATTTTTGCTCGCCAGCTCCGGCGCTATCGCCCTGCTGGGTTATTTGGTCAGCGCCGTTTCAAAGTTGTCGTTGCGCAAAATTCTGCGAGCAGGAGGAAGCGAAATTCGCTTg                                                                                            <  2:16274/147‑1 (MQ=255)                                                   gAACTATTACGCGCCTGCGAAGGTATTTAAATTTCTGATCGACAGCTCCGGCGCTATCGCCCTGCTGGTTTATTTAGTCATCGCCGTTTCAAAGTTGCGGATGCGCAAAATTCTGCGAGCAGAAGGAAGCGAAATTCGCTTGCGAAtgt                                                                                     >  2:316533/1‑149 (MQ=255)                                                    gatcatcaCGCGCCTGGGAAGGGATTTAAAGTTCGGATCGACAGCTGCGGCGCTAGCGCCTTGCTGGTTTATTTAGTCATCGCCGTTTCAAAGTTGCGGATGCGCAAAATTCTGCGAGCAGAAGGAAGCGAAATTCGCTTGCGAAtgtg                                                                                    <  2:272056/142‑1 (MQ=255)                                                    aaCTATTACGCGCCTGCGAAGGTATTTAAATTTCTGATCGACAGCTCCGGCGCTATGGGCTGGTGGGTGTATTCACACGTCTGCAGTTCAAAGTTGCGACAGCGGAAAGTTCTGCGAGCTGGAGGAAGCGAAAGTCGCTTTCGAAggtg                                                                                    >  1:579151/1‑149 (MQ=255)                                                    aaCTATTACGCGCCTGCGAAGGTATTTAAATTTCTGATCGACAGCTCCGGCGCTATCGCCCTGCTGGTTTATTTAGTCATCGCCGTTTCAAAGTTGCGGATGCGCAAAATTCTGCGAGCAGAAGGAAGCGAAATTCGCTTGCGAAtgtg                                                                                    >  2:111066/1‑149 (MQ=255)                                                     aCTATTACGCGCCTGCGAAGGTATTTAAATTTCTGATCGACAGCTCCGGCGCTATCGCCCTGCTGGTTTATTTAGTCATCGCCGTTTCAAAGTTGCGGATGCGCAAAATTCTGCGAGCAGAAGGAAGCGAAATTCGCTTGCGAATGTgg                                                                                   <  1:127577/149‑1 (MQ=255)                                                       tATTACGCGCCTGCGAAGGTATTTAAATTTCTGATCGACAGCTCCGGCGCTATCGCCCTGCTGGTTTATTTAGTCATCGCCGTTTCAAAGTTGCGGATGCGCAAAATTCTGCGAGCAGAAGGAAGCGAAATTCGCTTGCGAATGTGGCt                                                                                 >  1:193987/1‑149 (MQ=255)                                                               gcCTGCGAAGGTATTTAAATTTCTGATCGACAGCTCCGGCGCTATCGCCCTGCTGGTTTATTTAGTCATCGCCGTTTCAAAGTTGCGGATGCGCAAAATTCTGCGAGCAGAAGGAAGCGAAATTCGCTTGCGAATGTGGCTTTATCCGt                                                                         <  1:108622/149‑1 (MQ=255)                                                                   gCGAAGGTATTTAAATTTCTGATCGACAGCTAAGCCGAAATCGCCCTGCTGGTTTATTTAGTCATCGCCGTTTCAAAGTTGCGGATGCGCAAAATTCTGCGAGCAGAAGGAAGCGAAATTCGCTTGCGAATGTGGCTTTATCCGTGGCt                                                                     >  1:188216/1‑149 (MQ=255)                                                                    cGAAGGTATTTAAATTTCTGATCGACAGCTCCGGCGCTATCGCCCTGCTGGTTTATTTAGTCATCGCCGTTTCAAAGTTGCGGCTGCGCAAAATTCTGCGAGCAGAAGGAAGCGAAATTCGCTTGCGAATGTGGCTTTCTCCGTGGCTc                                                                    >  2:101801/1‑149 (MQ=255)                                                                        ggTATTTAAATTTCTGATCGACAGCTCCGGCGCTATCGCCCTGCTGGTTTATTTAGTCATCGCCGTTTCAAAGTTGCGGATGCGCAAAATTCTGCGAGCAGAAGGAAGCGAAATTCGCTTGCGAATGTGGCTTTATCCGTGGCTCACCt                                                                <  2:425971/149‑1 (MQ=255)                                                                            tttAAATTTCTGATCGACAGCTCCGGCGCTAGCGCCCTGCTGGTTTATTTAGTCATCGCCGTTTCAAAGTTGCGGATGCGCAAAATTCTGCGAGCAGAAGGAAGCGAAATTAGCTTGCGAATGTGGCTTTATCCGTGGCTAACCTGGCt                                                            >  2:133031/1‑149 (MQ=255)                                                                             ttAAATTTCTGATCGACAGCTCCGGCGCTATCGCCCTGCTGGTTTATTTAGTCATCGCCGTTTCAAAGTTGCGGATGCGCAAAATTCTGCGAGCAGAAGGAAGCGAAATTCGCTTGCGAATGTGGCTTTATCCGTGGCTCACCTGGCTa                                                           >  2:190864/1‑149 (MQ=255)                                                                             ttAAATTTCTGATCGACAGCTCCGGCGCTATCGCCCTGCTGGTTTATTTAGTCATCGCCGTTTCAAAGTTGCGGATGCGCAAAATTCTGCGAGCAGAAGGAAGCGAAATTCGCTTGCGAATGTGGCTTTATCCGTGGCTCACCTGGCTa                                                           <  1:47709/149‑1 (MQ=255)                                                                             ttAAATTTCTGATCGACAGCTCCGGCGCTAGCGCCCTGCGGGTTTATTTAGTCATCGCCGTTTCAAAGGTGCGGATGCGCAAAATTCTGCGAGCAGAAGGAAGCGAAATTCGCTTGCGAATGTGGCTTTATCCGTGGCTCACCTGGCTa                                                           <  2:613198/149‑1 (MQ=255)                                                                              tAAATTCCTGATCGACAGCTCCGGCGCTATCGCCCTGCTGGTTTATTTAGTCATCGCCGTTTCAAAGTTGCGGATGCGCAAAATTCTGCGAGCAGAAGGAAGCGAAATTCGCTTGCGAATGTGGCTTTATCCGTGGCTCACCTGGCTAg                                                          <  1:402101/149‑1 (MQ=255)                                                                                 aTTTCTCAGCGACAGCTCCGACGCTATCGCAGTGATGGTTGATTTAGTCATCACAGTATCAAAGTGGCGGCTGCGCAAAATTCTGCGAGCAGAAGGAAGCGAAATTCGCTTGCGAATGTGGCTTTATCCGTGGCTCACCTGGCTAGTcc                                                       >  2:14699/1‑148 (MQ=255)                                                                                  tttCTGATCGACAGCTCCGGCGCTATCGCCCTGCTGGTTTATTTAGTCATCGCCGTTTCAAAGTTGCGGATGCGCAAAATTCTGCGAGCAGAAGGAAGCGAAATTCGCTTGCGAATGTGGCTTTATCCGTGGCTCCCCTGGCTAGTCac                                                      >  2:531958/1‑148 (MQ=255)                                                                                    tCTGATCGACAGCCCCGGGGCTATCGCCCTGCTGGTTTATTTAGTCATCGCCGTTTCAAAGTTGCGGATGCGCAAAATTCTGCGAGCAGAAGGAAGCGAAATTCGCTTGCGAATGTGGCTTTATCCGTGGCTCACCTGGCTAGTCATCg                                                    <  2:421750/149‑1 (MQ=255)                                                                                     cTGATCGACAGCTCCGGCGCTATCGCCCTGCTGGTTTATTTAGTCATCGCCGTTTCAAAGTTGCGGATGCGCAAAATTCTGCGAGCAga                                                                                                               >  2:456824/1‑89 (MQ=255)                                                                                     cTGATCGACAGCTCCGGCGCTATCGCCCTGCTGGTTTATTTAGTCATCGCCGTTTCAAAGTTGCGGATGCGCAAAATTCTGCGAGCAga                                                                                                               <  1:456824/89‑1 (MQ=255)                                                                                         tCGACAGCTCCGGCGCTATCGCCCTGCTGGTTTATTTAGTCATCGCCGTTTCAAAGTTGCGGATGCGCAAAATTCTGCGAGCAGAAGGAAGCGAAATTCGCTTGCGAATGTGGCTTTATCCGTGGCTCACCTGGCTAGTCATCGGCttt                                               >  1:392366/1‑149 (MQ=255)                                                                                                    ggCGCTATCGCCCTGCTGGTTTATTTAGTCATCGCCGTTTCAAAGTTGCGGATGCGCAAAATTCTGCGAGCAGAAGGAAGCGAAATTCGCTTGCGAATGTGGCTTTATCCGTGGCCCACCTGGCTCGTCACCGGCTTTATTACCTTtgt                                    <  2:188216/149‑1 (MQ=255)                                                                                                       gcTATCGCCCTGCTGGTTTATTTAGTCATCGCCGTTTCAAAGTTGCGGATGCGCAAAATTCTGCGAGCAGAAGGAAGCGAAATTCGCTTGCGAATGTGGCTTTATCCGTGGCTCACCTGGCTAGTCATCGGCTTTATTACCTTTGTGtt                                 >  1:128962/1‑149 (MQ=255)                                                                                                           tCGCCCTGCTGGTTTATTTAGTCATCGCCGTTTCAAAGTTGCGGATGCGCaaaa                                                                                                                            >  1:637569/1‑54 (MQ=255)                                                                                                           tCGCCCTGCTGGTTTATTTAGTCATCGCCGTTTCAAAGTTGCGGATGCGCaaaa                                                                                                                            <  2:637569/54‑1 (MQ=255)                                                                                                               cggtcgggTTTATTTAGTCGTCGCCGTTTCAAAGTTGCGGATGCGCAAAATTCTGCGAGCAGAAGGAAGCGAAATTCGCTTGCGAATGTGGCTTTATCCGTGGCTCACCTGGCTAGTCATCGGCTTTATTACCTTTGTGTTGGTAGTGa                         <  2:193987/143‑1 (MQ=255)                                                                                                                     ggTTTATTTAGTCATCGCCGTTTCAAAGTTGCGGATGCGCAAAATTCTGCGAGCAGAAGGAAGCGAAATTCGCTTGCGAATGTGGCTTTATCCGTGGCTCACCTGGCTAGTCATCGGCTTTATTACCTTTGTGTTGGTAGTGATGCTAt                   <  1:316533/149‑1 (MQ=255)                                                                                                                            ttaGTCATCGCCGTTTCAAAGTTGCGGATGCGCAAAATTCTGCGCGCAGAAGGAAGCGAAATTCGCTTGCGAATGTGGCTTGATCGGTGGCTCACCTGGCTAGTCATCGGGTTTATTACCTTTGTGTGGGTGTTGATGCTATtacgtcc            >  2:425227/1‑149 (MQ=255)                                                                                                                            ttaGTCATCGCCGTTTCAAAGTTGCGGATGCGCAAAATTCTGCGAGCAGAAGGAAGCGAAATTCGCTTGCGAATGTGGCTTTATCCGTGGCTCACCTGGCTAGTCATCGGCTTTATTACCTTTGTGTTGGTAGTGATGCTATtccgtcc            <  2:392366/149‑1 (MQ=255)                                                                                                                                tCATCGCCGTTTCAAAGTTGCGGATGCGCAAAATTCTGCGAGCAGAAGGAAGCGAAATTCGCTTGCGAATGTGGCTTTATCCGTGGCTCACATGGCTAGTCATCGGCTTTCTTCCCTTTGTGTTGGTAGTGATGCTATTCCGTCCGgcg        <  2:582158/149‑1 (MQ=255)                                                                                                                                   tCGCCGTTTCAAAGTTGCGGATGCGCAAAATTCTGCGAGCAGAAGGAAGCGAAATTCGCTTGCGAATGTGGCTTTATCCGTGGCTCACCTGGCTAGTCATCGGCTTTATTACCTTTGTGTTGGTAGTGATGCTATTCCGTCCGGCGCaa     <  1:589330/149‑1 (MQ=255)                                                                                                                                      ccGTTTCAAAGTTGCGGATGCGCAAAATTCTGCGAGCAGAAGGAAGCGAAATTCGCTTGCGAATGTGGCTTTATCCGTGGCTCACCTGGCTAGTCATCGGCTTTATTACCTTTGTGTTGGTAGTGATGCTATTCCGTCCGGCGCAACAg  >  1:330346/1‑149 (MQ=255)                                                                                                                                      ccGTTTCAAAGTTGCGGATGCGCAAAATTCTGCGAGCAGAAGGAAGCGAAATTCGCTTGCGAATGTGGCTTTATCCGTGGCTCACCTGGCTAGTCATCGGCTTTATTACCTTTGTGTTGGTAGTGATGCTATTCCGTACGGCGCAAAAg  >  1:585262/1‑149 (MQ=255)                                                                                                                                              |                                                                                                                                             TACGTGGCGGTGTTACTCTCCACCGGCGCGGCATTCTTAACGGTGGTGGTGAACTATTACGCGCCTGCGAAGGTATTTAAATTTCTGATCGACAGCTCCGGCGCTATCGCCCTGCTGGTTTATTTAGTCATCGCCGTTTCACAGTTGCGGATGCGCAAAATTCTGCGAGCAGAAGGAAGCGAAATTCGCTTGCGAATGTGGCTTTATCCGTGGCTCACCTGGCTAGTCATCGGCTTTATTACCTTTGTGTTGGTAGTGATGCTATTCCGTCCGGCGCAACAG  >  NC_000913_3_vch_tpiA/2795255‑2795536

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