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

CCACCGGCGCGGCATTCTTAACGGTGGTGGTGAACTATTACGCGCCTGCGAAGGTATTTAAATTTCTGATCGACAGCTCCGGCGCTATCGCCCTGCTGGTTTATTTAGTCATCGCCGTTTCACAGTTGCGGATGCGCAAAATTCTGCGAGCAGAAGGAAGCGAAATTCGCTTGCGAATGTGGCTTTATCCGTGGCTCACCTGGCTAGTCATCGGCTTTATTACCTT  >  NC_000913_3_pae_tpiA/2795274‑2795499                                                                                                                           |                                                                                                        ccgccGGCGCGGCATTCTTACCGGTGGTGGTGAACTATTACGCGCCTGCGAAGGTATTTACATTTCTGATCGACAGCTCCGGCGCTATCGCCCTGCTGGTTTATTTAGTCATCGCCGTTTCAAAGTTGCGGATgcg                                                                                            <  2:131026/133‑1 (MQ=255) ccACCGGCGCGGCATTTGTAACGGTGGTGGTGAACTATTACGCGCCTGCGAAGGTATTTAAATTTCTGATCGACAGCTCCGGCGCTATCGCCCTGCTGGTTTATTTAGTCATCGCCGTTTCAAAGTTGCGGATgcg                                                                                            <  2:131282/136‑1 (MQ=255)                     aCGGTGGTGGTGAACTATTACGCGCCTGCGAAGGTATTTAAATTTCTGATCGACAGCTCCGGCGCTATCGCCCTGCTGGTTTATTTAGTCATCGCCGTTTCAAAGTTGCGGATGCGCAAAATTCTGCGAGCAgaag                                                                        <  2:245821/136‑1 (MQ=255)                     aCGGTGGTGGTGAACTATTACACGCCTGCGAAGATATTTAAATTTCTGATCGACAGCTCCGGCGCTATCGCCCTGCTGGTTTATTTAGTCATCGCCGTTTCAAAGTTGCGGATGCGCAAAATTCTGCGAGCCgaag                                                                        >  1:154063/1‑136 (MQ=255)                        gtggtggtgAACTATTACGCGCCTGCGAAGGTATTTAAATTTCTGATCGACAGCTCCGGCGCTATCGCCATGCTGGTTTATTTAGTCATCGCCGTTTCAAAGTTGCGGATGCGCAAAATTCTACGCGCAgaaggaa                                                                     >  2:187653/1‑136 (MQ=255)                         tggtggtgAACTATTACGCGCCTGCGAAGGTATTTAAATTTCTGATCGACAGCTCCGGCGCTATCGCCCTGCTGGTTTATTTAGTCATCGCCGTTTCAAAGTTGCGGATGCGCAAAATTCTGCGAGCAgaaggaag                                                                    >  2:147725/1‑136 (MQ=255)                              gtgAACTATTACGCGCCTGCGAAGGAATTTAAATTTCTGATCGACAGCTCCGGCGCTATCGCCCTGCTGGTTTATTTAGTCATCGCCGTTTCAAAGTTGCGGATGCGCAAAATTCTGCGAGCAGAAGGAAGCGaaa                                                               >  1:251884/1‑136 (MQ=255)                                   gtATTACCCGCCGTCGAAGGTATGTAAATTTCTGATCGACAGCTCCGGCGCTATCGCCCTGCTGGTTTATTTAGTCATCGCCGTTTCAAAGTTGCGGCTGCGCAAAATTCTGCGCGCAGAAGGAAGCGAAATTCGc                                                          <  2:247921/135‑1 (MQ=255)                                      ttACGCGCCTGCGAAAGTAATTAAATTTCTGAACGACAACTCCGGCGCTATCGCCCTGCTGGTTTATTTAGTCATCGCCGTTTCAAAGTTGCGGATGCGCAAAATTCTGTGAGCAGAAGGAAGCGAAATTCGCTTg                                                       >  1:236720/1‑136 (MQ=255)                                                             aaTTTCTGATCGACAGCTCCGGCGCTATCGCCCTGCTGGTTTATTTAGTCATCGCCGTTTCAAAGTTGCGGATGCGCAAAATTCTGCGAGCAGAAGGAAGCGAAATTCGCTTGCGAATGTGGCTTTATCCGTGGCt                                >  1:214498/1‑136 (MQ=255)                                                                       cGACCGCTCCGGCGCTATCGCCCTGCTGGTTTATTTAGTCATCGCCGTTTCAAAGTTGGGAATGCGAAAAATTCTGCGAGCAGAAGGAAGCGAAAATCGCTTGCGAATGTGGCTTTATACGTGGCTCACCTGGCTa                      >  1:86087/1‑136 (MQ=255)                                                                                           cccTGCTGGTTTATTTAGGCATAGCCGTTGCAAAGTTGCGGATGCGCAAAATTCTGCGAGCAGAAGGAAGCGAAATTCGCTTGCGAATGTGGCTTTATTCGTGGCTTATCTGGCTAGTCATCGGCTTTATTCCCtt  <  2:123830/136‑1 (MQ=255)                                                                                               gctgGTTTATTTAGTCATCGCCGTTTCAAAGTTGCGGATGCGCAAAATTCTGCGAGCAGAAGGAAGCGAAATTCGCTTGCGAATGTGGCTTTATCCGTGGCTCACCTGGCTAGTCATCGGCTTTATTACCt   >  1:84256/1‑131 (MQ=255)                                                                                               gctgGTTTATTTAGTCATCGCCGTTTCAAAGTTGCGGATGCGCAAAATTCTGCGAGCAGAAGGAAGCGAAATTCGCTTGCGAATGTGGCTTTATCCGTGGCTCACCTGGCTAGTCATCGGCTTTATTACCt   <  2:84256/131‑1 (MQ=255)                                                                                                                           |                                                                                                        CCACCGGCGCGGCATTCTTAACGGTGGTGGTGAACTATTACGCGCCTGCGAAGGTATTTAAATTTCTGATCGACAGCTCCGGCGCTATCGCCCTGCTGGTTTATTTAGTCATCGCCGTTTCACAGTTGCGGATGCGCAAAATTCTGCGAGCAGAAGGAAGCGAAATTCGCTTGCGAATGTGGCTTTATCCGTGGCTCACCTGGCTAGTCATCGGCTTTATTACCTT  >  NC_000913_3_pae_tpiA/2795274‑2795499

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