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

CGTACGTGGCGGTGTTACTCTCCACCGGCGCGGCATTCTTAACGGTGGTGGTGAACTATTACGCGCCTGCGAAGGTATTTAAATTTCTGATCGACAGCTCCGGCGCTATCGCCCTGCTGGTTTATTTAGTCATCGCCGTTTCACAGTTGCGGATGCGCAAAATTCTGCGAGCAGAAGGAAGCGAAATTCGCTTGCGAATGTGGCTTTATCCGTGGCTCACCTGGCTAGTCATCGGCTTTATTACCTTTGTGTTGGTAGTGATGCTATTCCGTCCG  >  NC_000913_3_vch_tpiA/2795253‑2795527                                                                                                                                                |                                                                                                                                    cGTACGTGGCGGGGTTACTCTCCGCCGGCGCGGCATTCTTGACGGTGGTGGTGAACTATTACGCGCCTGCGAAGGTCATTTAATTTCTTCTGTACATCTCCGGCGCTATCGCCCTGCTGGTTTATTTAGTCATCGCCGTTTCAAAGTTg                                                                                                                                <  2:234024/149‑1 (MQ=255)           ggTGTTACTCTCCACCGGCGCGGCATTATTAACGGTGGTGGTGAACAATTACTCCGCTGCGCACGTATTTAAATTTCTGATCGACAGCTCCGGCGCTATCGCCCTGCTGGTTTATTTAGTCATCGCCGCTTCAAAGTTGCGAATGcaca                                                                                                                      >  2:84164/1‑149 (MQ=255)                  ctctCCACCGGCGCGGCATTCTTAACGGTGGTGGTGAACTATTACGCGCCTGCGAAGGTATTTAAATTTCTGATCGACAGCTCCGGCGCTATCGCCCTGCTGGTTTATTTAGTCATCGCCGTTTCAAAGTTGCGGATGCGCAAAATTCt                                                                                                               >  1:159647/1‑149 (MQ=255)                   tctcCACCGGCGCGGCATTCTTAACGGTGGTGGTGAACTATTACGCGCCTGCGAAGGTCCACACAAAACTGATAGCAATCGCCATCGCTATCGCCCTGCTGGTTTATTTATTCATCGCCGTTTCAAAGTTGCGGATGCGCAAAATTCTg                                                                                                              >  1:245767/1‑149 (MQ=255)                      ccACCGGCGCGGCATTCTTAACGGTGGTGGTGAACTATTACGCGCCTGCGACGGTATTTAAATTTCTGATCGACAGCTCCGGCGCTATCGCCCTGCTGGTTTATTTAGTCATCGCCGTTTCAAAGTTGCGGATGCGCAAAATTCTGCGa                                                                                                           <  1:243076/149‑1 (MQ=255)                              gcgGCATTCTTAACGGTGGTGGTGAACTATTACGCGCCTGCGAAGGTATTTAAATTTCTGATCGACAGCTCCGGCGCTATCGCCCTGCTGGTTTATTTAGTCATCGCCGTTTCAAAGTTGCGGATGCGCAAAATTCTGCGAGCAgaagg                                                                                                   <  2:481167/149‑1 (MQ=255)                                ggCATTCTTAACGGTGGTGGTGAACTATTACGCGCCTGCGAAGGTATTTAAATTTCTGATCGACAGCTCCGGCGCTATCGCCCTGCTGGTTTATTTAGTCATCGCCGTTTCAAAGTTGCGGATGCGCAAAATTCTGCGAGCAgaaggaa                                                                                                 >  1:229535/1‑149 (MQ=255)                                ggCATTCTTAACGGTGGTGGTGAACTATTACGCGCCTGCGAAGGTATTAAAATTTCTGATCGACAGCTCCGGCGCTATCGCCCTGCTGGTTTATTTAGTCATCGCCGTTTCAAAGTTGCGGATGCGCAAAATTCTGCGAGCAgaaggaa                                                                                                 <  2:245767/149‑1 (MQ=255)                                 gCATTCTTAACGGTGGTGGTGAACTATTACGCGCCTGCGAAGGTACTTGCACTCACGATCGACAGCTCCAGCGCTATCGCCCTGCTGGTTTATTTAGTCATCGCCGTTTCAAAGTTGCGCATGCACAAAATTCTGCGCGCAgaaggaag                                                                                                >  2:456804/1‑149 (MQ=255)                                       ttAACGGTGGTGGTGAACTATTACGCGCCTGCGAAGGTATTTAAATTTCTGATCGACAGCTCCGGCGCTATGGCCCTGCTGGTTTATTTAGTCATCGCCGTTTCCAAGTTGCGGATGCGCAAAATTCTGCGAACAGAAGGAAGCGACAt                                                                                          >  2:430875/1‑149 (MQ=255)                                                     gAACTATTACGCGCCTGCGAAGGTATTTAAATTTCTGATCGACAGCTCCGGCGCTATCGCCCTGCTGGTTTATTTAGTCATCGCCGTTTCAAAGTTGCGGATGCGCAAAATTCTGCGAGCAGAAGGAAGCGAAATTCGCTTGCGAAtgt                                                                            <  1:154654/149‑1 (MQ=255)                                                                 gcCTGAGCAGGTATTTAAATTTCTGATCGACAGCTCCGGCGCTATCGCCCTGCTGGTTTATTTAGTCATCGCCGCTTCAAAGTTGCGGATGCGCAAAATTCTGCGAGCAGAAGGAAGCGAAATTCGCTTGCGAATGTGGCTTTATCCGt                                                                <  1:64710/149‑1 (MQ=255)                                                                  ccTGCGAAGGCATTTAAATTTCTGATCGACAGCCCCGGCGCTATCGCCCTGCTGGTTTATTTAGTCATCGCCGTTTCAAAGTTGCGGATGCGCAAAATTCTGCGAGCAGAAGGAAGCGAAATTCGCTTGCGAATGTGGCTTTATCCGTg                                                               <  2:394619/149‑1 (MQ=255)                                                                   cTGCGAAGGTATTTAAATTTCTGATCGACAGCTCCGGCGCTATCGCCCTGCTGGTTTATTTAGTCATCGCCGTTTCAAAGTTGCGGATGCGCAAAATTCTGCGAGCAGAAGGAAGCGAAATTCGCTTGCGAATGTGGCTTTATCCGTgg                                                              >  2:468679/1‑149 (MQ=255)                                                                                        tGATCGACAGCTCCGGCGCTATCGCCCTGCTGGTTTATTTAGTCATCGCCGTTTCAAAGTTGCGGATGCGCAAAATTCTGCGAGCAGAAGGAAGCGAAATTCGCTTGCGAATGTGGCTTTATCCGTGGCTCACCTGGCTAGTCATCGGc                                         >  2:462165/1‑149 (MQ=255)                                                                                          aTCGACAGCTCCGGCGCTATCGCCCTGCTGGTTTATTTAGTCATCGCCGTTTCAAAGTTGCGGATGCGCAAAATTCTGCGAGCAGAAGGAAGCGAAATTCGCTTGCGAATGTGGCTTTATCCGTGGCTCACCTGGCTAGTCATCGGCtt                                       <  1:364120/149‑1 (MQ=255)                                                                                           tCGACAGCTCCGGCGCTATCGCCCTGCTGGTTTATTTAGTCATCGCCGTGTCAAAGTTGCGGATGCGCAAAATTCTGCGAGCAGAAGGAAGCGAAATTCGCTTGCGAATGTGGCTTTATCCGTGGCTCACCTGGCTAGTCATAGGCttt                                      >  1:343061/1‑149 (MQ=255)                                                                                               cAGCTCCGGCGCTATCGCCCTGCTGGTTTATTTAGTCATCGCCGTTTCAAAGTTGCGGATGCGCAAAATTCTGCGAGCAGAAGGAAGCGAAATTCGCTTGCGAATGTGGCTTTATCCGTGGCTCACCTGGCTAGTCATCGGCTttatta                                  <  1:84164/149‑1 (MQ=255)                                                                                                   tCCGGCGCTATCGCCCTGCTGGTTTATTTAGTCATCGCCGTTTCAAAGTTGCGGATGCGCAAAATTCTGCGAGCAGAAGGAAGCGAAATTCGCTTGCGAAtgt                                                                            >  1:245116/1‑103 (MQ=255)                                                                                                   tCCGGCGCTATCGCCCTGCTGGTTTATTTAGTCATCGCCGTTTCAAAGTTGCGGATGCGCAAAATTCTGCGAGCAGAAGGAAGCGAAATTCGCTTGCGAAtgt                                                                            <  2:245116/103‑1 (MQ=255)                                                                                                                cccTGCTGGTTTATTTAGTCATCGCCGTTTCAAAGTTGCGGATGCGCAAAATTCTGCGAGCAGAAGGAAGCGAAATTCGCTTGCGAATGTGGCTTTATCCGTGGCTCACCTGGCTAGTCATCGGCTTTATTACCTTTGTGTTGGTAGTg                 >  1:336982/1‑149 (MQ=255)                                                                                                                 gcgggtgGTTTATTTAGTCATCGCCGTTTCAAAGTTGCGGAGGCGCAAAATTCTGCGAGCAGAAGGAAGCGAAATTCGCTTGCGAATGTGGCTTTATCCGGGGCTCACCTGGCTAGTCATCGGCTTTATTACCTTTGTGTTGGTAGTGa                <  2:186219/144‑1 (MQ=255)                                                                                                                             tttaGTCATCGCCGTTTCAAAGTTGCGGATGCGCAAAATTCTGCGAGCAGAAGGAAGCGAAATTCGATTGCGAATGTGGCTTTATCCGTGGCTCACCTGGCTAGTCATCGGCTTTATTACCTTTGAGTTGGTGGGTAGGCCATtccgtc    >  2:99492/1‑149 (MQ=255)                                                                                                                               taGTCATCGCCGTTTCAAAGTTGCGGATGCGCAAAATTCTGCGAGCAGAAGGAAGCGAAATTCGCTTGCGAATGTGGCTTTATCCGTGGCTCACCTGGCTAGTCATCGGCTTTATTACCTTTGTGTTGGTAGTGATGCTATtccgtccg  >  1:375955/1‑149 (MQ=255)                                                                                                                                                |                                                                                                                                    CGTACGTGGCGGTGTTACTCTCCACCGGCGCGGCATTCTTAACGGTGGTGGTGAACTATTACGCGCCTGCGAAGGTATTTAAATTTCTGATCGACAGCTCCGGCGCTATCGCCCTGCTGGTTTATTTAGTCATCGCCGTTTCACAGTTGCGGATGCGCAAAATTCTGCGAGCAGAAGGAAGCGAAATTCGCTTGCGAATGTGGCTTTATCCGTGGCTCACCTGGCTAGTCATCGGCTTTATTACCTTTGTGTTGGTAGTGATGCTATTCCGTCCG  >  NC_000913_3_vch_tpiA/2795253‑2795527

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