Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,695,997	C→T	*63* (TAG→TAA)	yhjR ←	DUF2629 family protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,695,997	0	C	T	100.0%	150.0 / NA	47	*63* (TAG→TAA)	yhjR	DUF2629 family protein
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (26/21);  total (26/21)

CATCGACCACCAGGACATTTTCTCCCAACATTTGTAATGACCAGGCTAATGCGGCGGTGATGGTTGTTGTCCCCACGCCTCCCCGCACCCCCTGCTATCCCAGTACGGCCATTCACCAATCCCCTACTTTTGTTGCGCAAACTCTGCCAGCAACGGCCAGCGTTTTAATGCCGCGGCCAACTGTTCGCGTTGGGAAATATCGGCATAATCAATATCAGGCAGTGAAAATGCCTGCTTTAACGCCACAATATCATTCTGG  >  NC_000913/3695874‑3696132                                                                                                                            |                                                                                                                                        caTCGACCACCAGGACATTTTCTCCCAACATTTGTAATGACCAGGCTAATGCGGCGGTGATGGTTGTTGTCCCCACGCCTCCCCGCACCCCCTGCTATCCCAGTACGGCCATTCACCAATCCCTTACTTTTGTTGCGCa                                                                                                                          >  2:1013699/1‑139 (MQ=255)   tCGACCACCAGGACATTTTCTCCCAACATTTGTAATGACCAGGCTAATGCGGCGGTGATGGTTGTTGTCCCCACGCCTCCCCGCACCCCCTGCTATCCCAGTACGGCCATTCACCAATCCCTTACTTTTGTTGCGCaaa                                                                                                                        <  2:313951/139‑1 (MQ=255)           caGGACATTTTCTCCCAACATTTGTAATGACCAGGCTAATGCGGCGGTGATGGTTGTTGTCCCCACGCCTCCCCGCACCCCCTGCTATCCCAGTACGGCCATTCACCAATCCCTTACTTTTGTTGCGCaa                                                                                                                         >  2:777833/1‑130 (MQ=255)           caGGACATTTTCTCCCAACATTTGTAATGACCAGGCTAATGCGGCGGTGATGGTTGTTGTCCCCACGCCTCCCCGCACCCCCTGCTATCCCAGTACGGCCATTCACCAATCCCTTACTTTTGTTGCGCaa                                                                                                                         <  1:777833/130‑1 (MQ=255)             ggACATTTTCTCCCAACATTTGTAATGACCAGGCTAATGCGGCGGTGATGGTTGTTGTCCCCACGCCTCCCCGCACCCCCTGCTATCCCAGTACGGCCATTCACCAATCCCTTACTTTTGTTGCGCAAACTCTGCCAGc                                                                                                              <  2:753622/139‑1 (MQ=255)                 aTTTTCTCCCAACATTTGTAATGACCAGGCTAATGCGGCGGTGATGGTTGTTGTCCCCACGCCTCCCCGCACCCCCTGCTATCCCAGTACGGCCATTCACCAATCCCTTACTTTTGTTGCGCAAACTCTGCCAGCAACg                                                                                                          >  2:424257/1‑139 (MQ=255)                    ttCTCCCAACATTTGTAATGACCAGGCTAATGCGGCGGTGATGGTTGTTGTCCCCACGCCTCCCCGCACCCCCTGCTATCCCAGTACGGCCATTCACCAATCCCTTACTTTTGTTGCGCAAACTCTGCCAGCAACGGcc                                                                                                       <  1:424257/139‑1 (MQ=255)                     tctcCCAACATTTGTAATGACCAGGCTAATGCGGCGGTGATGGTTGTTGTCCCCACGCCTCCCCGCACCCCCTGCTATCCCCGTACGGCCATTCACCAATCCCTTACTTTTGTTGCGCAAACTCTGCCAGCAACGGCCa                                                                                                      >  1:1107973/1‑139 (MQ=255)                                        aCCAGGCTAATGCGGCGGTGATGGTTGTTGTCCCCACGCCTCCCCGCACCCCCTGCTATCCCAGTACGGCCATTCACCAATCCCTTACTTTTGTTGCGCAAACTCTGCCAGCAACGGCCAGCGTTTTAATGCCGCGGcc                                                                                   >  1:810710/1‑139 (MQ=255)                                        aCCAGGCTAATGCGGCGGTGATGGTTGTTGTCCCCACGCCTCCCCGCACCCCCTGCTATCCCAGTACGGCCATTCACCAATCCCTTACTTTTGTTGCGCAAACTCTGCCAGCAACGGCCAGCGTTTTAATGCCGCGGcc                                                                                   >  1:838070/1‑139 (MQ=255)                                        aCCAGGCTAATGCAGCGGTGATGGTTGTTGTCCCCACGCCTCCCCGCACCCCCTGCTATCCCAGTACGGCCATTCACCAATCCCTTACTTTTGTTGCGCAAACTCTGCCAGCAACGGCCAGCGTTTTAATGCCGCGGcc                                                                                   >  2:71060/1‑139 (MQ=255)                                         ccAGGCTAATGCGGCGGTGATGGTTGTTGTCCCCACGCCTCCCCGCACCCCCTGCTATCCCAGTACGGCCATTCACCAATCCCTTACTTTTGTTGCGCAAACTCTGCCAGCAACGGCCAGCGTTTTAATGCCGCGGCCa                                                                                  >  2:686757/1‑139 (MQ=255)                                               tAATGCGGCGGTGATGGTTGTTGTCCCCACGCCTCCCCGCACCCCCTGCTATCCCAGTACGGCCATTCACCAATCCCTTACTTTTGTTGCGCAAACTCTGCCAGCAACGGCCAGCGTTTTAATGCCGCGGCCAACTGtt                                                                            <  2:623759/139‑1 (MQ=255)                                                    cggcggTGATGGTTGTTGTCCCCACGCCTCCCCGCACCCCCTGCTATCCCAGTACGGCCATTCACCAATCCCTTACTTTTGTTGCGCAAACTCTg                                                                                                                   >  1:182327/1‑95 (MQ=255)                                                    cggcggTGATGGTTGTTGTCCCCACGCCTCCCCGCACCCCCTGCTATCCCAGTACGGCCATTCACCAATCCCTTACTTTTGTTGCGCAAACTCTg                                                                                                                   <  2:182327/95‑1 (MQ=255)                                                       cggTGATGGTTGTTGTCCCCACGCCTCCCCGCACCCCCTGCTATCCCAGTACGGCCATTCACCAATCCCTTACTTTTGTTGCGCAAACTCTGCCAGCAACGGCCAGCGTTTTAATGCCGCGGCCAACTGTTCGCGTTgg                                                                    >  1:901653/1‑139 (MQ=255)                                                       cggTGATGGTTGTTGTCCCCACGCCTCCCCGCACCCCCTGCTATCCCAGTACGGCCATTCACCAATCCCTTACTTTTGTTGCGCAAACTCTGCCAGCAACGGCCAGCGTTTTAATGCCGCGGCCAACTGTTCGCGTTgg                                                                    >  1:903278/1‑139 (MQ=255)                                                              ggTTGTTGTCCCCACGCCTCCCCGCACCCCCTGCTATCCCAGTACGGCCATTCACCAATCCCTTACTTTTGTTGCGCAACCTCTGCCAGCAACGGCCAGCGTTTTAATGCCGCGGCCAACTGTTCGCGTTGTGAAatat                                                             <  2:1107973/139‑1 (MQ=255)                                                              ggTTGTTGTCCCCACGCCTCCCCGCACCCCCTGCTATCCCAGTACGGCCATTCACCAATCCCTTACTTTTGTTGCGCAAACTCTGCCAGCAACGGCCAGCGTTTTAATGCCGCGGCCAACTGTTCGCGTTGGGAAatat                                                             <  1:71060/139‑1 (MQ=255)                                                                 tgttgtCCCCACGTCTCCCCGCACCCCCTGCTATCCCAGTACGGCCATTCACCAATCCCTTACTTTTGTTGCGCAAACTCTGCCAGCAACGGCCAGCGTTTTAATGCCGCGGCCAACTGTTCGCGTTGGGAAATATCgg                                                          >  2:719775/1‑139 (MQ=255)                                                                  gttgtCCCCACGCCTCCCCGCACCCCCTGCTATCCCAGTACGGCCATTCACCAATCCCTTACTTTTGTTGCGCAAACTCTGCCAGCAACGGCCAGCGTTTTAATGCCGCGGCCAACTGTTCGCGTTGGGAAATATCGGc                                                         <  2:265227/139‑1 (MQ=255)                                                                   ttgtCCCCACGCCTCCCCGCACCCCCTGCTATCCCAGTACGGCCATTCACCAATCCCTTACTTTTGTTGCGCAAACTCTGCCAGCAACGGCCAGCGTTTTAATGCCGCGGCCAACTGTTCGCGTTGGGAAATATCGGCa                                                        <  1:299399/139‑1 (MQ=255)                                                                        cccACGCCTCCCCGCACCCCCTGCTATCCCAGTACGGCCATTCACCAATCCCTTACTTTTGTTGCGCAAACTCTGCCAGCAACGGCCAGCGTTTTAATGCCGCGGCCAACTGTTCGCGTTGGGAAATATCGGCATAATc                                                   >  1:436623/1‑139 (MQ=255)                                                                          cACGCCTCCCCGCACCCCCTGCTATCCCAGTACGGCCATTCACCAATCCCTTACTTTTGTTGCGCAAACTCTGCCAGCAACGGCCAGCGTTTTAATGCCGCGGCCAACTGTTCGCGTTGGGAAATATCGGCATAATCaa                                                 >  2:900207/1‑139 (MQ=255)                                                                                           ccTGCTATCCCAGTACGGCCATTCACCAATCCCTTACTTTTGTTGCGCAAACTCTGCCAGCAACGGCCAGCGTTTTAATGCCGCGGCCAACTGTTCGCGTTGGGAAATATCGGCATAATCAATATCAGGCAGTGAAAAt                                >  2:1074059/1‑139 (MQ=255)                                                                                                   cccAGTACGGCCATTCACCAATCCCTTACTTTTGTTGCGCAAACTCTGCCAGCAACGGCCAGCGTTTTAATGCCGCGGCCAACTGTTCGCGTTGGGAAATATCGGCATAATCAATATCAGGCAGTGAAAATGCCTGCtt                        >  2:655754/1‑139 (MQ=255)                                                                                                     cAGTACGGCCATTCACCAATCCCTTACTTTTGTTGCGCAAACTCTGCCAGCAACGGCCAGCGTTTTAATGCCGCGGCCAACTGTTCGCGTTGGGAAATATCGGCATAATCAATATCAGGCAGTGAAAATGCCTGCTTTa                      >  1:565195/1‑139 (MQ=255)                                                                                                      aGTACGGCCATTCACCAATCCCTTACTTTTGTTGCGCAAACTCTGCCAGCAACGGCCAGCGTTTTAATGCCGCGGCCAACTGTTCGCGTTGGGAAATATCGGCATAATCAATATCAGGCAGTGAAAATGCCTGCTTTaa                     <  2:810710/139‑1 (MQ=255)                                                                                                      aGTACGGCCATTCACCAATCCCTTACTTTTGTTGCGCAAACTCTGCCAGCAACGGCCAGCGTTTTAATGCCGCGGCCAACTGTTCGCGTTGGGAAATATCGGCATAATCAATATCAGGCAGTGAAAATGCCTGCTTTaa                     <  2:436623/139‑1 (MQ=255)                                                                                                       gTACGGCCATTCACCAATCCGTTACTTTTGTTGCGCAAACTCTGCCAGCAACGGCCAGCGTTTTAATGCCGCGGCCAACTGTTCGCGTTGGGAAATATCGGCATAATCAATATCAGGCAGTGAAAATGCCTGCTTTAAc                    >  2:653057/1‑139 (MQ=255)                                                                                                       gTACGGCCATTCACCAATCCCTTACTTTTGTTGCGCAAACTCTGCCAGCAACGGCCAGCGTTTTAATGCCGCGGCCAACTGTTCGCGTTGGGAAATATCGGCATAATCAATATCAGGCAGTGAAAATGCCTGCTTTAAc                    >  2:1028595/1‑139 (MQ=255)                                                                                                        tACGGCCATTCACCAATCCCTTACTTTTGTTGCGCAAACTCTGc                                                                                                                  >  2:656739/1‑44 (MQ=255)                                                                                                        tACGGCCATTCACCAATCCCTTACTTTTGTTGCGCAAACTCTGc                                                                                                                  <  1:656739/44‑1 (MQ=255)                                                                                                         aCGGCCATTCACCAATCCCTTACTTTTGTTGCGCAAACTCTGCCAGCAACGGCCAGCGTTTTAATGCCGCGGCCAACTGTTCGCGTTGGGAAATATCGGCATAATCAATATCAGGCAGTGAAAATGCCTGCTTTaa                     <  1:910503/136‑1 (MQ=255)                                                                                                         aCGGCCATTCACCAATCCCTTACTTTTGTTGCGCAAACTCTGCCAGCAACGGCCAGCGTTTTAATGCCGCGGCCAACTGTTCGCGTTGGGAAATATCGGCATAATCAATATCAGGCAGTGAAAATGCCTGCTTTaa                     >  2:910503/1‑136 (MQ=255)                                                                                                         aCGGCCATTCACCAATCCCTTACTTTTGTTGCGCAAACTCTGCCAGCAACGGCCAGCGTTTTAATGCCGCGGCCAACTGTTCGCGTTGGGAAATATCGGCATAATCAATATCAGGCAGTGAAAATGCCTGCTTTAACGc                  >  2:1008126/1‑139 (MQ=255)                                                                                                             ccATTCACCAATCCCTTACTTTTGTTGCGCAAACTCTGCCAGCAACGGCCAGCGTTTTAATGCCGCGGCCAACTGTTCGCGTTGGGAAATATCGGc                                                         >  2:704302/1‑96 (MQ=255)                                                                                                             ccATTCACCAATCCCTTACTTTTGTTGCGCAAACTCTGCCAGCAACGGCCAGCGTTTTAATGCCGCGGCCAACTGTTCGCGTTGGGAAATATCGGc                                                         <  1:704302/96‑1 (MQ=255)                                                                                                                ttCACCAATCCCTTACTTTTGTTGCGCAGACTCTGCCAGCAACGGCCAGCGTTTTAATGCCGCGGCCAACTGTTCGCGTTGGGAAATATCGGCATAATCAATATCAGGCAGTGAAAATGCCTGCTTTAACGCCACAata           <  1:667192/139‑1 (MQ=255)                                                                                                                 tCACCAATCCCTTACTTTTGTTGCGCAAACTCTGCCAGCAACGGCCAGCGTTTTAATGCCGCGGCCAACTGTTCGCGTTGGGAAATATCGGCATAATCAATATCAGGCAGTGAAAATGCCTGCTTTAACGCCACAatat          <  1:1028595/139‑1 (MQ=255)                                                                                                                  cACCAATCCCTTACTTTTGTTGCGCAAACTCTGCCAGCAACGGCCAGCGTTTTAATGCCGCGGCCAACTGTTCGCGTTGGGAAATATCGGCATAATCAATATCAGGCAGTGAAAATGCCTGCTTTAACGc                  <  2:799589/130‑1 (MQ=255)                                                                                                                  cACCAATCCCTTACTTTTGTTGCGCAAACTCTGCCAGCAACGGCCAGCGTTTTAATGCCGCGGCCAACTGTTCGCGTTGGGAAATATCGGCATAATCAATATCAGGCAGTGAAAATGCCTGCTTTAACGc                  >  1:799589/1‑130 (MQ=255)                                                                                                                     cAATCCCTTACTTTTGTTGCGCAAACTCTGCCAGCAACGGCCAGCGTTTTAATGCCGCGGCCAACTGTTCGCGTTGGGAAATATCGGCATAATCAATATCAGGCAGTGAAAATGCCTGCTTTAACGCCACAATATCAtt      <  1:1074059/139‑1 (MQ=255)                                                                                                                     cAATCCCTTACTTTTGTTGCGCAAACTCTGCCAGCAACGGCCAGCGTTTTAATGCCGCGGCCAACTGTTCGCGTTGGGAAATATCGGCATAATCAATATCAGGCAGTGAAAATGCCTGCTTTAACGCCACAATATCAtt      <  1:655754/139‑1 (MQ=255)                                                                                                                       aTCCCTTACTTTTGTTGCGCAAACTCTGCCAGCAACGGCCAGCGTTTTAATGCCGCGGCCAACTGTTCGCGTTGGGAAATATCGGc                                                         >  2:532876/1‑86 (MQ=255)                                                                                                                       aTCCCTTACTTTTGTTGCGCAAACTCTGCCAGCAACGGCCAGCGTTTTAATGCCGCGGCCAACTGTTCGCGTTGGGAAATATCGGc                                                         <  1:532876/86‑1 (MQ=255)                                                                                                                       aTCCCTTACTTTTGTTGCGCAAACTCTGCCAGCAACGGCCAGCGTTTTAATGCCGCGGCCAACTGTTCGCGTTGGGAAATATCGGCATAATCAATATCAGGCAGTGAAAATGCCTGCTTTAACGCCACAATATCATTCt    >  1:231822/1‑139 (MQ=255)                                                                                                                         ccctTACTTTTGTTGCGCAAACTCTGCCAGCAACGGCCAGCGTTTTAATGCCGCGGCCAACTGTTCGCGTTGGGAAATATCGGCATAATCAATATCAGGCAGTGAAAATGCCTGCTTTAACGCCACAATATCATTCTgg  <  2:947713/139‑1 (MQ=255)                                                                                                                            |                                                                                                                                        CATCGACCACCAGGACATTTTCTCCCAACATTTGTAATGACCAGGCTAATGCGGCGGTGATGGTTGTTGTCCCCACGCCTCCCCGCACCCCCTGCTATCCCAGTACGGCCATTCACCAATCCCCTACTTTTGTTGCGCAAACTCTGCCAGCAACGGCCAGCGTTTTAATGCCGCGGCCAACTGTTCGCGTTGGGAAATATCGGCATAATCAATATCAGGCAGTGAAAATGCCTGCTTTAACGCCACAATATCATTCTGG  >  NC_000913/3695874‑3696132

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