Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,681,841	G→A	Q34* (CAA→TAA)	yhjJ ←	peptidase M16 family protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,681,841	0	G	A	100.0%	129.5 / NA	36	Q34* (CAA→TAA)	yhjJ	peptidase M16 family protein
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (16/20);  total (16/20)

AATTTCAACACGATCGCTGGGACGCTGGGGGGTGGTCAGCACTTGCCACTGTAAACCGTTGGAAAGCGTCCCCTGTTGCCATGCTGGATCAGGCTGGAGCGCATCTGCCTGCACATAGCCAGCAGTGGCCATCATCAGCAAACCGCCCGCTAA  >  NC_000913/3681764‑3681916                                                                              |                                                                            aaTTTCAACACGATCGCTGGGACGCTGGGGGGTGGTCAGCACTTGCCACTGTAAACCGTTGGAAAGCGTCCCCTGTTACCa                                                                          >  1:716403/1‑81 (MQ=255)  aTTTCAACACGATCGCTGGGACGCTGGGGGGTGGTCAGCACTTGCCACTGTAAACCGTTGGAAAGCGTCCCCTGTTACCAt                                                                         >  2:716406/1‑81 (MQ=255)   tttCAACACGATCGCTGGGACGCTGGGGGGTGGTCAGCACTTGCCACTGTAAACCGTTGGAAAGCGTCCCCTGTTACCATg                                                                        >  1:716399/1‑81 (MQ=255)      cAACACGATCGCTGGGACGCTGGGGGGTGGTCAGCACTTGCCACTGTAAACCGTTGGAAAGCGTCCCCTGTTACCATGCTg                                                                     >  1:716409/1‑81 (MQ=255)            gATCGCTGGGACGCTGGGGGGTGGTCAGCACTTGCCACTGTAAACCGTTGGAAAGCGTCCCCTGTTACCATGCTGGATCAg                                                               >  2:716412/1‑81 (MQ=255)              tCGCTGGGACGCTGGGGGGTGGTCAGCACTTGCCACTGTAAACCGTTGGAAAGCGTCCCCTGTTACCATGCTGGATCAGGc                                                             <  1:716343/81‑1 (MQ=255)                    ggACGCTGGGGGGTGGTCAGCACTTGCCACTGTAAACCGTTGGAAAGCGTCCCCTGTTACCATGCTGGATCAGGCTGGa                                                         <  1:716344/79‑1 (MQ=255)                    ggACGCTGGGGGGTGGTCAGCACTTGCCACTGTAAACCGTTGGAAAGCGTCCCCTGTTACCATGCTGGATCAGGCTGGAg                                                        <  1:716345/80‑1 (MQ=255)                     gACGCTGGGGGGTGGTCAGCACTTGCCACTGTAAACCGTTGGAAAGCGTCCCCTGTTACCATGCTGGATCAGGCTGGAgcg                                                      >  2:716400/1‑81 (MQ=255)                        gCTGGGGGGTGGTCAGCACTTGCCACTGTAAACCGTTGGAAAGCGTCCCCTGTTACCATGCTGGATCAGGCTGGAGCGCAt                                                   <  2:716346/81‑1 (MQ=255)                               ggtggtCAGCACTTGCCACTGTAAACCGTTGGAAAGCGTCCCCTGTTACCATGCTGGATCAGGCTGGAGCGCATctgcctg                                            <  1:716347/81‑1 (MQ=255)                                   gtCAGCACTTGCCACTGTAAACCGTTGGAAAGCGTCCCCTGTTACCATGCTGGATCAGGCTGGAGCGCATCTGCCTGcac                                         <  2:716348/80‑1 (MQ=255)                                    tCAGCACTTGCCACTGTAAACCGTTGGAAAGCGTCCCCTGTTACCATGCTGGATCAGGCTGGAGCGCATCTGCCTGCACAt                                       <  2:716350/81‑1 (MQ=255)                                    tCAGCACTTGCCACTGTAAACCGTTGGAAAGCGTCCCCTGTTACCATGCTGGATCAGGCTGGAGCGCATCTGCCTGCACAt                                       <  1:716349/81‑1 (MQ=255)                                     cAGCACTTGCCACTGTAAACCGTTGGAAAGCGTCCCCTGTTACCATGCTGGATCAGGCTGGAGCGCATCTGCCTGcaca                                        >  1:716419/1‑79 (MQ=255)                                      aGCACTTGCCACTGTAAACCGTTGGAAAGCGTCCCCTGTTACCATGCTGGATCAGGCTGGAGCGCATCTGCCTGCACATa                                      <  1:716351/80‑1 (MQ=255)                                       gCACTTGCCACTGTAAACCGTTGGAAAGCGTCCCCTGTTACCATGCTGGATCAGGCTGGAGCGCATCTGCCTGCACATAGc                                    <  1:716352/81‑1 (MQ=255)                                                 cTGTAAACCGTTGGAAAGCGTCCCCTGTTACCATGCTGGATCAGGCTGGAGCGCATCTGCCTGCACATAGCCAGCAGTGGc                          >  2:716420/1‑81 (MQ=255)                                                        ccGTTGGAAAGCGTCCCCTGTTACCATGCTGGATCAGGCTGGAGCGCATCTGCCTGCACATAGCCAGCAGTGGCCatcatc                   <  2:716353/81‑1 (MQ=255)                                                          gTTGGAAAGCGTCCCCTGTTACCATGCTGGATCAGGCTGGAGCGCATCTGCCTGCACATAGCCAGCAGTGGCCATCATCAg                 >  1:716416/1‑81 (MQ=255)                                                          gTTGGAAAGCGTCCCCTGTTACCATGCTGGATCAGGCTGGAGCGCATCTGCCTGCACATAGCCAGCAGTGGCCATCATCAg                 >  1:716408/1‑81 (MQ=255)                                                             ggAAAGCGTCCCCTGTTACCATGCTGGATCAGGCTGGAGCGCATCTGCCTGCACATAGCCAGCAGTGGCCATCATCAGCa               >  2:716423/1‑80 (MQ=255)                                                                aaGCGTCCCCTGTTACCATGCTGGATCAGGCTGGAGCGCATCTGCCTGCACATAGCCAGCAGTGGCCATCATCAGCAAAcc           <  2:716354/81‑1 (MQ=255)                                                                 aGCGTCCCCTGTTACCATGCTGGATCAGGCTGGAGCGCATCTGCCTGCACATAGCCAGCAGTGGCCATCATCAGCAAAccg          <  1:716355/81‑1 (MQ=255)                                                                  gCGTCCCCTGTTACCATGCTGGATCAGGCTGGAGCGCATCTGCCTGCACATAGCCAGCAGTGGCCATCATCAGCAAAccgc         <  1:716356/81‑1 (MQ=255)                                                                   cGTCCCCTGTTACCATGCTGGATCAGGCTGGAGCGCATCTGCCTGCACATAGCCAGCAGTGGCCATCATCAGCAAAccgcc        >  2:716411/1‑81 (MQ=255)                                                                    gTCCCCTGTTACCATGCTGGATCAGGCTGGAGCGCATCTGCCTGCACATAGCCAGCAGTGGCCATCATCAGCAAAccgccc       >  1:716417/1‑81 (MQ=255)                                                                    gTCCCCTGTTACCATGCTGGATCAGGCTGGAGCGCATCTGCCTGCACATAGCCAGCAGTGGCCATCATCAGCAAAccgccc       <  1:716357/81‑1 (MQ=255)                                                                       cccTGTTACCATGCTGGATCAGGCTGGAGCGCATCTGCCTGCACATAGCCAGCAGTGGCCATCATCAGCAAAccgcccgc     >  1:716413/1‑80 (MQ=255)                                                                       cccTGTTACCATGCTGGATCAGGCTGGAGCGCATCTGCCTGCACATAGCCAGCAGTGGCCATCATCAGCAAACCGCCCGCt    <  2:716358/81‑1 (MQ=255)                                                                        ccTGTTACCATGCTGGATCAGGCTGGAGCGCATCTGCCTGCACATAGCCAGCAGTGGCCATCATCAGCAAACCGCCCGCTa   <  1:716360/81‑1 (MQ=255)                                                                        ccTGTTACCATGCTGGATCAGGCTGGAGCGCATCTGCCTGCACATAGCCAGCAGTGGCCATCATCAGCAAACCGCCCGCTa   <  1:716359/81‑1 (MQ=255)                                                                         cTGTTACCATGCTGGATCAGGCTGGAGCGCATCTGCCTGCACATAGCCAGCAGTGGCCATCATCAGCAAACCGCCCGCTaa  <  1:716361/81‑1 (MQ=255)                                                                         cTGTTACCATGCTGGATCAGGCTGGAGCGCATCTGCCTGCACATAGCCAGCAGTGGCCATCATCAGCAAACCGCCCGCTaa  <  2:716362/81‑1 (MQ=255)                                                                          tGTTACCATGCTGGATCAGGCTGGAGCGCATCTGCCTGCACATAGCCAGCAGTGGCCATCATCAGCAAACCGCCCGCTaa  >  1:716433/1‑80 (MQ=255)                                                                           gTTACCATGCTGGATCAGGCTGGAGCGCATCTGCCTGCACATAGCCAGCAGTGGCCATCATCAGCAAACCGCCCGCTaa  >  2:716414/1‑79 (MQ=255)                                                                              |                                                                            AATTTCAACACGATCGCTGGGACGCTGGGGGGTGGTCAGCACTTGCCACTGTAAACCGTTGGAAAGCGTCCCCTGTTGCCATGCTGGATCAGGCTGGAGCGCATCTGCCTGCACATAGCCAGCAGTGGCCATCATCAGCAAACCGCCCGCTAA  >  NC_000913/3681764‑3681916

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