Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,681,795	(G)6→7	coding (146/1497 nt)	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,789	1	.	G	100.0%	128.7 / NA	38	Q51P (CAG→CCG)	yhjJ	peptidase M16 family protein
Reads supporting (aligned to +/- strand):  ref base . (0/0);  new base G (16/22);  total (16/22)

CCTGTGCTGCGTCAAGGCCACCGCTTTGCGTTAGCGCAATACGAGGGATGGCGTGACTGTAACCGCTCTGTTGTGTACTTTCGGCGAGCGAACCGGTATTAACCAGCAGGCGAATTTCAACACGATCGCTGGGACGCT‑GGGGGGTGGTCAGCACTTGCCACTGTAAACCGTTGGAAAGCGTCCCCTGTTGCCATGCTGGATCAGGCTGGAGCGCATCTGCCTGCACATAGCCAGCAGTGGCCATCATCAGCAAACCGCCCGCTAAAAGTCGA  >  NC_000913/3681652‑3681923                                                                                                                                           |                                                                                                                                       ccTGTGCTGCGTCAAGGCCACCGCTTTGCGTTAGCGCAATACGAGGGATGGCGTGACTGTAACCGCTCTGTTGTGTACTTTCGGCGAGCGAACCGGTATTAACCAGCAGGCGAATTTCAACACGATCGCTGGGACGCT‑gg                                                                                                                                      >  1:754358/1‑140 (MQ=255) ccTGTGCTGCGTCAAGGCCACCGCTTTGCGTTAGCGCAATACGAGGGATGGCGTGACTGTAACCGCTCTGTTGTGTACTTTCGGCGAGCGAACCGGTATTAACCAGCAGGCGAATTTCAACACGATCGCTGGGACGCT‑gg                                                                                                                                      >  1:777415/1‑140 (MQ=255) ccTGTGCTGCGTCAAGGCCACCGCTTTGCGTTAGCGCAATACGAGGGATGGCGTGACTGTAACCGCTCTGTTGTGTACTTTCGGCGAGCGAACCGGTATTAACCAGCAGGCGAATTTCAACACGATCGCTGGGACGCT‑gg                                                                                                                                      >  2:129591/1‑140 (MQ=255)  cTGTGCTGCGTCAAGGCCACCGCTTTGCGTTAGCGCAATACGAGGGATGGCGTGACTGTAACCGCTCTGTTGTGTACTTTCGGCGAGCGAACCGGTATTAACCAGCAGGCGAATTTCAACACGATCGCTGGGACGCT‑ggg                                                                                                                                     <  1:201697/140‑1 (MQ=255)     tgctgcGTCAAGGCCCCCGCTTTGCGTTATCGCAATACGAGGGATGGCGTGACTGTAACCGCTCTGTTGTGTACTTTCGGCGAGCGAACCGGTATTAACCAGCAGGCGAATTTCAACACGATCGCTGGGACGCT‑gggggg                                                                                                                                  <  1:374923/140‑1 (MQ=255)     tgctgcGTCAAGGCCCCCGCTTAGCATTAACGCAATACGAGGGATGGCGTGACTGTAACCGCTCTGTTGTGTACTTTCGGCGAGCGAACCGGTATTAACCAGCAGGCGAATTTCAACACGATCGCTGGGACGCT‑gggggg                                                                                                                                  <  1:400071/140‑1 (MQ=255)     tgctgcGTCAAGGCCACCGCTTTGCGTTAGCGCAATACGAGGGATGGCGTGACTGTAACCGCTCTGTTGTGTACTTTCGGCGAGCGAACCGGTATTAACCAGCAGGCGAATTTCAACACGATCGCTGGGACGCT‑gggggg                                                                                                                                  >  2:149153/1‑140 (MQ=255)          cGTCAAGGCCACCGCTTTGCGTTAGCGCAATACGAGGGATGGCGTGACTGTAACCGCTCTGTTGTGTACTTTCGGCGAGCGAACCGGTATTAACCAGCAGGCGAATTTCAACACGATCGCTGGGACGCT‑gggggggtggt                                                                                                                             <  1:632867/140‑6 (MQ=255)          cGTCAAGGCCACCGCTTTGCGTTAGCGCAATACGAGGGATGGCGTGACTGTAACCGCTCTGTTGTGTACTTTCGGCGAGCGAACCGGTATTAACCAGCAGGCGAATTTCAACACGATCGCTGGGACGCT‑gggggggtggc                                                                                                                             >  2:75458/1‑135 (MQ=255)               aGGCCACCGCTTTGCGTTAGCGCAATACGAGGGATGGCGTGACTGTAACCGCTCTGTTGTGTACTTTCGGCGAGCGAACCGGTATTAACCAGCAGGCGAATTTCAACACGATCGCTGGGACGCTGGGGGGGTGGTCAGCa                                                                                                                         >  2:614008/1‑140 (MQ=255)                 gCCACCGCTTTGCGTTAGCGCAATACGAGGGATGGCGTGACTGTAACCGCTCTGTTGTGTACTTTCGGCGAGCGAACCGGTATTAACCAGCAGGCGAATTTCAACACGATCGCTGGGACGCTGGGGGGGTGGTCAGCACt                                                                                                                       <  1:505978/140‑1 (MQ=255)                 gCCACCGCTTTGCGTTAGCGCAATACGAGGGATGGCGTGACTGTAACCGCTCTGTTGTGTACTTTCGGCGAGCGAACCGGTATTAACCAGCAGGCGAATTTCAACACGATCGCTGGGACGCTGGGGGGGTGGTCAGCACt                                                                                                                       <  1:614008/140‑1 (MQ=255)                 gCCACCGCTTTGCGTTAGCGCAATACGAGGGATGGCGTGACTGTAACCGCTCTGTTGTGTACTTTCGGCGAGCGAACCGGTATTAACCAGCAGGCGAATTTCAACACGATCGCTGGGACGCTGGGGGGGTGGTCAGCACt                                                                                                                       <  1:641351/140‑1 (MQ=255)                   cACCGCTTTGCGTTAGCGCAATACGAGGGATGGCGTGACTGTAACCGCTCTGTTGTGTACTTTCGGCGAGCGAACCGGTATTAACCAGCAGGCGAATTTCAACACGATCGCTGGGACGCTGGGGGGGTGGTCAGCACTTg                                                                                                                     <  2:808817/140‑1 (MQ=255)                     ccGCTTTGCGTTAGCGCAATACGAGGGATGGCGTGACTGTAACCGCTCTGTTGTGTACTTTCGGCGAGCGAACCGGTATTAACCAGCAGGCGAATTTCAACACGATCGCTGGGACGCTGGGGGGGTGGTCAGCACTTGcc                                                                                                                   <  1:360214/140‑1 (MQ=255)                     ccGCTTTGCGTTAGCGCAATACGAGGGATGGCGTGACTGTAACCGCTCTGTTGTGTACTTTCGGCGAGCGAACCGGTATTAACCAGCAGGCGAATTTCAACACGATCGCTGGGACGCTGGGGGGGTGGTCAGCACTTGcc                                                                                                                   <  1:129591/140‑1 (MQ=255)                           tGCGTTAGCGCAATACGAGGGATGGCGTGACTGTAACCGCTCTGTTGTGTACTTTCGGCGAGCGAACCGGTATTAACCAGCAGGCGAATTTCAACACGATCGCTGGGACGCTGGGGGGGTGGTCAGCACTTGCCACTGTa                                                                                                             <  1:539589/140‑1 (MQ=255)                           tGCGTTAGCGCAATACGAGGGATGGCGTGACTGTAACCGCTCTGTTGTGTACTTTCGGCGAGCGAACCGGTATTAACCAGCAGGCGAATTTCAACACGATCGCTGGGACGCTGGGGGGGTGGTCAGCACTTGCCACTGTa                                                                                                             <  1:22997/140‑1 (MQ=255)                                tAGCGCAATACGAGGGATGGCGTGACTGTAACCGCTCTGTTGTGTACTTTCGGCGAGCGAACCGGTATTAACCAGCAGGCGAATTTCAGCACGATCGCTGGGACGCTGGGGGGGTGGTCAGCACTTGCCACTGTAAACCg                                                                                                        <  1:423542/140‑1 (MQ=255)                                 aGCGCAATACGAGGGATGGCGTGACTGTAACCGCTCTGTTGTGTACTTTCGGCGAGCGAACCGGTATTAACCAGCAGGCGAATTTCAACACGATCGCTGGGACGGTGGGGGGGTGGTAAGCAGTTGCCACTGTAAGCCGt                                                                                                       >  2:275743/1‑140 (MQ=255)                                            aGGGATGGCGTGACTGTAACCGCTCTGTTGTGTACTTTCGGCGAGCGAACCGGTATTAACCAGCAGGCGAATTTCAACACGATCGCTGGGACGCTGGGGGGGTGGTCAGCGCTTGCCGCTGTAAGCCGTTGGAAAGCGTc                                                                                            >  1:221562/1‑140 (MQ=255)                                               gATGGCGTGACTGTAATCGCTCTGTTGTGTAGTTTCGGCGAGCGAACCCGTATTAACCAGCAGGCGAATTTCAACACGATCGCTGGGACGCTGGGGGGGTGGTCAGCACTTGCCACTGTAAACCGTTGGAAAGCGTcccc                                                                                         <  2:366777/140‑1 (MQ=255)                                                aTGGCGTGACTGTAACCGCTCTGTTGTGTACTTTCGGCGAGCGAACCGGTATTAACCAGCAGGCGAATTTCAACACGATCGCTGGGACGCTGGGGGGGTGGTCAGCACTTGCCACTGTAAACCGTTGGAAAGCGTCCCCt                                                                                        >  2:748190/1‑140 (MQ=255)                                                  ggCGTGACTGTAACCGCTCTGTTGTGTACTTTCGGCGAGCGAACCGGTATTAACCAGCAGGCGAATTTCAACACGATCGCTGGGACGCTGGGGGGGTGGTCAGCACTTGCCACTGTAAACCGTTGGAAAGCGTCCCCTGt                                                                                      >  1:159090/1‑140 (MQ=255)                                                  ggCGTGACTGTAACCGCTCTGTTGTGTACTTTCGGCGAGCGAACCGGTATTAACCAGCAGGCGAATTTCAACACGATCGCTGGGACGCTGGGGGGGTGGTCAGCACTTGCCACTGTAAACCGTTGGAAAGCGTCCCCTGt                                                                                      >  2:7984/1‑140 (MQ=255)                                                  ggCGTGACTGTAACCGCTCTGTTGTGTACTTTCGGCGAGCGAACCGGTATTAACCAGCAGGCGAATTTCAACACGATCGCTGGGACGCTGGGGGGGTGGTCAGCACTTGCCACTGTAAACCGTTGGAAAGCGTCCCCTGt                                                                                      >  1:65088/1‑140 (MQ=255)                                                  ggCGTGACTGTAACCGCTCTGTTGTGTACTTTCGGCGAGCGAACCGGTATTAACCAGCAGGCGAATTTCAACACGATCGCTGGGACGCTGGGGGGGTGGTCAGCACTTGCCACTGTAAACCGTTGGAAAGCGTCCCCTGt                                                                                      >  1:690445/1‑140 (MQ=255)                                                   gCGTGACTGTAACCGCTCTGTTGTGTACTTTCGGCGAGCGAACCGGTATTAACCAGCAGGCGAATTTCAACACGATCGCTGGGACGCTGGGGGGGTGGTCAGTACTTGCCACTGGAAGCCGTTGGAAAGCGTCCCCTGtt                                                                                     >  1:366777/1‑140 (MQ=255)                                                    cGTGACTGTAACCGCTCTGTTGTGTACTTTCGGCGAGCGAACCGGTATTAACCAGCAGGCGAATTTCAACACGATCGCTGGGACGCTGGGGGGGTGGTCAGCACTTGCCACTGTAAACCGTTGGAAAGCGTCCCCTGTTg                                                                                    <  1:443207/140‑1 (MQ=255)                                                         cTGTAACCGCTCTGTTGTGTACTTTCGGCGAGCGAACCGGTATTAACCAGCAGGCGAATTTCAACACGATCGCTGGGACGCTGGGGGGGTGGTCAGCACTTGCCACTGTAAACCGTTGGAAAGCGTCCCCTGTTGCCATg                                                                               <  2:65088/140‑1 (MQ=255)                                                            tAACCGCTCTGTTGTGTACTTTCGGCGAGCGAACCGGTATTAACCAGCAGGCGAATTTCAACACGATCGCTGGGACGCTGGGGGGGTGGTCAGCACTTGCCACTGTAAACCGTTGGAAAGCGTCCCCTGTTGCCATGCTg                                                                            <  1:616472/140‑1 (MQ=255)                                                                 gCTCTGTTGTGTACTTTCGGCGAGCGAACCGGTATTAACCAGCAGGCGAATTTCAACACGATCGCTGGGACGCTGGGGGGGTGGTCAGCACTTGCCACTGTAAACCGTTGGAAAGCGTCCCCTGTTGCCATGCTGGATCa                                                                       <  2:151530/140‑1 (MQ=255)                                                                    ctGTTGTGTACTTTCGGCGAGCGAACCGGTATTAACCAGCAGGCGAATTTCAACACGATCGCTGGGACGCTGGGGGGGTGGTCAGCACTTGCCACTGTAAACCGTTGGAAAGCGTCCCCTGTTGCCATGCTGGATCAGGc                                                                    <  1:427441/140‑1 (MQ=255)                                                                      gttgtGTACTTTCGGCGAGCGAACCGGTATTAACCAGCAGGCGAATTTCAACACGATCGCTGGGACGCTGGGGGGGTGGTCAGCACTTGCCACTGTAAACCGTTGGAAAGCGTCCCCTGTTGCCATGCTGGATCAGGCTg                                                                  <  1:7984/140‑1 (MQ=255)                                                                      gttgtGTACTTTCGGCGAGCGAACCGGTATTAACCAGCAGGCGAATTTCAACACGATCGCTGGGACGCTGGGGGGGTGGTCAGCACTTGCCACTGTAAACCGTTGGAAAGCGTCCCCTGTTGCCATGCTGGATCAGGCTg                                                                  >  1:419416/1‑140 (MQ=255)                                                                        tgtgTACTTTCGGCGAGCGAACCGGTATTAACCAGCAGGCGAATTTCAACACGATCGCTGGGACGCTGGGGGGGTGGTCAGCACTTGCCACTGTAAACCGTTGGAAAGCGTCCCCTGTTGCCATGCTGGATCAGGCTGGa                                                                <  1:664416/140‑1 (MQ=255)                                                                                 tCGGCGAGCGAACCGGTATTAACCAGCAGGCGAATTTCAACACGATCGCTGGGACGCTGGGGGGGTGGTCAGCACTTGCCACTGTAAACCGTTGGAAAGCTTCCCCTGATGCCATGCCGGATCAGGCTGGAGCGCATctg                                                       >  2:369076/1‑140 (MQ=255)                                                                                 tCGGCGAGCGAACCGGTATTAACCAGCAGGCGAATTTCAACACGATCGCTGGGACGCTGGGGGGGTGGTCAGCACTTGCCACTGTAAACCGTTGGAAAGCGTCCCCTGTTGCCATGCTGGATCAGGCTGGAGCGCATctg                                                       >  2:416256/1‑140 (MQ=255)                                                                                        gcgaACCGGTATTAACCAGCAGTCGAATTTCAACACGCTCGCTGGGACGCTGGGGGGGTGGTCAGCACTTGCCACTGTAAACCGTTGGAAAGCGTCCCCTGTTGCCATGCTGGATCAGGCTGGAGCGCATCTGCCTGcac                                                <  1:265150/140‑1 (MQ=255)                                                                                         cgaACCGGTATTAACCAGCAGGCGAATTTCAACACGATCGCTGGGACGCTGGGGGGGTGGTCAGCACTTGCCACTGTAAACCGTTGGAAAGCGTCCCCTGTTGCCATGCTGGATCAGGCTGGAGCGCATCTGCCTGcaca                                               >  2:796718/1‑140 (MQ=255)                                                                                                  aTTAACCAGCAGGCGAATTTCAACACGATCGCTGGGACGCTGGGGGGGTGGTCAGCACTTGCCACTGTAAACCGTTGGAAAGCGTCCCCTGTTGCCATGCTGGATCAGGCTGGAGCGCATCTGCCTGCACATAGCcagca                                      <  2:750747/140‑1 (MQ=255)                                                                                                   ttAACCAGCAGGCGAATTTCAACACGATCGCTGGGACGCTGGGGGGGTGGTCAGCACTTGCCACTGTAAACCGTTGGAAAGCGTCCCCTGTTGCCATGCTGGATCAGGCTGGAGCGCATCTGCCTGCACATAGCcagcag                                     <  2:754358/140‑1 (MQ=255)                                                                                                           cagGCGAATTTCAACACGATCGCTGGGACGCTGGGGGGGTGGTCAGCACTTGCCACTGTAAACCGTTGGAAAGCGTCCCCTGTTGCCATGCTGGATCAGGCTTGAGCGCAACTGGCTCCACATATCCAGCAGTGGCcatc                             >  1:244903/1‑140 (MQ=255)                                                                                                           cagGCGAATTTCAACACGATCGCTGGGACGCTGGGGGGGTGGTCAGCACTTGCCACTGTAAACCGTTGGAAAGCGTCCCCTGTTGCCATGCTGGATCAGGCTGGAGCGCATCTGCCTGCACATAGCCAGCAGTGGCcatc                             <  2:159090/140‑1 (MQ=255)                                                                                                           cagGCGAATTTCAACACGATCGCTGGGACGCTGGGGGGGTGGTCAGCACTTGCCACTGTAAACCGTTGGAAAGCGTCCCCTGTTGCCATGCTGGATCAGGCTGGAGCGCATCTGCCTGCACATAGCCAGCAGTGGCcatc                             >  1:604344/1‑140 (MQ=255)                                                                                                           cagGCGAATTTCAACACGATCGCTGGGACGCTGGGGGGGTGGTCAGCACTTGCCACTGTAAACCGTTGGAAAGCGTCCCCTGTTGCCATGCTGGATCAGGCTGGAGCGCATCTGCCTGCACATAGCCAGCAGTGGCcatc                             >  1:774406/1‑140 (MQ=255)                                                                                                                        acacGATCGCTGGGACGCTGGGGGGGTGGTCAGCACTTGCCACTGTAAACCGTTGGAAAGCGTCCCCTGTTGCCATGCTGGATCAGGCTGGAGCGCATCTGCCTGCACATAGCCAGCAGTGGCCATCATCAGCAAAccgc                <  2:737789/140‑1 (MQ=255)                                                                                                                                      aCGCTGGGGGGGTGGTCAGCACTTGCCACTGTAAACCGTTGGAAAGCGTCCCCTGTTGCCATGCTGGATCAGGCTGGAGCGCATCTGCCTGCACATAGCCAGCAGTGGCCATCATCAGCAAACCGCCCGCTAAAAGTCGa  >  1:311192/1‑140 (MQ=255)                                                                                                                                           |                                                                                                                                       CCTGTGCTGCGTCAAGGCCACCGCTTTGCGTTAGCGCAATACGAGGGATGGCGTGACTGTAACCGCTCTGTTGTGTACTTTCGGCGAGCGAACCGGTATTAACCAGCAGGCGAATTTCAACACGATCGCTGGGACGCT‑GGGGGGTGGTCAGCACTTGCCACTGTAAACCGTTGGAAAGCGTCCCCTGTTGCCATGCTGGATCAGGCTGGAGCGCATCTGCCTGCACATAGCCAGCAGTGGCCATCATCAGCAAACCGCCCGCTAAAAGTCGA  >  NC_000913/3681652‑3681923

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