Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,165,858	T→A	V64D (GTC→GAC)	yegQ →	putative peptidase YegQ

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,165,858	0	T	A	100.0%	101.7 / NA	30	V64D (GTC→GAC)	yegQ	putative peptidase YegQ
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base A (15/15);  total (15/15)

CGCAACAACGAATTCAACCACGAAAATCTTCAGCTCGGCATCAATGAAGCCCACGCGCTGGGGAAAAAGTTTTATGTCGTGGTCAACATTGCACCGCACAACGCCAAGCTGAAAACCTTTATCCGTGACCTGAAACCGGTGGTGGAAATGGGGC  >  NC_000913/2165782‑2165935                                                                             |                                                                              cacaacatCGAATTCAACCACGAAAATCATCAGCTCGGGATCAATGAAGCCCACGCGCTGGGGAAAAAGTTTTATGACgg                                                                            <  2:368459/78‑2 (MQ=255)   caacaaCGAATTCAACCACGAAAATCTTCAGCTCGGCATCAATGAAGCCCACGCGCTGGGGAAAAAGTTTTATGACGTGGt                                                                         >  2:368495/1‑81 (MQ=255)       aacGAATTCAACCACGAAAATCTTCAGCTCGGCATCAATGAAGCCCACGCGCTGGGGAAAAAGTTTTATGACGTGGTCaa                                                                      >  2:368507/1‑80 (MQ=255)       aacGAATTCAACCACGAAAATCTTCAGCTCGGCATCAATGAAGCCCACGCGCTGGGGAAAAAGTTTTATGACGTGGTCaa                                                                      >  1:368512/1‑80 (MQ=255)           aaTTCAACCACGAAAATCTTCAGCTCGGCATCAATGAAGCCCACGCGCTGGGGAAAAAGTTTTATGACGTGGTCAACATTg                                                                 <  1:368460/81‑1 (MQ=255)             ttCAACCACGAAAATCTTCAGCTCGGCATCAATGAAGCCCACGCGCTGGGGAAAAAGTTTTATGACGTGGTCAACATTGCa                                                               <  1:368461/81‑1 (MQ=255)               cAACCACGAAAATCTTCAGCTCGGCATCAATGAAGCCCACGCGCTGGGGAAAAAGTTTTATGACGTGGTCAACATTGCAc                                                              <  1:368462/80‑1 (MQ=255)                aaCCACGAAAATCTTCAGCTCGGCATCAATGAAGCCCACGCGCTGGGGAAAAAGTTTTATGACGTGGTCAACATTGCACCg                                                            >  1:368500/1‑81 (MQ=255)                  ccACGAAAATCTTCAGCTCGGCATCAATGAAGCCCACGCGCTGGGGAAAAAGTTTTATGACGTGGTCAACATTGCACCGca                                                          <  1:368463/81‑1 (MQ=255)                     cGAAAATCTTCAGCTCGGCATCAATGAAGCCCACGCGCTGGGGAAAAAGTTTTATGACGTGGTCAACATTGCACCGCACaa                                                       <  2:368464/81‑1 (MQ=255)                      gAAAATCTTCAGCTCGGCATCAATGAAGCCCACGCGCTGGGGAAAAAGTTTTATGACGTGGTCAACATTGCACCGCACAAc                                                      >  2:368503/1‑81 (MQ=255)                      gAAAATCTTCAGCTCGGCATCAATGAAGCCCACGCGCTGGGGAAAAAGTTTTATGACGTGGTCAACATTGCACCGCACAAc                                                      >  1:368509/1‑81 (MQ=255)                               cAGCTCGGCATCAATGAAGCCCACGCGCTGGGGAAAAAGTTTTATGACGTGGTCAACATTGCACCGCACAACGCCaa                                                 <  1:368466/77‑1 (MQ=255)                               cAGCTCGGCATCAATGAAGCCCACGCGCTGGGGAAAAAGTTTTATGACGTGGTCAACATTGCACCGCACAACGCCAAGc                                               <  1:368465/79‑1 (MQ=255)                               cAGCTCGGCATCAATGAAGCCCACGCGCTGGGGAAAAAGTTTTATGACGTGGTCAACATTGCACCGCACAACGCCAAGCt                                              <  2:368467/80‑1 (MQ=255)                                    cGGCATCAATGAAGCCCACGCGCTGGGGAAAAAGTTTTATGACGTGGTCTACATTGCACCGCACAACGCCAAGCTGAAAAc                                        >  2:368532/1‑81 (MQ=255)                                            aTGAAGCCCACGCGCTGGGGAAAAAGTTTTATGACGTGGTCAACATTGCACCGCACAACGCCAAGCTGAAAACCTTTATcc                                <  2:368468/81‑1 (MQ=255)                                              gAAGCCCACGCGCTGGGGAAAAAGTTTTATGACGTGGTCAACATTGCACCGCACAACGCCAAGCTGAAAACCTTTATCCGt                              >  1:368510/1‑81 (MQ=255)                                              gAAGCCCACGCGCTGGGGAAAAAGTTTTATGACGTGGTCAACATTGCACCGCACAACGCCAAGCTGAAAACCTTTATCCGt                              >  2:368516/1‑81 (MQ=255)                                                aGCCCACGCGCTGGGGAAAAAGTTTTATGACGTGGTCAACATTGCACCGCACAACGCCAAGCTGAAAACCTTTATCCGTGa                            >  2:368525/1‑81 (MQ=255)                                                  cccACGCGCTGGGGAAAAAGTTTTATGACGTGGTCAACATTGCACCGCACAACGCCAAGCTGAAAACCTTTATCCGTGAcc                          <  1:368469/81‑1 (MQ=255)                                                       gcgcTGGGGAAAAAGTTTTATGACGTGGTCAACATTGCACCGCACAACGCCAAGCTGAAAACCTTTATCCGTGACCTGaaa                     <  1:368470/81‑1 (MQ=255)                                                          cTGGGGAAAAAGTTTTATGACGTGGTCAACATTGCACCGCACAACGCCAAGCTGAAAACCTTTATCCGTGACCTGAAACCg                  <  2:368471/81‑1 (MQ=255)                                                            ggggAAAAAGTTTTATGACGTGGTCAACATTGCACCGCACAACGCCAAGCTGAAAACCTTTATCCGTGACCTGAAACCggt                >  1:368529/1‑81 (MQ=255)                                                                aaaaaGTTTTATGACGTGGTCAACATTGCACCGCACAACGCCAAGCTGAAAACCTTTATCCGTGACCTGAAACCGgtggtg            >  2:368508/1‑81 (MQ=255)                                                                   aaGTTTTATGACGTGGTCAACATTGCACCGCACAACGCCAAGCTGAAAACCTTTATCCGTGACCTGAAACCGGTGGTGGa          >  1:368519/1‑80 (MQ=255)                                                                   aaGTTTTATGACGTGGTCAACATTGCACCGCACAACGCCAAGCTGAAAACCTTTATCCGTGACCTGAAACCGGTGGTGGa          >  2:368526/1‑80 (MQ=255)                                                                     gTTTTATGACGTGGTCAACATTGCACCGCACAACGCCAAGCTGAAAACCTTTATCCGTGACCTGAAACCGGTGGTGGAAAt       <  1:368472/81‑1 (MQ=255)                                                                        ttATGACGTGGTCAACATTGCACCGCACAACGCCAAGCTGAAAACCTTTATCCGTGACCTGAAACCGGTGGTGGAAATgg     <  1:368473/80‑1 (MQ=255)                                                                          aTGACGTGGTCAACATTGCACCGCACAACGCCAAGCTGAAAACCTTTATCCGTGACCTGAAACCGGTGGTGGAAATGGGGc  >  1:368521/1‑81 (MQ=255)                                                                             |                                                                              CGCAACAACGAATTCAACCACGAAAATCTTCAGCTCGGCATCAATGAAGCCCACGCGCTGGGGAAAAAGTTTTATGTCGTGGTCAACATTGCACCGCACAACGCCAAGCTGAAAACCTTTATCCGTGACCTGAAACCGGTGGTGGAAATGGGGC  >  NC_000913/2165782‑2165935

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 26 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑