Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,425,505	Δ1 bp	coding (412/687 nt)	hisQ ←	lysine/arginine/ornithine ABC transporter/histidine ABC transporter, membrane subunit HisQ

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,425,503	0	A	.	100.0%	107.6 / NA	26	coding (414/687 nt)	hisQ	lysine/arginine/ornithine ABC transporter/histidine ABC transporter, membrane subunit HisQ
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base . (11/15);  total (11/15)

CCTGCCAGTTGGGTGGCTTTGACCACATCTTCCAGGCCGAGTAACGAAACCAGTGCGGTAGATTTGAGGATCACCTGCCAGTTGTTGCCAATGCCTGGCAGGGCGTAACGCATCATCGACGGAAACATGATCCGCCGAAACACTTGCCCACGAGTAAAACCGAACGCCGTCGCCGCCTCTATATGTCCTTTCGGCACTGCCATAAAAGCACCACGGAACGTTTCGGTAAAATAAGCACCGTAGATAAAACCGAGA  >  NC_000913/2425366‑2425620                                                                                                                                          |                                                                                                                      cctgccAGTTGGGTGGCTTTGACCACATCTTCCAGGCCGAGTAACGAAACCAGTGCGGTAGATTTGAGGATCACCTGCCAGTTGTTGCCAATGCCTGGCAGGGCGTAACGCATCATCGACGGAAACATGATCCGCCGaac                                                                                                                     >  2:48043/1‑139 (MQ=255)         ttGGGTGGCTTTGACCACATCTTCCAGGCCGAGTAACGAAACCAGTGCGGTAGATTTGAGGATCACCTGCCAGTTGTTGCCAATGCCTGGCAGGGCGTAACGCATCATCGACGGAAACATGATCCGCCG‑AACACTTGccc                                                                                                            <  2:87467/140‑1 (MQ=255)                           aTCTTCCAGGCCGAGTAACGAAACCAGTGCGGTAGATTTGAGGATCACCTGCCAGTTGTTGCCAATGCCTGGCAGGGCGTAACGCATCATCGACGGAAACATGATCCGCCG‑AACACTTGCCCACGAGTa                                                                                                     <  2:349166/129‑1 (MQ=255)                           aTCTTCCAGGCCGAGTAACGAAACCAGTGCGGTAGATTTGAGGATCACCTGCCAGTTGTTGCCAATGCCTGGCAGGGCGTAACGCATCATCGACGGAAACATGATCCGCCG‑AACACTTGCCCACGAGTa                                                                                                     >  1:349166/1‑129 (MQ=255)                             cTTCCAGGCCGAGTAACGAAACCAGTGCGGTAGATTTGAGGATCACCTGCCAGTTGTTGCCAATGCCTGGCAGGGCGTAACGCATCATCGACGGAAACATGATCCGCCG‑AACACTTGCCCACGAGTAAAACCGAACGCCg                                                                                        >  2:34787/1‑140 (MQ=255)                                ccAGGCCGAGTAACGAAACCAGTGCGGTAGATTTGAGGATCACCTGCCAGTTGTTGCCAATGCCTGGCAGGGCGTAACGCATCATCGACGGAAACATGATCCGCCG‑AACACTTGCCCACGAGTAAAACCGaa                                                                                             >  1:156522/1‑132 (MQ=255)                                ccAGGCCGAGTAACGAAACCAGTGCGGTAGATTTGAGGATCACCTGCCAGTTGTTGCCAATGCCTGGCAGGGCGTAACGCATCATCGACGGAAACATGATCCGCCG‑AACACTTGCCCACGAGTAAAACCGaa                                                                                             <  2:156522/132‑1 (MQ=255)                                    gCCGAGTAACGAAACCAGTGCGGTAGATTTGAGGATCACCTGCCAGTTGTTGCCAATGCCTGGCAGGGCGTAACGCATCATCGACGGAAACATGATCCGCCG‑AACACTTGCCCACGAGTAAAACCGAACGCCGTcgccgc                                                                                 <  1:235648/140‑1 (MQ=255)                                     ccGAGTAACGAAACCAGTGCGGTAGATTTGAGGATCACCTGCCAGTTGTTGCCAATGCCTGGCAGGGCGTAACGCATCATCGACGGAAACATGATCCGCCG‑AACACTTGCCCACGAGTAAAACCGAACGCCGTcgc                                                                                    <  1:124273/136‑1 (MQ=255)                                     ccGAGTAACGAAACCAGTGCGGTAGATTTGAGGATCACCTGCCAGTTGTTGCCAATGCCTGGCAGGGCGTAACGCATCATCGACGGAAACATGATCCGCCG‑AACACTTGCCCACGAGTAAAACCGAACGCCGTcgc                                                                                    >  2:124273/1‑136 (MQ=255)                                                    aGTGCGGTAGATTTGAGGATCACCTGCCAGTTGTTGCCAATGCCTGGCAGGGCGTAACGCATCATCGACGGAAACATGATCCGCCG‑AACACTTGCCCACGAGTAAAACCGAACGCCGTCGCCGCCTCTATATGTCCTTTc                                                                 >  1:92552/1‑140 (MQ=255)                                                      tGCGGTAGATTTGAGGATCACCTGCCAGTTGTTGCCAATGCCTGGCAGGGCGTAACGCATCATCGACGGAAACATGATCCGCCG‑AACACTTGCCCACGAGTAAAACCGAACGCCGTCGCCGCCTCTatat                                                                         <  1:36424/130‑1 (MQ=255)                                                      tGCGGTAGATTTGAGGATCACCTGCCAGTTGTTGCCAATGCCTGGCAGGGCGTAACGCATCATCGACGGAAACATGATCCGCCG‑AACACTTGCCCACGAGTAAAACCGAACGCCGTCGCCGCCTCTatat                                                                         >  2:36424/1‑130 (MQ=255)                                                            aGATTTGAGGATCACCTGCCAGTTGTTGCCAATGCCTGGCAGGGCGTAACGCATCATCGACGGAAACATGATCCGCCG‑AACACTTGCCCACGAGTAAAACCGAACGCCGTCGCCGCCTCTATATGTCCTTTCGGCACTGc                                                         <  1:62621/140‑1 (MQ=255)                                                               tttGAGGATCACCTGCCAGTTGTTGCCAATGCCTGGCAGGGCGTAACGCATCATCGACGGAAACATGATCCGCCG‑AACACTTGCCCACGAGTAAAACCGAACGCCGTCGCCGCCTCTATATGTCCTTTCGGCACTGCCAt                                                      <  1:48043/140‑1 (MQ=255)                                                                         aCCTGCCAGTTGTTGCCAATGCCTGGCAGGGCGTAACGCATCATCGACGGAAACATGATCCGCCG‑AACACTTGCCCACGAGTAAAACCGAACGCCGTCGCCGCCTCTATATGTCCTTTCGGCACTGCCATAAAAGcacca                                            >  2:134945/1‑140 (MQ=255)                                                                                       gCCAATGCCTGGCAGGGCGTAACGCATCATCGACGGAAACATGATCCGCCG‑AACACTTGCCCACGAGTAAAACCGAACGCCGTCGCCGCCTCTATATGTCCTTTCGGCACTGCCATAAAAGCACCACGGAACGTTTCGGt                              <  2:92552/140‑1 (MQ=255)                                                                                               cTGGCAGGGCGTAACGCATCATCGACGGAAACATGATCCGCCG‑AACACTTGCCCACGAGTAAAACCGAACGCCGTCGCCGCCTCTATATGTCCTTTCGGCACTGCCATAAAAGCACCACGGAACGTTTCGGTAAAATAAg                      <  1:134945/140‑1 (MQ=255)                                                                                                   cAGGGCGTAACGCATCATCGACGGAAACATGATCCGCCG‑AACACTTGCCCACGAGTAAAACCGAACGCCGTCGCCGCCTCTATATGTCCTTTCGGCACTGCCATAAAAGCACCACGGAACGTTTCGGTAAAATAAGCAcc                  <  1:34787/140‑1 (MQ=255)                                                                                                     gggCGTAACGCATCATCGACGGAAACATGATCCGCCG‑AACACTTGCCCACGAGTAAAACCGAACGCCGTCGCCGCCTCTATATGTCCTTTCGGCACTGCCATAAAAGCACCACGGAACGTTTCGGTAAAATAAGCACCGt                <  1:27422/140‑1 (MQ=255)                                                                                                        cGTAACGCATCATCGACGGAAACATGATCCGCCG‑AACACTTGCCCACGAGTAAAACCGAACGCCGTCGCCGCCTCTATATGTCCTTTCGGCACTGCCATAAAAGCACCACGGAACGTTTCGGTAAAATAAGCACCGTAGa             >  1:297745/1‑140 (MQ=255)                                                                                                         gTAACGCATCATCGACGGAAACATGATCCGCCG‑AACACTTGCCCACGAGTAAAACCGAACGCCGTCGCCGCCTCTATATGTCCTTTCGGCACTGCCATAAAAGCACCACGGAACGTTTCGGTAAAATAAGCACCGTAGAt            >  2:136846/1‑140 (MQ=255)                                                                                                             cGCATCATCGACGGAAACATGATCCGCCG‑AACACTTGCCCACGAGTAAAACCGAACGCCGTCGCCGCCTCTATATGTCCTTTCGGCACTGCCATAAAAGCACCACGGAACGTTTCGGTAAAATAAGCACCGTAGATaaaa        >  2:164934/1‑140 (MQ=255)                                                                                                                  catcGACGGAAACATGATCCGCCG‑AACACTTGCCCACGAGTAAAACCGAACGCCGTCGCCGCCTCTATATGTCCTTTCGGCACTGCCATAAAAGCACCACGGAACGTTTCGGTAAAATAAGCACCGTAGATAAAACCgag   <  1:293679/140‑1 (MQ=255)                                                                                                                   atcGACGGAAACATGATCCGCCG‑AACACTTGCCCACGAGTAAAACCGAACGCCGTCGCCGCCTCTATATGTCCTTTCGGCACTGCCATAAAAGCACCACGGAACGTTTCGGTAAAATAAGCACCGTAGATAAAACCgaga  <  2:297745/140‑1 (MQ=255)                                                                                                                          gAAACATGATCCGCCG‑AACACTTGCCCACGAGTAAAACCGAACGCCGTCGCCGCCTCTATATGTCCTTTCGGCACTGCCATAAAAGCACCACGGAACGTTTCGGTAAAATAAGCACCGTAGATAAAACCgaga  >  1:111317/1‑133 (MQ=255)                                                                                                                          gAAACATGATCCGCCG‑AACACTTGCCCACGAGTAAAACCGAACGCCGTCGCCGCCTCTATATGTCCTTTCGGCACTGCCATAAAAGCACCACGGAACGTTTCGGTAAAATAAGCACCGTAGATAAAACCgaga  <  2:111317/133‑1 (MQ=255)                                                                                                                                          |                                                                                                                      CCTGCCAGTTGGGTGGCTTTGACCACATCTTCCAGGCCGAGTAACGAAACCAGTGCGGTAGATTTGAGGATCACCTGCCAGTTGTTGCCAATGCCTGGCAGGGCGTAACGCATCATCGACGGAAACATGATCCGCCGAAACACTTGCCCACGAGTAAAACCGAACGCCGTCGCCGCCTCTATATGTCCTTTCGGCACTGCCATAAAAGCACCACGGAACGTTTCGGTAAAATAAGCACCGTAGATAAAACCGAGA  >  NC_000913/2425366‑2425620

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