Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,302,518	Δ1 bp	coding (30/879 nt)	yhbV →	putative peptidase YhbV

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,302,517	0	G	.	100.0%	160.9 / NA	38	coding (29/879 nt)	yhbV	putative peptidase YhbV
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base . (23/15);  total (23/15)

CACGCAAACCACTCTTGGCGCATATCACCGTAAATGGCAGTGAGAAAAGCAATGAAATATTCCTTAGGGCCAGTGCTGTGGTACTGGCCAAAAGAGACGCTGGAAGAATTTTATCAGCAGGCCGCCACCAGCAGCGCCGACGTGATTTATCTTGGT  >  NC_000913/3302438‑3302593                                                                                |                                                                             cACGCAAACCACTCTTGGCGCATATCACCGTAAATGGCAGTGAGAAAAGCAATGAAATATTCCTTAGGGCCAGTGCtgtg                                                                              >  2:557056/1‑80 (MQ=255)    gCAAACCACTCTTGGCGCATATCACCGTAAATGGCAGTGAGAAAAGCAATGAAATATTCCTTAGGGCCAGTGCtgtgtact                                                                          >  1:557059/1‑77 (MQ=255)    gCAAACCACTCTTGGCGCATATCACCGTAAATGGCAGTGAGAAAAGCAATGAAATATTCCTTAGGGCCAGTGCtgtgtact                                                                          >  1:557041/1‑77 (MQ=255)        aCCACTCTTGGCGCATATCACCGTAAATGGCAGTGAGAAAAGCAATGAAATATTCCTTAGGGCCAGTGCTGT‑GTACTGGcc                                                                     <  1:556998/81‑1 (MQ=255)        aCCACTCTTGGCGCATATCACCGTAAATGGCAGTGAGAAAAGCAATGAAATATTCCTTAGGGCCAGTGCTGT‑GTACTGGc                                                                      >  2:557039/1‑80 (MQ=255)          cACTCTTGGCGCATATCACCGTAAATGGCAGTGAGAAAAGCAATGAAATATTCCTTAGGGCCAGTGCTGT‑GTACTGGCCa                                                                    >  2:557035/1‑80 (MQ=255)                  gcgcATATCACCGTAAATGGCAGTGAGAAAAGCAATGAAATATTCCTTAGGGCCAGTGCTGT‑GTACTGGCCAAAAGAGAc                                                            <  1:556999/80‑1 (MQ=255)                  gcgcATATCACCGTAAATGGCAGTGAGAAAAGCAATGAAATATTCCTTAGGGCCAGTGCTGT‑GTACTGGCCAAAAGAGACg                                                           >  1:557043/1‑81 (MQ=255)                   cgcATATCACCGTAAATGGCAGTGAGAAAAGCAATGAAATATTCCTTAGGGCCAGTGCTGT‑GTACTGGCCAAAAGAGACGc                                                          >  2:557047/1‑81 (MQ=255)                       tatCACCGTAAATGGCAGTGAGAAAAGCAATGAAATATTCCTTAGGGCCAGTGCTGT‑GTACTGGCCAAAAGAGACGCTGga                                                      <  2:556999/81‑1 (MQ=255)                                 aaTGGCAGTGAGAAAAGCAATGAAATATTCCTTAGGGCCAGTGCTGT‑GTACTGGCCAAAAGAGACGCTGGAAGAATTTTAt                                            >  1:557047/1‑81 (MQ=255)                                 aaTGGCAGTGAGAAAAGCAATGAAATATTCCTTAGGGCCAGTGCTGT‑GTACTGGCCAAAAGAGACGCTGGAAGAATTTTAt                                            <  1:557001/81‑1 (MQ=255)                                  aTGGCAGTGAGAAAAGCAATGAAATATTCCTTAGGGCCAGTGCTGT‑GTACTGGCCAAAAGAGACGCTGGAAGAATTTTAt                                            >  1:557054/1‑80 (MQ=255)                                   tGGCAGTGAGAAAAGCAATGAAATATTCCTTAGGGCCAGTGCTGT‑GTACTGGCCAAAAGAGACGCTGGAAGAATTTTAt                                            >  2:557050/1‑79 (MQ=255)                                    ggCAGTGAGAAAAGCAATGAAATATTCCTTAGGGCCAGTGCTGT‑GTACTGGCCAAAAGAGACGCTGGAAGAATTTTATcag                                         >  1:557050/1‑81 (MQ=255)                                     gCAGTGAGAAAAGCAATGAAATATTCCTTAGGGCCAGTGCTGT‑GTACTGGCCAAAAGAGACGCTGGAAGAATTTTATcagc                                        <  2:557001/81‑1 (MQ=255)                                       aGTGAGAAAAGCAATGAAATATTCCTTAGGGCCAGTGCTGT‑GTACTGGCCAAAAGAGACGCTGGAAGAATTTTATcagcag                                      <  1:557003/81‑1 (MQ=255)                                             aaaaGCAATGAAATATTCCTTAGGGCCAGTGCTGT‑GTACTGGCCAAAAGAGACGCTGGAAGAATTTTATCAGCAGgccgc                                 >  2:557055/1‑80 (MQ=255)                                             aaaaGCAATGAAATATTCCTTAGGGCCAGTGCTGT‑GTACTGGCCAAAAGAGACGCTGGAAGAATTTTATCAGCAGgccgc                                 >  1:557064/1‑80 (MQ=255)                                                aGCAATGAAATATTCCTTAGGGCCAGTGCTGT‑GTACTGGCCAAAAGAGACGCTGGAAGAATTTTATCAGCAGGCCGccacc                             <  1:557004/81‑1 (MQ=255)                                                   aaTGAAATATTCCTTAGGGCCAGTGCTGT‑GTACTGGCCAAAAGAGACGCTGGAAGAATTTTATCAGCAGGCCGCCACcagc                          <  1:557005/81‑1 (MQ=255)                                                     tGAAATATTCCTTAGGGCCAGTGCTGT‑GTACTGGCCAAAAGAGACGCTGGAAGAATTTTATCAGCAGGCCGCCACcagcag                        >  1:557060/1‑81 (MQ=255)                                                      gAAATATTCCTTAGGGCCAGTGCTGT‑GTACTGGCCAAAAGAGACGCTGGAAGAATTTTATCAGCAGGCCGCCACCagcagc                       >  1:557055/1‑81 (MQ=255)                                                         atatTCCTTAGGGCCAGTGCTGT‑GTACTGGCCAAAAGAGACGCTGGAAGAATTTTATCAGCAGGCCGCCACCAGCAGCGcc                    <  2:557005/81‑1 (MQ=255)                                                         atatTCCTTAGGGCCAGTGCTGT‑GTACTGGCCAAAAGAGACGCTGGAAGAATTTTATCAGCAGGCCGCCACCAGCAGCGcc                    >  2:557033/1‑81 (MQ=255)                                                           atTCCTTAGGGCCAGTGCTGT‑GTACTGGCCAAAAGAGACGCTGGAAGAATTTTATCAGCAGGCCGCCACCAGCAGCGCCGa                  >  2:557062/1‑81 (MQ=255)                                                            ttCCTTAGGGCCAGTGCTGT‑GTACTGGCCAAAAGAGACGCTGGAAGAATTTTATCAGCAGGCCGCCACCAGCAGCGCCGAc                 <  1:557007/81‑1 (MQ=255)                                                              ccTTAGGGCCAGTGCTGT‑GTACTGGCCAAAAGAGACGCTGGAAGAATTTTATCAGCAGGCCGCCACCAGCAGCGCCGACg                <  1:557008/80‑1 (MQ=255)                                                              ccTTAGGGCCAGTGCTGT‑GTACTGGCCAAAAGAGACGCTGGAAGAATTTTATCAGCAGGCCGCCACCAGCAGCGCCGACGt               >  1:557075/1‑81 (MQ=255)                                                              ccTTAGGGCCAGTGCTGT‑GTACTGGCCAAAAGAGACGCTGGAAGAATTTTATCAGCAGGCCGCCACCAGCAGCGCCGACGt               >  1:557053/1‑81 (MQ=255)                                                               cTTAGGGCCAGTGCTGT‑GTACTGGCCAAAAGAGACGCTGGAAGAATTTTATCAGCAGGCCGCCACCAGCAGCGCCGACGt               >  1:557029/1‑80 (MQ=255)                                                               cTTAGGGCCAGTGCTGT‑GTACTGGCCAAAAGAGACGCTGGAAGAATTTTATCAGCAGGCCGCCACCAGCAGCGCCGACGt               <  1:557009/80‑1 (MQ=255)                                                                ttAGGGCCAGTGCTGT‑GTACTGGCCAAAAGAGACGCTGGAAGAATTTTATCAGCAGGCCGCCACCAGCAGCGCCGACGTGa             >  2:557051/1‑81 (MQ=255)                                                                  aGGGCCAGTGCTGT‑GTACTGGCCAAAAGAGACGCTGGAAGAATTTTATCAGCAGGCCGCCACCAGCAGCGCCGACGTGAtt           >  2:557068/1‑81 (MQ=255)                                                                   gggCCAGTGCTGT‑GTACTGGCCAAAAGAGACGCTGGAAGAATTTTATCAGCAGGCCGCCACCAGCAGCGCCGACGTGAttt          <  1:557010/81‑1 (MQ=255)                                                                    ggCCAGTGCTGT‑GTACTGGCCAAAAGAGACGCTGGAAGAATTTTATCAGCAGGCCGCCACCAGCAGCGCCGACGTGATTTa         <  1:557011/81‑1 (MQ=255)                                                                     gCCAGTGCTGT‑GTACTGGCCAAAAGAGACGCTGGAAGAATTTTATCAGCAGGCCGCCACCAGCAGCGCCGACGTGATTTAt        >  2:557070/1‑81 (MQ=255)                                                                        aGTGCTGT‑GTACTGGCCAAAAGAGACGCTGGAAGAATTTTATCAGCAGGCCGCCACCAGCAGCGCCGACGTGATTTATCtt     <  2:557011/81‑1 (MQ=255)                                                                        aGTGCTGT‑GTACTGGCCAAAAGAGACGCTGGAAGAATTTTATCAGCAGGCCGCCACCAGCAGCGCCGACGTGATTTATCtt     >  1:557074/1‑81 (MQ=255)                                                                        aGTGCTGT‑GTACTGGCCAAAAGAGACGCTGGAAGAATTTTATCAGCAGGCCGCCACCAGCAGCGCCGACGTGATTTATCtt     >  2:557069/1‑81 (MQ=255)                                                                           gCTGT‑GTACTGGCCAAAAGAGACGCTGGAAGAATTTTATCAGCAGGCCGCCACCAGCAGCGCCGACGTGATTTATCTTGGt  >  1:557061/1‑81 (MQ=255)                                                                                |                                                                             CACGCAAACCACTCTTGGCGCATATCACCGTAAATGGCAGTGAGAAAAGCAATGAAATATTCCTTAGGGCCAGTGCTGTGGTACTGGCCAAAAGAGACGCTGGAAGAATTTTATCAGCAGGCCGCCACCAGCAGCGCCGACGTGATTTATCTTGGT  >  NC_000913/3302438‑3302593

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