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

CCCCGTACGTGGCGGTGTTACTCTCCACCGGCGCGGCATTCTTAACGGTGGTGGTGAACTATTACGCGCCTGCGAAGGTATTTAAATTTCTGATCGACAGCTCCGGCGCTATCGCCCTGCTGGTTTATTTAGTCATCGCCGTTTCACAGTTGCGGATGCGCAAAATTCTGCGAGCAGAAGGAAGCGAAATTCGCTTGCGAATGTGGCTTTATCCGTGGCTCACCTGGCTAGTCATCGGCTTTATTACCTTTGTGTTGGTAGTGATGCTATTCCGTCCGGCGCAACAGT  >  NC_000913_3_vch_tpiA/2795250‑2795537                                                                                                                                                   |                                                                                                                                              ccccGTACGTGGCGGTGTTACTCTCCACCGGCGCGGCATTCTTAACGGTGGTGGTGAACTATTACGCGCCTGCGAAGGTATTTAAATTTCTGATCGACAGCTCCGGCGCTATCGCCCTGCTGGTTTATTTAGTCATCGCCGTTTCAAAg                                                                                                                                             <  2:19399/149‑1 (MQ=255)                 gtTACTCTCCACCGGCGCGGCATTCTTAACGGTGGTGGTGAACTATTACGCGCCTGCGAAGGTATTTAAATTTCTGATCGACAGCTCCGGTGCTATCGCCCTGCTGGTTTATTTAGTCATCGCCGTTTCAAAGTTGCGGATGCGCaaaa                                                                                                                             <  1:41024/149‑1 (MQ=255)                         ccACCGGCTCGGGATTCTTAACGGGGGGGGTGAACTATTACGCGCCTTGGAAGGTATTTAAATTTCTGATCGACATCTCCGTCGCTATCTCCCTGCTTGTTTATTTAGTCATCGCCGTTTCAAAGTTGCGGATGCGCAAAATTCTGCGa                                                                                                                     <  1:554931/149‑1 (MQ=255)                          cACCGGCGCGGGATTCTTAACGGTGGTGGTGAACTATTACGCGCCTGCGAAGGTATTTAAATTTCTGATCGACAGCTCCGGCGCTATCGCCCTGCTGGTTTATTTAGTCATCGCCGTTTCAAAGTTGCGGATGCGCAAAATTCTGCGAg                                                                                                                    >  1:340675/1‑149 (MQ=255)                          cACCGGCGCGGCATTCTTAACGGTGGTGGTGAACTATTACGCGCGTGCGAAGGTGTTTAAATTTCTGATCTACAGCTCCGGCGCTATCGCCCTGCTGGTTTATTTAGTCATCGCCGTTTCAAAGTCGCGGATGCGCAAAATTCTGCGAg                                                                                                                    >  1:484999/1‑149 (MQ=255)                               gcgcTGCATTCTTAACGGTGGTGGTTAACTATTACGCGCCTGCGAAGGTATTTAAATTTCTGATCGACAGCTCCGGCGCTATCGCCCTGCTGGTTTATTTAGTCATCGCCGTTTCAAAGTTGCGGATGCGCAAAATTCTGCGAGCAgaa                                                                                                               <  1:301216/149‑1 (MQ=255)                                  cgGCATTCTTAACGGTGGTGGTGAACTATTACGCGCCTGCGAAGGTATTTAAATTTCTGATCGAGAGCTCCGGCGCTATCGCCCTACTGGATTAAATAGTAATCCCCGTTTCAAAGTTGCGGATACGCAAAATTATGCGGGCaggaggc                                                                                                            >  2:162834/1‑148 (MQ=255)                                  cgGCATTCTTAACGGTGGTGGTGAACTATTACGCGCCTGCGAAGGTATTTAAATTTCTGATCGACAGCTCCGGCGCTATCGCCCTGCTGGTTTATTTAGTCATCGCCGTTTCAAAGTTGCGGATGCGCAAAATTCTGCGAGCAgaagga                                                                                                            >  2:584642/1‑149 (MQ=255)                                       ttCTTAACGGTGGTGGTGAACTATTACGCGCCTGCGAAGGTATTTAAATTTCTGATCGACAGCTCCGGCGCTATCGCCCTGCTGGTTTATTTAGTCATCGCCGTTTCAAAGTTGCGGATGCGCAAAATTCTGCGAGCAGAAGGAAGCGa                                                                                                       >  2:301754/1‑149 (MQ=255)                                                gtggtggtgAACTATTACGCGCCTGCGAAGGTATTTAAATTTCTGATCGACAGCTCCGGCGCTATCGCCCTGCTGGTTTATTTAGTCATCGCCGTTTCAAAGTTGCGGATGCGCAAAATTCTGCGAGCAGAAGGAAGCGAAATTCGCtt                                                                                              >  1:298573/1‑149 (MQ=255)                                                gtggtggtgAACGATTACGCGCCTCCGAAGGTATTTAAATTTCTGATCGACAGCTCCGGCGCTATCGCCCTGCTGGTTTATTTAGTCATCGCCGTTTCAAAGTTGCGGATGCGCAAAATTCTGCGAGCAGAAGGAAGCGAAATTCGCtt                                                                                              <  2:208770/149‑1 (MQ=255)                                                      gtgAACTATTACGCGCCTGCGAAGGTATTTAAATTTCTGATCGACAGCTCCGGCGCTATCGCCCTGCTGGTTTATTTAGTCATCGCCGTTTAAAAGTTGGGGATGCGCAAAGTTCTGCGAGCAGAAGTAAGCAAAATTCGCTTTCGAAt                                                                                        >  2:292026/1‑149 (MQ=255)                                                               tACGCGCCTGCGAAGGTATTTAAATTTCTGATCGACAGCTCCGGCGCTATCGCCCTGCTGGTTTATTTAGTCATCGCCGTTTCAAAGTTGCGGATGCGCAAAATTCTGCGAGCAGAAGGAAGCGAAATTCGCTTGCGAATGTGGCTTTa                                                                               <  1:584642/149‑1 (MQ=255)                                                                tcgcgCTTGCGGAGGTATTTAAATTTCTGTTCGACAGCTCCTGCGCGATCGTCCTGTTGGTTTATTTAGTCATCGCCGTTTCAAAGTTGCGGATGCGCAAAATTCTGAGAGCAGAATGAAGCGAAATTCGCTTGCGAATGTGGCTTTAt                                                                              <  2:547433/148‑1 (MQ=255)                                                                 cgcgCCTGCGAAGGTATTTAAATTTCTGATCGACAGCTCCGGCGCTATCGCCCAGCTGGTTTATTTAGTCATCGCCGTTTCAAAGTTGCGGATGCGCAAAATTCTGCGAGCAGAAGGAAGCGAAATTCGCTTGCGAATGTGGCTTTATc                                                                             >  1:364909/1‑149 (MQ=255)                                                                      cTGCGAAGGTATTTAAATTTCTGATCGACAGCTCCGGCGCTATCGCCCTGCTGGTTTATTTAGTCATCGCCGTTTCAAAGTTGCGGATGCGCAAAATTCTGCGAGCAGAAGGAAGCGAAATTCGCTTGCGa                                                                                          >  1:87540/1‑131 (MQ=255)                                                                      cTGCGAAGGTATTTAAATTTCTGATCGACAGCTCCGGCGCTATCGCCCTGCTGGTTTATTTAGTCATCGCCGTTTCAAAGTTGCGGATGCGCAAAATTCTGCGAGCAGAAGGAAGCGAAATTCGCTTGCGa                                                                                          <  2:87540/131‑1 (MQ=255)                                                                      cTGCGAAGCTATTTAAATTTCTGATCGACAGCTCCGGCGCTATCGCCCTGCTGGTTTATTTAGTCATCGCCGTTTCAAAGTTGCGGATGCGCAAAATTCTGCGAGTAGAAGGAAGCGAAATTCGCTTGCGAATGTGGCTTTATCCGTgg                                                                        <  1:510189/149‑1 (MQ=255)                                                                              gTATTTAAATTTCTGATCGACAGCTCCGGCGCTATCGCCCTTCTGGTTTATTTAGTCATCGCCGTTTCAAAGTTGCTGATGCGCAAAATTCTTCGAGCAGAAGGAAGCGAAATTCGCTTGCGAATGTGGCTTTATCCGTGGCTCACCTg                                                                <  2:397857/149‑1 (MQ=255)                                                                                 tttAAATTTCTGATCGACAGCTCCGGCGCTATCGCCCTGCCGGTTTATTTAGTCATCGCCGTTTCAAAGTTGCGGATGCGCAAACCTCTGCGAGCCGAATGAATCGAAATTCGCTTGCGAATGTGGCTTTATCCGTGGCTCACCTGGCt                                                             <  2:274784/149‑1 (MQ=255)                                                                                             gtCGACAGCTCCGGCGCTATCGCCCTGCTGGTTTATTTAGTCATCGCCGTTTCAAAGTTGCGGATGCGCAAAATTCTGCGAGCAGAAGGAAGCGAAATTCGCTTGCGAATGTGGCTTTATCCGTGGCTCACCTGGCTCGTCATCGGCtc                                                 >  2:259770/2‑148 (MQ=255)                                                                                             aTCGACAGCTCCGGCGCTATCGCCCTGCTGGTTTATTTAGTCATCGCCGTTTCAAAGTTGCGGATGCGCAAAATTCTGCGAGCAGAAGGAAGCGAAATTCGCTTGCGAATGTGGc                                                                                   <  2:130651/115‑1 (MQ=255)                                                                                             aTCGACAGCTCCGGCGCTATCGCCCTGCTGGTTTATTTAGTCATCGCCGTTTCAAAGTTGCGGATGCGCAAAATTCTGCGAGCAGAAGGAAGCGAAATTCGCTTGCGAATGTGGc                                                                                   >  1:130651/1‑115 (MQ=255)                                                                                               cGACAGCTCCGGCGCTATCGCCCTGCTGGTTTATTTAGTCATCGCCGTTTCAAAGTTGCGGATGCGCAAAATTCTGCGAGCAGAAGGAAGCGAAATTCGCTTGCGAATGTGGCTTTATCCGTGGCTCACCTGGCTAGTCATCGGCTtta                                               <  1:301754/149‑1 (MQ=255)                                                                                                    gCTCCGGCGCTATCGCCCTGCTGGTTTATTTAGTCATCGCCGTTTCAAAGTTGCGGATGCGCAAAATTCTGCGAGCAGAAGGAAGCGAAATTCGCTTGCGAATGTGGCTTTATCCGTGGCTCACCTGGCTAGTCATCGGCTTTATTAcc                                          <  1:289935/149‑1 (MQ=255)                                                                                                         ggCGCTATCGCCCTGCTGGTTTATTTAGTCATCGCCGTTTCAAAGTTGCGGATGCACAAAATTCTGCGAGCAGAAGGAAGCGAAATTCGCTTGCGAATGTGGCTTTATCCGTGGCTCACCTGGCTAGTCATCGGCTTTATTACCTTtgt                                     >  1:337960/1‑149 (MQ=255)                                                                                                          gcgcTATCGCCCTGCTGGTTTATTTAGTCATCGCCGTTTCAAAGTTGCGGATACGCAAAATTCTGCGAGCAGAAGGAAGCGAAATTCGCTTGCGAATGTGGCTTTATCCGTGGCTCACCTGGCTAGTCATCGGCTTTATTACCTTtgtg                                    >  1:154165/1‑149 (MQ=255)                                                                                                              tATCGCCCTGCTGGTTCACTAAGTCAATGCCGATTCAAAGTTGCGGATGCGCAAAATTCTGCGAGCAGAAGGAAGCGAAATTCGCTTGCGAATGTGGCTTTATCCGTGGCTCACCTGGCTAGTCAACGGCTTTATTACCTTTGTGTTgg                                <  1:317596/149‑1 (MQ=255)                                                                                                                     ctgctgGTTTATTTAGTCATCGCCGTTTCAAAGTTGCGGATGCGCAAAATACTGCGAGCAGAAAGAAGCGAAATTCGCATGCGAATGTGGCTTTATCCGTGGCTCACCTGGCTAGTCATCGGCTTTATTACCTTTGTGTTTGTAGTGAt                         >  1:262715/1‑149 (MQ=255)                                                                                                                         tgGTTTATTTAGTCATCGCCGTTTCAACGTTGCGGATGCGCAAAATTCTGCAAGCAGAAGGAAGCGAAATTCGCTTGCGAATGTGGCTTTATCCGTGGCTCACCTGGCTAGTCATCGGCTTTATTACCTTTGTGTTGGTAGTGATGCTa                     >  1:564971/1‑149 (MQ=255)                                                                                                                           gTTTATTTAGTCATCGCCGTTTCAAAGTTGCGGATGCGCAAAATTCAGCGAGCAGAAAGAAGCGAAATTCGCTTGCGAATGTGGCGTTATCCGCGGCTCACCTGGCTAGTCATCAGCTTTATTACCTTTGTGTTGGTAGTGATGCTAtt                   >  1:514292/1‑149 (MQ=255)                                                                                                                                 ttaGTCATCGCCGTTTCAAAGTTGCGGATGCGCAAAATTCTGCGAGCAGAAGGAAGCGAAATTCGCTTGCGAATGTGGCTTTATCCGTGGCTCACCTGGCTAGTCATCGGCTttatta                                            >  2:324307/1‑118 (MQ=255)                                                                                                                                 ttaGTCATCGCCGTTTCAAAGTTGCGGATGCGCAAAATTCTGCGAGCAGAAGGAAGCGAAATTCGCTTGCGAATGTGGCTTTATCCGTGGCTCACCTGGCTAGTCATCGGCTttatta                                            <  1:324307/118‑1 (MQ=255)                                                                                                                                       aTCGCCGTTTCAAAGTTGCGGATGCGCAAAATTCTGCGAGCAga                                                                                                                <  2:22773/44‑1 (MQ=255)                                                                                                                                       aTCGCCGTTTCAAAGTTGCGGATGCGCAAAATTCTGCGAGCAga                                                                                                                >  1:22773/1‑44 (MQ=255)                                                                                                                                           ccGTTTCAAAGTTGCGGATGCGCAAAATTCTGCGAGCAGAAGGAAGCGAAATTCGCTTGCGAATGTGGCTTTATCCGTGGCTCACCTGGCTAGTCATCGGCTTTATTACCTTTGTGTTGGTAGTGATGCTATTCCGTCCGGCGCAACAg   >  2:350842/1‑149 (MQ=255)                                                                                                                                            cGTTTCAAAGTTGCGGATGCGCAAAATTCTGCGAGCAGAAGGAAGCGAAATTCGCTTGCGAATGTGGCTTTATCCGTGGCTCACCTGGCTAGTCATCGGCTTTATTACCTTTGTGTTGGTAGTGATGCTATTCCGTCCGGCGCAACAGt  <  2:154165/149‑1 (MQ=255)                                                                                                                                                   |                                                                                                                                              CCCCGTACGTGGCGGTGTTACTCTCCACCGGCGCGGCATTCTTAACGGTGGTGGTGAACTATTACGCGCCTGCGAAGGTATTTAAATTTCTGATCGACAGCTCCGGCGCTATCGCCCTGCTGGTTTATTTAGTCATCGCCGTTTCACAGTTGCGGATGCGCAAAATTCTGCGAGCAGAAGGAAGCGAAATTCGCTTGCGAATGTGGCTTTATCCGTGGCTCACCTGGCTAGTCATCGGCTTTATTACCTTTGTGTTGGTAGTGATGCTATTCCGTCCGGCGCAACAGT  >  NC_000913_3_vch_tpiA/2795250‑2795537

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